template<class Physics, class Problem>
CompFlow class
CompFlow used polymorphically with tk::DGPDE.
Contents
The template arguments specify policies and are used to configure the behavior of the class. The policies are:
- Physics - physics configuration, see PDE/CompFlow/Physics.h
Problem - problem configuration, see PDE/CompFlow/Problem.h
Constructors, destructors, conversion operators
- CompFlow() explicit
- Constructor.
Public functions
- auto nprim() const -> std::size_t
- auto nmat() const -> std::size_t
- void numEquationDofs(std::vector<std::size_t>& numEqDof) const
- void IcBoxElems(const tk::Fields& geoElem, std::size_t nielem, std::vector<std::unordered_set<std::size_t>>& inbox) const
- auto nstiffeq() const -> std::size_t
- auto nnonstiffeq() const -> std::size_t
- void setStiffEqIdx(std::vector<std::size_t>& stiffEqIdx) const
- void setNonStiffEqIdx(std::vector<std::size_t>& nonStiffEqIdx) const
-
void initialize(const tk::Fields& L,
const std::vector<std::size_t>& inpoel,
const tk::
UnsMesh:: Coords& coord, const std::vector<std::unordered_set<std::size_t>>& inbox, const std::unordered_map<std::size_t, std::set<std::size_t>>&, tk::Fields& unk, tk:: real t, const std::size_t nielem) const -
void computeDensityConstr(std::size_t,
tk::Fields&,
std::vector<tk::
real>& densityConstr) const - Compute density constraint for a given material.
- void lhs(const tk::Fields& geoElem, tk::Fields& l) const
- void updateInterfaceCells(tk::Fields&, std::size_t, std::vector<std::size_t>&, std::vector<std::size_t>&) const
- void updatePrimitives(const tk::Fields&, const tk::Fields&, const tk::Fields&, tk::Fields&, std::size_t, std::vector<std::size_t>&) const
-
void cleanTraceMaterial(tk::
real, const tk::Fields&, tk::Fields&, tk::Fields&, std::size_t) const -
void reconstruct(tk::
real, const tk::Fields&, const tk::Fields&, const inciter:: FaceData&, const std::map<std::size_t, std::vector<std::size_t>>&, const std::vector<std::size_t>&, const tk:: UnsMesh:: Coords&, tk::Fields&, tk::Fields&, const bool, const std::vector<std::size_t>&) const - Reconstruct second-order solution from first-order using least-squares.
-
void limit(] tk::
real t, ] const bool pref, ] const tk::Fields& geoFace, const tk::Fields& geoElem, const inciter:: FaceData& fd, const std::map<std::size_t, std::vector<std::size_t>>& esup, const std::vector<std::size_t>& inpoel, const tk:: UnsMesh:: Coords& coord, const std::vector<std::size_t>& ndofel, const std::vector<std::size_t>& gid, const std::unordered_map<std::size_t, std::size_t>& bid, const std::vector<std::vector<tk:: real>>& uNodalExtrm, const std::vector<std::vector<tk:: real>>&, const std::vector<std::vector<tk:: real>>& mtInv, tk::Fields& U, tk::Fields&, std::vector<std::size_t>& shockmarker) const -
void CPL(const tk::Fields&,
const tk::Fields&,
const std::vector<std::size_t>&,
const tk::
UnsMesh:: Coords&, tk::Fields&, std::size_t) const -
auto cellAvgDeformGrad(const tk::Fields&,
std::size_t) const -> std::array<std::vector<tk::
real>, 9> -
void resetAdapSol(const inciter::
FaceData& fd, tk::Fields& unk, tk::Fields&, const std::vector<std::size_t>& ndofel) const -
void rhs(tk::
real t, const bool pref, const tk::Fields& geoFace, const tk::Fields& geoElem, const inciter:: FaceData& fd, const std::vector<std::size_t>& inpoel, const std::vector<std::unordered_set<std::size_t>>& boxelems, const tk:: UnsMesh:: Coords& coord, const tk::Fields& U, const tk::Fields& P, const std::vector<std::size_t>& ndofel, const tk:: real dt, tk::Fields& R) const -
void eval_ndof(std::size_t nunk,
const tk::
UnsMesh:: Coords& coord, const std::vector<std::size_t>& inpoel, const inciter:: FaceData& fd, const tk::Fields& unk, ] const tk::Fields& prim, inciter:: ctr:: PrefIndicatorType indicator, std::size_t ndof, std::size_t ndofmax, tk:: real tolref, std::vector<std::size_t>& ndofel) const -
auto dt(const std::array<std::vector<tk::
real>, 3>& coord, const std::vector<std::size_t>& inpoel, const inciter:: FaceData& fd, const tk::Fields& geoFace, const tk::Fields& geoElem, const std::vector<std::size_t>& ndofel, const tk::Fields& U, const tk::Fields&, const std::size_t) const -> tk:: real -
void stiff_rhs(std::size_t,
const tk::Fields&,
const std::vector<std::size_t>&,
const tk::
UnsMesh:: Coords&, const tk::Fields&, const tk::Fields&, const std::vector<std::size_t>&, tk::Fields&) const - Compute stiff terms for a single element, not implemented here.
-
auto velocity(const tk::Fields& U,
const std::array<std::vector<tk::
real>, 3>&, const std::array<std::size_t, 4>& N) const -> std::array<std::array<tk:: real, 4>, 3> -
auto OutVarFn() const -> std::map<std::string, tk::
GetVarFn> - auto analyticFieldNames() const -> std::vector<std::string>
- auto histNames() const -> std::vector<std::string>
-
auto surfOutput(const std::map<int, std::vector<std::size_t>>&,
tk::Fields&) const -> std::vector<std::vector<tk::
real>> - Return surface field output going to file.
-
auto histOutput(const std::vector<HistData>& h,
const std::vector<std::size_t>& inpoel,
const tk::
UnsMesh:: Coords& coord, const tk::Fields& U, const tk::Fields&) const -> std::vector<std::vector<tk:: real>> - auto names() const -> std::vector<std::string>
-
auto analyticSolution(tk::
real xi, tk:: real yi, tk:: real zi, tk:: real t) const -> std::vector<tk:: real> -
auto solution(tk::
real xi, tk:: real yi, tk:: real zi, tk:: real t) const -> std::vector<tk:: real> -
auto sp_totalenergy(std::size_t e,
const tk::Fields& unk) const -> tk::
real
Function documentation
template<class Physics, class Problem>
std::size_t inciter::
| Returns | The number of primitive quantities required to be stored for this PDE system |
|---|
Find the number of primitive quantities required for this PDE system
template<class Physics, class Problem>
std::size_t inciter::
| Returns | The number of materials set up for this PDE system |
|---|
Find the number of materials set up for this PDE system
template<class Physics, class Problem>
void inciter::
| Parameters | |
|---|---|
| numEqDof in/out | Array storing number of Dofs for each PDE equation |
Assign number of DOFs per equation in the PDE system
template<class Physics, class Problem>
void inciter::
| Parameters | |
|---|---|
| geoElem in | Element geometry array |
| nielem in | Number of internal elements |
| inbox in/out | List of nodes at which box user ICs are set for each IC box |
Determine elements that lie inside the user-defined IC box
template<class Physics, class Problem>
std::size_t inciter::
| Returns | number of stiff equations |
|---|
Find how 'stiff equations', which we currently have none for Compflow
template<class Physics, class Problem>
std::size_t inciter::
| Returns | number of non-stiff equations |
|---|
Find how 'nonstiff equations', which we currently don't use for Compflow
template<class Physics, class Problem>
void inciter::
| Parameters | |
|---|---|
| stiffEqIdx out | list |
Locate the stiff equations. Unused for compflow.
template<class Physics, class Problem>
void inciter::
| Parameters | |
|---|---|
| nonStiffEqIdx out | list |
Locate the nonstiff equations. Unused for compflow.
template<class Physics, class Problem>
void inciter::
| Parameters | |
|---|---|
| L in | Block diagonal mass matrix |
| inpoel in | Element-node connectivity |
| coord in | Array of nodal coordinates |
| inbox in | List of elements at which box user ICs are set for each IC box |
| unk in/out | Array of unknowns |
| t in | Physical time |
| nielem in | Number of internal elements |
Initalize the compressible flow equations, prepare for time integration
template<class Physics, class Problem>
void inciter::
Compute density constraint for a given material.
| Parameters | |
|---|---|
| densityConstr out | Density Constraint: rho/(rho0*det(g)) |
template<class Physics, class Problem>
void inciter::
| Parameters | |
|---|---|
| geoElem in | Element geometry array |
| l in/out | Block diagonal mass matrix |
Compute the left hand side block-diagonal mass matrix
template<class Physics, class Problem>
void inciter::
Update the interface cells to first order dofs
This function resets the high-order terms in interface cells, and is currently not used in compflow.
template<class Physics, class Problem>
void inciter::
Update the primitives for this PDE system
This function computes and stores the dofs for primitive quantities, which is currently unused for compflow. But if a limiter requires primitive variables for example, this would be the place to add the computation of the primitive variables.
template<class Physics, class Problem>
void inciter::
Clean up the state of trace materials for this PDE system
This function cleans up the state of materials present in trace quantities in each cell. This is unused for compflow.
template<class Physics, class Problem>
void inciter::
| Parameters | |
|---|---|
| t in | Physical time |
| pref | |
| geoFace in | Face geometry array |
| geoElem in | Element geometry array |
| fd in | Face connectivity and boundary conditions object |
| esup in | Elements surrounding points |
| inpoel in | Element-node connectivity |
| coord in | Array of nodal coordinates |
| ndofel in | Vector of local number of degrees of freedome |
| gid in | Local->global node id map |
| bid in | Local chare-boundary node ids (value) associated to global node ids (key) |
| uNodalExtrm in | Chare-boundary nodal extrema for conservative variables |
| mtInv in | Inverse of Taylor mass matrix |
| U in/out | Solution vector at recent time step |
| shockmarker in/out | Vector of shock-marker values |
Limit second-order solution
template<class Physics, class Problem>
void inciter::
Update the conservative variable solution for this PDE system
This function computes the updated dofs for conservative quantities based on the limited solution and is currently not used in compflow.
template<class Physics, class Problem>
std::array<std::vector<tk::
Return cell-average deformation gradient tensor (no-op for compflow)
This function is a no-op in compflow.
template<class Physics, class Problem>
void inciter::
| Parameters | |
|---|---|
| fd in | Face connectivity and boundary conditions object |
| unk in/out | Solution vector at recent time step |
| ndofel in | Vector of local number of degrees of freedome |
Reset the high order solution for p-adaptive scheme This function reset the high order coefficient for p-adaptive solution polynomials.
template<class Physics, class Problem>
void inciter::
| Parameters | |
|---|---|
| t in | Physical time |
| pref in | Indicator for p-adaptive algorithm |
| geoFace in | Face geometry array |
| geoElem in | Element geometry array |
| fd in | Face connectivity and boundary conditions object |
| inpoel in | Element-node connectivity |
| boxelems in | Mesh node ids within user-defined IC boxes |
| coord in | Array of nodal coordinates |
| U in | Solution vector at recent time step |
| P in | Primitive vector at recent time step |
| ndofel in | Vector of local number of degrees of freedom |
| dt in | Delta time |
| R in/out | Right-hand side vector computed |
Compute right hand side
template<class Physics, class Problem>
void inciter::
| Parameters | |
|---|---|
| nunk in | Number of unknowns |
| coord in | Array of nodal coordinates |
| inpoel in | Element-node connectivity |
| fd in | Face connectivity and boundary conditions object |
| unk in | Array of unknowns |
| prim in | Array of primitive quantities |
| indicator in | p-refinement indicator type |
| ndof in | Number of degrees of freedom in the solution |
| ndofmax in | Max number of degrees of freedom for p-refinement |
| tolref in | Tolerance for p-refinement |
| ndofel in/out | Vector of local number of degrees of freedome |
Evaluate the adaptive indicator and mark the ndof for each element
template<class Physics, class Problem>
tk::
| Parameters | |
|---|---|
| coord in | Mesh node coordinates |
| inpoel in | Mesh element connectivity |
| fd in | Face connectivity and boundary conditions object |
| geoFace in | Face geometry array |
| geoElem in | Element geometry array |
| ndofel in | Vector of local number of degrees of freedom |
| U in | Solution vector at recent time step |
| Returns | Minimum time step size |
Compute the minimum time step size
template<class Physics, class Problem>
std::array<std::array<tk::
| Parameters | |
|---|---|
| U in | Solution vector at recent time step |
| N in | Element node indices |
| Returns | Array of the four values of the velocity field |
Extract the velocity field at cell nodes. Currently unused.
template<class Physics, class Problem>
std::map<std::string, tk::
| Returns | Map that associates user-specified strings to functions that compute relevant quantities to be output to file |
|---|
Return a map that associates user-specified strings to functions
template<class Physics, class Problem>
std::vector<std::string> inciter::
| Returns | Vector of strings labelling analytic fields output in file |
|---|
Return analytic field names to be output to file
template<class Physics, class Problem>
std::vector<std::string> inciter::
| Returns | Vector of strings labeling time history fields output in file |
|---|
Return time history field names to be output to file
template<class Physics, class Problem>
std::vector<std::vector<tk::
| Parameters | |
|---|---|
| h in | History point data |
| inpoel in | Element-node connectivity |
| coord in | Array of nodal coordinates |
| U in | Array of unknowns |
Return time history field output evaluated at time history points
template<class Physics, class Problem>
std::vector<std::string> inciter::
| Returns | Vector of strings labelling integral variables output |
|---|
Return names of integral variables to be output to diagnostics file
template<class Physics, class Problem>
std::vector<tk::
| Parameters | |
|---|---|
| xi in | X-coordinate at which to evaluate the analytic solution |
| yi in | Y-coordinate at which to evaluate the analytic solution |
| zi in | Z-coordinate at which to evaluate the analytic solution |
| t in | Physical time at which to evaluate the analytic solution |
| Returns | Vector of analytic solution at given location and time |
Return analytic solution (if defined by Problem) at xi, yi, zi, t
template<class Physics, class Problem>
std::vector<tk::
| Parameters | |
|---|---|
| xi in | X-coordinate at which to evaluate the analytic solution |
| yi in | Y-coordinate at which to evaluate the analytic solution |
| zi in | Z-coordinate at which to evaluate the analytic solution |
| t in | Physical time at which to evaluate the analytic solution |
| Returns | Vector of analytic solution at given location and time |
Return analytic solution for conserved variables
template<class Physics, class Problem>
tk::
| Parameters | |
|---|---|
| e in | Element id for which total energy is required |
| unk in | Vector of conserved quantities |
| Returns | Cell-averaged specific total energy for given element |
Return cell-averaged specific total energy for an element