From 82f35001ff03da11b3ff17123a8d6281b80e6a2e Mon Sep 17 00:00:00 2001 From: Patrick R Oare Date: Thu, 25 Sep 2025 15:36:42 -0400 Subject: [PATCH] small bug fix for wilson spectrum since we're actually running DWF --- examples/Example_wilson_spectrum.cc | 12 +++++++++--- 1 file changed, 9 insertions(+), 3 deletions(-) diff --git a/examples/Example_wilson_spectrum.cc b/examples/Example_wilson_spectrum.cc index 6709cce7..bd356bc5 100644 --- a/examples/Example_wilson_spectrum.cc +++ b/examples/Example_wilson_spectrum.cc @@ -288,8 +288,8 @@ int main (int argc, char ** argv) GridDefaultSimd(Nd,vComplex::Nsimd()), GridDefaultMpi()); GridRedBlackCartesian * UrbGrid = SpaceTimeGrid::makeFourDimRedBlackGrid(UGrid); - GridCartesian* FGrid = UGrid; - GridRedBlackCartesian* FrbGrid = UrbGrid; + GridCartesian * FGrid = SpaceTimeGrid::makeFiveDimGrid(Ls,UGrid); + GridRedBlackCartesian * FrbGrid = SpaceTimeGrid::makeFiveDimRedBlackGrid(Ls,UGrid); std::vector seeds4({1,2,3,4}); GridParallelRNG RNG4(UGrid); @@ -301,13 +301,19 @@ int main (int argc, char ** argv) FieldMetaData header; NerscIO::readConfiguration(Umu, header, file); + std::cout << GridLogMessage << "Loaded configuration" << std::endl; + RealD mass = 0.01; 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); NonHermitianLinearOperator DLinOp (Ddwf); + std::cout << GridLogMessage << "DWF operator defined" << std::endl; + // Define PV^dag D (if we want) DomainWallFermionD Dpv(Umu, *FGrid, *FrbGrid, *UGrid, *UrbGrid, 1.0, M5); typedef PVdagMLinearOperator PVdagM_t;