#include <F_GcrLS.hpp>
Inheritance diagram for F_GcrLS:
Purpose: GCR(m) krylov solver with restart each m iterations, for non symmetric/indefinite linear system
Used to solve a set of linearized differencial equations,....
GCR is a Krylov solver working on residuals. It approximates a solution via generation of the suitable
partial Krylov orthonormal base, especially in the restart case.
We hope that a small number of Krylov sub-spaces are sufficient to lead to a good solution.
GCR is competitive with GMRES, but it may breakdown.
Papers: "Modified Generalized Conjugate Residuals Method for Nonsymmetric Systems of Linear Equations."
Zbigniew Leyk, report 1993, Australian National University.
"Conjugate Direction Waveform Methods for Transient 2 Dimensional Simulation of MOS Devices.",
A. Lumsdaine, M. Reichelt, J. White, Research Lab of Electronics,
Massachusetts Institut of Technology, Cmbridge, MA 02139, IEEE 1991.
Public Methods | |
| F_GcrLS (int maxIterations, int m=10, double epsilon=1.0e-4) | |
| virtual | ~F_GcrLS () |
| virtual FloatVector * | solve (AFloatMatrix *m, FloatVector *b, FloatVector *dest) |
Protected Methods | |
| virtual FloatVector * | singleSolve (AFloatMatrix *A, FloatVector *b, FloatVector *x=0) |
Private Attributes | |
| int | m |
| int | k |
| FloatVector * | r |
| FloatVector * | ar |
| FloatVector ** | p |
| FloatVector ** | ap |
| double * | apap |
| AFSymMatrix * | h |
 |
