Hi everyone. We wanted to change the code to do some more calculations for nearest neighbor in spinmc. So far I've tracked down the problem to the file alps/applications/mc/spins/localupdate.h. We've also added the edge type for nearest neighbour in xml file defining the unit lattice. The problem is that getting edge type for the graph always returns zero. Can you please help me find our how should we track down the problem?
Getting the edge type would be like this:
int edgeType = boost::get(alps::edge_type_t(), graph, *iter); // iter is out_edge_iterator
but the return value is always zero. (executing the exact same code for when the library parses the xml file produces the correct results. namely in file alps/src/alps/lattice/unitcell.c)
Yours, Kourosh