Chombo + EB  3.0
Public Member Functions | Public Attributes | Protected Member Functions | Protected Attributes | List of all members
PetscSolver< T > Class Template Referenceabstract

#include <PetscSolver.H>

Inheritance diagram for PetscSolver< T >:
Inheritance graph
[legend]

Public Member Functions

 PetscSolver ()
 
virtual ~PetscSolver ()
 
void destroy ()
 
virtual void setHomogeneous (bool a_homogeneous)
 
virtual void define (LinearOp< T > *a_operator, bool a_homogeneous=false)
 
virtual void solve (T &a_phi, const T &a_rhs)
 
int solve_private (T &a_phi, const T &a_rhs)
 
Real computeResidual ()
 
int applyOp (T &a_phi, const T &a_rhs)
 
void setInitialGuessNonzero (bool b=true)
 
virtual int getNNZPerRow () const
 
virtual BaseFab< Real > & getRegFab (T &a_fab, DataIterator &dit)=0
 
virtual const BaseFab< Real > & getRegFab (const T &a_fab, DataIterator &dit) const =0
 
virtual const BaseFab< Real > & getRegFab (const T &a_fab, DataIterator &dit, Box &a_box) const =0
 
virtual void defineData (T &a_data, const T &a_template)=0
 
void setNull (bool n=true)
 
- Public Member Functions inherited from LinearSolver< T >
virtual ~LinearSolver ()
 
virtual void setConvergenceMetrics (Real a_metric, Real a_tolerance)
 Set a convergence metric, along with solver tolerance, if desired. More...
 

Public Attributes

bool m_homogeneous
 
Real m_dx
 

Protected Member Functions

virtual Real addBCdiagValue (const IntVect &a_iv, const IntVect &a_jv, const T &a_rhs, DataIterator dit, const Real coeff=1)
 
int resetOperator ()
 
virtual Real addBCrhsValue (const IntVect &a_iv, const T &a_phi, DataIterator dit, const Real &coeff=1)
 
int create_mat_vec (const T &a_phi)
 
int setup_solver (const T &a_phi)
 

Protected Attributes

bool m_null
 
bool m_nz_init_guess
 
m_gids
 
LevelData< BaseFab< bool > > m_bccode
 
int m_gid0
 

Detailed Description

template<class T>
class PetscSolver< T >

Interface to PETSC linear solvers BC implementation hooks: addBCdiagValue (formMatrix) and addBCrhsValue (solveprivate) m_bccode should be encoded true for cells that are at the domain boundaries at formMatrix time.

Constructor & Destructor Documentation

◆ PetscSolver()

template<class T >
PetscSolver< T >::PetscSolver ( )

◆ ~PetscSolver()

template<class T >
PetscSolver< T >::~PetscSolver ( )
virtual

Member Function Documentation

◆ destroy()

template<class T >
void PetscSolver< T >::destroy ( )

◆ setHomogeneous()

template<class T>
virtual void PetscSolver< T >::setHomogeneous ( bool  a_homogeneous)
inlinevirtual

Set whether the solver uses only homogeneous boundary conditions

Implements LinearSolver< T >.

◆ define()

template<class T>
void PetscSolver< T >::define ( LinearOp< T > *  a_operator,
bool  a_homogeneous = false 
)
virtual

Set Function F(u) = 0 and Jacobian dF(u)/du for nonlinear solves

Implements LinearSolver< T >.

Reimplemented in PetscSolverViscousTensor< T >, and PetscSolverPoisson< T >.

Referenced by PetscSolver< LevelData< FArrayBox > >::setHomogeneous().

◆ solve()

template<class T>
void PetscSolver< T >::solve ( T &  a_phi,
const T &  a_rhs 
)
virtual

◆ solve_private()

template<class T>
int PetscSolver< T >::solve_private ( T &  a_phi,
const T &  a_rhs 
)

◆ computeResidual()

template<class T >
Real PetscSolver< T >::computeResidual ( )

◆ applyOp()

template<class T>
int PetscSolver< T >::applyOp ( T &  a_phi,
const T &  a_rhs 
)

◆ setInitialGuessNonzero()

template<class T>
void PetscSolver< T >::setInitialGuessNonzero ( bool  b = true)
inline

set initial guess non-zero

◆ getNNZPerRow()

template<class T>
virtual int PetscSolver< T >::getNNZPerRow ( ) const
inlinevirtual

get an estimate of the number of nnz/row for matrix allocation

Reimplemented in PetscSolverViscousTensor< T >.

Referenced by PetscSolver< LevelData< FArrayBox > >::create_mat_vec().

◆ addBCdiagValue()

template<class T>
virtual Real PetscSolver< T >::addBCdiagValue ( const IntVect a_iv,
const IntVect a_jv,
const T &  a_rhs,
DataIterator  dit,
const Real  coeff = 1 
)
inlineprotectedvirtual

handling boundary conditions, turn it into a term to be added to the diag term this function coordinates with addBCrhsValue for Dirichlet BC for Neumann BC no RHS modification is required

◆ resetOperator()

template<class T>
int PetscSolver< T >::resetOperator ( )
inlineprotected

◆ addBCrhsValue()

template<class T>
virtual Real PetscSolver< T >::addBCrhsValue ( const IntVect a_iv,
const T &  a_phi,
DataIterator  dit,
const Real coeff = 1 
)
inlineprotectedvirtual

handling boundary conditions, turn it into a term to be added to the RHS this function coordinates with addBCdiagValue for Dirichlet BC, should return a term that is to be added to RHS

Referenced by PetscSolver< LevelData< FArrayBox > >::solve_private().

◆ create_mat_vec()

template<class T>
int PetscSolver< T >::create_mat_vec ( const T &  a_phi)
protected

◆ setup_solver()

template<class T>
int PetscSolver< T >::setup_solver ( const T &  a_phi)
protected

◆ getRegFab() [1/3]

template<class T>
virtual BaseFab<Real>& PetscSolver< T >::getRegFab ( T &  a_fab,
DataIterator dit 
)
pure virtual

◆ getRegFab() [2/3]

template<class T>
virtual const BaseFab<Real>& PetscSolver< T >::getRegFab ( const T &  a_fab,
DataIterator dit 
) const
pure virtual

◆ getRegFab() [3/3]

template<class T>
virtual const BaseFab<Real>& PetscSolver< T >::getRegFab ( const T &  a_fab,
DataIterator dit,
Box a_box 
) const
pure virtual

◆ defineData()

template<class T>
virtual void PetscSolver< T >::defineData ( T &  a_data,
const T &  a_template 
)
pure virtual

◆ setNull()

template<class T >
void PetscSolver< T >::setNull ( bool  n = true)

Member Data Documentation

◆ m_homogeneous

template<class T>
bool PetscSolver< T >::m_homogeneous

public member data: whether or not to use inhomogeneous boundary conditions.

Referenced by PetscSolver< LevelData< FArrayBox > >::define(), and PetscSolver< LevelData< FArrayBox > >::setHomogeneous().

◆ m_dx

template<class T>
Real PetscSolver< T >::m_dx

◆ m_null

template<class T>
bool PetscSolver< T >::m_null
protected

◆ m_nz_init_guess

template<class T>
bool PetscSolver< T >::m_nz_init_guess
protected

◆ m_gids

template<class T>
T PetscSolver< T >::m_gids
protected

◆ m_bccode

template<class T>
LevelData<BaseFab<bool> > PetscSolver< T >::m_bccode
protected

◆ m_gid0

template<class T>
int PetscSolver< T >::m_gid0
protected

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