ML
Version of the Day
|
#include <ml_RefMaxwell.h>
Public Member Functions | |
ML_RMP (const Epetra_CrsMatrix &SM_Matrix, const Teuchos::ParameterList &List, const bool ComputePrec=true) | |
ML_RMP (const Epetra_CrsMatrix &SM_Matrix, const Epetra_CrsMatrix &D0_Clean_Matrix, const Epetra_CrsMatrix &Ms_Matrix, const Epetra_CrsMatrix &M0inv_Matrix, const Epetra_CrsMatrix &M1_Matrix, const Teuchos::ParameterList &List, const bool ComputePrec=true) | |
Constructs a RefMaxwellPreconditioner. | |
Destructor | |
~ML_RMP () | |
Destructor. | |
int | Apply (const Epetra_MultiVector &, Epetra_MultiVector &) const |
Apply the inverse of the preconditioner to an Epetra_MultiVector (NOT AVAILABLE) | |
int | ApplyInverse (const Epetra_MultiVector &B, Epetra_MultiVector &X) const |
Apply the preconditioner w/ RHS B and get result X. | |
int | ComputePreconditioner (const bool CheckFiltering=false) |
Computes the multilevel hierarchy. More... | |
int | ReComputePreconditioner () |
Recomputes the preconditioner. | |
void | Print (int whichHierarchy=11) |
Print the individual operators in the multigrid hierarchy. | |
int | DestroyPreconditioner () |
Destroys all structures allocated in ComputePreconditioner() if the preconditioner has been computed. | |
int | SetUseTranspose (bool) |
Sets use transpose (not implemented). | |
double | NormInf () const |
Returns the infinity norm (not implemented). | |
bool | UseTranspose () const |
Returns the current UseTranspose setting. | |
bool | HasNormInf () const |
Returns true if the this object can provide an approximate Inf-norm, false otherwise. | |
const Epetra_Comm & | Comm () const |
Returns a pointer to the Epetra_Comm communicator associated with this operator. | |
const Epetra_Map & | OperatorDomainMap () const |
Returns the Epetra_Map object associated with the domain of this operator. | |
const Epetra_Map & | OperatorRangeMap () const |
Returns the Epetra_Map object associated with the range of this operator. | |
void | Complexities (double &complexity, double &fineNnz) |
Return operator complexity and #nonzeros in fine grid matrix. | |
The preconditioner(s) for the Eddy Current Maxwell's equations using compatible discretizations (edge elements). Since the preconditioners involve both an edge and nodal component, different combinations of additive and/or multiplicative preconditioners can be used. These options are controlled with the "refmaxwell: mode" option in the Techos::ParameterList. The sublists "refmaxwell: 11list" and "refmaxwell: 22list" will be passed on to the edge and nodal sub-solvers appropriately.
Detail on this preconditioner can be found in Bochev, Hu, Siefert and Tuminaro, 2007.
int ML_Epetra::ML_RMP::ComputePreconditioner | ( | const bool | CheckFiltering = false | ) |
Computes the multilevel hierarchy.
Computes the multilevel hierarchy. This function retrives the user's defines parameters (as specified in the input ParameterList), or takes default values otherwise, and creates the ML objects for aggregation and hierarchy. Allocated data can be freed used DestroyPreconditioner(), or by the destructor,
In a Newton-type procedure, several linear systems have to be solved, Often, these systems are not too different. In this case, it might be convenient to keep the already computed preconditioner (with hierarchy, coarse solver, smoothers), and use it to precondition the next linear system. ML offers a way to determine whether the already available preconditioner is "good enough" for the next linear system. The user should proceed as follows:
"reuse: enable"
== true
ComputePreconditioner(true)
It is supposed that the pointer to the Epetra_RowMatrix remains constant. Currently, it is not possible to modify this pointer (other than creating a new preconditioner)Computes the preconditioner