Go to the documentation of this file.
18 #ifndef __DUDLEY_NODEFILE_H__
19 #define __DUDLEY_NODEFILE_H__
24 #include <escript/Distribution.h>
27 #include <paso/Coupler.h>
29 #ifdef ESYS_HAVE_TRILINOS
30 #include <trilinoswrap/types.h>
116 #ifdef ESYS_HAVE_TRILINOS
118 esys_trilinos::const_TrilinosGraph_ptr getTrilinosGraph()
const {
137 #ifdef ESYS_HAVE_TRILINOS
138 esys_trilinos::const_TrilinosGraph_ptr m_graph;
169 #ifdef ESYS_HAVE_PASO
247 #endif // __DUDLEY_NODEFILE_H__
index_t * globalDegreesOfFreedom
Definition: dudley/src/NodeFile.h:157
void updateTagList()
Definition: finley/src/NodeFile.h:326
IndexPair getFlaggedMinMaxInt(dim_t N, const index_t *values, index_t ignore)
Definition: finley/src/Util.cpp:305
int MPI_Status
Definition: EsysMPI.h:45
int getNumDataPointsPerSample() const
Return the number of data points per sample.
Definition: Data.h:544
NodeFile(int nDim, escript::JMPI MPIInfo)
Definition: finley/src/NodeFile.cpp:104
bool isComplex() const
True if components of this data are stored as complex.
Definition: Data.cpp:1160
void gather(const index_t *index, const NodeFile *in)
Definition: NodeFile_gather.cpp:77
index_t * globalReducedNodesIndex
Definition: finley/src/NodeFile.h:181
double real_t
type of all real-valued scalars in escript
Definition: DataTypes.h:51
Definition: finley/src/NodeMapping.h:30
dim_t numTargets
size of map (number of target nodes, e.g. DOF, reduced DOF, etc.)
Definition: dudley/src/NodeMapping.h:96
escript::Distribution_ptr nodesDistribution
MPI distribution of nodes.
Definition: dudley/src/NodeFile.h:164
int numDim
number of spatial dimensions
Definition: dudley/src/NodeFile.h:146
dim_t createDenseDOFLabeling()
Definition: finley/src/NodeFile.cpp:525
void assign(const std::vector< index_t > &theTarget, index_t unused)
Definition: finley/src/NodeMapping.h:41
#define DUDLEY_INITIAL_STATUS
Definition: Dudley.h:67
Definition: dudley/src/NodeFile.h:40
std::pair< index_t, index_t > getGlobalIdRange() const
Definition: finley/src/NodeFile.cpp:226
index_t * globalDegreesOfFreedom
Definition: finley/src/NodeFile.h:173
boost::shared_ptr< Distribution > Distribution_ptr
Definition: Distribution.h:25
void createNodeMappings(const IndexVector &dofDistribution, const IndexVector &nodeDistribution)
Definition: NodeFile_createMappings.cpp:224
void setCoordinates(const escript::Data &newX)
Definition: dudley/src/NodeFile.cpp:172
void createDOFMappingAndCoupling()
Definition: NodeFile_createMappings.cpp:29
std::pair< index_t, index_t > getGlobalDOFRange() const
Definition: finley/src/NodeFile.cpp:231
#define INDEX2(_X1_, _X2_, _N1_)
Definition: index.h:23
const index_t * borrowTargetDegreesOfFreedom() const
returns the mapping from local degrees of freedom to a target
Definition: dudley/src/NodeFile.h:234
void assignMPIRankToDOFs(std::vector< int > &mpiRankOfDOF, const IndexVector &distribution)
Definition: finley/src/NodeFile.cpp:438
int * Tag
Tag[i] is the tag of node i.
Definition: finley/src/NodeFile.h:166
const DataTypes::real_t * getSampleDataRO(DataTypes::RealVectorType::size_type sampleNo, DataTypes::real_t dummy=0) const
Return the sample data for the given sample no. Please do not use this unless you NEED to access samp...
Definition: Data.h:1989
IndexPair getMinMaxInt(int dim, dim_t N, const index_t *values)
Definition: dudley/src/Util.cpp:159
Definition: FinleyException.h:27
void clear()
resets both map and target.
Definition: finley/src/NodeMapping.h:32
escript::Distribution_ptr nodesDistribution
MPI distribution of nodes.
Definition: finley/src/NodeFile.h:186
#define MPI_SUM
Definition: EsysMPI.h:52
dim_t getNumNodes() const
returns the number of FEM nodes (on this rank)
Definition: dudley/src/NodeFile.h:204
dim_t createDenseDOFLabeling()
Definition: NodeFile_createDenseLabelings.cpp:22
index_t * reducedNodesId
Definition: finley/src/NodeFile.h:203
index_t * globalReducedDOFIndex
Definition: finley/src/NodeFile.h:178
dim_t numNodes
number of nodes
Definition: finley/src/NodeFile.h:151
void print() const
Definition: dudley/src/NodeFile.cpp:108
Definition: SharedComponents.h:42
escript::JMPI MPIInfo
MPI information.
Definition: dudley/src/NodeFile.h:144
NodeMapping reducedNodesMapping
Definition: finley/src/NodeFile.h:157
std::pair< index_t, index_t > getGlobalIdRange() const
Definition: dudley/src/NodeFile.cpp:157
void updateTagList()
Definition: dudley/src/NodeFile.h:239
void setTags(int newTag, const escript::Data &mask)
set tags to newTag where mask > 0
Definition: dudley/src/NodeFile.cpp:201
bool checkResult(int res, int &mres, const JMPI &info)
Everyone puts in their error code and everyone gets the largest one.
Definition: EsysMPI.cpp:110
int getDataPointSize() const
Return the size of the data point. It is the product of the data point shape dimensions.
Definition: Data.cpp:1360
index_t getLastNode() const
Definition: dudley/src/NodeFile.h:189
void gather(const index_t *index, const NodeFile *in)
Definition: finley/src/NodeFile.cpp:321
void createDOFMappingAndCoupling(bool reduced)
Definition: finley/src/NodeFile.cpp:772
NodeMapping nodesMapping
Definition: dudley/src/NodeFile.h:131
std::pair< index_t, index_t > getGlobalNodeIDIndexRange() const
Definition: finley/src/NodeFile.cpp:236
std::pair< index_t, index_t > getGlobalDOFRange() const
Definition: dudley/src/NodeFile.cpp:162
Definition: dudley/src/NodeMapping.h:28
int status
Definition: dudley/src/NodeFile.h:177
std::pair< index_t, index_t > getGlobalNodeIDIndexRange() const
Definition: dudley/src/NodeFile.cpp:167
std::vector< index_t > target
target[i] defines the target of FEM node i=0,...,numNodes-1
Definition: finley/src/NodeMapping.h:82
static std::pair< index_t, index_t > getGlobalRange(dim_t n, const index_t *id, escript::JMPI mpiInfo)
Definition: finley/src/NodeFile.cpp:30
Definition: Distribution.h:33
#define MPI_INT
Definition: EsysMPI.h:46
NodeMapping reducedDegreesOfFreedomMapping
Definition: finley/src/NodeFile.h:148
void copyTable(index_t offset, index_t idOffset, index_t dofOffset, const NodeFile *in)
Definition: dudley/src/NodeFile.cpp:125
#define MPI_MAX
Definition: EsysMPI.h:54
const index_t * borrowDegreesOfFreedomTarget() const
returns the mapping from target to the local degrees of freedom
Definition: dudley/src/NodeFile.h:224
void freeTable()
empties the node table and frees all memory
Definition: dudley/src/NodeFile.cpp:90
dim_t getNumNodes() const
returns the number of FEM nodes (on this rank)
Definition: finley/src/NodeFile.h:256
const index_t * borrowNodesTarget() const
returns the mapping from target to the local nodes
Definition: dudley/src/NodeFile.h:219
void createNodeMappings(const IndexVector &indexReducedNodes, const IndexVector &dofDistribution, const IndexVector &nodeDistribution)
Definition: finley/src/NodeFile.cpp:1004
int getNumSamples() const
Return the number of samples.
Definition: Data.h:533
dim_t getNumDegreesOfFreedomTargets() const
returns the number of degrees of freedom targets (own and shared)
Definition: dudley/src/NodeFile.h:214
double * Coordinates
Coordinates[INDEX2(k,i,numDim)] is the k-th coordinate of node i.
Definition: finley/src/NodeFile.h:175
A suite of factory methods for creating 2D and 3D dudley domains.
Definition: dudley/src/Assemble.h:32
Data represents a collection of datapoints.
Definition: Data.h:64
void assignMPIRankToDOFs(int *mpiRankOfDOF, const IndexVector &distribution)
Definition: dudley/src/NodeFile.cpp:220
index_t dim_t
Definition: DataTypes.h:65
boost::shared_ptr< JMPI_ > JMPI
Definition: EsysMPI.h:72
NodeMapping degreesOfFreedomMapping
Definition: dudley/src/NodeFile.h:132
const index_t * borrowTargetNodes() const
returns the mapping from local nodes to a target
Definition: dudley/src/NodeFile.h:229
dim_t createDenseReducedLabeling(const std::vector< short > &reducedMask, bool useNodes)
Definition: finley/src/NodeFile.cpp:702
std::vector< int > tagsInUse
vector of tags which are actually used
Definition: dudley/src/NodeFile.h:152
IndexPair getMinMaxInt(int dim, dim_t N, const index_t *values)
Definition: finley/src/Util.cpp:278
void gather_global(const index_t *index, const NodeFile *in)
Definition: NodeFile_gather.cpp:85
~NodeFile()
destructor
Definition: finley/src/NodeFile.cpp:122
escript::JMPI MPIInfo
MPI information.
Definition: finley/src/NodeFile.h:160
dim_t prepareLabeling(const std::vector< short > &mask, IndexVector &buffer, IndexVector &distribution, bool useNodes)
Definition: finley/src/NodeFile.cpp:464
void freeTable()
empties the node table and frees all memory
Definition: finley/src/NodeFile.cpp:162
std::vector< index_t > map
maps the target nodes back to the FEM nodes: target[map[i]]=i
Definition: finley/src/NodeMapping.h:84
dim_t createDenseNodeLabeling(IndexVector &nodeDistribution, const IndexVector &dofDistribution)
An exception class for assertions within escript.
Definition: EsysException.h:60
index_t * target
target[i] defines the target of FEM node i=0,...,numNodes
Definition: dudley/src/NodeMapping.h:93
index_t * Id
Id[i] is the unique ID number of FEM node i.
Definition: finley/src/NodeFile.h:164
dim_t createDenseNodeLabeling(IndexVector &nodeDistribution, const IndexVector &dofDistribution)
Definition: finley/src/NodeFile.cpp:588
NodeMapping nodesMapping
Definition: finley/src/NodeFile.h:146
static void gatherEntries(dim_t n, const index_t *index, index_t min_index, index_t max_index, index_t *Id_out, const index_t *Id_in, int *Tag_out, const int *Tag_in, index_t *globalDegreesOfFreedom_out, const index_t *globalDegreesOfFreedom_in, int numDim, double *Coordinates_out, const double *Coordinates_in)
Definition: finley/src/NodeFile.cpp:77
int status
Definition: finley/src/NodeFile.h:209
void setCoordinates(const escript::Data &newX)
copies the array newX into this->coordinates
Definition: finley/src/NodeFile.cpp:242
boost::shared_ptr< SharedComponents > SharedComponents_ptr
Definition: SharedComponents.h:36
void print() const
Definition: finley/src/NodeFile.cpp:196
static std::pair< index_t, index_t > getGlobalRange(dim_t n, const index_t *id, escript::JMPI mpiInfo)
Definition: dudley/src/NodeFile.cpp:24
NodeFile(int nDim, escript::JMPI MPIInfo)
Definition: dudley/src/NodeFile.cpp:44
index_t getFirstNode() const
Definition: dudley/src/NodeFile.h:184
dim_t getGlobalNumNodes() const
Definition: dudley/src/NodeFile.h:194
escript::Distribution_ptr reducedDegreesOfFreedomDistribution
Definition: finley/src/NodeFile.h:189
void setValuesInUse(const int *values, dim_t numValues, std::vector< int > &valuesInUse, escript::JMPI mpiinfo)
Definition: dudley/src/Util.cpp:224
#define MPI_DOUBLE
Definition: EsysMPI.h:47
std::pair< index_t, index_t > getDOFRange() const
Definition: dudley/src/NodeFile.cpp:146
void allocTable(dim_t numNodes)
allocates the node table within this node file to hold numNodes nodes.
Definition: dudley/src/NodeFile.cpp:63
dim_t prepareLabeling(const std::vector< short > &mask, IndexVector &buffer, IndexVector &distribution, bool useNodes)
escript::Distribution_ptr degreesOfFreedomDistribution
Definition: finley/src/NodeFile.h:188
index_t * map
maps the target nodes back to the FEM nodes: target[map[i]]=i
Definition: dudley/src/NodeMapping.h:99
std::pair< index_t, index_t > getDOFRange() const
Definition: finley/src/NodeFile.cpp:215
int index_t
type for array/matrix indices used both globally and on each rank
Definition: DataTypes.h:60
escript::Distribution_ptr dofDistribution
MPI distribution of degrees of freedom.
Definition: dudley/src/NodeFile.h:167
Definition: finley/src/NodeFile.h:42
Definition: AbstractContinuousDomain.cpp:23
index_t * reducedDegreesOfFreedomId
Definition: finley/src/NodeFile.h:205
NodeMapping degreesOfFreedomMapping
Definition: finley/src/NodeFile.h:147
index_t * globalNodesIndex
assigns each local node a global unique ID in a dense labeling
Definition: finley/src/NodeFile.h:183
index_t * Id
Id[i] is the unique ID number of FEM node i.
Definition: dudley/src/NodeFile.h:148
int * Tag
Tag[i] is the tag of node i.
Definition: dudley/src/NodeFile.h:150
escript::Distribution_ptr reducedNodesDistribution
Definition: finley/src/NodeFile.h:187
boost::shared_ptr< Connector > Connector_ptr
Definition: Coupler.h:38
dim_t numNodes
number of nodes
Definition: dudley/src/NodeFile.h:135
static void scatterEntries(dim_t n, const index_t *index, index_t min_index, index_t max_index, index_t *Id_out, const index_t *Id_in, int *Tag_out, const int *Tag_in, index_t *globalDegreesOfFreedom_out, const index_t *globalDegreesOfFreedom_in, int numDim, double *Coordinates_out, const double *Coordinates_in)
Definition: finley/src/NodeFile.cpp:51
index_t * degreesOfFreedomId
Definition: finley/src/NodeFile.h:204
Definition: escriptcore/src/IndexList.h:29
const index_t * borrowGlobalNodesIndex() const
Definition: dudley/src/NodeFile.h:199
void scatter(const index_t *index, const NodeFile *in)
Definition: finley/src/NodeFile.cpp:310
void gather_global(const index_t *index, const NodeFile *in)
Definition: finley/src/NodeFile.cpp:329
bool shipString(const char *src, char **dest, MPI_Comm &comm)
Definition: EsysMPI.cpp:164
std::vector< int > tagsInUse
vector of tags which are actually used
Definition: finley/src/NodeFile.h:168
index_t * degreesOfFreedomId
Definition: dudley/src/NodeFile.h:173
~NodeFile()
destructor
Definition: dudley/src/NodeFile.cpp:58
A suite of factory methods for creating various finley domains.
Definition: finley/src/Assemble.h:32
void setTags(int newTag, const escript::Data &mask)
set tags to newTag where mask > 0
Definition: finley/src/NodeFile.cpp:267
#define FINLEY_INITIAL_STATUS
Definition: Finley.h:77
int numDim
number of spatial dimensions
Definition: finley/src/NodeFile.h:162
#define ESYS_ASSERT(a, b)
EsysAssert is a MACRO that will throw an exception if the boolean condition specified is false.
Definition: Assert.h:79
std::vector< index_t > packMask(const std::vector< short > &mask)
Definition: finley/src/Util.cpp:332
dim_t getNumDegreesOfFreedom() const
returns the number of degrees of freedom (on this rank)
Definition: dudley/src/NodeFile.h:209
void allocTable(dim_t numNodes)
allocates the node table within this node file to hold numNodes nodes.
Definition: finley/src/NodeFile.cpp:127
double * Coordinates
Coordinates[INDEX2(k,i,numDim)] is the k-th coordinate of node i.
Definition: dudley/src/NodeFile.h:159
An exception class that signals an invalid argument value.
Definition: EsysException.h:90
void clear()
resets both map and target
Definition: dudley/src/NodeMapping.h:32
void copyTable(index_t offset, index_t idOffset, index_t dofOffset, const NodeFile *in)
Definition: finley/src/NodeFile.cpp:285
index_t * globalNodesIndex
assigns each local node a global unique ID in a dense labeling
Definition: dudley/src/NodeFile.h:161
std::vector< index_t > IndexVector
Definition: DataTypes.h:63
dim_t getNumTargets() const
returns the number of target nodes (number of items in the map array)
Definition: finley/src/NodeMapping.h:79