Public Types | Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
oomph::TrilinosAztecOOSolver Class Reference

An interface to the Trilinos AztecOO classes allowing it to be used as an Oomph-lib LinearSolver. The AztecOO solver is a Krylov Subspace solver; the solver type (either CG, GMRES or BiCGStab) can be set using solver_type(). This solver can be preconditioned with Trilinos Preconditioners (derived from TrilinosPreconditionerBase) or Oomph-lib preconditioners (derived from Preconditioner). Preconditioners are set using preconditioner_pt(). More...

#include <trilinos_solver.h>

+ Inheritance diagram for oomph::TrilinosAztecOOSolver:

Public Types

enum  AztecOO_solver_types { CG , GMRES , BiCGStab }
 if this solver is using an oomph-lib preconditioner then the vectors passed to preconditioner_solve(...) should be using the values in the epetra vector as external data. If the vectors are using external data then rebuild(...) methods cannot be used be used in the preconditioner. More...
 

Public Member Functions

 TrilinosAztecOOSolver ()
 Constructor. More...
 
 ~TrilinosAztecOOSolver ()
 Destructor - delete the solver and the matrices. More...
 
 TrilinosAztecOOSolver (const TrilinosAztecOOSolver &)=delete
 Broken copy constructor. More...
 
void operator= (const TrilinosAztecOOSolver &)=delete
 Broken assignment operator. More...
 
void enable_aztecoo_workaround_for_epetra_matrix_setup ()
 Enable workaround for creating of epetra matrix that respects aztecoo's ordering requirements. More...
 
void disable_aztecoo_workaround_for_epetra_matrix_setup ()
 Disable workaround for creating of epetra matrix that respects aztecoo's ordering requirements. More...
 
bool is_aztecoo_workaround_for_epetra_matrix_setup_enabled ()
 Is workaround for creating of epetra matrix that respects aztecoo's ordering requirements enabled? More...
 
void clean_up_memory ()
 Clean up method - deletes the solver, the matrices and the preconditioner. More...
 
void solve (Problem *const &problem_pt, DoubleVector &solution)
 Function which uses problem_pt's get_jacobian(...) function to generate a linear system which is then solved. This function deletes any existing internal data and then generates a new AztecOO solver. More...
 
void solve (DoubleMatrixBase *const &matrix_pt, const DoubleVector &rhs, DoubleVector &solution)
 Function to solve the linear system defined by matrix_pt and rhs. NOTE 1. The matrix has to be of type CRDoubleMatrix or DistributedCRDoubleMatrix. NOTE 2. This function will delete any existing internal data and generate a new AztecOO solver. More...
 
void resolve (const DoubleVector &rhs, DoubleVector &solution)
 Function to resolve a linear system using the existing solver data, allowing a solve with a new right hand side vector. This function must be used after a call to solve(...) with enable_resolve set to true. More...
 
void disable_resolve ()
 Disable resolve function (overloads the LinearSolver disable_resolve function). More...
 
void enable_delete_matrix ()
 Call if the matrix can be deleted. More...
 
void disable_delete_matrix ()
 Call if the matrix can not be deleted (default) More...
 
unsigned & max_iter ()
 Access function to Max_iter. More...
 
unsigned iterations () const
 Acess function to Iterations. More...
 
double & tolerance ()
 Access function to Tolerance. More...
 
unsigned & solver_type ()
 Access function to Solver_type. More...
 
double jacobian_setup_time ()
 Function to return Jacobian_setup_time;. More...
 
double linear_solver_solution_time ()
 Function to return Linear_solver_solution_time. More...
 
void enable_assemble_serial_jacobian ()
 Set the assembly of the serial jacobian when performing a problem-based solve. More...
 
void disable_assemble_serial_jacobian ()
 Unset the assembly of the serial jacobian. More...
 
- Public Member Functions inherited from oomph::IterativeLinearSolver
 IterativeLinearSolver ()
 Constructor: Set (default) trivial preconditioner and set defaults for tolerance and max. number of iterations. More...
 
 IterativeLinearSolver (const IterativeLinearSolver &)=delete
 Broken copy constructor. More...
 
void operator= (const IterativeLinearSolver &)=delete
 Broken assignment operator. More...
 
virtual ~IterativeLinearSolver ()
 Destructor (empty) More...
 
Preconditioner *& preconditioner_pt ()
 Access function to preconditioner. More...
 
Preconditioner *const & preconditioner_pt () const
 Access function to preconditioner (const version) More...
 
double & tolerance ()
 Access to convergence tolerance. More...
 
unsigned & max_iter ()
 Access to max. number of iterations. More...
 
void enable_doc_convergence_history ()
 Enable documentation of the convergence history. More...
 
void disable_doc_convergence_history ()
 Disable documentation of the convergence history. More...
 
void open_convergence_history_file_stream (const std::string &file_name, const std::string &zone_title="")
 Write convergence history into file with specified filename (automatically switches on doc). Optional second argument is a string that can be used (as a zone title) to identify what case we're running (e.g. what combination of linear solver and preconditioner or parameter values are used). More...
 
void close_convergence_history_file_stream ()
 Close convergence history output stream. More...
 
double jacobian_setup_time () const
 returns the time taken to assemble the jacobian matrix and residual vector More...
 
double linear_solver_solution_time () const
 return the time taken to solve the linear system More...
 
virtual double preconditioner_setup_time () const
 returns the the time taken to setup the preconditioner More...
 
void enable_setup_preconditioner_before_solve ()
 Setup the preconditioner before the solve. More...
 
void disable_setup_preconditioner_before_solve ()
 Don't set up the preconditioner before the solve. More...
 
void enable_error_after_max_iter ()
 Throw an error if we don't converge within max_iter. More...
 
void disable_error_after_max_iter ()
 Don't throw an error if we don't converge within max_iter (default). More...
 
void enable_iterative_solver_as_preconditioner ()
 Enables the iterative solver be used as preconditioner (when calling the solve method it bypass the setup solver method – currently only used by Trilinos solver —) More...
 
void disable_iterative_solver_as_preconditioner ()
 Disables the iterative solver be used as preconditioner (when calling the solve method it bypass the setup solver method – currently only used by Trilinos solver —) More...
 
- Public Member Functions inherited from oomph::LinearSolver
 LinearSolver ()
 Empty constructor, initialise the member data. More...
 
 LinearSolver (const LinearSolver &dummy)=delete
 Broken copy constructor. More...
 
void operator= (const LinearSolver &)=delete
 Broken assignment operator. More...
 
virtual ~LinearSolver ()
 Empty virtual destructor. More...
 
void enable_doc_time ()
 Enable documentation of solve times. More...
 
void disable_doc_time ()
 Disable documentation of solve times. More...
 
bool is_doc_time_enabled () const
 Is documentation of solve times enabled? More...
 
bool is_resolve_enabled () const
 Boolean flag indicating if resolves are enabled. More...
 
virtual void enable_resolve ()
 Enable resolve (i.e. store matrix and/or LU decomposition, say) Virtual so it can be overloaded to perform additional tasks. More...
 
virtual void solve (DoubleMatrixBase *const &matrix_pt, const Vector< double > &rhs, Vector< double > &result)
 Linear-algebra-type solver: Takes pointer to a matrix and rhs vector and returns the solution of the linear system. More...
 
virtual void solve_transpose (Problem *const &problem_pt, DoubleVector &result)
 Solver: Takes pointer to problem and returns the results vector which contains the solution of the linear system defined by the problem's fully assembled Jacobian and residual vector (broken virtual). More...
 
virtual void solve_transpose (DoubleMatrixBase *const &matrix_pt, const DoubleVector &rhs, DoubleVector &result)
 Linear-algebra-type solver: Takes pointer to a matrix and rhs vector and returns the solution of the linear system. More...
 
virtual void solve_transpose (DoubleMatrixBase *const &matrix_pt, const Vector< double > &rhs, Vector< double > &result)
 Linear-algebra-type solver: Takes pointer to a matrix and rhs vector and returns the solution of the linear system. More...
 
virtual void resolve_transpose (const DoubleVector &rhs, DoubleVector &result)
 Solver: Resolve the system defined by the last assembled jacobian and the rhs vector. Solution is returned in the vector result. (broken virtual) More...
 
virtual void enable_computation_of_gradient ()
 function to enable the computation of the gradient required for the globally convergent Newton method More...
 
void disable_computation_of_gradient ()
 function to disable the computation of the gradient required for the globally convergent Newton method More...
 
void reset_gradient ()
 function to reset the size of the gradient before each Newton solve More...
 
void get_gradient (DoubleVector &gradient)
 function to access the gradient, provided it has been computed More...
 
- Public Member Functions inherited from oomph::DistributableLinearAlgebraObject
 DistributableLinearAlgebraObject ()
 Default constructor - create a distribution. More...
 
 DistributableLinearAlgebraObject (const DistributableLinearAlgebraObject &matrix)=delete
 Broken copy constructor. More...
 
void operator= (const DistributableLinearAlgebraObject &)=delete
 Broken assignment operator. More...
 
virtual ~DistributableLinearAlgebraObject ()
 Destructor. More...
 
LinearAlgebraDistributiondistribution_pt () const
 access to the LinearAlgebraDistribution More...
 
unsigned nrow () const
 access function to the number of global rows. More...
 
unsigned nrow_local () const
 access function for the num of local rows on this processor. More...
 
unsigned nrow_local (const unsigned &p) const
 access function for the num of local rows on this processor. More...
 
unsigned first_row () const
 access function for the first row on this processor More...
 
unsigned first_row (const unsigned &p) const
 access function for the first row on this processor More...
 
bool distributed () const
 distribution is serial or distributed More...
 
bool distribution_built () const
 if the communicator_pt is null then the distribution is not setup then false is returned, otherwise return true More...
 
void build_distribution (const LinearAlgebraDistribution *const dist_pt)
 setup the distribution of this distributable linear algebra object More...
 
void build_distribution (const LinearAlgebraDistribution &dist)
 setup the distribution of this distributable linear algebra object More...
 

Protected Member Functions

void solve_using_AztecOO (Epetra_Vector *&rhs_pt, Epetra_Vector *&soln_pt)
 Helper function performs the actual solve once the AztecOO solver is set up. More...
 
void solver_setup (DoubleMatrixBase *const &matrix_pt)
 Helper function for setting up the solver. Converts the oomph-lib matrices to Epetra matrices, sets up the preconditioner, creates the Trilinos Aztec00 solver and passes in the matrix, preconditioner and parameters. More...
 
- Protected Member Functions inherited from oomph::DistributableLinearAlgebraObject
void clear_distribution ()
 clear the distribution of this distributable linear algebra object More...
 

Protected Attributes

bool Use_aztecoo_workaround_for_epetra_matrix_setup
 Use workaround for creating of epetra matrix that respects aztecoo's ordering requirements. More...
 
unsigned Iterations
 Stores number of iterations used. More...
 
AztecOO * AztecOO_solver_pt
 Pointer to the AztecOO solver. More...
 
double Jacobian_setup_time
 Stores set up time for Jacobian. More...
 
double Linear_solver_solution_time
 Stores time for the solution (excludes time to set up preconditioner) More...
 
bool Delete_matrix
 Trilinos copies matrix data from oomph-lib's own CRDoubleMatrix or DistributedCRDoubleMatrix to Trilinos's Epetra format - the Trilinos solver no longer requires the oomph-lib matrices and therefore they could be deleted to save memory. This must be requested explicitly by setting this flag to TRUE. NOTE: The matrix is deleted after the preconditioner is setup. More...
 
bool Assemble_serial_jacobian
 If true, when performing a problem based solve a serial matrix will be requested from Problem::get_jacobian(...). Defaults to true. More...
 
unsigned Solver_type
 Defines which solver is set up - available types are defined in AztecOO_solver_types. More...
 
bool Using_problem_based_solve
 Helper flag keeping track of whether we called the linear algebra or problem-based solve function. More...
 
Epetra_CrsMatrix * Epetra_matrix_pt
 A pointer for the linear system matrix in Epetra_CrsMatrix format. More...
 
Epetra_Operator * Epetra_preconditioner_pt
 A pointer to the Epetra_Operator for the preconditioner. This is only used if the preconditioner NOT a Trilinos preconditioner. More...
 
DoubleMatrixBaseOomph_matrix_pt
 Oomph lib matrix pointer. More...
 
ProblemProblem_pt
 A pointer to the underlying problem (NULL if MATRIX based solve) The problem_pt is stored here in a problem based solve for the preconditioner. More...
 
bool If_oomphlib_preconditioner_use_epetra_values
 if this solver is using an oomph-lib preconditioner then the vectors passed to preconditioner_solve(...) should be using the values in the epetra vector as external data. If the vectors are using external data then rebuild(...) methods cannot be used. More...
 
- Protected Attributes inherited from oomph::IterativeLinearSolver
bool Doc_convergence_history
 Flag indicating if the convergence history is to be documented. More...
 
std::ofstream Output_file_stream
 Output file stream for convergence history. More...
 
double Tolerance
 Convergence tolerance. More...
 
unsigned Max_iter
 Maximum number of iterations. More...
 
PreconditionerPreconditioner_pt
 Pointer to the preconditioner. More...
 
double Jacobian_setup_time
 Jacobian setup time. More...
 
double Solution_time
 linear solver solution time More...
 
double Preconditioner_setup_time
 Preconditioner setup time. More...
 
bool Setup_preconditioner_before_solve
 indicates whether the preconditioner should be setup before solve. Default = true; More...
 
bool Throw_error_after_max_iter
 Should we throw an error instead of just returning when we hit the max iterations? More...
 
bool Use_iterative_solver_as_preconditioner
 Use the iterative solver as preconditioner. More...
 
bool First_time_solve_when_used_as_preconditioner
 When the iterative solver is used a preconditioner then we call the setup of solver method only once (the first time the solve method is called) More...
 
- Protected Attributes inherited from oomph::LinearSolver
bool Enable_resolve
 Boolean that indicates whether the matrix (or its factors, in the case of direct solver) should be stored so that the resolve function can be used. More...
 
bool Doc_time
 Boolean flag that indicates whether the time taken. More...
 
bool Compute_gradient
 flag that indicates whether the gradient required for the globally convergent Newton method should be computed or not More...
 
bool Gradient_has_been_computed
 flag that indicates whether the gradient was computed or not More...
 
DoubleVector Gradient_for_glob_conv_newton_solve
 DoubleVector storing the gradient for the globally convergent Newton method. More...
 

Additional Inherited Members

- Static Protected Attributes inherited from oomph::IterativeLinearSolver
static IdentityPreconditioner Default_preconditioner
 Default preconditioner: The base class for preconditioners is a fully functional (if trivial!) preconditioner. More...
 

Detailed Description

An interface to the Trilinos AztecOO classes allowing it to be used as an Oomph-lib LinearSolver. The AztecOO solver is a Krylov Subspace solver; the solver type (either CG, GMRES or BiCGStab) can be set using solver_type(). This solver can be preconditioned with Trilinos Preconditioners (derived from TrilinosPreconditionerBase) or Oomph-lib preconditioners (derived from Preconditioner). Preconditioners are set using preconditioner_pt().

Definition at line 266 of file trilinos_solver.h.

Member Enumeration Documentation

◆ AztecOO_solver_types

if this solver is using an oomph-lib preconditioner then the vectors passed to preconditioner_solve(...) should be using the values in the epetra vector as external data. If the vectors are using external data then rebuild(...) methods cannot be used be used in the preconditioner.

Enumerated list to define which AztecOO solver is used

Enumerator
CG 
GMRES 
BiCGStab 

Definition at line 483 of file trilinos_solver.h.

Constructor & Destructor Documentation

◆ TrilinosAztecOOSolver() [1/2]

oomph::TrilinosAztecOOSolver::TrilinosAztecOOSolver ( )
inline

◆ ~TrilinosAztecOOSolver()

oomph::TrilinosAztecOOSolver::~TrilinosAztecOOSolver ( )
inline

Destructor - delete the solver and the matrices.

Definition at line 304 of file trilinos_solver.h.

References clean_up_memory(), Oomph_matrix_pt, and Using_problem_based_solve.

◆ TrilinosAztecOOSolver() [2/2]

oomph::TrilinosAztecOOSolver::TrilinosAztecOOSolver ( const TrilinosAztecOOSolver )
delete

Broken copy constructor.

Member Function Documentation

◆ clean_up_memory()

void oomph::TrilinosAztecOOSolver::clean_up_memory ( )
inlinevirtual

Clean up method - deletes the solver, the matrices and the preconditioner.

Reimplemented from oomph::LinearSolver.

Definition at line 347 of file trilinos_solver.h.

References AztecOO_solver_pt, oomph::Preconditioner::clean_up_memory(), Epetra_matrix_pt, Epetra_preconditioner_pt, and oomph::IterativeLinearSolver::preconditioner_pt().

Referenced by disable_resolve(), solve(), solver_setup(), and ~TrilinosAztecOOSolver().

◆ disable_assemble_serial_jacobian()

void oomph::TrilinosAztecOOSolver::disable_assemble_serial_jacobian ( )
inline

Unset the assembly of the serial jacobian.

Definition at line 467 of file trilinos_solver.h.

References Assemble_serial_jacobian.

◆ disable_aztecoo_workaround_for_epetra_matrix_setup()

void oomph::TrilinosAztecOOSolver::disable_aztecoo_workaround_for_epetra_matrix_setup ( )
inline

Disable workaround for creating of epetra matrix that respects aztecoo's ordering requirements.

Definition at line 333 of file trilinos_solver.h.

References Use_aztecoo_workaround_for_epetra_matrix_setup.

◆ disable_delete_matrix()

void oomph::TrilinosAztecOOSolver::disable_delete_matrix ( )
inline

Call if the matrix can not be deleted (default)

Definition at line 418 of file trilinos_solver.h.

References Delete_matrix.

◆ disable_resolve()

void oomph::TrilinosAztecOOSolver::disable_resolve ( )
inlinevirtual

Disable resolve function (overloads the LinearSolver disable_resolve function).

Reimplemented from oomph::LinearSolver.

Definition at line 405 of file trilinos_solver.h.

References clean_up_memory(), and oomph::LinearSolver::Enable_resolve.

◆ enable_assemble_serial_jacobian()

void oomph::TrilinosAztecOOSolver::enable_assemble_serial_jacobian ( )
inline

Set the assembly of the serial jacobian when performing a problem-based solve.

Definition at line 461 of file trilinos_solver.h.

References Assemble_serial_jacobian.

◆ enable_aztecoo_workaround_for_epetra_matrix_setup()

void oomph::TrilinosAztecOOSolver::enable_aztecoo_workaround_for_epetra_matrix_setup ( )
inline

Enable workaround for creating of epetra matrix that respects aztecoo's ordering requirements.

Definition at line 326 of file trilinos_solver.h.

References Use_aztecoo_workaround_for_epetra_matrix_setup.

◆ enable_delete_matrix()

void oomph::TrilinosAztecOOSolver::enable_delete_matrix ( )
inline

Call if the matrix can be deleted.

Definition at line 412 of file trilinos_solver.h.

References Delete_matrix.

◆ is_aztecoo_workaround_for_epetra_matrix_setup_enabled()

bool oomph::TrilinosAztecOOSolver::is_aztecoo_workaround_for_epetra_matrix_setup_enabled ( )
inline

Is workaround for creating of epetra matrix that respects aztecoo's ordering requirements enabled?

Definition at line 340 of file trilinos_solver.h.

References Use_aztecoo_workaround_for_epetra_matrix_setup.

◆ iterations()

unsigned oomph::TrilinosAztecOOSolver::iterations ( ) const
inlinevirtual

Acess function to Iterations.

Implements oomph::IterativeLinearSolver.

Definition at line 430 of file trilinos_solver.h.

References Iterations.

◆ jacobian_setup_time()

double oomph::TrilinosAztecOOSolver::jacobian_setup_time ( )
inline

Function to return Jacobian_setup_time;.

Definition at line 448 of file trilinos_solver.h.

References Jacobian_setup_time.

◆ linear_solver_solution_time()

double oomph::TrilinosAztecOOSolver::linear_solver_solution_time ( )
inline

Function to return Linear_solver_solution_time.

Definition at line 454 of file trilinos_solver.h.

References Linear_solver_solution_time.

◆ max_iter()

unsigned& oomph::TrilinosAztecOOSolver::max_iter ( )
inline

Access function to Max_iter.

Definition at line 424 of file trilinos_solver.h.

References oomph::IterativeLinearSolver::Max_iter.

◆ operator=()

void oomph::TrilinosAztecOOSolver::operator= ( const TrilinosAztecOOSolver )
delete

Broken assignment operator.

◆ resolve()

void oomph::TrilinosAztecOOSolver::resolve ( const DoubleVector rhs,
DoubleVector solution 
)
virtual

◆ solve() [1/2]

void oomph::TrilinosAztecOOSolver::solve ( DoubleMatrixBase *const &  matrix_pt,
const DoubleVector rhs,
DoubleVector solution 
)
virtual

◆ solve() [2/2]

void oomph::TrilinosAztecOOSolver::solve ( Problem *const &  problem_pt,
DoubleVector solution 
)
virtual

Function which uses problem_pt's get_jacobian(...) function to generate a linear system which is then solved. This function deletes any existing internal data and then generates a new AztecOO solver.

////////////////////////////////////////////////////////////////////////////

//////////////////////////////////////////////////////////////////////////// Function which uses problem_pt's get_jacobian(...) function to generate a linear system which is then solved. This function deletes any existing internal data and then generates a new AztecOO solver.

Implements oomph::LinearSolver.

Definition at line 43 of file trilinos_solver.cc.

References oomph::DoubleVector::build(), oomph::DistributableLinearAlgebraObject::build_distribution(), oomph::DoubleVector::built(), clean_up_memory(), oomph::DistributableLinearAlgebraObject::distribution_pt(), oomph::LinearSolver::Doc_time, oomph::Problem::get_jacobian(), Jacobian_setup_time, oomph::oomph_info, Oomph_matrix_pt, Problem_pt, oomph::DoubleVector::redistribute(), oomph::TimingHelpers::timer(), and Using_problem_based_solve.

◆ solve_using_AztecOO()

void oomph::TrilinosAztecOOSolver::solve_using_AztecOO ( Epetra_Vector *&  rhs_pt,
Epetra_Vector *&  soln_pt 
)
protected

Helper function performs the actual solve once the AztecOO solver is set up.

Helper function performs the actual solve once the AztecOO solver is set up (i.e. solver_setup() is called)

Definition at line 570 of file trilinos_solver.cc.

References AztecOO_solver_pt, oomph::LinearSolver::Doc_time, Iterations, oomph::IterativeLinearSolver::Max_iter, oomph::oomph_info, and oomph::IterativeLinearSolver::Tolerance.

Referenced by resolve(), and solve().

◆ solver_setup()

void oomph::TrilinosAztecOOSolver::solver_setup ( DoubleMatrixBase *const &  matrix_pt)
protected

◆ solver_type()

unsigned& oomph::TrilinosAztecOOSolver::solver_type ( )
inline

Access function to Solver_type.

Definition at line 442 of file trilinos_solver.h.

References Solver_type.

◆ tolerance()

double& oomph::TrilinosAztecOOSolver::tolerance ( )
inline

Access function to Tolerance.

Definition at line 436 of file trilinos_solver.h.

References oomph::IterativeLinearSolver::Tolerance.

Member Data Documentation

◆ Assemble_serial_jacobian

bool oomph::TrilinosAztecOOSolver::Assemble_serial_jacobian
protected

If true, when performing a problem based solve a serial matrix will be requested from Problem::get_jacobian(...). Defaults to true.

Definition at line 527 of file trilinos_solver.h.

Referenced by disable_assemble_serial_jacobian(), enable_assemble_serial_jacobian(), and TrilinosAztecOOSolver().

◆ AztecOO_solver_pt

AztecOO* oomph::TrilinosAztecOOSolver::AztecOO_solver_pt
protected

Pointer to the AztecOO solver.

Definition at line 509 of file trilinos_solver.h.

Referenced by clean_up_memory(), solve_using_AztecOO(), solver_setup(), and TrilinosAztecOOSolver().

◆ Delete_matrix

bool oomph::TrilinosAztecOOSolver::Delete_matrix
protected

Trilinos copies matrix data from oomph-lib's own CRDoubleMatrix or DistributedCRDoubleMatrix to Trilinos's Epetra format - the Trilinos solver no longer requires the oomph-lib matrices and therefore they could be deleted to save memory. This must be requested explicitly by setting this flag to TRUE. NOTE: The matrix is deleted after the preconditioner is setup.

Definition at line 523 of file trilinos_solver.h.

Referenced by disable_delete_matrix(), enable_delete_matrix(), solver_setup(), and TrilinosAztecOOSolver().

◆ Epetra_matrix_pt

Epetra_CrsMatrix* oomph::TrilinosAztecOOSolver::Epetra_matrix_pt
protected

A pointer for the linear system matrix in Epetra_CrsMatrix format.

Definition at line 538 of file trilinos_solver.h.

Referenced by clean_up_memory(), resolve(), solver_setup(), and TrilinosAztecOOSolver().

◆ Epetra_preconditioner_pt

Epetra_Operator* oomph::TrilinosAztecOOSolver::Epetra_preconditioner_pt
protected

A pointer to the Epetra_Operator for the preconditioner. This is only used if the preconditioner NOT a Trilinos preconditioner.

Definition at line 542 of file trilinos_solver.h.

Referenced by clean_up_memory(), solver_setup(), and TrilinosAztecOOSolver().

◆ If_oomphlib_preconditioner_use_epetra_values

bool oomph::TrilinosAztecOOSolver::If_oomphlib_preconditioner_use_epetra_values
protected

if this solver is using an oomph-lib preconditioner then the vectors passed to preconditioner_solve(...) should be using the values in the epetra vector as external data. If the vectors are using external data then rebuild(...) methods cannot be used.

Definition at line 556 of file trilinos_solver.h.

Referenced by TrilinosAztecOOSolver().

◆ Iterations

unsigned oomph::TrilinosAztecOOSolver::Iterations
protected

Stores number of iterations used.

Definition at line 506 of file trilinos_solver.h.

Referenced by iterations(), and solve_using_AztecOO().

◆ Jacobian_setup_time

double oomph::TrilinosAztecOOSolver::Jacobian_setup_time
protected

Stores set up time for Jacobian.

Definition at line 512 of file trilinos_solver.h.

Referenced by jacobian_setup_time(), and solve().

◆ Linear_solver_solution_time

double oomph::TrilinosAztecOOSolver::Linear_solver_solution_time
protected

Stores time for the solution (excludes time to set up preconditioner)

Definition at line 515 of file trilinos_solver.h.

Referenced by linear_solver_solution_time(), resolve(), and solve().

◆ Oomph_matrix_pt

DoubleMatrixBase* oomph::TrilinosAztecOOSolver::Oomph_matrix_pt
protected

Oomph lib matrix pointer.

Definition at line 545 of file trilinos_solver.h.

Referenced by solve(), solver_setup(), and ~TrilinosAztecOOSolver().

◆ Problem_pt

Problem* oomph::TrilinosAztecOOSolver::Problem_pt
protected

A pointer to the underlying problem (NULL if MATRIX based solve) The problem_pt is stored here in a problem based solve for the preconditioner.

Definition at line 550 of file trilinos_solver.h.

Referenced by solve(), and TrilinosAztecOOSolver().

◆ Solver_type

unsigned oomph::TrilinosAztecOOSolver::Solver_type
protected

Defines which solver is set up - available types are defined in AztecOO_solver_types.

Definition at line 531 of file trilinos_solver.h.

Referenced by solver_setup(), solver_type(), and TrilinosAztecOOSolver().

◆ Use_aztecoo_workaround_for_epetra_matrix_setup

bool oomph::TrilinosAztecOOSolver::Use_aztecoo_workaround_for_epetra_matrix_setup
protected

Use workaround for creating of epetra matrix that respects aztecoo's ordering requirements.

Definition at line 503 of file trilinos_solver.h.

Referenced by disable_aztecoo_workaround_for_epetra_matrix_setup(), enable_aztecoo_workaround_for_epetra_matrix_setup(), is_aztecoo_workaround_for_epetra_matrix_setup_enabled(), solver_setup(), and TrilinosAztecOOSolver().

◆ Using_problem_based_solve

bool oomph::TrilinosAztecOOSolver::Using_problem_based_solve
protected

Helper flag keeping track of whether we called the linear algebra or problem-based solve function.

Definition at line 535 of file trilinos_solver.h.

Referenced by solve(), solver_setup(), TrilinosAztecOOSolver(), and ~TrilinosAztecOOSolver().


The documentation for this class was generated from the following files: