portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-06-10 19:36:56 +01:00
Code Releases Activity

Bug in two flav pseudofermion corrected to reimport gauge field upon rejection.

Browse Source 
exp(-DeltaH) = 1 now, and plaquette is sensible. Will reproduce an old Wilson Gauge
Wilson Fermion SCRI plaquette with precision in mass matching shortly.
This commit is contained in:
Peter Boyle
2015-07-29 21:02:07 +09:00
parent 4fe110bd07
commit 9ff0b2987c
4 changed files with 10 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
lib/qcd/hmc/HMC.h
View File

@ -94,6 +94,7 @@ namespace Grid{
void evolve(LatticeLorentzColourMatrix& Uin){
Real DeltaH;
// Thermalizations
@ -106,14 +107,16 @@ namespace Grid{
// Actual updates (evolve a copy Ucopy then copy back eventually)
LatticeLorentzColourMatrix Ucopy(Uin._grid);
for(int iter=Params.StartingConfig;
iter < Params.Nsweeps+Params.StartingConfig; ++iter){
for(int iter=Params.StartingConfig; iter < Params.Nsweeps+Params.StartingConfig; ++iter){
std::cout<<GridLogMessage << "-- # Sweep = "<< iter << "\n";
Ucopy = Uin;
DeltaH = evolve_step(Ucopy);
if(metropolis_test(DeltaH)) Uin = Ucopy;
//need sync? // Query Guido on this comment
}

2
lib/qcd/hmc/integrators/Integrator_base.h
View File

@ -60,10 +60,10 @@ namespace Grid{
IntegratorPolicy TheIntegrator;// contains parameters too
void update_P(LatticeLorentzColourMatrix&U, int level,double ep){
for(int a=0; a<as[level].size(); ++a){
LatticeLorentzColourMatrix force(U._grid);
as[level].at(a)->deriv(U,force);
Complex dSdt=0.0;