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

Making tests compile

Browse Source
This commit is contained in:
Chulwoo Jung
2018-03-08 23:02:19 -05:00
parent 3686827df5
commit 8bb11e5039
5 changed files with 28 additions and 24 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
tests/lanczos/Test_dwf_block_lanczos.cc
View File

@ -46,7 +46,7 @@ class CmdJobParams
double M5;
double mob_b;
std::vector<ComplexD> omega;
std::vector<ComplexD> boundary_phase;
std::vector<Complex> boundary_phase;
LanczosType Impl;
int Nu;

1
tests/lanczos/Test_dwf_compressed_lanczos.cc
View File

@ -22,6 +22,7 @@
*/
#include <Grid/Grid.h>
#include <Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h>
#include <Grid/parallelIO/IldgIO.h>
/////////////////////////////////////////////////////////////////////////////
// The following are now decoupled from the Lanczos and deal with grids.
// Safe to replace functionality

4
tests/lanczos/Test_dwf_compressed_lanczos_reorg.cc
View File

@ -33,6 +33,8 @@ Author: Peter Boyle <paboyle@ph.ed.ac.uk>
#include <Grid/Grid.h>
#include <Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h>
#include <Grid/algorithms/iterative/LocalCoherenceLanczos.h>
#include <Grid/parallelIO/IldgIOtypes.h>
#include <Grid/parallelIO/IldgIO.h>
using namespace std;
using namespace Grid;
@ -57,7 +59,7 @@ public:
void checkpointFine(std::string evecs_file,std::string evals_file)
{
assert(this->subspace.size()==nbasis);
emptyUserRecord record;
Grid::emptyUserRecord record;
Grid::QCD::ScidacWriter WR;
WR.open(evecs_file);
for(int k=0;k<nbasis;k++) {

2
tests/lanczos/Test_dwf_lanczos.cc
View File

@ -84,7 +84,7 @@ int main (int argc, char ** argv)
std::vector<double> Coeffs { 0.,-1.};
Polynomial<FermionField> PolyX(Coeffs);
Chebyshev<FermionField> Cheb(0.2,5.5,11);
Chebyshev<FermionField> Cheby(0.2,5.5,11);
FunctionHermOp<FermionField> OpCheby(Cheby,HermOp);
PlainHermOp<FermionField> Op (HermOp);

43
tests/solver/Test_dwf_mrhs_cg.cc
View File

@ -26,17 +26,18 @@ Author: Peter Boyle <paboyle@ph.ed.ac.uk>
*************************************************************************************/
/* END LEGAL */
#include <Grid/Grid.h>
#include <Grid/parallelIO/IldgIOtypes.h>
#include <Grid/algorithms/iterative/BlockConjugateGradient.h>
using namespace std;
//using namespace std;
using namespace Grid;
using namespace Grid::QCD;
int main (int argc, char ** argv)
{
typedef typename DomainWallFermionR::FermionField FermionField;
typedef typename DomainWallFermionR::ComplexField ComplexField;
typename DomainWallFermionR::ImplParams params;
typedef typename Grid::QCD::DomainWallFermionR::FermionField FermionField;
typedef typename Grid::QCD::DomainWallFermionR::ComplexField ComplexField;
typename Grid::QCD::DomainWallFermionR::ImplParams params;
const int Ls=4;
@ -102,7 +103,7 @@ int main (int argc, char ** argv)
std::cout << GridLogMessage << "****************************************************************** "<<std::endl;
std::cout << GridLogMessage << " Writing out in parallel view "<<std::endl;
std::cout << GridLogMessage << "****************************************************************** "<<std::endl;
emptyUserRecord record;
// Grid::emptyUserRecord record;
std::string file("./scratch.scidac");
std::string filef("./scratch.scidac.ferm");
@ -114,21 +115,21 @@ int main (int argc, char ** argv)
{
FGrid->Barrier();
ScidacWriter _ScidacWriter;
_ScidacWriter.open(file);
// ScidacWriter _ScidacWriter;
// _ScidacWriter.open(file);
std::cout << GridLogMessage << "****************************************************************** "<<std::endl;
std::cout << GridLogMessage << " Writing out gauge field "<<std::endl;
std::cout << GridLogMessage << "****************************************************************** "<<std::endl;
_ScidacWriter.writeScidacFieldRecord(Umu,record);
_ScidacWriter.close();
// _ScidacWriter.writeScidacFieldRecord(Umu,record);
// _ScidacWriter.close();
FGrid->Barrier();
std::cout << GridLogMessage << "****************************************************************** "<<std::endl;
std::cout << GridLogMessage << " Reading in gauge field "<<std::endl;
std::cout << GridLogMessage << "****************************************************************** "<<std::endl;
ScidacReader _ScidacReader;
_ScidacReader.open(file);
_ScidacReader.readScidacFieldRecord(s_Umu,record);
_ScidacReader.close();
// ScidacReader _ScidacReader;
// _ScidacReader.open(file);
// _ScidacReader.readScidacFieldRecord(s_Umu,record);
// _ScidacReader.close();
FGrid->Barrier();
std::cout << GridLogMessage << "****************************************************************** "<<std::endl;
std::cout << GridLogMessage << " Read in gauge field "<<std::endl;
@ -144,10 +145,10 @@ int main (int argc, char ** argv)
std::cout << GridLogMessage << "****************************************************************** "<<std::endl;
std::stringstream filefn; filefn << filef << "."<< n;
ScidacWriter _ScidacWriter;
_ScidacWriter.open(filefn.str());
_ScidacWriter.writeScidacFieldRecord(src[n],record);
_ScidacWriter.close();
// ScidacWriter _ScidacWriter;
// _ScidacWriter.open(filefn.str());
// _ScidacWriter.writeScidacFieldRecord(src[n],record);
// _ScidacWriter.close();
}
FGrid->Barrier();
@ -159,10 +160,10 @@ int main (int argc, char ** argv)
for(int n=0;n<nrhs;n++){
if ( n==me ) {
std::stringstream filefn; filefn << filef << "."<< n;
ScidacReader _ScidacReader;
_ScidacReader.open(filefn.str());
_ScidacReader.readScidacFieldRecord(s_src,record);
_ScidacReader.close();
// ScidacReader _ScidacReader;
// _ScidacReader.open(filefn.str());
// _ScidacReader.readScidacFieldRecord(s_src,record);
// _ScidacReader.close();
}
}
FGrid->Barrier();