|
| | CellConvectiveOperator (std::string object_name, SAMRAI::tbox::Pointer< SAMRAI::pdat::CellVariable< NDIM, double > > Q_cell_var, int Q_min_ghost_cell_width, SAMRAI::tbox::Pointer< SAMRAI::tbox::Database > input_db, ConvectiveDifferencingType difference_form, std::vector< SAMRAI::solv::RobinBcCoefStrategy< NDIM > * > bc_coefs) |
| | Class constructor. More...
|
| |
| | ~CellConvectiveOperator ()=default |
| | Default destructor. More...
|
| |
| virtual void | interpolateToFaceOnHierarchy (int q_interp_idx, int Q_cell_idx, int u_idx, bool synch_cf_bdry=true) |
| | Interpolate a cell-centered field Q to a face-centered field q, possibly using the provided advection velocity field. More...
|
| |
| virtual void | evaluateAdvectiveFluxOnHierarchy (int q_flux_idx, int Q_cell_idx, int u_idx, const bool synch_cf_bdry=true) |
| | Evaluate the face-centered flux Q to a face-centered field q using the provided advection velocity field. More...
|
| |
| virtual void | computeAdvectiveDerivativeOnHierarchy (int N_cell_idx, int q_interp_idx, int u_idx, bool synch_cf_bdry=true) |
| | Compute the advective derivative N = u * grad Q on the patch hierarchy. More...
|
| |
| virtual void | computeConservativeDerivativeOnHierarchy (int N_cell_idx, int q_flux_idx, bool synch_cf_bdry=true) |
| | Compute the conservative derivative N = div(Q u) on the patch hierarchy. More...
|
| |
| virtual void | computeSkewSymmetricDerivativeOnHierarchy (int N_cell_idx, int q_interp_idx, int q_flux_idx, int u_idx, bool synch_cf_bdry=true) |
| | Compute the skew-symmetric derivative N = 0.5[u * grad Q + div(Q u)] on the patch hierarchy. More...
|
| |
| virtual void | interpolateToFaceOnPatch (SAMRAI::pdat::FaceData< NDIM, double > &q_interp_data, const SAMRAI::pdat::CellData< NDIM, double > &Q_cell_data, const SAMRAI::pdat::FaceData< NDIM, double > &u_data, const SAMRAI::hier::Patch< NDIM > &patch)=0 |
| | Interpolate a cell-centered field Q to a face-centered field q on a single grid patch. More...
|
| |
| virtual void | evaluateAdvectiveFluxOnPatch (SAMRAI::pdat::FaceData< NDIM, double > &q_flux_data, const SAMRAI::pdat::CellData< NDIM, double > &Q_cell_data, const SAMRAI::pdat::FaceData< NDIM, double > &u_data, const SAMRAI::hier::Patch< NDIM > &patch) |
| | Evaluate the face-centered flux Q to a face-centered field q using the provided advection velocity field. More...
|
| |
| void | applyConvectiveOperator (int Q_idx, int N_idx) override |
| | Compute N = u * grad Q or N = div(Q u). More...
|
| |
| void | setAdvectionVelocity (int u_idx) |
| | Set the patch data index corresponding to the advection velocity to be used when computing the convective derivative. More...
|
| |
| int | getAdvectionVelocity () const |
| | Get the patch data index corresponding to the advection velocity used when computing the convective derivative. More...
|
| |
| void | setConvectiveDifferencingType (ConvectiveDifferencingType difference_form) |
| | Set the convective differencing form to be used by the operator. More...
|
| |
| ConvectiveDifferencingType | getConvectiveDifferencingType () const |
| | Get the convective differencing form used by the operator. More...
|
| |
The vector arguments for apply(), applyAdd(), etc, need not match those for initializeOperatorState(). However, there must be a certain degree of similarity, including
- hierarchy configuration (hierarchy pointer and level range)
- number, type and alignment of vector component data
- ghost cell widths of data in the input and output vectors
- Note
- It is generally necessary to reinitialize the operator state when the hierarchy configuration changes.
It is safe to call initializeOperatorState() when the state is already initialized. In this case, the operator state is first deallocated and then reinitialized.
Conditions on arguments:
- input and output vectors must have same hierarchy
- input and output vectors must have same structure, depth, etc.
Call deallocateOperatorState() to remove any data allocated by this method.
- See also
- deallocateOperatorState
- Parameters
-
| in | input vector |
| out | output vector |
Reimplemented from IBTK::GeneralOperator.
Before calling apply(), the form of the vectors x and y should be set properly by the user on all patch interiors on the specified range of levels in the patch hierarchy. The user is responsible for all data management for the quantities associated with the vectors. In particular, patch data in these vectors must be allocated prior to calling this method.
- Parameters
-
| x | input vector |
| y | output vector, i.e., \(y=F[x]\) |
Conditions on Parameters:
- vectors x and y must have same hierarchy
- vectors x and y must have same structure, depth, etc.
In general, the vectors x and y cannot be the same.
- See also
- initializeOperatorState
Implements IBTK::GeneralOperator.
Before calling applyAdd(), the form of the vectors x, y, and z should be set properly by the user on all patch interiors on the specified range of levels in the patch hierarchy. The user is responsible for all data management for the quantities associated with the vectors. In particular, patch data in these vectors must be allocated prior to calling this method.
- Parameters
-
| x | input vector |
| y | input vector |
| z | output vector, i.e., \(z=F[x]+y\) |
Conditions on Parameters:
- vectors x, y, and z must have same hierarchy
- vectors x, y, and z must have same structure, depth, etc.
In general, the vectors x and y cannot be the same.
- Note
- Subclasses may require that the operator be initialized prior to calling applyAdd().
- See also
- initializeOperatorState
- Note
- A default implementation is provided which employs apply() and SAMRAI::solv::SAMRAIVectorReal::add().
Reimplemented in IBTK::PETScSNESJacobianJOWrapper, and IBTK::PETScMatLOWrapper.