Public Types | Public Member Functions | Static Public Attributes | Private Member Functions | Private Attributes | List of all members
oomph::SuperLUSolver Class Reference

. SuperLU Project Solver class. This is a combined wrapper for both SuperLU and SuperLU Dist. See http://crd.lbl.gov/~xiaoye/SuperLU/ Default Behaviour: If this solver is distributed over more than one processor then SuperLU Dist is used. Member data naming convention: member data associated with the SuperLU Dist solver begins Dist_... and member data associated with the serial SuperLU solver begins Serial_... . More...

#include <linear_solver.h>

+ Inheritance diagram for oomph::SuperLUSolver:

Public Types

enum  Type { Default, Serial, Distributed }
 enum type to specify the solver behaviour. Default - will employ superlu dist if more than 1 processor. Serial - will always try use superlu (serial). Distributed - will always try to use superlu dist. More...
 

Public Member Functions

 SuperLUSolver ()
 Constructor. Set the defaults. More...
 
 SuperLUSolver (const SuperLUSolver &dummy)
 Broken copy constructor. More...
 
void operator= (const SuperLUSolver &)
 Broken assignment operator. More...
 
 ~SuperLUSolver ()
 Destructor, clean up the stored matrices. More...
 
void enable_computation_of_gradient ()
 function to enable the computation of the gradient More...
 
void disable_resolve ()
 Overload disable resolve so that it cleans up memory too. More...
 
void solve (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. More...
 
void solve (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. The function returns the global result Vector. Note: if Delete_matrix_data is true the function matrix_pt->clean_up_memory() will be used to wipe the matrix data. More...
 
void resolve (const DoubleVector &rhs, DoubleVector &result)
 Resolve the system defined by the last assembled jacobian and the specified rhs vector if resolve has been enabled. Note: returns the global result Vector. More...
 
void enable_doc_stats ()
 Enable documentation of solver statistics. More...
 
void disable_doc_stats ()
 Disable documentation of solver statistics. More...
 
double jacobian_setup_time () const
 returns the time taken to assemble the jacobian matrix and residual vector More...
 
virtual double linear_solver_solution_time () const
 return the time taken to solve the linear system (needs to be overloaded for each linear solver) More...
 
void factorise (DoubleMatrixBase *const &matrix_pt)
 Do the factorisation stage Note: if Delete_matrix_data is true the function matrix_pt->clean_up_memory() will be used to wipe the matrix data. More...
 
void backsub (const DoubleVector &rhs, DoubleVector &result)
 Do the backsubstitution for SuperLU solver Note: returns the global result Vector. More...
 
void clean_up_memory ()
 Clean up the memory allocated by the solver. More...
 
void set_solver_type (const Type &t)
 Specify the solve type. Either default, serial or distributed. See enum SuperLU_solver_type for more details. More...
 
void use_compressed_row_for_superlu_serial ()
 Use the compressed row format in superlu serial. More...
 
void use_compressed_column_for_superlu_serial ()
 Use the compressed column format in superlu serial. More...
 
void enable_delete_matrix_data_in_superlu_dist ()
 Set Delete_matrix_data flag. SuperLU_dist needs its own copy of the input matrix, therefore a copy must be made if any matrix used with this solver is to be preserved. If the input matrix can be deleted the flag can be set to true to reduce the amount of memory required, and the matrix data will be wiped using its clean_up_memory() function. Default value is false. More...
 
void disable_delete_matrix_data_in_superlu_dist ()
 Unset Delete_matrix_data flag. SuperLU_dist needs its own copy of the input matrix, therefore a copy must be made if any matrix used with this solver is to be preserved. If the input matrix can be. More...
 
void enable_row_and_col_permutations_in_superlu_dist ()
 Set flag so that SuperLU_DIST is allowed to permute matrix rows and columns during factorisation. This is the default for SuperLU_DIST, and can lead to significantly faster solves. More...
 
void disable_row_and_col_permutations_in_superlu_dist ()
 Set flag so that SuperLU_DIST is not allowed to permute matrix rows and columns during factorisation. More...
 
void use_global_solve_in_superlu_dist ()
 Calling this method will ensure that when the problem based solve interface is used, a global (serial) jacobian will be assembled. Note: calling this function will delete any distributed solve data. More...
 
void use_distributed_solve_in_superlu_dist ()
 Calling this method will ensure that when the problem based solve interface is used, a distributed jacobian will be assembled. Note: calling this function will delete any global solve data. More...
 
- Public Member Functions inherited from oomph::LinearSolver
 LinearSolver ()
 Empty constructor, initialise the member data. More...
 
 LinearSolver (const LinearSolver &dummy)
 Broken copy constructor. More...
 
void operator= (const LinearSolver &)
 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...
 
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)
 Broken copy constructor. More...
 
void operator= (const DistributableLinearAlgebraObject &)
 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
 
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...
 

Static Public Attributes

static bool Suppress_incorrect_rhs_distribution_warning_in_resolve =false
 Static flag that determines whether the warning about incorrect distribution of RHSs will be printed or not. More...
 

Private Member Functions

void factorise_serial (DoubleMatrixBase *const &matrix_pt)
 factorise method for SuperLU (serial) More...
 
void backsub_serial (const DoubleVector &rhs, DoubleVector &result)
 backsub method for SuperLU (serial) More...
 
void factorise_distributed (DoubleMatrixBase *const &matrix_pt)
 factorise method for SuperLU Dist More...
 
void backsub_distributed (const DoubleVector &rhs, DoubleVector &result)
 backsub method for SuperLU Dist More...
 

Private Attributes

double Jacobian_setup_time
 Jacobian setup time. More...
 
double Solution_time
 Solution time. More...
 
bool Suppress_solve
 Suppress solve? More...
 
bool Doc_stats
 Set to true to output statistics (false by default). More...
 
Type Solver_type
 the solver type. see SuperLU_solver_type for details. More...
 
bool Using_dist
 boolean flag indicating whether superlu dist is being used More...
 
void * Serial_f_factors
 Storage for the LU factors as required by SuperLU. More...
 
int Serial_info
 Info flag for the SuperLU solver. More...
 
unsigned long Serial_n_dof
 The number of unknowns in the linear system. More...
 
int Serial_sign_of_determinant_of_matrix
 Sign of the determinant of the matrix. More...
 
bool Serial_compressed_row_flag
 Use compressed row version? More...
 
bool Dist_use_global_solver
 Flag that determines whether the MPIProblem based solve function uses the global or distributed version of SuperLU_DIST (default value is false). More...
 
bool Dist_global_solve_data_allocated
 Flag is true if solve data has been generated for a global matrix. More...
 
bool Dist_distributed_solve_data_allocated
 Flag is true if solve data has been generated for distributed matrix. More...
 
void * Dist_solver_data_pt
 Storage for the LU factors and other data required by SuperLU. More...
 
int Dist_nprow
 Number of rows for the process grid. More...
 
int Dist_npcol
 Number of columns for the process grid. More...
 
int Dist_info
 Info flag for the SuperLU solver. More...
 
bool Dist_allow_row_and_col_permutations
 If true then SuperLU_DIST is allowed to permute matrix rows and columns during factorisation. This is the default for SuperLU_DIST, and can lead to significantly faster solves, but has been known to fail, hence the default value is 0. More...
 
bool Dist_delete_matrix_data
 Delete_matrix_data flag. SuperLU_dist needs its own copy of the input matrix, therefore a copy must be made if any matrix used with this solver is to be preserved. If the input matrix can be deleted the flag can be set to true to reduce the amount of memory required, and the matrix data will be wiped using its clean_up_memory() function. Default value is false. More...
 
double * Dist_value_pt
 Pointer for storage of the matrix values required by SuperLU_DIST. More...
 
int * Dist_index_pt
 Pointer for storage of matrix rows or column indices required by SuperLU_DIST. More...
 
int * Dist_start_pt
 Pointers for storage of matrix column or row starts. More...
 

Additional Inherited Members

- Protected Member Functions inherited from oomph::DistributableLinearAlgebraObject
void clear_distribution ()
 clear the distribution of this distributable linear algebra object 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...
 

Detailed Description

. SuperLU Project Solver class. This is a combined wrapper for both SuperLU and SuperLU Dist. See http://crd.lbl.gov/~xiaoye/SuperLU/ Default Behaviour: If this solver is distributed over more than one processor then SuperLU Dist is used. Member data naming convention: member data associated with the SuperLU Dist solver begins Dist_... and member data associated with the serial SuperLU solver begins Serial_... .

Definition at line 424 of file linear_solver.h.

Member Enumeration Documentation

enum type to specify the solver behaviour. Default - will employ superlu dist if more than 1 processor. Serial - will always try use superlu (serial). Distributed - will always try to use superlu dist.

Enumerator
Default 
Serial 
Distributed 

Definition at line 433 of file linear_solver.h.

Constructor & Destructor Documentation

oomph::SuperLUSolver::SuperLUSolver ( )
inline
oomph::SuperLUSolver::SuperLUSolver ( const SuperLUSolver dummy)
inline

Broken copy constructor.

Definition at line 465 of file linear_solver.h.

References oomph::BrokenCopy::broken_copy().

oomph::SuperLUSolver::~SuperLUSolver ( )
inline

Destructor, clean up the stored matrices.

Definition at line 477 of file linear_solver.h.

References clean_up_memory().

Member Function Documentation

void oomph::SuperLUSolver::backsub ( const DoubleVector rhs,
DoubleVector result 
)

Do the backsubstitution for SuperLU solver Note: returns the global result Vector.

Do the backsubstitution for SuperLUSolver. Note - this method performs no paranoid checks - these are all performed in solve(...) and resolve(...)

Definition at line 1539 of file linear_solver.cc.

References backsub_distributed(), backsub_serial(), and Using_dist.

Referenced by resolve(), and solve().

void oomph::SuperLUSolver::backsub_distributed ( const DoubleVector rhs,
DoubleVector result 
)
private
void oomph::SuperLUSolver::backsub_serial ( const DoubleVector rhs,
DoubleVector result 
)
private
void oomph::SuperLUSolver::clean_up_memory ( )
virtual
void oomph::SuperLUSolver::disable_delete_matrix_data_in_superlu_dist ( )
inline

Unset Delete_matrix_data flag. SuperLU_dist needs its own copy of the input matrix, therefore a copy must be made if any matrix used with this solver is to be preserved. If the input matrix can be.

Definition at line 597 of file linear_solver.h.

References Dist_delete_matrix_data.

void oomph::SuperLUSolver::disable_doc_stats ( )
inline

Disable documentation of solver statistics.

Definition at line 532 of file linear_solver.h.

References Doc_stats.

Referenced by oomph::SuperLUPreconditioner::SuperLUPreconditioner().

void oomph::SuperLUSolver::disable_resolve ( )
inlinevirtual

Overload disable resolve so that it cleans up memory too.

Reimplemented from oomph::LinearSolver.

Definition at line 492 of file linear_solver.h.

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

void oomph::SuperLUSolver::disable_row_and_col_permutations_in_superlu_dist ( )
inline

Set flag so that SuperLU_DIST is not allowed to permute matrix rows and columns during factorisation.

Definition at line 608 of file linear_solver.h.

References Dist_allow_row_and_col_permutations.

void oomph::SuperLUSolver::enable_computation_of_gradient ( )
inlinevirtual

function to enable the computation of the gradient

Reimplemented from oomph::LinearSolver.

Definition at line 483 of file linear_solver.h.

References oomph::LinearSolver::Compute_gradient.

void oomph::SuperLUSolver::enable_delete_matrix_data_in_superlu_dist ( )
inline

Set Delete_matrix_data flag. SuperLU_dist needs its own copy of the input matrix, therefore a copy must be made if any matrix used with this solver is to be preserved. If the input matrix can be deleted the flag can be set to true to reduce the amount of memory required, and the matrix data will be wiped using its clean_up_memory() function. Default value is false.

Definition at line 591 of file linear_solver.h.

References Dist_delete_matrix_data.

void oomph::SuperLUSolver::enable_doc_stats ( )
inline

Enable documentation of solver statistics.

Definition at line 529 of file linear_solver.h.

References Doc_stats.

void oomph::SuperLUSolver::enable_row_and_col_permutations_in_superlu_dist ( )
inline

Set flag so that SuperLU_DIST is allowed to permute matrix rows and columns during factorisation. This is the default for SuperLU_DIST, and can lead to significantly faster solves.

Definition at line 603 of file linear_solver.h.

References Dist_allow_row_and_col_permutations.

void oomph::SuperLUSolver::factorise ( DoubleMatrixBase *const &  matrix_pt)

Do the factorisation stage Note: if Delete_matrix_data is true the function matrix_pt->clean_up_memory() will be used to wipe the matrix data.

LU decompose the matrix addressed by matrix_pt by using the SuperLU solver. The resulting matrix factors are stored internally.

Definition at line 1104 of file linear_solver.cc.

References clean_up_memory(), oomph::LinearAlgebraDistribution::communicator_pt(), Default, Distributed, oomph::DistributableLinearAlgebraObject::distribution_pt(), factorise_distributed(), factorise_serial(), oomph::MPI_Helpers::mpi_has_been_initialised(), Solver_type, and Using_dist.

Referenced by oomph::SuperLUPreconditioner::setup(), and solve().

void oomph::SuperLUSolver::factorise_distributed ( DoubleMatrixBase *const &  matrix_pt)
private

factorise method for SuperLU Dist

LU decompose the matrix addressed by matrix_pt using the SuperLU_DIST solver. The resulting matrix factors are stored internally.

Definition at line 1151 of file linear_solver.cc.

References oomph::DistributableLinearAlgebraObject::build_distribution(), oomph::CRDoubleMatrix::built(), clean_up_memory(), oomph::CCMatrix< T >::clean_up_memory(), oomph::CRDoubleMatrix::clear(), oomph::CRDoubleMatrix::column_index(), oomph::CCMatrix< T >::column_start(), oomph::LinearAlgebraDistribution::communicator_pt(), oomph::MPI_Helpers::communicator_pt(), Dist_allow_row_and_col_permutations, Dist_delete_matrix_data, Dist_distributed_solve_data_allocated, Dist_global_solve_data_allocated, Dist_index_pt, Dist_info, Dist_npcol, Dist_nprow, Dist_solver_data_pt, Dist_start_pt, Dist_value_pt, oomph::DistributableLinearAlgebraObject::distributed(), oomph::DistributableLinearAlgebraObject::distribution_pt(), oomph::LeakCheckNames::doc(), Doc_stats, oomph::DistributableLinearAlgebraObject::first_row(), i, oomph::DoubleMatrixBase::ncol(), oomph::SparseMatrix< T, MATRIX_TYPE >::nnz(), oomph::CRDoubleMatrix::nnz(), oomph::LinearAlgebraDistribution::nrow(), oomph::DoubleMatrixBase::nrow(), oomph::DistributableLinearAlgebraObject::nrow(), oomph::CRDoubleMatrix::nrow(), oomph::CCDoubleMatrix::nrow(), oomph::DistributableLinearAlgebraObject::nrow_local(), oomph::CCMatrix< T >::row_index(), oomph::CRDoubleMatrix::row_start(), oomph::superlu_cr_to_cc(), oomph::superlu_dist_distributed_matrix(), oomph::superlu_dist_global_matrix(), oomph::SparseMatrix< T, MATRIX_TYPE >::value(), and oomph::CRDoubleMatrix::value().

Referenced by factorise().

void oomph::SuperLUSolver::factorise_serial ( DoubleMatrixBase *const &  matrix_pt)
private
double oomph::SuperLUSolver::jacobian_setup_time ( ) const
inlinevirtual

returns the time taken to assemble the jacobian matrix and residual vector

Reimplemented from oomph::LinearSolver.

Definition at line 536 of file linear_solver.h.

References Jacobian_setup_time.

virtual double oomph::SuperLUSolver::linear_solver_solution_time ( ) const
inlinevirtual

return the time taken to solve the linear system (needs to be overloaded for each linear solver)

Reimplemented from oomph::LinearSolver.

Definition at line 543 of file linear_solver.h.

References Solution_time.

void oomph::SuperLUSolver::operator= ( const SuperLUSolver )
inline

Broken assignment operator.

Definition at line 471 of file linear_solver.h.

References oomph::BrokenCopy::broken_assign().

void oomph::SuperLUSolver::resolve ( const DoubleVector rhs,
DoubleVector result 
)
virtual

Resolve the system defined by the last assembled jacobian and the specified rhs vector if resolve has been enabled. Note: returns the global result Vector.

Resolve the system for a given RHS.

Reimplemented from oomph::LinearSolver.

Definition at line 1080 of file linear_solver.cc.

References backsub(), oomph::LinearSolver::Doc_time, oomph::oomph_info, Solution_time, and oomph::TimingHelpers::timer().

Referenced by oomph::SuperLUPreconditioner::preconditioner_solve().

void oomph::SuperLUSolver::set_solver_type ( const Type t)
inline

Specify the solve type. Either default, serial or distributed. See enum SuperLU_solver_type for more details.

Definition at line 563 of file linear_solver.h.

References clean_up_memory(), Solver_type, and t.

void oomph::SuperLUSolver::solve ( Problem *const &  problem_pt,
DoubleVector result 
)
virtual
void oomph::SuperLUSolver::solve ( DoubleMatrixBase *const &  matrix_pt,
const DoubleVector rhs,
DoubleVector result 
)
virtual

Linear-algebra-type solver: Takes pointer to a matrix and rhs vector and returns the solution of the linear system. The function returns the global result Vector. Note: if Delete_matrix_data is true the function matrix_pt->clean_up_memory() will be used to wipe the matrix data.

Linear-algebra-type solver: Takes pointer to a matrix and rhs vector and returns the solution of the linear system. Problem pointer defaults to NULL and can be omitted. The function returns the global result Vector. Note: if Delete_matrix_data is true the function matrix_pt->clean_up_memory() will be used to wipe the matrix data.

Reimplemented from oomph::LinearSolver.

Definition at line 933 of file linear_solver.cc.

References backsub(), oomph::DistributableLinearAlgebraObject::build_distribution(), oomph::DoubleVector::built(), clean_up_memory(), oomph::LinearAlgebraDistribution::distributed(), oomph::DistributableLinearAlgebraObject::distribution_pt(), oomph::LinearSolver::Doc_time, oomph::LinearSolver::Enable_resolve, factorise(), oomph::DoubleMatrixBase::ncol(), oomph::CRDoubleMatrix::nnz(), oomph::DoubleMatrixBase::nrow(), oomph::DistributableLinearAlgebraObject::nrow(), oomph::oomph_info, Solution_time, and oomph::TimingHelpers::timer().

void oomph::SuperLUSolver::use_compressed_column_for_superlu_serial ( )
inline

Use the compressed column format in superlu serial.

Definition at line 577 of file linear_solver.h.

References Serial_compressed_row_flag.

void oomph::SuperLUSolver::use_compressed_row_for_superlu_serial ( )
inline

Use the compressed row format in superlu serial.

Definition at line 573 of file linear_solver.h.

References Serial_compressed_row_flag.

void oomph::SuperLUSolver::use_distributed_solve_in_superlu_dist ( )
inline

Calling this method will ensure that when the problem based solve interface is used, a distributed jacobian will be assembled. Note: calling this function will delete any global solve data.

Definition at line 628 of file linear_solver.h.

References clean_up_memory(), and Dist_use_global_solver.

void oomph::SuperLUSolver::use_global_solve_in_superlu_dist ( )
inline

Calling this method will ensure that when the problem based solve interface is used, a global (serial) jacobian will be assembled. Note: calling this function will delete any distributed solve data.

Definition at line 615 of file linear_solver.h.

References clean_up_memory(), and Dist_use_global_solver.

Member Data Documentation

bool oomph::SuperLUSolver::Dist_allow_row_and_col_permutations
private

If true then SuperLU_DIST is allowed to permute matrix rows and columns during factorisation. This is the default for SuperLU_DIST, and can lead to significantly faster solves, but has been known to fail, hence the default value is 0.

Definition at line 733 of file linear_solver.h.

Referenced by disable_row_and_col_permutations_in_superlu_dist(), enable_row_and_col_permutations_in_superlu_dist(), factorise_distributed(), and SuperLUSolver().

bool oomph::SuperLUSolver::Dist_delete_matrix_data
private

Delete_matrix_data flag. SuperLU_dist needs its own copy of the input matrix, therefore a copy must be made if any matrix used with this solver is to be preserved. If the input matrix can be deleted the flag can be set to true to reduce the amount of memory required, and the matrix data will be wiped using its clean_up_memory() function. Default value is false.

Definition at line 741 of file linear_solver.h.

Referenced by disable_delete_matrix_data_in_superlu_dist(), enable_delete_matrix_data_in_superlu_dist(), factorise_distributed(), solve(), and SuperLUSolver().

bool oomph::SuperLUSolver::Dist_distributed_solve_data_allocated
private

Flag is true if solve data has been generated for distributed matrix.

Definition at line 715 of file linear_solver.h.

Referenced by backsub_distributed(), clean_up_memory(), factorise_distributed(), and SuperLUSolver().

bool oomph::SuperLUSolver::Dist_global_solve_data_allocated
private

Flag is true if solve data has been generated for a global matrix.

Definition at line 712 of file linear_solver.h.

Referenced by backsub_distributed(), clean_up_memory(), factorise_distributed(), and SuperLUSolver().

int* oomph::SuperLUSolver::Dist_index_pt
private

Pointer for storage of matrix rows or column indices required by SuperLU_DIST.

Definition at line 748 of file linear_solver.h.

Referenced by clean_up_memory(), factorise_distributed(), and SuperLUSolver().

int oomph::SuperLUSolver::Dist_info
private

Info flag for the SuperLU solver.

Definition at line 727 of file linear_solver.h.

Referenced by backsub_distributed(), clean_up_memory(), and factorise_distributed().

int oomph::SuperLUSolver::Dist_npcol
private

Number of columns for the process grid.

Definition at line 724 of file linear_solver.h.

Referenced by backsub_distributed(), clean_up_memory(), and factorise_distributed().

int oomph::SuperLUSolver::Dist_nprow
private

Number of rows for the process grid.

Definition at line 721 of file linear_solver.h.

Referenced by backsub_distributed(), clean_up_memory(), and factorise_distributed().

void* oomph::SuperLUSolver::Dist_solver_data_pt
private

Storage for the LU factors and other data required by SuperLU.

Definition at line 718 of file linear_solver.h.

Referenced by backsub_distributed(), clean_up_memory(), factorise_distributed(), and SuperLUSolver().

int* oomph::SuperLUSolver::Dist_start_pt
private

Pointers for storage of matrix column or row starts.

Definition at line 752 of file linear_solver.h.

Referenced by clean_up_memory(), factorise_distributed(), and SuperLUSolver().

bool oomph::SuperLUSolver::Dist_use_global_solver
private

Flag that determines whether the MPIProblem based solve function uses the global or distributed version of SuperLU_DIST (default value is false).

Definition at line 709 of file linear_solver.h.

Referenced by solve(), SuperLUSolver(), use_distributed_solve_in_superlu_dist(), and use_global_solve_in_superlu_dist().

double* oomph::SuperLUSolver::Dist_value_pt
private

Pointer for storage of the matrix values required by SuperLU_DIST.

Definition at line 744 of file linear_solver.h.

Referenced by clean_up_memory(), factorise_distributed(), and SuperLUSolver().

bool oomph::SuperLUSolver::Doc_stats
private

Set to true to output statistics (false by default).

Definition at line 670 of file linear_solver.h.

Referenced by backsub_distributed(), backsub_serial(), clean_up_memory(), disable_doc_stats(), enable_doc_stats(), factorise_distributed(), factorise_serial(), and SuperLUSolver().

double oomph::SuperLUSolver::Jacobian_setup_time
private

Jacobian setup time.

Definition at line 661 of file linear_solver.h.

Referenced by jacobian_setup_time(), and solve().

bool oomph::SuperLUSolver::Serial_compressed_row_flag
private
void* oomph::SuperLUSolver::Serial_f_factors
private

Storage for the LU factors as required by SuperLU.

Definition at line 682 of file linear_solver.h.

Referenced by backsub_serial(), clean_up_memory(), and factorise_serial().

int oomph::SuperLUSolver::Serial_info
private

Info flag for the SuperLU solver.

Definition at line 685 of file linear_solver.h.

Referenced by backsub_serial(), clean_up_memory(), and factorise_serial().

unsigned long oomph::SuperLUSolver::Serial_n_dof
private

The number of unknowns in the linear system.

Definition at line 688 of file linear_solver.h.

Referenced by backsub_serial(), clean_up_memory(), factorise_serial(), and SuperLUSolver().

int oomph::SuperLUSolver::Serial_sign_of_determinant_of_matrix
private

Sign of the determinant of the matrix.

Definition at line 691 of file linear_solver.h.

Referenced by factorise_serial(), solve(), and SuperLUSolver().

double oomph::SuperLUSolver::Solution_time
private

Solution time.

Definition at line 664 of file linear_solver.h.

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

Type oomph::SuperLUSolver::Solver_type
private

the solver type. see SuperLU_solver_type for details.

Definition at line 673 of file linear_solver.h.

Referenced by factorise(), set_solver_type(), solve(), and SuperLUSolver().

bool oomph::SuperLUSolver::Suppress_incorrect_rhs_distribution_warning_in_resolve =false
static

Static flag that determines whether the warning about incorrect distribution of RHSs will be printed or not.

Static warning to suppress warnings about incorrect distribution of RHS vector. Default is false

Definition at line 703 of file linear_solver.h.

Referenced by backsub_distributed().

bool oomph::SuperLUSolver::Suppress_solve
private

Suppress solve?

Definition at line 667 of file linear_solver.h.

Referenced by solve(), and SuperLUSolver().

bool oomph::SuperLUSolver::Using_dist
private

boolean flag indicating whether superlu dist is being used

Definition at line 676 of file linear_solver.h.

Referenced by backsub(), factorise(), and SuperLUSolver().


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