From eb7cf239d93b2ec5576871737a3f9bf8fb4f980f Mon Sep 17 00:00:00 2001 From: Daniel Richtmann Date: Thu, 8 Feb 2018 17:43:47 +0100 Subject: [PATCH] Print warning messages in CAGMRES solvers Currently, the implementation of these algorithms doesn't differ from their non communication-avoiding versions. --- .../iterative/CommunicationAvoidingGeneralisedMinimalResidual.h | 2 ++ .../FlexibleCommunicationAvoidingGeneralisedMinimalResidual.h | 2 ++ 2 files changed, 4 insertions(+) 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 &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 &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);