Go to the documentation of this file.
18 #ifndef __PASO_PRECONDITIONER_H__
19 #define __PASO_PRECONDITIONER_H__
27 struct Preconditioner;
31 struct Preconditioner_Smoother;
81 dim_t sweeps,
bool x_is_initial);
85 dim_t sweeps,
bool x_is_initial);
89 double atol,
dim_t* sweeps,
bool x_is_initial);
146 #endif // __PASO_PRECONDITIONER_H__
index_t * mask_C
Definition: Preconditioner.h:136
dim_t n_block
Definition: Preconditioner.h:126
dim_t n
Definition: Preconditioner.h:125
dim_t n_C
Definition: Preconditioner.h:128
Solver_ILU * Solver_getILU(SparseMatrix_ptr A, bool verbose)
constructs the incomplete block factorization
Definition: ILU.cpp:57
Preconditioner_LocalSmoother * Preconditioner_LocalSmoother_alloc(SparseMatrix_ptr A, bool jacobi, bool verbose)
Definition: Smoother.cpp:77
void Preconditioner_Smoother_free(Preconditioner_Smoother *in)
Definition: Smoother.cpp:47
bool Jacobi
Definition: Preconditioner.h:70
dim_t type
Definition: Preconditioner.h:50
void Preconditioner_LocalSmoother_Sweep_tiled(SparseMatrix_ptr A, Preconditioner_LocalSmoother *gs, double *x)
#define PASO_RILU
Definition: Options.h:63
double * b_F
Definition: Preconditioner.h:138
boost::shared_ptr< const Preconditioner > const_Preconditioner_ptr
Definition: Preconditioner.h:41
#define PASO_ILU0
Definition: Options.h:44
double * b_C
Definition: Preconditioner.h:140
boost::shared_ptr< SystemMatrix > SystemMatrix_ptr
Definition: SystemMatrix.h:53
ILU preconditioner.
Definition: Preconditioner.h:117
Solver_RILU * rilu
RILU preconditioner.
Definition: Preconditioner.h:59
Preconditioner_Smoother * jacobi
Jacobi preconditioner.
Definition: Preconditioner.h:53
Preconditioner_Smoother * gs
Gauss-Seidel preconditioner.
Definition: Preconditioner.h:55
#define PASO_NO_PRECONDITIONER
Definition: Options.h:65
void Solver_RILU_free(Solver_RILU *in)
Definition: RILU.cpp:49
double * diag
Definition: Preconditioner.h:71
Preconditioner_LocalSmoother * localSmoother
Definition: Preconditioner.h:78
dim_t n_F
Definition: Preconditioner.h:127
void Solver_updateIncompleteSchurComplement(SparseMatrix_ptr A_CC, SparseMatrix_ptr A_CF, double *invA_FF, index_t *A_FF_pivot, SparseMatrix_ptr A_FC)
Definition: SchurComplement.cpp:41
void Preconditioner_LocalSmoother_Sweep_colored(SparseMatrix_ptr A, Preconditioner_LocalSmoother *gs, double *x)
Definition: Smoother.cpp:351
void copy(dim_t N, double *out, const double *in)
out = in
Definition: PasoUtil.h:113
void Preconditioner_LocalSmoother_free(Preconditioner_LocalSmoother *in)
Definition: Smoother.cpp:55
void Preconditioner_Smoother_solve(SystemMatrix_ptr A, Preconditioner_Smoother *gs, double *x, const double *b, dim_t sweeps, bool x_is_initial)
Definition: Smoother.cpp:113
SparseMatrix_ptr A_FC
Definition: Preconditioner.h:131
SparseMatrix_ptr A_CF
Definition: Preconditioner.h:132
index_t dim_t
Definition: DataTypes.h:90
Solver_RILU * RILU_of_Schur
Definition: Preconditioner.h:141
RILU preconditioner.
Definition: Preconditioner.h:123
index_t * rows_in_C
Definition: Preconditioner.h:134
Definition: Preconditioner.h:76
SolverResult
Definition: Paso.h:55
index_t * mask_F
Definition: Preconditioner.h:135
#define PASO_GS
Definition: Options.h:62
bool is_local
Definition: Preconditioner.h:79
void Solver_solveRILU(Solver_RILU *rilu, double *x, double *b)
Definition: RILU.cpp:303
boost::shared_ptr< SparseMatrix > SparseMatrix_ptr
Definition: SparseMatrix.h:48
void Solver_ILU_free(Solver_ILU *in)
Definition: ILU.cpp:48
double * inv_A_FF
Definition: Preconditioner.h:129
Definition: Preconditioner.h:48
dim_t sweeps
Definition: Preconditioner.h:51
#define PASO_JACOBI
Definition: Options.h:46
int index_t
type for array/matrix indices used both globally and on each rank
Definition: DataTypes.h:85
double * x_C
Definition: Preconditioner.h:139
Preconditioner_Smoother * Preconditioner_Smoother_alloc(SystemMatrix_ptr A, bool jacobi, bool is_local, bool verbose)
constructs the symmetric Gauss-Seidel preconditioner
Definition: Smoother.cpp:67
double * buffer
Definition: Preconditioner.h:72
SolverResult Preconditioner_Smoother_solve_byTolerance(SystemMatrix_ptr A, Preconditioner_Smoother *gs, double *x, const double *b, double atol, dim_t *sweeps, bool x_is_initial)
Definition: Smoother.cpp:140
void Preconditioner_LocalSmoother_Sweep_sequential(SparseMatrix_ptr A, Preconditioner_LocalSmoother *gs, double *x)
inplace Gauss-Seidel sweep in sequential mode
Definition: Smoother.cpp:234
index_t * pivot
Definition: Preconditioner.h:73
double * x_F
Definition: Preconditioner.h:137
Preconditioner * Preconditioner_alloc(SystemMatrix_ptr A, Options *options)
Definition: Preconditioner.cpp:61
index_t * A_FF_pivot
Definition: Preconditioner.h:130
void Preconditioner_LocalSmoother_solve(SparseMatrix_ptr A, Preconditioner_LocalSmoother *gs, double *x, const double *b, dim_t sweeps, bool x_is_initial)
Definition: Smoother.cpp:174
double * factors
Definition: Preconditioner.h:119
Solver_RILU * Solver_getRILU(SparseMatrix_ptr A, bool verbose)
Definition: RILU.cpp:85
Solver_ILU * ilu
ILU preconditioner.
Definition: Preconditioner.h:57
Definition: BiCGStab.cpp:26
void Preconditioner_free(Preconditioner *in)
Definition: Preconditioner.cpp:50
void Preconditioner_solve(Preconditioner *prec, SystemMatrix_ptr A, double *x, double *b)
Definition: Preconditioner.cpp:126
void Solver_solveILU(SparseMatrix_ptr A, Solver_ILU *ilu, double *x, const double *b)
Definition: ILU.cpp:329
index_t * rows_in_F
Definition: Preconditioner.h:133
Definition: Preconditioner.h:68
void Preconditioner_LocalSmoother_Sweep(SparseMatrix_ptr A, Preconditioner_LocalSmoother *gs, double *x)
Definition: Smoother.cpp:214
boost::shared_ptr< Preconditioner > Preconditioner_ptr
Definition: Preconditioner.h:39