Discretization class
Discretization Charm++ chare array holding common functinoality to all discretization schemes.
Contents
Public types
-
template<typename BCType>struct BCNodes
- Function object for querying the node ids at which a particular BCType BC is configured by the user, called for each PDE type.
Public static functions
- static void registerReducers()
- Configure Charm++ reduction types.
Constructors, destructors, conversion operators
- Discretization(CkMigrateMessage*) explicit
- Migrate constructor.
Public functions
-
auto Discretization(const CProxy_DistFCT& fctproxy,
const CProxy_Transporter& transporter,
const tk::CProxy_MeshWriter& meshwriter,
const std::vector<std::size_t>& ginpoel,
const tk::
UnsMesh:: CoordMap& coordmap, const tk:: CommMaps& msum, int nc) -> Discretization_SDAG_CODE explicit - Constructor.
-
void resizePostAMR(const tk::
UnsMesh:: Chunk& chunk, const tk:: UnsMesh:: Coords& coord, const tk:: NodeCommMap& nodeCommMap) - Resize mesh data structures (e.g., after mesh refinement)
- void startvol()
- Get ready for (re-)computing/communicating nodal volumes.
- void vol()
- Sum mesh volumes to nodes, start communicating them on chare-boundaries.
- void setRefiner(const CProxy_Refiner& ref)
- Set Refiner Charm++ proxy.
-
void comvol(const std::vector<std::size_t>& gid,
const std::vector<tk::
real>& nodevol) - Collect nodal volumes across chare boundaries.
- void totalvol()
- Sum mesh volumes and contribute own mesh volume to total volume.
-
void stat(tk::
real mesh_volume) - Compute mesh cell statistics.
- void boxvol(const std::vector<std::size_t>& boxnodes)
- Compute total box IC volume.
-
void setdt(tk::
real newdt) - Set time step size.
- void next()
- Prepare for next step.
- void status()
- Otput one-liner status report.
- auto histfilename(const std::string& id, kw::precision::info::expect::type precision) -> std::string
- Construct history output filename.
- void histheader(std::vector<std::string>&& names)
- Output headers for time history files (one for each point)
-
void history(std::vector<std::vector<tk::
real>>&& data) - Output time history for a time step.
-
void write(const std::vector<std::size_t>& inpoel,
const tk::
UnsMesh:: Coords& coord, const std::map<int, std::vector<std::size_t>>& bface, const std::map<int, std::vector<std::size_t>>& bnode, const std::vector<std::size_t>& triinpoel, const std::vector<std::string>& elemfieldnames, const std::vector<std::string>& nodefieldnames, const std::vector<std::string>& nodesurfnames, const std::vector<std::vector<tk:: real>>& elemfields, const std::vector<std::vector<tk:: real>>& nodefields, const std::vector<std::vector<tk:: real>>& nodesurfs, CkCallback c) - Output mesh and fields data (solution dump) to file(s)
- void grindZero()
- Zero grind-timer.
- auto restarted(int nrestart) -> bool
- Detect if just returned from a checkpoint and if so, zero timers.
- void remap(const std::unordered_map<std::size_t, std::size_t>& map)
- Remap mesh data due to new local ids.
-
template<typename BCType>auto bcnodes(const std::map<int, std::vector<std::size_t>>& bface, const std::vector<std::size_t>& triinpoel) const -> std::unordered_map<int, std::unordered_set<std::size_t>>
- Query nodes at which BCType boundary conditions are set for all PDE types.
Accessors
-
auto Coord() const -> const tk::
UnsMesh:: Coords& -
auto Coord() -> tk::
UnsMesh:: Coords& - Coordinates accessors as non-const ref.
- auto Gid() const -> const std::vector<std::size_t>&
- Global ids accessors as const-ref.
- auto Gid() -> std::vector<std::size_t>&
- Global ids accessors as non-const-ref.
- auto Lid() const -> const std::unordered_map<std::size_t, std::size_t>&
- Local ids accessors as const-ref.
- auto Lid() -> std::unordered_map<std::size_t, std::size_t>&
- Local ids accessors as non-const-ref.
- auto Inpoel() const -> const std::vector<std::size_t>&
- Tetrahedron element connectivity (with local ids) accessors as const-ref.
- auto Inpoel() -> std::vector<std::size_t>&
- Tetrahedron element connectivity (with local ids) accessors as non-const-ref.
-
auto meshvol() const -> tk::
real - Total mesh volume accessor.
-
auto V() const -> const std::vector<tk::
real>& - Nodal mesh volume accessors const-ref.
-
auto V() -> std::vector<tk::
real>& - Nodal mesh volume accessors non-const-ref.
-
auto Vol() const -> const std::vector<tk::
real>& - Nodal mesh volumes accessors as const-ref.
-
auto Vol() -> std::vector<tk::
real>& - Nodal mesh volumes accessors as non-const-ref.
- auto Hist() const -> const std::vector<HistData>&
- History points data accessor as const-ref.
-
auto Boxvol() -> tk::
real& - Box volume accessor.
-
auto Dt() const -> tk::
real - Time step size accessor.
-
auto T() const -> tk::
real - Physical time accessor.
- auto It() const -> uint64_t
- Iteration count accessor.
- auto Itr() -> uint64_t&
- Non-const-ref refinement iteration count accessor.
- auto Itf() -> uint64_t&
- Non-const-ref field-output iteration count accessor.
- auto Nrestart() -> int&
- Non-const-ref number of restarts accessor.
-
auto Timer() const -> const tk::
Timer& - Timer accessor as const-ref.
-
auto Timer() -> tk::
Timer& - Timer accessor as non-const-ref.
- auto refined() const -> int
- Accessor to flag indicating if the mesh was refined as a value.
- auto refined() -> int&
- Accessor to flag indicating if the mesh was refined as non-const-ref.
- auto Tr() const -> const CProxy_Transporter&
- Transporter proxy accessor as const-ref.
- auto Tr() -> CProxy_Transporter&
- Transporter proxy accessor as non-const-ref.
- auto Ref() const -> Refiner*
- Access bound Refiner class pointer.
- auto FCT() const -> DistFCT*
- Access bound DistFCT class pointer.
- auto Bid() const -> const std::unordered_map<std::size_t, std::size_t>&
- Boundary node ids accessor as const-ref.
- auto Bid() -> std::unordered_map<std::size_t, std::size_t>&
- Boundary node ids accessor as non-const-ref.
-
auto NodeCommMap() const -> const tk::
NodeCommMap& - Node communication map accessor as const-ref.
-
auto NodeCommMap() -> tk::
NodeCommMap& - Node communication map accessor as non-const-ref.
-
auto EdgeCommMap() const -> const tk::
EdgeCommMap& - Edge communication map accessor as const-ref.
-
auto EdgeCommMap() -> tk::
EdgeCommMap& - Edge communication map accessor as non-const-ref.
Charm++ pack/unpack serializer member functions
- void pup(PUP::er& p) override
- Pack/Unpack serialize member function.
- void operator|(PUP::er& p, Discretization& i)
- Pack/Unpack serialize operator|.
Function documentation
static void inciter::
Configure Charm++ reduction types.
Since this is a [initnode] routine, see the .ci file, the Charm++ 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:/
Discretization_SDAG_CODE inciter::
Constructor.
| Parameters | |
|---|---|
| fctproxy in | Distributed FCT proxy |
| transporter in | Host (Transporter) proxy |
| meshwriter in | Mesh writer proxy |
| ginpoel in | Vector of mesh element connectivity owned (global IDs) |
| coordmap in | Coordinates of mesh nodes and their global IDs |
| msum in | Communication maps associated to chare IDs bordering the mesh chunk we operate on |
| nc in | Total number of Discretization chares |
void inciter::
Resize mesh data structures (e.g., after mesh refinement)
| Parameters | |
|---|---|
| chunk in | New mesh chunk (connectivity and global<->local id maps) |
| coord in | New mesh node coordinates |
| nodeCommMap in | New node communication map |
void inciter::
Set Refiner Charm++ proxy.
| Parameters | |
|---|---|
| ref in | Incoming refiner proxy to store |
void inciter::
Collect nodal volumes across chare boundaries.
| Parameters | |
|---|---|
| gid in | Global mesh node IDs at which we receive volume contributions |
| nodevol in | Partial sums of nodal volume contributions to chare-boundary nodes |
This function receives contributions to m_vol, which stores the nodal volumes. While m_vol stores own contributions, m_volc collects the neighbor chare contributions during communication. This way work on m_vol and m_volc is overlapped. The contributions are applied in totalvol().
std::string inciter::
Construct history output filename.
| Parameters | |
|---|---|
| id in | History point id |
| precision in | Floating point precision to use for output |
| Returns | History file name |
void inciter::
Output headers for time history files (one for each point)
| Parameters | |
|---|---|
| names in | History output variable names |
void inciter::
Output mesh and fields data (solution dump) to file(s)
| Parameters | |
|---|---|
| inpoel in | Mesh connectivity for the mesh chunk to be written |
| coord in | Node coordinates of the mesh chunk to be written |
| bface in | Map of boundary-face lists mapped to corresponding side set ids for this mesh chunk |
| bnode in | Map of boundary-node lists mapped to corresponding side set ids for this mesh chunk |
| triinpoel in | Interconnectivity of points and boundary-face in this mesh chunk |
| elemfieldnames in | Names of element fields to be output to file |
| nodefieldnames in | Names of node fields to be output to file |
| nodesurfnames in | Names of node surface fields to be output to file |
| elemfields in | Field data in mesh elements to output to file |
| nodefields in | Field data in mesh nodes to output to file |
| nodesurfs in | Surface field data in mesh nodes to output to file |
| c in | Function to continue with after the write |
Since m_meshwriter is a Charm++ chare group, it never migrates and an instance is guaranteed on every PE. We index the first PE on every logical compute node. In Charm++'s non-SMP mode, a node is the same as a PE, so the index is the same as CkMyPe(). In SMP mode the index is the first PE on every logical node. In non-SMP mode this yields one or more output files per PE with zero or non-zero virtualization, respectively. If there are multiple chares on a PE, the writes are serialized per PE, since only a single entry method call can be executed at any given time. In SMP mode, still the same number of files are output (one per chare), but the output is serialized through the first PE of each compute node. In SMP mode, channeling multiple files via a single PE on each node is required by NetCDF and HDF5, as well as ExodusII, since none of these libraries are thread-safe.
bool inciter::
Detect if just returned from a checkpoint and if so, zero timers.
| Parameters | |
|---|---|
| nrestart in | Number of times restarted |
| Returns | True if restart detected |
void inciter::
Remap mesh data due to new local ids.
| Parameters | |
|---|---|
| map in | Mapping of old->new local ids |
template<typename BCType>
std::unordered_map<int, std::unordered_set<std::size_t>> inciter::
Query nodes at which BCType boundary conditions are set for all PDE types.
| Template parameters | |
|---|---|
| BCType | Boundary condition type |
| Parameters | |
| bface in | Boundary-faces mapped to side set ids |
| triinpoel in | Boundary-face connectivity |
| Returns | Node ids at which BCType BCs are set (value), associated to sides set id on which the BC is set (key) |
const tk::
Coordinates accessors as const-ref
void inciter::
Pack/Unpack serialize member function.
| Parameters | |
|---|---|
| p in/out | Charm++'s PUP::er serializer object reference |
void inciter::
Pack/Unpack serialize operator|.
| Parameters | |
|---|---|
| p in/out | Charm++'s PUP::er serializer object reference |
| i in/out | Discretization object reference |