diff --git a/lib/algorithms/iterative/CommunicationAvoidingGeneralisedMinimalResidual.h b/lib/algorithms/iterative/CommunicationAvoidingGeneralisedMinimalResidual.h index b56c8afc..4469c7ed 100644 --- a/lib/algorithms/iterative/CommunicationAvoidingGeneralisedMinimalResidual.h +++ b/lib/algorithms/iterative/CommunicationAvoidingGeneralisedMinimalResidual.h @@ -72,6 +72,8 @@ class CommunicationAvoidingGeneralisedMinimalResidual : public OperatorFunction< 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; conformable(psi, src); diff --git a/lib/algorithms/iterative/FlexibleCommunicationAvoidingGeneralisedMinimalResidual.h b/lib/algorithms/iterative/FlexibleCommunicationAvoidingGeneralisedMinimalResidual.h index f49eb7dc..dc4885d6 100644 --- a/lib/algorithms/iterative/FlexibleCommunicationAvoidingGeneralisedMinimalResidual.h +++ b/lib/algorithms/iterative/FlexibleCommunicationAvoidingGeneralisedMinimalResidual.h @@ -77,6 +77,8 @@ class FlexibleCommunicationAvoidingGeneralisedMinimalResidual : public OperatorF 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; conformable(psi, src);