Commit 51217eee by Francois Gygi

Fix BasisMapping single-task map for np0!=np1

parent 712c6f87
...@@ -208,8 +208,8 @@ BasisMapping::BasisMapping (const Basis &basis, int np0, int np1, int np2) : ...@@ -208,8 +208,8 @@ BasisMapping::BasisMapping (const Basis &basis, int np0, int np1, int np2) :
if (hm < 0) hm += np0_; if (hm < 0) hm += np0_;
if (km < 0) km += np1_; if (km < 0) km += np1_;
zvec_to_val_.push_back(hp+np1_*kp); zvec_to_val_.push_back(hp+np0_*kp);
zvec_to_val_.push_back(hm+np1_*km); zvec_to_val_.push_back(hm+np0_*km);
} }
} }
else else
...@@ -222,7 +222,7 @@ BasisMapping::BasisMapping (const Basis &basis, int np0, int np1, int np2) : ...@@ -222,7 +222,7 @@ BasisMapping::BasisMapping (const Basis &basis, int np0, int np1, int np2) :
k = basis_.rod_k(0,irod); k = basis_.rod_k(0,irod);
if (h < 0) h += np0_; if (h < 0) h += np0_;
if (k < 0) k += np1_; if (k < 0) k += np1_;
zvec_to_val_.push_back(h+np1_*k); zvec_to_val_.push_back(h+np0_*k);
} }
} }
} }
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment