DG class
DG Charm++ chare array used to advance PDEs in time with DG+RK.
Contents
Public static functions
- static void registerReducers()
- Configure Charm++ reduction types for concatenating BC nodelists.
Constructors, destructors, conversion operators
- DG(CkMigrateMessage* msg) explicit
- Migrate constructor.
Public functions
- auto DG(const CProxy_Discretization& disc, const CProxy_Ghosts& ghostsproxy, const std::map<int, std::vector<std::size_t>>& bface, const std::map<int, std::vector<std::size_t>>&, const std::vector<std::size_t>& triinpoel) -> DG_SDAG_CODE explicit
- Constructor.
- void ResumeFromSync() override
- Return from migration.
- void resizeSolVectors()
- Resize solution vectors after extension due to Ghosts and continue setup.
- void setup()
- Setup: query boundary conditions, output mesh, etc.
-
void box(tk::
real v) - Receive total box IC volume and set conditions in box.
- void evalLB(int nrestart)
- void start()
- Start time stepping.
- void next()
- Continue to next time step.
-
void comlim(int fromch,
const std::vector<std::size_t>& tetid,
const std::vector<std::vector<tk::
real>>& u, const std::vector<std::vector<tk:: real>>& prim, const std::vector<std::size_t>& ndof) - Receive chare-boundary limiter function data from neighboring chares.
-
void comreco(int fromch,
const std::vector<std::size_t>& tetid,
const std::vector<std::vector<tk::
real>>& u, const std::vector<std::vector<tk:: real>>& prim, const std::vector<std::size_t>& ndof) - Receive chare-boundary reconstructed data from neighboring chares.
-
void comsol(int fromch,
std::size_t fromstage,
const std::vector<std::size_t>& tetid,
const std::vector<std::vector<tk::
real>>& u, const std::vector<std::vector<tk:: real>>& prim, const std::vector<std::size_t>& ndof) - Receive chare-boundary ghost data from neighboring chares.
-
void comnodalExtrema(const std::vector<std::size_t>& gid,
const std::vector<std::vector<tk::
real>>& G1, const std::vector<std::vector<tk:: real>>& G2) - Receive contributions to nodal gradients on chare-boundaries.
- void resizeNodalExtremac()
- Initialize the vector of nodal extrema.
-
void evalNodalExtrm(const std::size_t ncomp,
const std::size_t nprim,
const std::size_t ndof_NodalExtrm,
const std::vector<std::size_t>& bndel,
const std::vector<std::size_t>& inpoel,
const tk::
UnsMesh:: Coords& coord, const std::vector<std::size_t>& gid, const std::unordered_map<std::size_t, std::size_t>& bid, const tk::Fields& U, const tk::Fields& P, std::vector<std::vector<tk:: real>>& uNodalExtrm, std::vector<std::vector<tk:: real>>& pNodalExtrm) - Compute the nodal extrema for chare-boundary nodes.
-
void comnodeout(const std::vector<std::size_t>& gid,
const std::vector<std::size_t>& nesup,
const std::vector<std::vector<tk::
real>>& L) - Receive nodal solution (ofor field output) contributions from neighboring chares.
-
void refine(const std::vector<tk::
real>& l2res) - Optionally refine/derefine mesh.
-
void resizePostAMR(const std::vector<std::size_t>&,
const tk::
UnsMesh:: Chunk& chunk, const tk:: UnsMesh:: Coords& coord, const std::unordered_map<std::size_t, tk:: UnsMesh:: Edge>&, const std::unordered_map<std::size_t, std::size_t>& addedTets, const std::set<std::size_t>&, const tk:: NodeCommMap& nodeCommMap, const std::map<int, std::vector<std::size_t>>& bface, const std::map<int, std::vector<std::size_t>>&, const std::vector<std::size_t>& triinpoel) - Receive new mesh from Refiner.
-
void extractFieldOutput(const std::vector<std::size_t>&,
const tk::
UnsMesh:: Chunk& chunk, const tk:: UnsMesh:: Coords& coord, const std::unordered_map<std::size_t, tk:: UnsMesh:: Edge>&, const std::unordered_map<std::size_t, std::size_t>& addedTets, const tk:: NodeCommMap& nodeCommMap, const std::map<int, std::vector<std::size_t>>& bface, const std::map<int, std::vector<std::size_t>>&, const std::vector<std::size_t>& triinpoel, CkCallback c) - Extract field output going to file.
- auto solution() const -> const tk::Fields&
- void lhs()
- Compute left hand side.
- void resized()
- Unused in DG.
-
void solve(tk::
real newdt) - Compute right hand side and solve system.
- void step()
- Evaluate whether to continue with next time step.
Charm++ pack/unpack serializer member functions
- CProxy_Discretization m_disc
- Discretization proxy.
- CProxy_Ghosts m_ghosts
- Distributed Ghosts proxy.
- std::size_t m_ndof_NodalExtrm
- Degree of freedom for nodal extrema vector. When DGP1 is applied, there is one degree of freedom for cell average variable. When DGP2 is applied, the degree of freedom is 4 which refers to cell average and gradients in three directions.
- std::size_t m_nsol
- Counter signaling that we have received all our solution ghost data.
- std::size_t m_ninitsol
- Counter signaling that we have received all our solution ghost data during setup.
- std::size_t m_nlim
- Counter signaling that we have received all our limiter function ghost data.
- std::size_t m_nnod
- Counter signaling that we have received all our node solution contributions.
- std::size_t m_nreco
- Counter signaling that we have received all our reconstructed ghost data.
- std::size_t m_nnodalExtrema
- Counter signaling that we have received all our nodal extrema from ghost chare partitions.
- tk::Fields m_u
- Vector of unknown/solution average over each mesh element.
- tk::Fields m_un
- Vector of unknown at previous time-step.
- tk::Fields m_p
- Vector of primitive quantities over each mesh element.
- tk::Fields m_geoElem
- Element geometry.
- tk::Fields m_lhs
- Left-hand side mass-matrix which is a diagonal matrix.
- tk::Fields m_rhs
- Vector of right-hand side.
- std::vector<std::vector<tk::
real>> m_uNodalExtrm - Vector of nodal extrema for conservative variables.
- std::vector<std::vector<tk::
real>> m_pNodalExtrm - Vector of nodal extrema for primitive variables.
- std::unordered_map<std::size_t, std::vector<tk::
real>> m_uNodalExtrmc - Buffer for vector of nodal extrema for conservative variables.
- std::unordered_map<std::size_t, std::vector<tk::
real>> m_pNodalExtrmc - Buffer for vector of nodal extrema for primitive variables.
- std::size_t m_npoin
- Counter for number of nodes on this chare excluding ghosts.
- ElemDiagnostics m_diag
- Diagnostics object.
- std::size_t m_stage
- Runge-Kutta stage counter.
- std::vector<std::size_t> m_ndof
- Vector of local number of degrees of freedom for each element.
- std::vector<std::size_t> m_numEqDof
- Vector of number of degrees of freedom for each PDE equation/component.
- std::array<std::vector<std::vector<tk::
real>>, 3> m_uc - Solution receive buffers for ghosts only.
- std::array<std::vector<std::vector<tk::
real>>, 3> m_pc - Primitive-variable receive buffers for ghosts only.
- std::array<std::vector<std::size_t>, 3> m_ndofc
- Number of degrees of freedom (for p-adaptive) receive buffers for ghosts only.
- std::size_t m_initial
- 1 if starting time stepping, 0 if during time stepping
- std::vector<std::vector<tk::
real>> m_elemfields - Elem output fields.
- std::vector<std::vector<tk::
real>> m_nodefields - Node output fields.
- std::unordered_map<std::size_t, std::pair<std::vector<tk::
real>, std::size_t>> m_nodefieldsc - Ghosts::
OutMesh m_outmesh - Storage for refined mesh used for field output.
- std::vector<std::unordered_set<std::size_t>> m_boxelems
- Element ids at which box ICs are defined by user (multiple boxes)
- std::vector<std::size_t> m_shockmarker
- Shock detection marker for field output.
- void pup(PUP::er& p) override
- Pack/Unpack serialize member function.
- void operator|(PUP::er& p, DG& i)
- Pack/Unpack serialize operator|.
Function documentation
static void inciter::
Configure Charm++ reduction types for concatenating BC nodelists.
Since this is a [initnode] routine, the runtime system executes the routine exactly once on every logical node early on in the Charm++ init sequence. Must be static as it is called without an object. See also: Section "Initializations at Program Startup" at in the Charm++ manual http:/
DG_SDAG_CODE inciter::
Constructor.
| Parameters | |
|---|---|
| disc in | Discretization proxy |
| ghostsproxy | |
| bface in | Boundary-faces mapped to side set ids |
| triinpoel in | Boundary-face connectivity |
void inciter::
Return from migration.
This is called when load balancing (LB) completes. The presence of this function does not affect whether or not we block on LB.
void inciter::
| Parameters | |
|---|---|
| nrestart in | Number of times restarted |
void inciter::
Receive chare-boundary limiter function data from neighboring chares.
| Parameters | |
|---|---|
| fromch in | Sender chare id |
| tetid in | Ghost tet ids we receive solution data for |
| u in | Limited high-order solution |
| prim in | Limited high-order primitive quantities |
| ndof in | Number of degrees of freedom for chare-boundary elements |
This function receives contributions to the limited solution from fellow chares.
void inciter::
Receive chare-boundary reconstructed data from neighboring chares.
| Parameters | |
|---|---|
| fromch in | Sender chare id |
| tetid in | Ghost tet ids we receive solution data for |
| u in | Reconstructed high-order solution |
| prim in | Limited high-order primitive quantities |
| ndof in | Number of degrees of freedom for chare-boundary elements |
This function receives contributions to the reconstructed solution from fellow chares.
void inciter::
Receive chare-boundary ghost data from neighboring chares.
| Parameters | |
|---|---|
| fromch in | Sender chare id |
| fromstage in | Sender chare time step stage |
| tetid in | Ghost tet ids we receive solution data for |
| u in | Solution ghost data |
| prim in | Primitive variables in ghost cells |
| ndof in | Number of degrees of freedom for chare-boundary elements |
This function receives contributions to the unlimited solution from fellow chares.
void inciter::
Receive contributions to nodal gradients on chare-boundaries.
| Parameters | |
|---|---|
| gid in | Global mesh node IDs at which we receive grad contributions |
| G1 in | Partial contributions of extrema for conservative variables to chare-boundary nodes |
| G2 in | Partial contributions of extrema for primitive variables to chare-boundary nodes |
This function receives contributions to m_uNodalExtrm/m_pNodalExtrm , which stores nodal extrems at mesh chare-boundary nodes. While m_uNodalExtrm/m_pNodalExtrm stores own contributions, m_uNodalExtrmc /m_pNodalExtrmc collects the neighbor chare contributions during communication.
void inciter::
Compute the nodal extrema for chare-boundary nodes.
| Parameters | |
|---|---|
| ncomp in | Number of conservative variables |
| nprim in | Number of primitive variables |
| ndof_NodalExtrm in | Degree of freedom for nodal extrema |
| bndel in | List of elements contributing to chare-boundary nodes |
| inpoel in | Element-node connectivity for element e |
| coord in | Array of nodal coordinates |
| gid in | Local->global node id map |
| bid in | Local chare-boundary node ids (value) associated to global node ids (key) |
| U in | Vector of conservative variables |
| P in | Vector of primitive variables |
| uNodalExtrm in/out | Chare-boundary nodal extrema for conservative variables |
| pNodalExtrm in/out | Chare-boundary nodal extrema for primitive variables |
void inciter::
Receive nodal solution (ofor field output) contributions from neighboring chares.
| Parameters | |
|---|---|
| gid in | Global mesh node IDs at which we receive contributions |
| nesup in | Number of elements surrounding points |
| L in | Partial contributions of node fields to chare-boundary nodes |
void inciter::
Receive new mesh from Refiner.
| Parameters | |
|---|---|
| chunk in | New mesh chunk (connectivity and global<->local id maps) |
| coord in | New mesh node coordinates |
| addedTets in | Newly added mesh cells and their parents (local ids) |
| nodeCommMap in | New node communication map |
| bface in | Boundary-faces mapped to side set ids |
| triinpoel in | Boundary-face connectivity |
void inciter::
Extract field output going to file.
| Parameters | |
|---|---|
| chunk in | Field-output mesh chunk (connectivity and global<->local id maps) |
| coord in | Field-output mesh node coordinates |
| addedTets in | Field-output mesh cells and their parents (local ids) |
| nodeCommMap in | Field-output mesh node communication map |
| bface in | Field-output meshndary-faces mapped to side set ids |
| triinpoel in | Field-output mesh boundary-face connectivity |
| c in | Function to continue with after the write |
const tk::Fields& inciter::
| Returns | Const-ref to current solution |
|---|
Const-ref access to current solution
void inciter::
Pack/Unpack serialize member function.
| Parameters | |
|---|---|
| p in/out | Charm++'s PUP::er serializer object reference |
Variable documentation
std::unordered_map<std::size_t, std::pair<std::vector<tk::
Receive buffer for communication of node output fields
Key: global node id, value: output fields and number of elements surrounding the node