escript  Revision_
RipleyElements.h
Go to the documentation of this file.
1 
2 /*****************************************************************************
3 *
4 * Copyright (c) 2003-2020 by The University of Queensland
5 * http://www.uq.edu.au
6 *
7 * Primary Business: Queensland, Australia
8 * Licensed under the Apache License, version 2.0
9 * http://www.apache.org/licenses/LICENSE-2.0
10 *
11 * Development until 2012 by Earth Systems Science Computational Center (ESSCC)
12 * Development 2012-2013 by School of Earth Sciences
13 * Development from 2014-2017 by Centre for Geoscience Computing (GeoComp)
14 * Development from 2019 by School of Earth and Environmental Sciences
15 **
16 *****************************************************************************/
17 
18 #ifndef __WEIPA_RIPLEYELEMENTS_H__
19 #define __WEIPA_RIPLEYELEMENTS_H__
20 
21 #include <weipa/ElementData.h>
22 #include <weipa/RipleyNodes.h>
23 
24 class DBfile;
25 
26 namespace ripley {
27 class RipleyDomain;
28 }
29 
30 
31 namespace weipa {
32 
34 typedef boost::shared_ptr<RipleyElements> RipleyElements_ptr;
35 
46 {
47 public:
48 
50  RipleyElements(const std::string& elementName, RipleyNodes_ptr nodes);
51 
54 
56  virtual ~RipleyElements() {}
57 
59  bool initFromRipley(const ripley::RipleyDomain* ripleyDomain, int fsType);
60 
63  void reorderGhostZones(int ownIndex);
64 
66  void removeGhostZones(int ownIndex);
67 
69  virtual void writeConnectivityVTK(std::ostream& os);
70 
76  bool writeToSilo(DBfile* dbfile, const std::string& siloPath,
77  const StringVec& labels, const StringVec& units,
78  bool writeMeshData);
79 
81  virtual StringVec getMeshNames() const;
82 
84  virtual StringVec getVarNames() const;
85 
87  virtual int getNumElements() const { return numElements; }
88 
90  virtual int getNodesPerElement() const { return nodesPerElement; }
91 
93  virtual int getGhostCount() const { return numGhostElements; }
94 
96  virtual ZoneType getType() const { return type; }
97 
99  virtual const IntVec& getNodeList() const { return nodes; }
100 
102  virtual const IntVec& getIDs() const { return ID; }
103 
107  virtual const IntVec& getVarDataByName(const std::string varName) const;
108 
110  virtual NodeData_ptr getNodes() const { return nodeMesh; }
111 
113  virtual ElementData_ptr getReducedElements() const { return ElementData_ptr(); }
114 
116  virtual const QuadMaskInfo& getQuadMask(int functionSpace) const { return quadMask; }
117 
121  virtual int getElementFactor() const { return 1; }
122 
123 private:
125  void buildMeshes();
126  IntVec prepareGhostIndices(int ownIndex);
127  void reorderArray(IntVec& v, const IntVec& idx, int elementsPerIndex);
128 
131  std::string name;
139 };
140 
141 } // namespace weipa
142 
143 #endif // __WEIPA_RIPLEYELEMENTS_H__
144 
weipa::RipleyElements::getType
virtual ZoneType getType() const
Returns the type of the elements.
Definition: RipleyElements.h:96
weipa::RipleyElements::ID
IntVec ID
Definition: RipleyElements.h:137
ripley
Definition: ripley/src/AbstractAssembler.h:26
weipa::RipleyElements::numGhostElements
int numGhostElements
Definition: RipleyElements.h:133
ripley::RipleyDomain::getOwnerVector
virtual RankVector getOwnerVector(int fsType) const =0
returns a vector of rank numbers where vec[i]=n means that rank n 'owns' element/face element i.
weipa::RipleyElements::removeGhostZones
void removeGhostZones(int ownIndex)
Removes "ghost" elements.
Definition: RipleyElements.cpp:328
weipa::RipleyElements::reorderArray
void reorderArray(IntVec &v, const IntVec &idx, int elementsPerIndex)
Definition: RipleyElements.cpp:273
weipa::ZONETYPE_HEX
@ ZONETYPE_HEX
Definition: ElementData.h:58
weipa::RipleyElements::nodesPerElement
int nodesPerElement
Definition: RipleyElements.h:134
weipa::RipleyElements
Stores and manipulates one type of ripley mesh elements (cells, faces).
Definition: RipleyElements.h:45
ripley::RipleyDomain::getNumNodesPerDim
virtual const dim_t * getNumNodesPerDim() const =0
returns the number of nodes per MPI rank in each dimension
weipa::RipleyElements::getNodes
virtual NodeData_ptr getNodes() const
Returns the node mesh instance used by the elements.
Definition: RipleyElements.h:110
weipa::RipleyElements::getGhostCount
virtual int getGhostCount() const
Returns the number of "ghost" elements.
Definition: RipleyElements.h:93
ripley::Elements
@ Elements
Definition: Ripley.h:71
weipa::QuadMaskInfo
This struct holds a mask (0's and 1's) that indicates which quad nodes contribute to a sub-element wh...
Definition: ElementData.h:54
weipa::RipleyElements::getIDs
virtual const IntVec & getIDs() const
Returns a vector of element IDs.
Definition: RipleyElements.h:102
weipa::ZONETYPE_UNKNOWN
@ ZONETYPE_UNKNOWN
Definition: ElementData.h:56
weipa::RipleyElements::getElementFactor
virtual int getElementFactor() const
If the original element type is not supported they are subdivided into N smaller elements (e....
Definition: RipleyElements.h:121
weipa
Definition: DataVar.cpp:51
weipa::ZoneType
ZoneType
Definition: ElementData.h:41
weipa::RipleyElements_ptr
boost::shared_ptr< RipleyElements > RipleyElements_ptr
Definition: RipleyElements.h:33
weipa::RipleyElements::getVarDataByName
virtual const IntVec & getVarDataByName(const std::string varName) const
Returns an array of data values for the name provided.
Definition: RipleyElements.cpp:261
weipa::RipleyNodes
Stores and manipulates ripley mesh nodes.
Definition: RipleyNodes.h:39
weipa::RipleyElements::prepareGhostIndices
IntVec prepareGhostIndices(int ownIndex)
Definition: RipleyElements.cpp:294
paso::util::copy
void copy(dim_t N, double *out, const double *in)
out = in
Definition: PasoUtil.h:116
weipa::RipleyElements::getMeshNames
virtual StringVec getMeshNames() const
Returns the names of the meshes associated with the elements.
Definition: RipleyElements.cpp:244
ripley::RipleyDomain
RipleyDomain extends the AbstractContinuousDomain interface for the Ripley library and is the base cl...
Definition: ripley/src/RipleyDomain.h:116
escript::DataTypes::dim_t
index_t dim_t
Definition: DataTypes.h:93
weipa::CoordArray
std::vector< float * > CoordArray
Definition: weipa.h:62
weipa::RipleyElements::getQuadMask
virtual const QuadMaskInfo & getQuadMask(int functionSpace) const
Returns a QuadMaskInfo structure for given functionspace code.
Definition: RipleyElements.h:116
weipa::ElementData_ptr
boost::shared_ptr< ElementData > ElementData_ptr
Definition: weipa.h:67
weipa::ElementData
Stores and manipulates one type of domain elements.
Definition: ElementData.h:63
weipa::RipleyElements::reorderGhostZones
void reorderGhostZones(int ownIndex)
Moves "ghost" elements (whose owner does not match ownIndex) and the corresponding data to the end of...
Definition: RipleyElements.cpp:315
weipa::RipleyElements::writeToSilo
bool writeToSilo(DBfile *dbfile, const std::string &siloPath, const StringVec &labels, const StringVec &units, bool writeMeshData)
Writes element data into given directory in given Silo file.
Definition: RipleyElements.cpp:388
ripley::RipleyDomain::getDim
virtual int getDim() const
returns the number of spatial dimensions of the domain
Definition: ripley/src/RipleyDomain.h:191
weipa::RipleyElements::name
std::string name
Definition: RipleyElements.h:131
ripley::FaceElements
@ FaceElements
Definition: Ripley.h:73
weipa::RipleyElements::getVarNames
virtual StringVec getVarNames() const
Returns a vector with the mesh variable names.
Definition: RipleyElements.cpp:252
weipa::ZONETYPE_QUAD
@ ZONETYPE_QUAD
Definition: ElementData.h:60
ripley::RipleyDomain::borrowSampleReferenceIDs
const dim_t * borrowSampleReferenceIDs(int fsType) const =0
returns the array of reference numbers for a function space type
weipa::RipleyElements::nodes
IntVec nodes
Definition: RipleyElements.h:136
weipa::RipleyElements::owner
IntVec owner
Definition: RipleyElements.h:137
weipa::RipleyElements::writeConnectivityVTK
virtual void writeConnectivityVTK(std::ostream &os)
Writes connectivity data to a stream in VTK text format.
Definition: RipleyElements.cpp:359
weipa::RipleyElements::initFromRipley
bool initFromRipley(const ripley::RipleyDomain *ripleyDomain, int fsType)
Initialises with data from a RipleyDomain instance.
Definition: RipleyElements.cpp:75
weipa::NodeData_ptr
boost::shared_ptr< NodeData > NodeData_ptr
Definition: weipa.h:69
weipa::RipleyElements::originalMesh
RipleyNodes_ptr originalMesh
Definition: RipleyElements.h:130
weipa::StringVec
std::vector< std::string > StringVec
Definition: weipa.h:61
weipa::RipleyElements::RipleyElements
RipleyElements()
Definition: RipleyElements.h:124
weipa::RipleyElements::getNodeList
virtual const IntVec & getNodeList() const
Returns a vector of the node IDs used by the elements.
Definition: RipleyElements.h:99
weipa::RipleyElements::tag
IntVec tag
Definition: RipleyElements.h:137
ripley::RipleyDomain::getNumFacesPerBoundary
virtual const dim_t * getNumFacesPerBoundary() const =0
returns the number of face elements in the order (left,right,bottom,top,[front,back]) on current MPI ...
weipa::ZONETYPE_POLYGON
@ ZONETYPE_POLYGON
Definition: ElementData.h:59
weipa::RipleyElements::~RipleyElements
virtual ~RipleyElements()
Destructor.
Definition: RipleyElements.h:56
weipa::RipleyElements::buildMeshes
void buildMeshes()
Definition: RipleyElements.cpp:342
weipa::RipleyElements::getNodesPerElement
virtual int getNodesPerElement() const
Returns the number of nodes per element.
Definition: RipleyElements.h:90
ripley::RipleyDomain::getDataShape
virtual std::pair< int, dim_t > getDataShape(int fsType) const
returns the number of data points per sample, and the number of samples as a pair.
Definition: ripley/src/RipleyDomain.cpp:152
weipa::RipleyElements::quadMask
QuadMaskInfo quadMask
Definition: RipleyElements.h:138
weipa::IntVec
std::vector< int > IntVec
Definition: weipa.h:60
ripley::RipleyDomain::getNumElementsPerDim
virtual const dim_t * getNumElementsPerDim() const =0
returns the number of elements per MPI rank in each dimension
weipa::RipleyElements::type
ZoneType type
Definition: RipleyElements.h:135
weipa::RipleyElements::nodeMesh
RipleyNodes_ptr nodeMesh
Definition: RipleyElements.h:129
weipa::ZONETYPE_BEAM
@ ZONETYPE_BEAM
Definition: ElementData.h:57
weipa::RipleyElements::numElements
int numElements
Definition: RipleyElements.h:132
weipa::RipleyElements::getNumElements
virtual int getNumElements() const
Returns the number of elements.
Definition: RipleyElements.h:87
weipa::RipleyNodes_ptr
boost::shared_ptr< RipleyNodes > RipleyNodes_ptr
Definition: RipleyNodes.h:31
weipa::RipleyElements::getReducedElements
virtual ElementData_ptr getReducedElements() const
Returns the reduced elements.
Definition: RipleyElements.h:113