1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
|
/*****************************************************************************
*
* Copyright (c) 2003-2020 by The University of Queensland
* http://www.uq.edu.au
*
* Primary Business: Queensland, Australia
* Licensed under the Apache License, version 2.0
* http://www.apache.org/licenses/LICENSE-2.0
*
* Development until 2012 by Earth Systems Science Computational Center (ESSCC)
* Development 2012-2013 by School of Earth Sciences
* Development from 2014-2017 by Centre for Geoscience Computing (GeoComp)
* Development from 2019 by School of Earth and Environmental Sciences
**
*****************************************************************************/
#ifndef __WEIPA_FINLEYDOMAIN_H__
#define __WEIPA_FINLEYDOMAIN_H__
#include <weipa/DomainChunk.h>
#include <weipa/FinleyElements.h>
#include <boost/enable_shared_from_this.hpp>
namespace finley {
class Mesh;
}
namespace weipa {
/// \brief Represents a full Finley or Dudley domain including nodes and
/// elements.
///
/// This class represents a Finley or Dudley domain including nodes, cells,
/// and face elements (plus contact elements for Finley). It provides
/// functionality to read a domain from a netCDF file (generated by the
/// domain's dump() method) or directly through an instance of finley::Mesh or
/// Dudley_Mesh.
///
/// Once initialised, the domain can be saved to a Silo file or its nodes
/// and elements accessed through the respective methods.
class FinleyDomain : public DomainChunk, public boost::enable_shared_from_this<FinleyDomain>
{
public:
FinleyDomain();
FinleyDomain(const FinleyDomain& m);
virtual ~FinleyDomain();
virtual bool initFromEscript(const escript::AbstractDomain* domain);
virtual bool initFromFile(const std::string& filename);
virtual bool writeToSilo(DBfile* dbfile, const std::string& pathInSilo,
const StringVec& labels, const StringVec& units,
bool writeMeshData);
virtual void reorderGhostZones(int ownIndex);
virtual void removeGhostZones(int ownIndex);
virtual StringVec getMeshNames() const;
virtual StringVec getVarNames() const;
virtual ElementData_ptr getElementsByName(const std::string& name) const;
virtual NodeData_ptr getMeshByName(const std::string& name) const;
virtual DataVar_ptr getDataVarByName(const std::string& name) const;
virtual Centering getCenteringForFunctionSpace(int fsCode) const;
virtual NodeData_ptr getMeshForFunctionSpace(int fsCode) const;
virtual ElementData_ptr getElementsForFunctionSpace(int fsCode) const;
virtual NodeData_ptr getNodes() const { return nodes; }
virtual std::string getSiloPath() const { return siloPath; }
virtual void setSiloPath(const std::string& path) { siloPath = path; }
private:
void cleanup();
bool initialized;
FinleyNodes_ptr nodes;
FinleyElements_ptr cells;
FinleyElements_ptr faces;
FinleyElements_ptr contacts;
std::string siloPath;
};
} // namespace weipa
#endif // __WEIPA_FINLEYDOMAIN_H__
|
