Go to the documentation of this file.
17 #ifndef __WEIPA_ESCRIPTDATASET_H__
18 #define __WEIPA_ESCRIPTDATASET_H__
20 #include <weipa/weipa.h>
88 const std::string units =
"");
100 bool loadNetCDF(
const std::string domainFilePattern,
110 void setCycleAndTime(
int c,
double t) { cycle=c; time=t; }
113 int getCycle()
const {
return cycle; }
116 double getTime()
const {
return time; }
120 void setMeshLabels(
const std::string x,
const std::string y,
const std::string z=
"");
124 void setMeshUnits(
const std::string x,
const std::string y,
const std::string z=
"");
128 void setMetadataSchemaString(
const std::string schema,
129 const std::string metadata)
130 { mdSchema=schema; mdString=metadata; }
133 void setSaveMeshData(
bool flag) { wantsMeshVars=flag; }
136 bool saveSilo(
const std::string fileName,
bool useMultiMesh=
true);
139 void saveVTK(
const std::string fileName);
142 DomainChunks getConvertedDomain() {
return domainChunks; }
145 const VarVector& getVariables()
const {
return variables; }
148 const VarVector& getMeshVariables()
const {
return meshVariables; }
155 getMPIComm() {
return mpiComm; }
158 bool loadDomain(
const std::string filePattern,
int nChunks);
160 bool loadData(
const std::string filePattern,
const std::string name,
161 const std::string units);
163 void convertMeshVariables();
164 void updateSampleDistribution(
VarInfo& vi);
165 void putSiloMultiMesh(DBfile* dbfile,
const std::string& meshName);
166 void putSiloMultiTensor(DBfile* dbfile,
const VarInfo& vi);
167 void putSiloMultiVar(DBfile* dbfile,
const VarInfo& vi,
168 bool useMeshFile =
false);
169 void saveVTKsingle(
const std::string& fileName,
170 const std::string& meshName,
const VarVector& vars);
171 void writeVarToVTK(
const VarInfo& varInfo, std::ostream& os);
175 std::string mdSchema, mdString;
177 bool externalDomain, wantsMeshVars;
180 int mpiRank, mpiSize;
190 #endif // __WEIPA_ESCRIPTDATASET_H__
IntVec sampleDistribution
Definition: EscriptDataset.h:40
Data interpolateFromTable3DP(boost::python::object table, DataTypes::real_t Amin, DataTypes::real_t Astep, Data &B, DataTypes::real_t Bmin, DataTypes::real_t Bstep, Data &C, DataTypes::real_t Cmin, DataTypes::real_t Cstep, DataTypes::real_t undef, bool check_boundaries)
Definition: Data.cpp:4740
void setAcceptanceConvergenceFailure(bool acceptance)
Definition: SolverOptions.cpp:770
Definition: SolverOptions.h:155
DataTypes::dim_t getReferenceIDFromDataPointNo(DataTypes::dim_t dataPointNo) const
Returns the reference number associated with the given data-point number.
Definition: FunctionSpace.cpp:177
double getDropTolerance() const
Definition: SolverOptions.cpp:624
bool isReady() const
Return true if this data is ready.
Definition: Data.cpp:1156
BOOST_PYTHON_MODULE(escriptcpp)
Definition: escriptcpp.cpp:116
Definition: SolverOptions.h:161
virtual void saveMM(const std::string &filename) const
writes the matrix to a file using the Matrix Market file format
Definition: AbstractSystemMatrix.cpp:130
std::string mdSchema
Definition: EscriptDataset.h:174
bool isComplex() const
True if components of this data are stored as complex.
Definition: Data.cpp:1163
void MPIBarrierWorld()
performs a barrier synchronization across all processors.
Definition: Utils.cpp:376
escript::Data getSize() const
Returns the sample size (e.g. the diameter of elements, radius of particles).
Definition: FunctionSpace.cpp:245
Data erf() const
Return the error function erf of each data point of this Data object.
Definition: Data.cpp:2114
Domain_ptr getDomainPython() const
Return the function space domain. Internal use only! This gets around some python difficulties by cas...
Definition: FunctionSpace.cpp:107
int getSvnVersion()
some functions
Definition: Utils.cpp:42
Data sign() const
Return the sign of each data point of this Data object. -1 for negative values, zero for zero values,...
Definition: Data.cpp:2160
void saveDataCSV(const std::string &filename, bp::dict arg, const std::string &sep, const std::string &csep, bool refid, bool append)
Definition: Utils.cpp:387
Reducer_ptr makeDataReducer(std::string type)
Definition: MPIDataReducer.cpp:45
virtual int getTransportTypeId(int solver, int preconditioner, int package, bool symmetry) const
return the identifier of the matrix type to be used for the global stiffness matrix when a particular...
Definition: AbstractContinuousDomain.cpp:141
void setInnerToleranceAdaptionOff()
Definition: SolverOptions.cpp:742
std::string mdString
Definition: EscriptDataset.h:174
void getNumpy(bp::dict arg)
Definition: Utils.cpp:919
void setDropStorage(double drop)
Definition: SolverOptions.cpp:629
DataTypes::real_t getScalarVariable(const std::string &name)
Definition: SplitWorld.cpp:28
Data atanh() const
Return the atanh of each data point of this Data object.
Definition: Data.cpp:2139
Reducer_ptr makeScalarReducer(std::string type)
Definition: MPIScalarReducer.cpp:30
double getMaxFloat()
Definition: Utils.cpp:371
Data atan() const
Return the atan of each data point of this Data object.
Definition: Data.cpp:2085
#define REGISTER_ESCRIPT_EXCEPTION_TRANSLATORS
Definition: ExceptionTranslators.h:24
Definition: SolverOptions.h:127
Data rtruedivO(const boost::python::object &left)
Newer style division operator for python.
Definition: Data.h:2191
FunctionSpace reducedFunctionOnBoundary(const AbstractDomain &domain)
Return a function on boundary FunctionSpace with reduced integration order.
Definition: FunctionSpaceFactory.cpp:71
Data antihermitian() const
Return the anti-hermitian part of a matrix which is half the matrix minus its hermitian.
Definition: Data.cpp:2811
FunctionSpace reducedContinuousFunction(const AbstractDomain &domain)
Return a continuous with reduced order FunctionSpace (overlapped node values on reduced element order...
Definition: FunctionSpaceFactory.cpp:47
void setToZero()
set all values to zero
Definition: Data.cpp:818
Data cos() const
Return the cos of each data point of this Data object.
Definition: Data.cpp:2056
void saveVTKsingle(const std::string &fileName, const std::string &meshName, const VarVector &vars)
Definition: EscriptDataset.cpp:581
DataTypes::real_t Lsup()
Return the maximum absolute value of this Data object.
Definition: Data.cpp:2218
Data grad() const
Definition: Data.cpp:1353
bool hasConverged() const
Definition: SolverOptions.cpp:328
double getDiagnostics(const std::string name) const
Definition: SolverOptions.cpp:304
bool useLocalPreconditioner() const
Definition: SolverOptions.cpp:778
Definition: SolverOptions.h:126
virtual int getTag(const std::string &name) const =0
Return the tag key for tag name.
Represents a full Speckley domain including nodes and elements.
Definition: weipa/src/SpeckleyDomain.h:46
Definition: SolverOptions.h:128
bool loadDomain(const std::string filePattern, int nChunks)
Definition: EscriptDataset.cpp:836
Data Tensor4FromObj(bp::object o, const FunctionSpace &what, bool expanded)
Definition: DataFactory.cpp:209
void removeVariable(std::string name)
Definition: SplitWorld.cpp:298
void setValueOfDataPoint(int dataPointNo, const DataTypes::real_t)
sets the values of a data-point on this process
Definition: Data.cpp:1658
int getDim()
Definition: SolverOptions.cpp:852
bool isEmpty() const
returns true if the matrix is empty
Definition: AbstractSystemMatrix.h:96
boost::python::object getVarPyInfo()
Definition: SplitWorld.cpp:330
void setTruncation(int truncation)
Definition: SolverOptions.cpp:533
Definition: SolverOptions.h:130
virtual ASM_ptr newSystemMatrix(const int row_blocksize, const escript::FunctionSpace &row_functionspace, const int column_blocksize, const escript::FunctionSpace &column_functionspace, const int type) const
creates a SystemMatrixAdapter stiffness matrix and initializes it with zeros:
Definition: AbstractContinuousDomain.cpp:177
void setVerbosityOff()
Definition: SolverOptions.cpp:719
double getTolerance() const
Definition: SolverOptions.cpp:588
Data whereNonNegative() const
Return a Data with a 1 for +ive or 0 values and a 0 for -ive values.
Definition: Data.cpp:1278
Definition: SolverOptions.h:125
bool loadConfigured()
returns true if the load funtion is configured.
Definition: DataFactory.cpp:704
virtual void saveHB(const std::string &filename) const
writes the matrix to a file using the Harwell-Boeing file format
Definition: AbstractSystemMatrix.cpp:135
Data tan() const
Return the tan of each data point of this Data object.
Definition: Data.cpp:2063
Definition: SolverOptions.h:122
std::string getSummary() const
Definition: SolverOptions.cpp:73
void tag()
If possible convert this Data to DataTagged. This will only allow Constant data to be converted to ta...
Definition: Data.cpp:1207
boost::python::object integrateToTuple()
Calculate the integral over the function space domain as a python tuple.
Definition: Data.cpp:1787
double getAbsoluteTolerance() const
Definition: SolverOptions.cpp:600
int getSubWorldCount()
Definition: SplitWorld.cpp:426
Definition: SolverOptions.h:149
void setAcceptanceConvergenceFailureOn()
Definition: SolverOptions.cpp:760
virtual bool isValidTagName(const std::string &name) const
Returns True if name is a defined tag name.
Definition: AbstractDomain.cpp:64
int getTruncation() const
Definition: SolverOptions.cpp:540
Data sinh() const
Return the sinh of each data point of this Data object.
Definition: Data.cpp:2092
Data abs() const
Return the absolute value of each data point of this Data object.
Definition: Data.cpp:2168
WeipaException exception class.
Definition: WeipaException.h:38
virtual DataTypes::dim_t getNumDataPointsGlobal() const
Return the number of data points summed across all MPI processes.
Definition: AbstractContinuousDomain.cpp:197
void setTaggedValueByName(std::string name, const boost::python::object &value)
Assign the given value to the tag assocciated with name. Implicitly converts this object to type Data...
Definition: Data.cpp:3580
Data getItem(const boost::python::object &key) const
Returns a slice from this Data object.
Definition: Data.cpp:3471
Definition: SolverOptions.h:141
Data swapaxes(const int axis0, const int axis1) const
swaps the components axis0 and axis1
Definition: Data.cpp:2668
SolverOptions getPackage() const
Definition: SolverOptions.cpp:487
void resetDiagnostics(bool all=false)
Definition: SolverOptions.cpp:179
VarVector meshVariables
Definition: EscriptDataset.h:178
FunctionSpace continuousFunction(const AbstractDomain &domain)
Create function space objects.
Definition: FunctionSpaceFactory.cpp:41
bool saveSilo(const std::string fileName, bool useMultiMesh=true)
Saves the dataset in the Silo file format.
Definition: EscriptDataset.cpp:307
int getNumSweeps() const
Definition: SolverOptions.cpp:576
void setReordering(int ordering)
Definition: SolverOptions.cpp:492
void resolveGroup(bp::object obj)
Definition: Utils.cpp:924
int getEscriptParamInt(const std::string &name, int sentinel=0)
get the value of a named parameter. See listEscriptParams() for available parameters.
Definition: EscriptParams.h:88
virtual void addPDEToSystem(AbstractSystemMatrix &mat, escript::Data &rhs, const escript::Data &A, const escript::Data &B, const escript::Data &C, const escript::Data &D, const escript::Data &X, const escript::Data &Y, const escript::Data &d, const escript::Data &y, const escript::Data &d_contact, const escript::Data &y_contact, const escript::Data &d_dirac, const escript::Data &y_dirac) const
Return the domain as const AbstractContinuousDomain& / static const AbstractContinuousDomain& asAbst...
Definition: AbstractContinuousDomain.cpp:146
Data interpolate(const FunctionSpace &functionspace) const
Interpolates this onto the given functionspace and returns the result as a Data object.
Definition: Data.cpp:1314
Data exp() const
Return the exponential function of each data point of this Data object.
Definition: Data.cpp:2194
boost::python::object __mul__(const boost::python::object &right)
wrapper for python multiply operation
Definition: Data.cpp:5988
void setSymmetryOn()
Definition: SolverOptions.cpp:668
Definition: FileWriter.h:39
Data interpolateFromTable2DP(boost::python::object table, DataTypes::real_t Amin, DataTypes::real_t Astep, Data &B, DataTypes::real_t Bmin, DataTypes::real_t Bstep, DataTypes::real_t undef, bool check_boundaries)
Definition: Data.cpp:4750
SolverOptions getSolverMethod() const
Definition: SolverOptions.cpp:424
bool isTagged() const
Return true if this Data is tagged.
Definition: Data.cpp:1128
Data ScalarFromObj(boost::python::object o, const FunctionSpace &what, bool expanded)
Python interface for the Scalar factory.
Definition: DataFactory.cpp:75
Data besselSecondKind(int order)
Return the Bessel function of the second kind for each data point of this Data object.
Definition: Data.cpp:1867
#define MPI_SUM
Definition: EsysMPI.h:50
boost::python::object raw_addJob(boost::python::tuple t, boost::python::dict kwargs)
Definition: SplitWorld.cpp:471
Data TensorFromObj(bp::object o, const FunctionSpace &what, bool expanded)
Definition: DataFactory.cpp:133
void writeVarToVTK(const VarInfo &varInfo, std::ostream &os)
Definition: EscriptDataset.cpp:808
bool loadNetCDF(const std::string domainFilePattern, const StringVec &varFiles, const StringVec &varNames, int nChunks)
Loads domain and variables from escript NetCDF files.
Definition: EscriptDataset.cpp:253
int getNumRefinements() const
Definition: SolverOptions.cpp:808
Data gradOn(const FunctionSpace &functionspace) const
Calculates the gradient of the data at the data points of functionspace. If functionspace is not pres...
Definition: Data.cpp:1326
StringVec meshUnits
Definition: EscriptDataset.h:175
virtual escript::Data getX() const =0
Returns locations in the domain. The function space is chosen appropriately.
void setRestart(int restart)
Definition: SolverOptions.cpp:512
std::string units
Definition: EscriptDataset.h:38
Data vectorMultiply(const Data &right) const
returns the matrix-vector product this*right
Definition: AbstractSystemMatrix.cpp:75
Definition: SolverOptions.h:151
std::vector< DomainChunk_ptr > DomainChunks
Definition: EscriptDataset.h:34
Data phase() const
Return the phase/arg/angular-part of complex values.
Definition: Data.cpp:2026
DataTypes::real_t sup()
Return the maximum value of this Data object.
Definition: Data.cpp:2267
void setSymmetryOff()
Definition: SolverOptions.cpp:673
Definition: SolverOptions.h:145
bool isConstant() const
Return true if this Data is constant.
Definition: Data.cpp:1142
Definition: SolverOptions.h:159
void putSiloMultiTensor(DBfile *dbfile, const VarInfo &vi)
Definition: EscriptDataset.cpp:1177
int getMPISizeWorld()
returns the total number of available MPI processes for MPI_COMM_WORLD
Definition: Utils.cpp:93
Definition: SolverOptions.h:142
bool addData(escript::Data &data, const std::string name, const std::string units="")
Adds an escript data instance to this dataset. You must ensure that the data is defined on the domain...
Definition: EscriptDataset.cpp:217
bool isHermitian() const
Definition: SolverOptions.cpp:686
boost::python::object __add__(const boost::python::object &right)
wrapper for python add operation
Definition: Data.cpp:5915
boost::python::list listFeatures()
returns a list of features escript was compiled with.
Definition: EscriptParams.h:115
bool isComplex() const
Definition: SolverOptions.cpp:653
void setHermitianOn()
Definition: SolverOptions.cpp:691
Data interpolateFromTable1DP(boost::python::object table, DataTypes::real_t Amin, DataTypes::real_t Astep, DataTypes::real_t undef, bool check_boundaries)
Definition: Data.cpp:4759
Definition: FunctionSpace.h:45
Data truedivO(const boost::python::object &right)
Newer style division operator for python.
Definition: Data.h:2178
Domain_ptr getDomainPython() const
Return the domain. TODO: For internal use only. This should be removed.
Definition: Data.h:487
void updateDiagnosticsPy(const std::string &key, const boost::python::object &value)
Definition: SolverOptions.cpp:245
Data tanh() const
Return the tanh of each data point of this Data object.
Definition: Data.cpp:2106
Definition: SolverOptions.h:132
Definition: SolverOptions.h:156
const FunctionSpace & getFunctionSpace() const
Return the function space.
Definition: Data.h:461
void complicate()
Definition: Data.cpp:6060
void setIterMax(int iter_max)
Definition: SolverOptions.cpp:557
void setNumSweeps(int sweeps)
Definition: SolverOptions.cpp:569
boost::shared_ptr< DataVar > DataVar_ptr
Definition: weipa.h:63
Definition: DataVar.cpp:49
VarVector variables
Definition: EscriptDataset.h:178
void insertConstraint(Data &source, Data &q, Data &r)
inserts constraint u_{,t}=r where q>0 into the problem
Definition: AbstractTransportProblem.cpp:83
#define MPI_INT
Definition: EsysMPI.h:44
void setNumRefinements(int refinements)
Definition: SolverOptions.cpp:801
void copyVariable(const std::string &src, const std::string &dest)
Definition: SplitWorld.cpp:436
Data whereZero(DataTypes::real_t tol=0.0) const
Return a Data with a 1 for 0 values and a 0 for +ive or -ive values.
Definition: Data.cpp:1300
Represents a full Ripley domain including nodes and elements.
Definition: weipa/src/RipleyDomain.h:46
#define MPI_MAX
Definition: EsysMPI.h:52
Data C_GeneralTensorProduct(Data &arg_0, Data &arg_1, int axis_offset=0, int transpose=0)
Compute a tensor product of two Data objects.
Definition: Data.cpp:3693
void setMeshLabels(const std::string x, const std::string y, const std::string z="")
Sets labels for the mesh axes.
Definition: EscriptDataset.cpp:557
void setDim(int dim)
Definition: SolverOptions.cpp:845
void setAbsoluteTolerance(double atol)
Definition: SolverOptions.cpp:593
void setLocalPreconditioner(bool local)
Definition: SolverOptions.cpp:793
virtual StatusType getStatus() const
Returns a status indicator of the domain. The status identifier should be unique over the live time i...
Definition: AbstractDomain.cpp:59
FunctionSpace reducedSolution(const AbstractDomain &domain)
Return a FunctionSpace with reduced integration order.
Definition: FunctionSpaceFactory.cpp:107
std::vector< DataVar_ptr > DataChunks
Definition: EscriptDataset.h:33
virtual void resetValues(bool preserveSolverData=false)
resets the matrix entries
Definition: AbstractSystemMatrix.cpp:140
boost::python::list listEscriptParams()
describe available parameters.
Definition: EscriptParams.h:97
bool isVerbose() const
Definition: SolverOptions.cpp:709
Definition: SolverOptions.h:112
Reducer_ptr makeNonReducedVariable()
Definition: NonReducedVariable.cpp:121
virtual void addPDEToRHS(escript::Data &rhs, const escript::Data &X, const escript::Data &Y, const escript::Data &y, const escript::Data &y_contact, const escript::Data &y_dirac) const
adds a PDE onto the lumped stiffness matrix matrix / virtual void addPDEToLumpedSystem( escript::Dat...
Definition: AbstractContinuousDomain.cpp:157
int getNumberOfThreads()
returns the number of threads
Definition: Utils.cpp:84
virtual void nullifyRowsAndCols(Data &row_q, Data &col_q, double mdv)
sets matrix entries to zero in specified rows and columns. The rows and columns are marked by positiv...
Definition: AbstractSystemMatrix.cpp:123
escript::Data getNormal() const
Returns the surface normal field.
Definition: FunctionSpace.cpp:236
Data powD(const Data &right) const
Return the given power of each data point of this Data object.
Definition: Data.cpp:3340
FunctionSpace functionOnContactZero(const AbstractDomain &domain)
Return a FunctionSpace on left side of contact.
Definition: FunctionSpaceFactory.cpp:77
void setRelaxationFactor(double factor)
Definition: SolverOptions.cpp:641
char NcFType(const std::string &name)
Definition: NCHelper.cpp:38
virtual bool onMasterProcessor() const =0
Return true if on MPI master, else false.
Data represents a collection of datapoints.
Definition: Data.h:62
std::string toString() const
Write the data as a string. For large amounts of data, a summary is printed.
Definition: Data.cpp:4646
DomainChunks domainChunks
Definition: EscriptDataset.h:177
FunctionSpace reducedFunctionOnContactZero(const AbstractDomain &domain)
Return a FunctionSpace on left side of contact with reduced integration order.
Definition: FunctionSpaceFactory.cpp:83
Data whereNonPositive() const
Return a Data with a 1 for -ive or 0 values and a 0 for +ive values.
Definition: Data.cpp:1289
bool valid
Definition: EscriptDataset.h:41
DataChunks dataChunks
Definition: EscriptDataset.h:39
void resolve()
If this data is lazy, then convert it to ready data. What type of ready data depends on the expressio...
Definition: Data.cpp:1233
escript::Data getX() const
Returns the spatial locations of the data points.
Definition: FunctionSpace.cpp:227
Data powO(const boost::python::object &right) const
Return the given power of each data point of this boost python object.
Definition: Data.cpp:3333
const boost::python::tuple getShapeTuple() const
Return the data point shape as a tuple of integers.
Definition: Data.cpp:727
bool externalDomain
Definition: EscriptDataset.h:176
int runMPIProgram(bp::list args)
Definition: Utils.cpp:311
Definition: SolverOptions.h:158
void putSiloMultiMesh(DBfile *dbfile, const std::string &meshName)
Definition: EscriptDataset.cpp:1053
Data symmetric() const
Return the symmetric part of a matrix which is half the matrix plus its transpose.
Definition: Data.cpp:2725
int cycle
Definition: EscriptDataset.h:172
bool acceptConvergenceFailure() const
Definition: SolverOptions.cpp:755
void setComplex(bool complex)
Definition: SolverOptions.cpp:658
virtual double getUnlimitedTimeStepSize() const
Definition: AbstractTransportProblem.cpp:128
void setTags(const int newTag, const escript::Data &mask) const
Assigns new tag newTag to all samples with a positive value of mask for any of its sample points.
Definition: FunctionSpace.cpp:254
const boost::python::object toListOfTuples(bool scalarastuple=true)
returns the values of the object as a list of tuples (one for each datapoint).
Definition: Data.cpp:1381
FunctionSpace diracDeltaFunctions(const AbstractDomain &domain)
Return a FunctionSpace.
Definition: FunctionSpaceFactory.cpp:113
virtual std::string getDescription() const
Return a description for this domain.
Definition: AbstractContinuousDomain.cpp:38
int getApproximationOrder() const
Returns the approximation order used for this function space.
Definition: FunctionSpace.cpp:309
virtual void dump(const std::string &filename) const =0
dumps the domain to an external file filename.
Data real() const
Definition: Data.cpp:1987
virtual void setTagMap(const std::string &name, int tag)=0
sets a map from a clear tag name to a tag key
Data sin() const
Return the sin of each data point of this Data object.
Definition: Data.cpp:2049
Data maxval() const
Return the maximum value of each data point of this Data object.
Definition: Data.cpp:2650
double getRelaxationFactor() const
Definition: SolverOptions.cpp:648
Data solve(const Data &in, boost::python::object &options) const
returns the solution u of the linear system this*u=in
Definition: AbstractSystemMatrix.cpp:99
boost::shared_ptr< ElementData > ElementData_ptr
Definition: weipa.h:65
Data besselFirstKind(int order)
Return the Bessel function of the first kind for each data point of this Data object.
Definition: Data.cpp:1860
void * mpiComm
Definition: EscriptDataset.h:183
int getInnerIterMax() const
Definition: SolverOptions.cpp:552
int mpiRank
Definition: EscriptDataset.h:179
virtual void resetTransport(bool preserveSolverData) const
resets the transport operator typically as they have been updated.
Definition: AbstractTransportProblem.cpp:120
Data acosh() const
Return the acosh of each data point of this Data object.
Definition: Data.cpp:2132
int getMPIWorldSum(int val)
returns sum of an integer over all processors with MPI_COMM_WORLD
Definition: Utils.cpp:123
double getInnerTolerance() const
Definition: SolverOptions.cpp:612
void setTolerance(double rtol)
Definition: SolverOptions.cpp:581
virtual bool supportsContactElements() const =0
Data cosh() const
Return the cosh of each data point of this Data object.
Definition: Data.cpp:2099
virtual escript::Data getNormal() const =0
Return boundary normals. The function space is chosen appropriately.
Definition: SolverOptions.h:150
boost::python::list getListOfTags() const
Returns a list of the tags used in this function space.
Definition: FunctionSpace.cpp:293
const boost::python::tuple eigenvalues_and_eigenvectors(const DataTypes::real_t tol=1.e-12) const
Return the eigenvalues and corresponding eigenvcetors of the symmetric part at each data point of thi...
Definition: Data.cpp:2957
Data trace(int axis_offset) const
Return the trace of a matrix.
Definition: Data.cpp:2851
void setTaggedValue(int tagKey, const boost::python::object &value)
Assign the given value to the tag. Implicitly converts this object to type DataTagged if it is consta...
Definition: Data.cpp:3597
const char * getName(int key) const
Definition: SolverOptions.cpp:117
Definition: SolverOptions.h:119
const boost::python::tuple minGlobalDataPoint() const
Definition: Data.cpp:2983
FunctionSpace reducedFunction(const AbstractDomain &domain)
Return a function FunctionSpace with reduced integration order.
Definition: FunctionSpaceFactory.cpp:59
Data copySelf() const
Return a pointer to a deep copy of this object.
Definition: Data.cpp:773
Definition: SolverOptions.h:118
bool isEmpty() const
Definition: Data.cpp:1135
void replaceInfPython(boost::python::object obj)
Definition: Data.cpp:2468
int getNumDataPoints() const
Return the number of data points.
Definition: Data.h:509
EscriptDataset()
Default constructor.
Definition: EscriptDataset.cpp:74
void runJobs()
Definition: SplitWorld.cpp:122
bool setExternalDomain(const DomainChunks &domain)
Definition: EscriptDataset.cpp:909
Definition: EscriptDataset.h:36
void dump(const std::string fileName) const
dumps the object into a netCDF file
Definition: Data.cpp:5420
void setPackage(int package)
Definition: SolverOptions.cpp:429
virtual int getDim() const =0
Returns the spatial dimension of the domain.
const boost::python::object getValueOfGlobalDataPointAsTuple(int procNo, int dataPointNo)
Return a data point across all processors as a python tuple.
Definition: Data.cpp:1702
int getIterMax() const
Definition: SolverOptions.cpp:564
Definition: SolverOptions.h:135
Definition: SolverOptions.h:139
FunctionSpace functionOnContactOne(const AbstractDomain &domain)
Return a FunctionSpace on right side of contact.
Definition: FunctionSpaceFactory.cpp:89
void setItemD(const boost::python::object &key, const Data &value)
Copies slice from value into this Data object.
Definition: Data.cpp:3501
Definition: SolverOptions.h:129
void setInnerToleranceAdaption(bool adaption)
Definition: SolverOptions.cpp:747
const boost::python::tuple maxGlobalDataPoint() const
Definition: Data.cpp:3086
virtual void setNewX(const escript::Data &arg)
assigns new location to the domain
Definition: AbstractContinuousDomain.cpp:209
Data Tensor3(double value, const FunctionSpace &what, bool expanded)
Return a Data object containing tensor3 datapoints. ie: rank 3 data-points.
Definition: DataFactory.cpp:158
void copyWithMask(const Data &other, const Data &mask)
Copy other Data object into this Data object where mask is positive.
Definition: Data.cpp:854
virtual int getSystemMatrixTypeId(const boost::python::object &options) const
return the identifier of the matrix type to be used for the global stiffness matrix when particular s...
Definition: AbstractContinuousDomain.cpp:136
SolverOptions getReordering() const
Definition: SolverOptions.cpp:507
FunctionSpace reducedFunctionOnContactOne(const AbstractDomain &domain)
Return a FunctionSpace on right side of contact with reduced integration order.
Definition: FunctionSpaceFactory.cpp:95
bool isSymmetric() const
Definition: SolverOptions.cpp:663
boost::python::object __rsub__(const boost::python::object &right)
wrapper for python reverse subtract operation
Definition: Data.cpp:5963
Data transpose(int axis_offset) const
Transpose each data point of this Data object around the given axis.
Definition: Data.cpp:2906
Data delay()
produce a delayed evaluation version of this Data.
Definition: Data.cpp:789
Definition: SolverOptions.h:160
void setEscriptParamInt(const std::string &name, int value)
Set the value of a named parameter. See listEscriptParams() for available parameters.
Definition: EscriptParams.h:79
Definition: SolverOptions.h:117
FunctionSpace getTestDomainFunctionSpace(int dpps, DataTypes::dim_t samples, int dpsize)
Definition: TestDomain.cpp:349
SolverOptions getPreconditioner() const
Definition: SolverOptions.cpp:355
void releaseUnusedMemory()
releases unused memory in the memory manager.
Definition: DataVectorTaipan.cpp:39
virtual escript::Data getSize() const =0
Returns the local size of samples. The function space is chosen appropriately.
Data hermitian() const
Return the hermitian part of a matrix which is half the matrix plus its adjoint.
Definition: Data.cpp:2784
virtual MPI_Comm getMPIComm() const =0
get the communicator for this domain. Returns an integer on non-MPI builds Routine must be implemente...
Data load(const std::string fileName, const AbstractDomain &domain)
reads Data on domain from file in netCDF format
Definition: DataFactory.cpp:494
Data solve(Data &u0, Data &source, const double dt, boost::python::object &options)
returns the solution u for a time step dt>0 with initial value u0 at time t=0
Definition: AbstractTransportProblem.cpp:60
Data whereNegative() const
Return a Data with a 1 for -ive values and a 0 for +ive or 0 values.
Definition: Data.cpp:1267
Data minval() const
Return the minimum value of each data point of this Data object.
Definition: Data.cpp:2659
int getSubWorldID()
Definition: SplitWorld.cpp:431
void setHermitianOff()
Definition: SolverOptions.cpp:696
virtual double getSafeTimeStepSize() const
Definition: AbstractTransportProblem.cpp:124
const boost::python::object getValueOfDataPointAsTuple(int dataPointNo)
Return the value of a data point as a python tuple.
Definition: Data.cpp:1523
void setInnerIterMax(int iter_max)
Definition: SolverOptions.cpp:545
void setInnerToleranceAdaptionOn()
Definition: SolverOptions.cpp:737
Definition: SolverOptions.h:107
Definition: SolverOptions.h:131
Data log() const
Return the natural log of each data point of this Data object.
Definition: Data.cpp:2153
Data truedivD(const Data &right)
Newer style division operator for python.
Definition: Data.h:2166
void printParallelThreadCnt()
print a message about how many MPI CPUs and OpenMP threads we're using
Definition: Utils.cpp:135
Data rpowO(const boost::python::object &left) const
Return the given power of each data point of this boost python object.
Definition: Data.cpp:3326
boost::python::object __div__(const boost::python::object &right)
wrapper for python divide operation
Definition: Data.cpp:6012
void updateSampleDistribution(VarInfo &vi)
Definition: EscriptDataset.cpp:1029
void setVerbosityOn()
Definition: SolverOptions.cpp:714
Definition: SolverOptions.h:124
void setSolverMethod(int method)
Definition: SolverOptions.cpp:360
bool canInterpolate(FunctionSpace src, FunctionSpace dest)
Definition: FunctionSpace.cpp:32
#define WEIPA_DLL_API
Definition: weipa.h:36
Definition: AbstractContinuousDomain.cpp:22
Data matrixInverse() const
return inverse of matricies.
Definition: Data.cpp:3298
DataTypes::real_t inf()
Return the minimum value of this Data object.
Definition: Data.cpp:2306
StringVec meshLabels
Definition: EscriptDataset.h:175
Data wherePositive() const
Return a Data with a 1 for +ive values and a 0 for 0 or -ive values.
Definition: Data.cpp:1256
double time
Definition: EscriptDataset.h:173
bool isExpanded() const
Return true if this Data is expanded.
Definition: Data.cpp:1114
Definition: SolverOptions.h:121
void setPreconditioner(int preconditioner)
Definition: SolverOptions.cpp:333
Definition: SolverOptions.h:123
bool hasInf()
Definition: Data.cpp:2437
void setSymmetry(bool symmetry)
Definition: SolverOptions.cpp:678
Data Tensor(double value, const FunctionSpace &what, bool expanded)
Return a Data object containing tensor datapoints. ie: rank 2 data-points.
Definition: DataFactory.cpp:120
Data nonuniformslope(boost::python::object in, boost::python::object out, bool check_boundaries)
Definition: Data.cpp:5311
int getTagNumber(int dpno)
Return the tag number associated with the given data-point.
Definition: Data.cpp:3676
boost::python::object __sub__(const boost::python::object &right)
wrapper for python subtract operation
Definition: Data.cpp:5939
Data antisymmetric() const
Return the antisymmetric part of a matrix which is half the matrix minus its transpose.
Definition: Data.cpp:2748
void saveVTK(const std::string fileName)
Saves the dataset in the VTK XML file format.
Definition: EscriptDataset.cpp:486
void setTagsByString(const std::string &name, const escript::Data &mask) const
Definition: FunctionSpace.cpp:264
Definition: SolverOptions.h:140
Definition: SolverOptions.h:104
int getMPIWorldMax(int val)
returns the maximum value of an integer over all processors within MPI_COMM_WORLD
Definition: Utils.cpp:111
int mpiSize
Definition: EscriptDataset.h:179
bool hasFeature(const std::string &name)
returns true if escript was compiled with the feature name, false otherwise.
Definition: EscriptParams.h:106
Definition: SolverOptions.h:144
std::vector< std::string > StringVec
Definition: weipa.h:59
Data asinh() const
Return the asinh of each data point of this Data object.
Definition: Data.cpp:2125
void setDropTolerance(double drop_tol)
Definition: SolverOptions.cpp:617
boost::python::object getVarPyList()
Definition: SplitWorld.cpp:315
int getDim() const
Return the number of spatial dimensions of the underlying domain.
Definition: FunctionSpace.h:186
std::string toString() const
Returns a text description of the function space.
Definition: FunctionSpace.cpp:116
virtual void MPIBarrier() const =0
If compiled for MPI then execute an MPI_Barrier, else do nothing.
bool isLazy() const
Return true if this Data is lazy.
Definition: Data.cpp:1149
int isEmpty() const
Definition: AbstractTransportProblem.cpp:54
FunctionSpace function(const AbstractDomain &domain)
Return a function FunctionSpace.
Definition: FunctionSpaceFactory.cpp:53
Data Vector(double value, const FunctionSpace &what, bool expanded)
Return a Data object containing vector data-points. ie: rank 1 data-points.
Definition: DataFactory.cpp:96
Data whereNonZero(DataTypes::real_t tol=0.0) const
Return a Data with a 0 for 0 values and a 1 for +ive or -ive values.
Definition: Data.cpp:1307
std::vector< VarInfo > VarVector
Definition: EscriptDataset.h:44
void setTupleForGlobalDataPoint(int id, int proc, boost::python::object)
Set the value of a global data point.
Definition: Data.cpp:1570
void setVerbosity(bool verbose)
Definition: SolverOptions.cpp:724
virtual int getMPISize() const =0
return the number of processors used for this domain
Represents an escript dataset including a domain and data variables for one timestep.
Definition: EscriptDataset.h:61
void setValueOfDataPointToArray(int dataPointNo, const boost::python::object &)
sets the values of a data-point from a array-like object on this process
Definition: Data.cpp:1615
FunctionSpace functionOnBoundary(const AbstractDomain &domain)
Return a function on boundary FunctionSpace.
Definition: FunctionSpaceFactory.cpp:65
void setItemO(const boost::python::object &key, const boost::python::object &value)
Definition: Data.cpp:3493
Data condEval(escript::Data &mask, escript::Data &trueval, escript::Data &falseval)
Definition: Data.cpp:5651
void setODESolver(int solver)
Definition: SolverOptions.cpp:813
int MPI_Comm
Definition: EsysMPI.h:40
Data log10() const
Return the log to base 10 of each data point of this Data object.
Definition: Data.cpp:2146
Data pos() const
Return the identity of each data point of this Data object. Simply returns this object unmodified.
Definition: Data.cpp:2182
boost::python::object __rdiv__(const boost::python::object &right)
wrapper for python reverse divide operation
Definition: Data.cpp:6036
void setHermitian(bool hermitian)
Definition: SolverOptions.cpp:701
void replaceNaNPython(boost::python::object obj)
replaces all NaN values with value
Definition: Data.cpp:2423
Data VectorFromObj(bp::object o, const FunctionSpace &what, bool expanded)
Definition: DataFactory.cpp:102
Data imag() const
Definition: Data.cpp:2006
bool hasNaN()
returns return true if data contains NaN.
Definition: Data.cpp:2392
void setAcceptanceConvergenceFailureOff()
Definition: SolverOptions.cpp:765
virtual std::pair< int, DataTypes::dim_t > getDataShape(int functionSpaceCode) const
Return the number of data points per sample, and the number of samples as a pair.
Definition: AbstractContinuousDomain.cpp:203
void setInnerTolerance(double rtol)
Definition: SolverOptions.cpp:605
std::string varName
Definition: EscriptDataset.h:37
Base class for all escript domains.
Definition: AbstractDomain.h:56
Data Tensor3FromObj(bp::object o, const FunctionSpace &what, bool expanded)
Definition: DataFactory.cpp:171
virtual ATP_ptr newTransportProblem(const int blocksize, const escript::FunctionSpace &functionspace, const int type) const
creates a TransportProblemAdapter
Definition: AbstractContinuousDomain.cpp:188
bool setDomain(const escript::AbstractDomain *domain)
Sets the escript domain for this dataset. This method can only be called once and should be the first...
Definition: EscriptDataset.cpp:110
const char * MESH_VARS
Definition: EscriptDataset.cpp:69
Definition: SolverOptions.h:108
void copy(const Data &other)
Make this object a deep copy of "other".
Definition: Data.cpp:780
int getMPIRankWorld()
returns the MPI processor number within MPI_COMM_WORLD
Definition: Utils.cpp:102
virtual std::string showTagNames() const =0
Returns all tag names in a single string sperated by commas.
boost::python::object raw_addJobPerWorld(boost::python::tuple t, boost::python::dict kwargs)
Definition: SplitWorld.cpp:490
bool loadData(const std::string filePattern, const std::string name, const std::string units)
Definition: EscriptDataset.cpp:945
boost::python::object raw_addVariable(boost::python::tuple t, boost::python::dict kwargs)
Definition: SplitWorld.cpp:511
boost::python::object getLocalObjectVariable(const std::string &name)
Definition: SplitWorld.cpp:34
Data neg() const
Return the negation of each data point of this Data object.
Definition: Data.cpp:2175
Definition: SolverOptions.h:146
Data nonuniforminterp(boost::python::object in, boost::python::object out, bool check_boundaries)
Definition: Data.cpp:5235
Definition: SolverOptions.h:113
void setTrilinosParameter(const std::string &name, const boost::python::object &value)
Definition: SolverOptions.cpp:832
boost::shared_ptr< DomainChunk > DomainChunk_ptr
Definition: weipa.h:64
bool adaptInnerTolerance() const
Definition: SolverOptions.cpp:732
void setValueOfDataPointToPyObject(int dataPointNo, const boost::python::object &py_object)
sets the values of a data-point from a python object on this process
Definition: Data.cpp:1562
Data acos() const
Return the acos of each data point of this Data object.
Definition: Data.cpp:2077
virtual void Print_Mesh_Info(const bool full=false) const
Definition: AbstractContinuousDomain.cpp:215
Definition: SolverOptions.h:111
Definition: SolverOptions.h:143
void expand()
Whatever the current Data type make this into a DataExpanded.
Definition: Data.cpp:1183
std::vector< int > IntVec
Definition: weipa.h:58
void setLocalPreconditionerOn()
Definition: SolverOptions.cpp:783
Definition: SolverOptions.h:134
~EscriptDataset()
Destructor.
Definition: EscriptDataset.cpp:103
bool wantsMeshVars
Definition: EscriptDataset.h:176
double getDropStorage() const
Definition: SolverOptions.cpp:636
Definition: SolverOptions.h:114
void setNumberOfThreads(int num_threads)
set the number of threads
Definition: Utils.cpp:77
void setLocalPreconditionerOff()
Definition: SolverOptions.cpp:788
Data randomData(const boost::python::tuple &shape, const FunctionSpace &what, long seed, const boost::python::tuple &filter)
Create a new Expanded Data object filled with pseudo-random data.
Definition: Data.cpp:5875
Data conjugate() const
For complex values return the conjugate values. For non-complex data return a copy.
Definition: Data.cpp:1967
Definition: SolverOptions.h:136
Definition: SolverOptions.h:133
A class that provides functionality to read an escript data object from a dump file or an escript::Da...
Definition: DataVar.h:34
void setProtection()
switches on update protection
Definition: Data.cpp:1169
Data Tensor4(double value, const FunctionSpace &what, bool expanded)
Return a Data object containing tensor4 datapoints. ie: rank 4 data-points.
Definition: DataFactory.cpp:196
Data asin() const
Return the asin of each data point of this Data object.
Definition: Data.cpp:2070
boost::python::dict getTrilinosParameters() const
Definition: SolverOptions.cpp:840
Data sqrt() const
Return the square root of each data point of this Data object.
Definition: Data.cpp:2201
boost::python::object raw_buildDomains(boost::python::tuple t, boost::python::dict kwargs)
Definition: SplitWorld.cpp:454
void setMeshUnits(const std::string x, const std::string y, const std::string z="")
Sets units for the mesh axes.
Definition: EscriptDataset.cpp:569
bool isProtected() const
Returns true, if the data object is protected against update.
Definition: Data.cpp:1175
SolverOptions getODESolver() const
Definition: SolverOptions.cpp:827
void clearVariable(std::string name)
Definition: SplitWorld.cpp:304
virtual int getMPIRank() const =0
return the number MPI rank of this processor
double getMachinePrecision()
returns the machine precision
Definition: Utils.cpp:366
Data eigenvalues() const
Return the eigenvalues of the symmetric part at each data point of this Data object in increasing val...
Definition: Data.cpp:2930
int getTagFromDataPointNo(DataTypes::dim_t dataPointNo) const
Returns the tag associated with the given data-point number.
Definition: FunctionSpace.cpp:148
void convertMeshVariables()
Definition: EscriptDataset.cpp:1002
Definition: SolverOptions.h:120
int getRestart() const
Definition: SolverOptions.cpp:520
FunctionSpace solution(const AbstractDomain &domain)
Return a FunctionSpace.
Definition: FunctionSpaceFactory.cpp:101
virtual void addPDEToTransportProblem(AbstractTransportProblem &tp, escript::Data &source, const escript::Data &M, const escript::Data &A, const escript::Data &B, const escript::Data &C, const escript::Data &D, const escript::Data &X, const escript::Data &Y, const escript::Data &d, const escript::Data &y, const escript::Data &d_contact, const escript::Data &y_contact, const escript::Data &d_dirac, const escript::Data &y_dirac) const
adds a PDE onto a transport problem
Definition: AbstractContinuousDomain.cpp:165
void putSiloMultiVar(DBfile *dbfile, const VarInfo &vi, bool useMeshFile=false)
Definition: EscriptDataset.cpp:1121
Represents a full Finley or Dudley domain including nodes and elements.
Definition: weipa/src/FinleyDomain.h:40
Definition: SolverOptions.h:154
unsigned int getDataPointRank() const
Return the rank of the point data.
Definition: Data.h:498