Chombo + EB
3.2
|
#include <BackwardEuler.H>
Public Member Functions | |
~BackwardEuler () | |
BackwardEuler (const RefCountedPtr< AMRMultiGrid< LevelData< FArrayBox > > > &a_solver, const AMRLevelOpFactory< LevelData< FArrayBox > > &a_factory, const ProblemDomain &a_level0Domain, const Vector< int > &a_refRat, int a_numLevels=-1, int a_verbosity=3) | |
void | oneStep (Vector< LevelData< FArrayBox > *> &a_phiNew, Vector< LevelData< FArrayBox > *> &a_phiOld, Vector< LevelData< FArrayBox > *> &a_source, const Real &a_dt, int a_lbase, int a_lmax, bool a_zeroPhi=true) |
void | resetAlphaAndBeta (const Real &a_alpha, const Real &a_beta) |
void | computeDiffusion (Vector< LevelData< FArrayBox > * > &a_diffusiveTerm, Vector< LevelData< FArrayBox > * > &a_phiOld, Vector< LevelData< FArrayBox > * > &a_src, Real a_oldTime, Real a_dt, int a_lbase, int a_lmax, bool a_zeroPhi) |
Protected Member Functions | |
void | solveHelm (Vector< LevelData< FArrayBox > * > &a_ans, Vector< LevelData< FArrayBox > * > &a_rhs, int a_lbase, int a_lmax, Real a_dt, bool a_zeroPhi) |
void | createEulerRHS (Vector< LevelData< FArrayBox > * > &a_ans, Vector< LevelData< FArrayBox > * > &a_source, Vector< LevelData< FArrayBox > * > &a_phiOld, int a_lbase, int a_lmax, Real a_dt) |
void | createData (Vector< LevelData< FArrayBox > * > &a_source, int a_lbase, int a_lmax) |
TGAHelmOp< LevelData< FArrayBox > > * | newOp (const ProblemDomain &a_indexSpace, const AMRLevelOpFactory< LevelData< FArrayBox > > &a_opFact) |
Private Member Functions | |
BackwardEuler (const BackwardEuler &a_opin) | |
void | operator= (const BackwardEuler &a_opin) |
BackwardEuler () | |
weak construction is bad. Ref Counted pointers are your friends. More... | |
Private Attributes | |
Vector< TGAHelmOp< LevelData< FArrayBox > > *> | m_ops |
Vector< LevelData< FArrayBox > *> | m_rhst |
ProblemDomain | m_level0Domain |
Vector< int > | m_refRat |
RefCountedPtr< AMRMultiGrid< LevelData< FArrayBox > > > | m_solver |
int | m_verbosity |
int | m_numLevels |
bool | m_dataCreated |
Solves dphi/dt = L phi + rho using backward euler. Uses TGAHelmOp functionaltity.
BackwardEuler::~BackwardEuler | ( | ) |
BackwardEuler::BackwardEuler | ( | const RefCountedPtr< AMRMultiGrid< LevelData< FArrayBox > > > & | a_solver, |
const AMRLevelOpFactory< LevelData< FArrayBox > > & | a_factory, | ||
const ProblemDomain & | a_level0Domain, | ||
const Vector< int > & | a_refRat, | ||
int | a_numLevels = -1 , |
||
int | a_verbosity = 3 |
||
) |
|
inlineprivate |
References MayDay::Error().
|
inlineprivate |
weak construction is bad. Ref Counted pointers are your friends.
References MayDay::Error().
void BackwardEuler::oneStep | ( | Vector< LevelData< FArrayBox > *> & | a_phiNew, |
Vector< LevelData< FArrayBox > *> & | a_phiOld, | ||
Vector< LevelData< FArrayBox > *> & | a_source, | ||
const Real & | a_dt, | ||
int | a_lbase, | ||
int | a_lmax, | ||
bool | a_zeroPhi = true |
||
) |
This advances a parabolic pde from a_phiOld to a_phiNew using TGA on a non-moving domain with source term a_source
void BackwardEuler::computeDiffusion | ( | Vector< LevelData< FArrayBox > * > & | a_diffusiveTerm, |
Vector< LevelData< FArrayBox > * > & | a_phiOld, | ||
Vector< LevelData< FArrayBox > * > & | a_src, | ||
Real | a_oldTime, | ||
Real | a_dt, | ||
int | a_lbase, | ||
int | a_lmax, | ||
bool | a_zeroPhi | ||
) |
|
protected |
|
protected |
|
protected |
|
protected |
|
inlineprivate |
References MayDay::Error().
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |