portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2026-05-16 07:04:32 +01:00
Code Releases Activity

small bug fix for wilson spectrum since we're actually running DWF

Browse Source
This commit is contained in:
Patrick R Oare
2025-09-25 15:36:42 -04:00
parent fa30c791aa
commit 82f35001ff
1 changed files with 9 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files
examples/Example_wilson_spectrum.cc
+9 -3
View File
@@ -288,8 +288,8 @@ int main (int argc, char ** argv)
GridDefaultSimd(Nd,vComplex::Nsimd()), GridDefaultSimd(Nd,vComplex::Nsimd()),
GridDefaultMpi()); GridDefaultMpi());
GridRedBlackCartesian * UrbGrid = SpaceTimeGrid::makeFourDimRedBlackGrid(UGrid); GridRedBlackCartesian * UrbGrid = SpaceTimeGrid::makeFourDimRedBlackGrid(UGrid);
GridCartesian* FGrid = UGrid; GridCartesian * FGrid = SpaceTimeGrid::makeFiveDimGrid(Ls,UGrid);
GridRedBlackCartesian* FrbGrid = UrbGrid; GridRedBlackCartesian * FrbGrid = SpaceTimeGrid::makeFiveDimRedBlackGrid(Ls,UGrid);
std::vector<int> seeds4({1,2,3,4}); std::vector<int> seeds4({1,2,3,4});
GridParallelRNG RNG4(UGrid); GridParallelRNG RNG4(UGrid);
@@ -301,13 +301,19 @@ int main (int argc, char ** argv)
FieldMetaData header; FieldMetaData header;
NerscIO::readConfiguration(Umu, header, file); NerscIO::readConfiguration(Umu, header, file);
std::cout << GridLogMessage << "Loaded configuration" << std::endl;
RealD mass = 0.01; RealD mass = 0.01;
RealD M5 = 1.8; RealD M5 = 1.8;
// Define domain wall D std::cout << GridLogMessage << "masses specified" << std::endl;
// Define domain wall D. This is giving a floating point exception?
DomainWallFermionD Ddwf(Umu, *FGrid, *FrbGrid, *UGrid, *UrbGrid, mass, M5); DomainWallFermionD Ddwf(Umu, *FGrid, *FrbGrid, *UGrid, *UrbGrid, mass, M5);
NonHermitianLinearOperator<DomainWallFermionD, LatticeFermionD> DLinOp (Ddwf); NonHermitianLinearOperator<DomainWallFermionD, LatticeFermionD> DLinOp (Ddwf);
std::cout << GridLogMessage << "DWF operator defined" << std::endl;
// Define PV^dag D (if we want) // Define PV^dag D (if we want)
DomainWallFermionD Dpv(Umu, *FGrid, *FrbGrid, *UGrid, *UrbGrid, 1.0, M5); DomainWallFermionD Dpv(Umu, *FGrid, *FrbGrid, *UGrid, *UrbGrid, 1.0, M5);
typedef PVdagMLinearOperator<DomainWallFermionD,LatticeFermionD> PVdagM_t; typedef PVdagMLinearOperator<DomainWallFermionD,LatticeFermionD> PVdagM_t;