portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-11 22:50:45 +01:00
Code Releases Activity

No compile babbage fix

Browse Source
This commit is contained in:
paboyle 2015-09-30 16:02:44 -07:00
parent f4b6d1dfea
commit a660ce716b
1 changed files with 7 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
lib/qcd/utils/SUn.h
View File

@ -524,16 +524,22 @@ Note that in step D setting B ~ X - A and using B in place of A in step E will g
// reunitarise?? // reunitarise??
static void LieRandomize(GridParallelRNG &pRNG,LatticeMatrix &out,double scale=1.0){ static void LieRandomize(GridParallelRNG &pRNG,LatticeMatrix &out,double scale=1.0){
GridBase *grid = out._grid; GridBase *grid = out._grid;
LatticeComplex ca (grid); LatticeComplex ca (grid);
LatticeMatrix lie(grid); LatticeMatrix lie(grid);
LatticeMatrix la (grid); LatticeMatrix la (grid);
Complex ci(0.0,scale); Complex ci(0.0,scale);
Complex cone(1.0,0.0);
Matrix ta; Matrix ta;
lie=zero; lie=zero;
for(int a=0;a<generators();a++){ for(int a=0;a<generators();a++){
random(pRNG,ca); ca=real(ca)-0.5; random(pRNG,ca);
ca = (ca+conjugate(ca))*0.5;
ca = ca - 0.5;
generator(a,ta); generator(a,ta);
la=ci*ca*ta; la=ci*ca*ta;