portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-06-14 13:57:07 +01:00
Code Releases Activity

Skeleton of HMC/Integrators

Browse Source
This commit is contained in:
neo
2015-07-03 16:51:41 +09:00
parent a666e66e36
commit 9655d43017
13 changed files with 263 additions and 15 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -543,6 +543,26 @@ Note that in step D setting B ~ X - A and using B in place of A in step E will g
taExp(lie,out);
}
static void GaussianLieAlgebraMatrix(GridParallelRNG &pRNG,LatticeMatrix &out,double scale=1.0){
GridBase *grid = out._grid;
LatticeComplex ca (grid);
LatticeMatrix lie(grid);
LatticeMatrix la (grid);
Complex ci(0.0,scale);
Matrix ta;
out=zero;
for(int a=0;a<generators();a++){
gaussian(pRNG,ca);
generator(a,ta);
la=ci*ca*ta;
out = lie+la; // e^{i la ta}
}
}
static void HotConfiguration(GridParallelRNG &pRNG,LatticeGaugeField &out){
LatticeMatrix Umu(out._grid);
for(int mu=0;mu<Nd;mu++){