PetscSolver< T > Class Template Reference

#include <PetscSolver.H>

Inheritance diagram for PetscSolver< T >:

Inheritance graph
[legend]

List of all members.


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.

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 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

Constructor & Destructor Documentation

template<class T>
PetscSolver< T >::PetscSolver (  )  [inline]

template<class T>
PetscSolver< T >::~PetscSolver (  )  [inline, virtual]


Member Function Documentation

template<class T>
void PetscSolver< T >::destroy (  )  [inline]

template<class T>
virtual void PetscSolver< T >::setHomogeneous ( bool  a_homogeneous  )  [inline, virtual]

Set whether the solver uses only homogeneous boundary conditions

Implements LinearSolver< T >.

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

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

Implements LinearSolver< T >.

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

References CH_assert, LinearOp< T >::dx(), PetscSolver< T >::m_dx, and PetscSolver< T >::m_homogeneous.

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

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

template<class T>
Real PetscSolver< T >::computeResidual (  )  [inline]

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

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

set initial guess non-zero

template<class T>
virtual int PetscSolver< T >::getNNZPerRow (  )  const [inline, virtual]

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

Reimplemented in PetscSolverViscousTensor< T >.

Referenced by PetscSolver< T >::create_mat_vec().

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 
) [inline, protected, virtual]

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

template<class T>
int PetscSolver< T >::resetOperator (  )  [inline, protected]

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

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< T >::solve_private().

template<class T>
int PetscSolver< T >::create_mat_vec ( const T &  a_phi  )  [inline, protected]

template<class T>
int PetscSolver< T >::setup_solver ( const T &  a_phi  )  [inline, protected]

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

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

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

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

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


Member Data Documentation

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

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

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

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

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

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

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

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

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


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

Generated on Sat Jun 24 03:43:08 2017 for Chombo + EB by  doxygen 1.5.5