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

Print warning messages in CAGMRES solvers

Browse Source 
Currently, the implementation of these algorithms doesn't differ from their non
communication-avoiding versions.
This commit is contained in:
Daniel Richtmann 2018-02-08 17:43:47 +01:00
parent 13ae371ef8
commit eb7cf239d9
No known key found for this signature in database
GPG Key ID: B33C490AF0772057
2 changed files with 4 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
lib/algorithms/iterative/CommunicationAvoidingGeneralisedMinimalResidual.h
View File

@ -72,6 +72,8 @@ class CommunicationAvoidingGeneralisedMinimalResidual : public OperatorFunction<
void operator()(LinearOperatorBase<Field> &LinOp, const Field &src, Field &psi) { void operator()(LinearOperatorBase<Field> &LinOp, const Field &src, Field &psi) {
std::cout << GridLogWarning << "This algorithm currently doesn't differ from regular GMRES" << std::endl;
psi.checkerboard = src.checkerboard; psi.checkerboard = src.checkerboard;
conformable(psi, src); conformable(psi, src);

2
lib/algorithms/iterative/FlexibleCommunicationAvoidingGeneralisedMinimalResidual.h
View File

@ -77,6 +77,8 @@ class FlexibleCommunicationAvoidingGeneralisedMinimalResidual : public OperatorF
void operator()(LinearOperatorBase<Field> &LinOp, const Field &src, Field &psi) { void operator()(LinearOperatorBase<Field> &LinOp, const Field &src, Field &psi) {
std::cout << GridLogWarning << "This algorithm currently doesn't differ from regular FGMRES" << std::endl;
psi.checkerboard = src.checkerboard; psi.checkerboard = src.checkerboard;
conformable(psi, src); conformable(psi, src);