Utility implementing Fischer's RHS / solution projection scheme. More...
#include "BelosConfigDefs.hpp"#include "BelosTypes.hpp"#include "BelosLinearProblem.hpp"#include "BelosMultiVecTraits.hpp"#include "BelosOperatorTraits.hpp"#include "BelosTeuchosDenseAdapter.hpp"#include "BelosKokkosDenseAdapter.hpp"#include "Teuchos_ParameterList.hpp"#include "Teuchos_RCP.hpp"#include "Teuchos_ScalarTraits.hpp"#include <string>#include <vector>
Go to the source code of this file.
Classes | |
| class | Belos::SolutionProjectionFailure |
| class | Belos::SolutionProjection< ScalarType, MV, OP, DM > |
| Fischer Method 1 style projection utility. More... | |
Namespaces | |
| namespace | Belos |
Detailed Description
Utility implementing Fischer's RHS / solution projection scheme.
This utility stores paired bases U and C satisfying
A U = C, C^H C = I,
and applies the projection
x <- x + U C^H r, r <- r - C C^H r,
where r = b - A x.
This is Fischer's Method 1 when U stores previous solution-like vectors and C stores their corresponding operator images.
The projection tolerance is relative: a candidate image vector c is rejected if, after orthogonalization against the current C basis,
||c_orth|| <= projectionTol * ||c_original||.
The default replacement strategy is "Restart": if the basis is full and a new vector is accepted, the previous basis is discarded and the new vector starts the next basis. Thus, when the history space is full, the next accepted vector "kicks off" the new space using the most recent information.
Definition in file BelosSolutionProjection.hpp.
