portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-09 21:50:45 +01:00
Code Releases Activity

Patch to broken assertion

Browse Source
This commit is contained in:
Peter Boyle 2018-11-06 22:18:17 +00:00
parent 8c3a599148
commit b62b9ac214
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
Grid/algorithms/iterative/SchurRedBlack.h
View File

@ -474,16 +474,17 @@ namespace Grid {
Field tmp(grid); Field tmp(grid);
guess(src_o,sol_o); guess(src_o,sol_o);
// std::cout<<GridLogMessage << "SchurRedBlack saving the guess" <<std::endl;
Field guess_save(grid); Field guess_save(grid);
// std::cout<<GridLogMessage << "SchurRedBlack operator =saving the guess" <<std::endl;
guess_save = sol_o; guess_save = sol_o;
std::cout<<GridLogMessage << "SchurRedBlack solver calling the MpcDagMp solver" <<std::endl; std::cout<<GridLogMessage << "SchurRedBlack solver calling the MpcDagMp solver" <<std::endl;
////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////
// Call the red-black solver // Call the red-black solver
////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////
SchurDiagTwoOperator<Matrix,Field> _HermOpEO(_Matrix); SchurDiagTwoOperator<Matrix,Field> _HermOpEO(_Matrix);
_HermitianRBSolver(_HermOpEO,src_o,sol_o); assert(tmp.checkerboard==Odd); _HermitianRBSolver(_HermOpEO,src_o,sol_o); assert(sol_o.checkerboard==Odd);
//////////////////////////////////////////////// ////////////////////////////////////////////////
// Fionn A2A boolean behavioural control // Fionn A2A boolean behavioural control