Hi everyone! I'm trying to write a code to do different calculations for different bond types, so I changed my lattice file and set edge type to 8888 for test, a part of lattice file is as follows:
<EDGE type="8888"><SOURCE vertex="1"/><TARGET vertex="2" offset=" 0 0 0"/></EDGE> <EDGE type="8888"><SOURCE vertex="3"/><TARGET vertex="4" offset=" 0 0 0"/></EDGE> <EDGE type="8888"><SOURCE vertex="2"/><TARGET vertex="1" offset=" 1 0 0"/></EDGE>
and I added a test code in alps/applications/mc/spins/spinsim.h as following:
typedef typename AbstractSpinSim<MIdMatrix<double, M::dim>
::bond_iterator bond_iterator;
for(bond_iterator bi = this->bonds().first; bi != this->bonds().second; bi++) { int et = this->edge_type(*bi);
std::cout << et << std::endl; }
but all I get from bond types is zero. Do you know how can I get the correct bond type here?
Yours, Kourosh