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

Updated python output

Browse Source
This commit is contained in:
Peter Boyle
2025-06-03 14:09:02 -04:00
parent d418f78352
commit 3737a24096
1 changed files with 17 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

19
tests/lanczos/Test_dwf_G5R5.cc
View File

@ -314,7 +314,7 @@ int main(int argc, char** argv) {
} }
} }
} }
for(int i = 0; i < Nconv; i++){ for(int i = 0; i < Nconv; i++){
G5R5Herm.HermOpAndNorm(finalevec[i], G5R5Mevec[i], eMe[i], eMMe[i]); G5R5Herm.HermOpAndNorm(finalevec[i], G5R5Mevec[i], eMe[i], eMMe[i]);
} }
cout << "Re<evec, G5R5M(evec)>: " << endl; cout << "Re<evec, G5R5M(evec)>: " << endl;
@ -322,6 +322,7 @@ int main(int argc, char** argv) {
cout << "<G5R5M(evec), G5R5M(evec)>" << endl; cout << "<G5R5M(evec), G5R5M(evec)>" << endl;
cout << eMMe << endl; cout << eMMe << endl;
// vector<LatticeFermion> finalevec(Nconv, FGrid); // vector<LatticeFermion> finalevec(Nconv, FGrid);
// temporary, until doing rotation // temporary, until doing rotation
@ -343,6 +344,20 @@ int main(int argc, char** argv) {
} }
} }
for(int i = 0; i < Nconv; i++){
Ddwf.M(finalevec[i], G5R5Mevec[i]);
for(int j = 0; j < Nconv; j++){
std::cout << "<"<<j<<"|Ddwf|"<<i<<"> = "<<innerProduct(finalevec[j],G5R5Mevec[i])<<std::endl;
}
}
for(int i = 0; i < Nconv; i++){
RealD t1,t2;
G5R5Herm.HermOpAndNorm(finalevec[i], G5R5Mevec[i], t1, t2);
for(int j = 0; j < Nconv; j++){
std::cout << "<"<<j<<"|G5R5 M|"<<i<<"> = "<<innerProduct(finalevec[j],G5R5Mevec[i])<<std::endl;
}
}
for(int i = 0; i < Nconv; i++){ for(int i = 0; i < Nconv; i++){
chiral_matrix_real[i].resize(Nconv); chiral_matrix_real[i].resize(Nconv);
chiral_matrix[i].resize(Nconv); chiral_matrix[i].resize(Nconv);
@ -380,7 +395,7 @@ int main(int argc, char** argv) {
PYTHON_LINE("ax = fig.add_subplot(projection='3d')"); PYTHON_LINE("ax = fig.add_subplot(projection='3d')");
PYTHON_LINE(""); PYTHON_LINE("");
PYTHON_LINE("x, y = np.random.rand(2, 100) * 4"); PYTHON_LINE("x, y = np.random.rand(2, 100) * 4");
PYTHON_LINE("hist, xedges, yedges = np.histogram2d(x, y, bins=10, range=[[0, 9], [0, 9]])"); fprintf(fp,"hist, xedges, yedges = np.histogram2d(x, y, bins=%d, range=[[0, %d], [0, %d]])\n",Nconv,Nconv-1,Nconv-1);
PYTHON_LINE(""); PYTHON_LINE("");
PYTHON_LINE("# Construct arrays for the anchor positions of the 16 bars"); PYTHON_LINE("# Construct arrays for the anchor positions of the 16 bars");
PYTHON_LINE("xpos, ypos = np.meshgrid(xedges[:-1] + 0.25, yedges[:-1] + 0.25, indexing=\"ij\")"); PYTHON_LINE("xpos, ypos = np.meshgrid(xedges[:-1] + 0.25, yedges[:-1] + 0.25, indexing=\"ij\")");