19 #ifndef __ESCRIPT_DATA_H__ 20 #define __ESCRIPT_DATA_H__ 33 #include <boost/python/object.hpp> 34 #include <boost/python/tuple.hpp> 35 #include <boost/math/special_functions/bessel.hpp> 37 #ifndef ESCRIPT_MAX_DATA_RANK 38 #define ESCRIPT_MAX_DATA_RANK 4 169 Data(
const boost::python::object& value,
192 Data(boost::python::object value,
193 boost::python::object par1=boost::python::object(),
194 boost::python::object par2=boost::python::object(),
195 boost::python::object par3=boost::python::object());
269 const boost::python::object
300 const boost::python::object
454 return m_data->getFunctionSpace();
512 return m_data->getNumSamples();
523 return m_data->getNumDPPSample();
559 return m_data->getNoValues();
568 dump(
const std::string fileName)
const;
576 const boost::python::object
633 return m_data->getSampleDataByTag(tag, dummy);
640 return m_data->getSampleDataByTag(tag, dummy);
673 return m_data->getPointOffset(sampleNo,dataPointNo);
684 return m_data->getShape();
691 const boost::python::tuple
716 return m_data->getNumSamples()==0;
729 const boost::python::object& value);
742 const boost::python::object& value);
824 nonuniforminterp(boost::python::object in, boost::python::object out,
bool check_boundaries);
827 nonuniformslope(boost::python::object in, boost::python::object out,
bool check_boundaries);
845 boost::python::object
853 boost::python::object
1008 const boost::python::tuple
1018 const boost::python::tuple
1071 trace(
int axis_offset)
const;
1099 const boost::python::tuple
1108 swapaxes(
const int axis0,
const int axis1)
const;
1322 powO(
const boost::python::object& right)
const;
1333 rpowO(
const boost::python::object& left)
const;
1395 boost::python::object
__add__(
const boost::python::object& right);
1402 boost::python::object
__sub__(
const boost::python::object& right);
1408 boost::python::object
__rsub__(
const boost::python::object& right);
1414 boost::python::object
__mul__(
const boost::python::object& right);
1420 boost::python::object
__div__(
const boost::python::object& right);
1426 boost::python::object
__rdiv__(
const boost::python::object& right);
1457 getItem(
const boost::python::object& key)
const;
1471 setItemD(
const boost::python::object& key,
1475 setItemO(
const boost::python::object& key,
1476 const boost::python::object& value);
1485 template <
class UnaryFunction>
1615 template <
typename S>
1619 template <
class BinaryOp>
1636 template<
typename Scalar>
1637 boost::python::object
1674 template <
class BinaryFunction>
1687 template <
class BinaryFunction>
1758 ESYS_ASSERT(dr!=0,
"error casting to DataReady.");
1766 ESYS_ASSERT(dr!=0,
"error casting to DataReady.");
1778 ESYS_ASSERT(dr.get()!=0,
"error casting to DataReady.");
1786 ESYS_ASSERT(dr.get()!=0,
"error casting to DataReady.");
1806 #ifdef SLOWSHARECHECK 1807 return m_data->isShared();
1809 return !m_data.unique();
1818 if (omp_in_parallel())
1820 throw DataException(
"Please do not call forceResolve() in a parallel region.");
1834 if (omp_in_parallel())
1836 throw DataException(
"Programming error. Please do not run exclusiveWrite() in multi-threaded sections.");
1846 m_data->exclusivewritecalled=
true;
1857 std::ostringstream oss;
1858 oss <<
"Programming error. ExclusiveWrite required - please call requireWrite() isLazy=" <<
isLazy() <<
" isShared()=" <<
isShared();
1876 #ifdef IKNOWWHATIMDOING 1879 template <
typename S>
1881 friend Data randomData(
const boost::python::tuple& shape,
const FunctionSpace& what,
long seed,
const boost::python::tuple& filter);
1886 #ifdef IKNOWWHATIMDOING 1901 long seed,
const boost::python::tuple& filter);
1927 throw DataException(
"Error, attempt to acquire RW access to lazy data. Please call requireWrite() first.");
1930 if (!
getReady()->exclusivewritecalled)
1932 throw DataException(
"Error, call to Data::getSampleDataRW without a preceeding call to requireWrite/exclusiveWrite.");
1944 throw DataException(
"Error, attempt to acquire RW access to lazy data. Please call requireWrite() first.");
1947 if (!
getReady()->exclusivewritecalled)
1949 throw DataException(
"Error, call to Data::getSampleDataRW without a preceeding call to requireWrite/exclusiveWrite.");
1965 return &((*res)[offset]);
1977 throw DataException(
"Programming error: complex lazy objects are not supported.");
1989 throw DataException(
"Programmer error - getDataRO must not be called on Lazy Data.");
2007 throw DataException(
"Programmer error - getDataRO must not be called on Lazy Data.");
2151 return *
this / right;
2189 template <
class BinaryFunction>
2196 ESYS_ASSERT(leftC!=0,
"Programming error - casting to DataExpanded.");
2204 const auto& vec=data.
getTypedVectorRO(
typename BinaryFunction::first_argument_type(0));
2208 #pragma omp parallel private(local_current_value) 2210 local_current_value=initial_value;
2211 #pragma omp for private(i,j) schedule(static) 2212 for (i=0;i<numSamples;i++) {
2213 for (j=0;j<numDPPSample;j++) {
2218 #pragma omp critical 2219 global_current_value=operation(global_current_value,local_current_value);
2221 return global_current_value;
2224 ESYS_ASSERT(leftC!=0,
"Programming error - casting to DataTagged.");
2229 const auto& vec=data.
getTypedVectorRO(
typename BinaryFunction::first_argument_type(0));
2233 for (std::list<int>::const_iterator i=used.begin();i!=used.end();++i)
2236 DataTagged::DataMapType::const_iterator it=lookup.find(tag);
2237 if ((tag==0) || (it==lookup.end()))
2246 return current_value;
2249 ESYS_ASSERT(leftC!=0,
"Programming error - casting to DataConstant.");
2252 throw DataException(
"Error - Operations (algorithm) not permitted on instances of DataEmpty.");
2254 throw DataException(
"Error - Operations not permitted on instances of DataLazy.");
2256 throw DataException(
"Error - Data encapsulates an unknown type.");
2268 template <
class BinaryFunction>
2274 throw DataException(
"Error - Operations (dp_algorithm) not permitted on instances of DataEmpty.");
2280 ESYS_ASSERT(dataE!=0,
"Programming error - casting data to DataExpanded.");
2281 ESYS_ASSERT(resultE!=0,
"Programming error - casting result to DataExpanded.");
2295 #pragma omp parallel for private(i,j) schedule(static) 2296 for (i=0;i<numSamples;i++) {
2297 for (j=0;j<numDPPSample;j++) {
2308 ESYS_ASSERT(dataT!=0,
"Programming error - casting data to DataTagged.");
2317 for (DataTagged::DataMapType::const_iterator i=lookup.begin(); i!=lookup.end(); i++) {
2327 return Data(resultT);
2333 ESYS_ASSERT(dataC!=0,
"Programming error - casting data to DataConstant.");
2334 ESYS_ASSERT(resultC!=0,
"Programming error - casting result to DataConstant.");
2343 throw DataException(
"Error - Operations not permitted on instances of DataLazy.");
2345 throw DataException(
"Error - Data encapsulates an unknown type.");
2370 #endif // __ESCRIPT_DATA_H__ Definition: FunctionSpace.h:34
std::ostream & operator<<(std::ostream &o, const Data &data)
Output operator.
Definition: Data.cpp:3602
void typeMatchLeft(Data &right) const
Convert the data type of the RHS to match this.
Definition: Data.cpp:3448
DataAbstract_ptr m_data
Definition: Data.h:1750
DataTypes::RealVectorType::size_type getLength() const
Return the number of doubles stored for this Data.
Definition: Data.cpp:1349
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:1360
DataTypes::real_t * getSampleDataByTag(int tag, DataTypes::real_t dummy=0)
Return the sample data for the given tag. If an attempt is made to access data that isn't tagged an e...
Definition: Data.h:631
Data sqrt() const
Return the square root of each data point of this Data object.
Definition: Data.cpp:2162
void setItemO(const boost::python::object &key, const boost::python::object &value)
Definition: Data.cpp:3409
Data antihermitian() const
Return the anti-hermitian part of a matrix which is half the matrix minus its hermitian.
Definition: Data.cpp:2727
Data interpolateFromTable1D(const WrappedArray &table, DataTypes::real_t Amin, DataTypes::real_t Astep, DataTypes::real_t undef, bool check_boundaries)
Definition: Data.cpp:4684
boost::python::object __add__(const boost::python::object &right)
wrapper for python add operation
Definition: Data.cpp:5831
const_DataReady_ptr getReadyPtr() const
Definition: Data.h:1783
void checkExclusiveWrite()
checks if caller can have exclusive write to the object
Definition: Data.h:1853
boost::python::object __rsub__(const boost::python::object &right)
wrapper for python reverse subtract operation
Definition: Data.cpp:5879
std::string toString() const
Write the data as a string. For large amounts of data, a summary is printed.
Definition: Data.cpp:4562
void setToZero()
set all values to zero
Definition: Data.cpp:797
DataTypes::real_t rpow(DataTypes::real_t x, DataTypes::real_t y)
Definition: Data.h:2023
ElementType & reference
Definition: DataVectorAlt.h:50
DataTypes::real_t inf()
Return the minimum value of this Data object.
Definition: Data.cpp:2267
void replaceNaN(DataTypes::real_t value)
replaces all NaN values with value
Definition: Data.cpp:2364
unsigned int getDataPointRank() const
Return the rank of the point data.
Definition: Data.h:489
int MPI_Op
Definition: EsysMPI.h:43
Data delay()
produce a delayed evaluation version of this Data.
Definition: Data.cpp:768
Data & operator*=(const Data &right)
Overloaded operator *=.
Definition: Data.cpp:3156
void tag()
If possible convert this Data to DataTagged. This will only allow Constant data to be converted to ta...
Definition: Data.cpp:1186
bool isDataPointShapeEqual(int rank, const int *dimensions) const
Returns true if the shape matches the vector (dimensions[0],..., dimensions[rank-1]). DataEmpty always returns true.
Definition: Data.h:544
void set_m_data(DataAbstract_ptr p)
Modify the data abstract hosted by this Data object For internal use only. Passing a pointer to null ...
Definition: Data.cpp:612
virtual DataTypes::RealVectorType & getTypedVectorRW(DataTypes::real_t dummy)
These versions use the type system rather than method name to determine return type.
Definition: DataExpanded.cpp:1438
boost::python::object integrateToTuple()
Calculate the integral over the function space domain as a python tuple.
Definition: Data.cpp:1748
Data getSlice(const DataTypes::RegionType ®ion) const
Return a Data object containing the specified slice of this Data object.
Definition: Data.cpp:3401
DataTypes::real_t supWorker() const
Definition: Data.cpp:2456
boost::shared_ptr< AbstractDomain > Domain_ptr
Definition: AbstractDomain.h:36
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:2121
Definition: DataReady.h:35
DataTypes::vec_size_type size_type
Definition: DataVectorAlt.h:49
std::vector< std::pair< int, int > > RegionType
Definition: DataTypes.h:43
int getDataPointSize() const
Return the size of the data point. It is the product of the data point shape dimensions.
Definition: Data.cpp:1342
Definition: AbstractContinuousDomain.cpp:22
Data imag() const
Definition: Data.cpp:1967
Data atan() const
Return the atan of each data point of this Data object.
Definition: Data.cpp:2046
Data log10() const
Return the log to base 10 of each data point of this Data object.
Definition: Data.cpp:2107
Data pos() const
Return the identity of each data point of this Data object. Simply returns this object unmodified...
Definition: Data.cpp:2143
boost::shared_ptr< const DataReady > const_DataReady_ptr
Definition: DataAbstract.h:59
void setItemD(const boost::python::object &key, const Data &value)
Copies slice from value into this Data object.
Definition: Data.cpp:3417
Data operator-(const Data &left, const Data &right)
Operator- Takes two Data objects.
Definition: Data.cpp:3278
MPI_Comm get_MPIComm(void) const
return the MPI rank number of the local data MPI_COMM_WORLD is assumed and returned.
Definition: Data.cpp:5382
DataTypes::real_t sup_const() const
Definition: Data.cpp:2214
Data besselSecondKind(int order)
Return the Bessel function of the second kind for each data point of this Data object.
Definition: Data.cpp:1828
void print(void)
print the data values to stdout. Used for debugging
Definition: Data.cpp:5306
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:4656
void typeMatchRight(const Data &right)
Convert the data type of this to match the RHS.
Definition: Data.cpp:3468
void setSlice(const Data &value, const DataTypes::RegionType ®ion)
Copy the specified slice from the given value into this Data object.
Definition: Data.cpp:3433
DataConstant stores a single data point which represents the entire function space.
Definition: DataConstant.h:37
static const ShapeType scalarShape
Use this instead of creating empty shape objects for scalars.
Definition: DataTypes.h:46
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:1541
DataTypes::RealVectorType::const_reference getDataPointRO(int sampleNo, int dataPointNo)
Return a reference into the DataVector which points to the specified data point.
Definition: Data.cpp:4632
bool isLazy() const
Return true if this Data is lazy.
Definition: Data.cpp:1128
Data C_TensorUnaryOperation(Data const &arg_0, escript::ES_optype operation, DataTypes::real_t tol=0)
Definition: Data.cpp:5998
DataAbstract_ptr borrowDataPtr(void) const
Definition: Data.cpp:4547
Data real() const
Definition: Data.cpp:1948
virtual const DataTypes::RealVectorType & getTypedVectorRO(DataTypes::real_t dummy) const
Definition: DataConstant.cpp:786
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:3496
Give a short description of what DataExpanded does.
Definition: DataExpanded.h:43
virtual const DataTypes::RealVectorType & getTypedVectorRO(DataTypes::real_t dummy) const
Definition: DataTagged.cpp:1591
Data sinh() const
Return the sinh of each data point of this Data object.
Definition: Data.cpp:2053
Data atanh() const
Return the atanh of each data point of this Data object.
Definition: Data.cpp:2100
DataReady_ptr getReadyPtr()
Definition: Data.h:1775
void TensorSelfUpdateBinaryOperation(const Data &right, escript::ES_optype operation)
Definition: Data.cpp:6450
DataTypes::RealVectorType::reference getDataByTagRW(int tag, DataTypes::RealVectorType::size_type i, DataTypes::real_t dummy=0)
getDataByTag
Definition: DataTagged.cpp:932
const_Domain_ptr getDomain() const
Returns the function space domain.
Definition: FunctionSpace.cpp:102
Data rpowO(const boost::python::object &left) const
Return the given power of each data point of this boost python object.
Definition: Data.cpp:3242
friend 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.
DataTypes::real_t lazyAlgWorker(DataTypes::real_t init)
Data interpolateFromTable1DP(boost::python::object table, DataTypes::real_t Amin, DataTypes::real_t Astep, DataTypes::real_t undef, bool check_boundaries)
Definition: Data.cpp:4675
std::list< int > getListOfTagsSTL() const
Returns an stl list of the tags used in this function space.
Definition: FunctionSpace.cpp:286
(Testing use only) Provides a domain to wrap a collection of values.
Definition: TestDomain.h:41
void setTaggedValueFromCPP(int tagKey, const DataTypes::ShapeType &pointshape, const DataTypes::RealVectorType &value, int dataOffset=0)
Assign the given value to the tag. Implicitly converts this object to type DataTagged if it is consta...
Definition: Data.cpp:3553
Data nonuniformslope(boost::python::object in, boost::python::object out, bool check_boundaries)
Definition: Data.cpp:5227
std::vector< int > ShapeType
The shape of a single datapoint.
Definition: DataTypes.h:42
boost::python::object __sub__(const boost::python::object &right)
wrapper for python subtract operation
Definition: Data.cpp:5855
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:4666
Data acos() const
Return the acos of each data point of this Data object.
Definition: Data.cpp:2038
Simulates a full dataset accessible via sampleNo and dataPointNo.
Definition: DataTagged.h:44
bool probeInterpolation(const FunctionSpace &functionspace) const
Returns true if this can be interpolated to functionspace.
Definition: Data.cpp:1299
boost::shared_ptr< DataAbstract > DataAbstract_ptr
Definition: DataAbstract.h:51
Data phase() const
Return the phase/arg/angular-part of complex values.
Definition: Data.cpp:1987
DataTypes::real_t * getSampleDataRW(DataTypes::RealVectorType::size_type sampleNo, DataTypes::real_t dummy=0)
Return the sample data for the given sample number.
Definition: DataReady.h:181
Data getItem(const boost::python::object &key) const
Returns a slice from this Data object.
Definition: Data.cpp:3387
void initialise(const DataTypes::RealVectorType &value, const DataTypes::ShapeType &shape, const FunctionSpace &what, bool expanded)
Construct a Data object of the appropriate type.
Definition: Data.cpp:641
Data maxval() const
Return the maximum value of each data point of this Data object.
Definition: Data.cpp:2566
void setValueOfDataPointC(int dataPointNo, const DataTypes::cplx_t)
Definition: Data.cpp:1658
Data cos() const
Return the cos of each data point of this Data object.
Definition: Data.cpp:2017
const DataTypes::real_t * getDataRO(DataTypes::real_t dummy=0) const
Return a pointer to the beginning of the underlying data.
Definition: Data.h:1985
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:3609
DataTypes::RealVectorType::const_reference getDataAtOffsetRO(DataTypes::RealVectorType::size_type i, DataTypes::real_t dummy)
Return a pointer to the beginning of the datapoint at the specified offset. TODO Eventually these sho...
Definition: Data.cpp:4609
bool m_protected
Definition: Data.h:1744
Data cosh() const
Return the cosh of each data point of this Data object.
Definition: Data.cpp:2060
void setValueOfDataPoint(int dataPointNo, const DataTypes::real_t)
sets the values of a data-point on this process
Definition: Data.cpp:1637
Data interpolate(const FunctionSpace &functionspace) const
Interpolates this onto the given functionspace and returns the result as a Data object.
Definition: Data.cpp:1293
Data matrixInverse() const
return inverse of matricies.
Definition: Data.cpp:3214
Data & operator/=(const Data &right)
Overloaded operator /=.
Definition: Data.cpp:3183
DataReady * getReady()
Definition: Data.h:1763
Data copySelf() const
Return a pointer to a deep copy of this object.
Definition: Data.cpp:752
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:1212
bool isEmpty() const
Return true if this Data holds an instance of DataEmpty. This is _not_ the same as asking if the obj...
Definition: Data.cpp:1114
Data tanh() const
Return the tanh of each data point of this Data object.
Definition: Data.cpp:2067
void calc_minGlobalDataPoint(int &ProcNo, int &DataPointNo) const
Definition: Data.cpp:2916
void operandCheck(const Data &right) const
Check *this and the right operand are compatible. Throws an exception if they aren't.
Definition: Data.h:1664
DataAbstract * borrowData(void) const
return the object produced by the factory, which is a DataConstant or DataExpanded TODO Ownership of ...
Definition: Data.cpp:4540
Data erf() const
Return the error function erf of each data point of this Data object.
Definition: Data.cpp:2075
const DataTypes::ShapeType & getShape() const
Return the shape information for the point data.
Definition: DataAbstract.h:641
int getNumSamples() const
Return the number of samples.
Definition: DataAbstract.h:622
int getNumDPPSample() const
Return the number of data points per sample.
Definition: DataAbstract.h:611
DataTypes::real_t infWorker() const
Definition: Data.cpp:2499
int getNoValues() const
Return the number of values in the shape for this object.
Definition: Data.h:557
Describes binary operations performed on DataVector.
virtual DataTypes::RealVectorType::size_type getPointOffset(int sampleNo, int dataPointNo) const
Return the offset for the given given data point. This returns the offset in bytes for the given poin...
Definition: DataExpanded.cpp:534
boost::python::object integrateWorker() const
Data rtruedivO(const boost::python::object &left)
Newer style division operator for python.
Definition: Data.h:2174
const boost::python::tuple maxGlobalDataPoint() const
Return the (sample number, data-point number) of the data point with the minimum component value in t...
Definition: Data.cpp:3002
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:1958
std::map< int, int > DataMapType
Definition: DataTagged.h:57
Data powD(const Data &right) const
Return the given power of each data point of this Data object.
Definition: Data.cpp:3256
void replaceNaNPython(boost::python::object obj)
replaces all NaN values with value
Definition: Data.cpp:2384
Data trace(int axis_offset) const
Return the trace of a matrix.
Definition: Data.cpp:2767
DataTypes::RealVectorType::size_type getOffsetForTag(int tag) const
getOffsetForTag
Definition: DataTagged.cpp:910
const boost::python::tuple minGlobalDataPoint() const
Return the (sample number, data-point number) of the data point with the minimum component value in t...
Definition: Data.cpp:2899
Data minval_nonlazy() const
Definition: Data.cpp:2544
int get_MPISize(void) const
return the MPI rank number of the local data MPI_COMM_WORLD is assumed and the result of MPI_Comm_ran...
Definition: Data.cpp:5358
const DataTypes::ShapeType & getDataPointShape() const
Return a reference to the data point shape.
Definition: Data.h:682
const FunctionSpace & getFunctionSpace() const
Return the function space.
Definition: Data.h:452
void delaySelf()
convert the current data into lazy data.
Definition: Data.cpp:779
Data transpose(int axis_offset) const
Transpose each data point of this Data object around the given axis.
Definition: Data.cpp:2822
int get_MPIRank(void) const
return the MPI rank number of the local data MPI_COMM_WORLD is assumed and the result of MPI_Comm_siz...
Definition: Data.cpp:5370
void expand()
Whatever the current Data type make this into a DataExpanded.
Definition: Data.cpp:1162
DataTypes::real_t reduction(BinaryFunction operation, DataTypes::real_t initial_value) const
Perform the specified reduction algorithm on every element of every data point in this Data object ac...
Definition: Data.h:2192
Data operator+(const Data &left, const Data &right)
Operator+ Takes two Data objects.
Definition: Data.cpp:3268
Data represents a collection of datapoints.
Definition: Data.h:63
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:1279
DataTypes::RealVectorType::reference getDataPointRW(int sampleNo, int dataPointNo)
Return a reference into the DataVector which points to the specified data point.
Definition: Data.cpp:4648
Data asinh() const
Return the asinh of each data point of this Data object.
Definition: Data.cpp:2086
void maskWorker(Data &other2, Data &mask2, S sentinel)
Definition: Data.cpp:936
virtual DataTypes::RealVectorType & getTypedVectorRW(DataTypes::real_t dummy)
These versions use the type system rather than method name to determine return type.
Definition: DataConstant.cpp:779
DataTypes::RealVectorType::reference getDataAtOffsetRW(DataTypes::RealVectorType::size_type i, DataTypes::real_t dummy)
Definition: Data.cpp:4601
ES_optype
Definition: ES_optype.h:26
Data asin() const
Return the asin of each data point of this Data object.
Definition: Data.cpp:2031
bool isConstant() const
Return true if this Data is constant.
Definition: Data.cpp:1121
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:1594
bool isComplex() const
True if components of this data are stored as complex.
Definition: Data.cpp:1142
bool isShared() const
Definition: Data.h:1804
bool numSamplesEqual(int numDataPointsPerSample, int numSamples) const
Returns true if the number of data points per sample and the number of samples match the respective a...
Definition: Data.h:532
const DataTypes::RealVectorType * resolveSample(int sampleNo, size_t &roffset) const
Compute the value of the expression for the given sample.
Definition: DataLazy.cpp:2334
DataReady_ptr borrowReadyPtr(void) const
Definition: Data.cpp:4554
Data acosh() const
Return the acosh of each data point of this Data object.
Definition: Data.cpp:2093
Data wherePositive() const
Return a Data with a 1 for +ive values and a 0 for 0 or -ive values.
Definition: Data.cpp:1235
bool actsExpanded() const
Return true if this Data is expanded or resolves to expanded. That is, if it has a separate value for...
Definition: Data.cpp:1100
Data & operator=(const Data &other)
Definition: Data.cpp:3121
void forceResolve()
Definition: Data.h:1813
DataTypes::cplx_t * getSampleDataByTag(int tag, DataTypes::cplx_t dummy)
Definition: Data.h:638
virtual DataAbstract * deepCopy() const =0
Return a deep copy of the current object.
void requireWrite()
Ensures data is ready for write access. This means that the data will be resolved if lazy and will be...
Definition: Data.cpp:1221
DataTypes::real_t Lsup_const() const
Definition: Data.cpp:2169
int getNumDataPointsPerSample() const
Return the number of data points per sample.
Definition: Data.h:521
const DataReady * getReady() const
Definition: Data.h:1755
DataTypes::RealVectorType::size_type getDefaultOffset() const
Returns the offset in the structure which stores the default value.
Definition: DataTagged.h:703
DataTypes::real_t * getSampleDataRW(DataTypes::RealVectorType::size_type sampleNo, DataTypes::real_t dummy=0)
Return the sample data for the given sample no. Please do not use this unless you NEED to access samp...
Definition: Data.h:1923
int MPI_Comm
Definition: EsysMPI.h:41
Data hermitian() const
Return the hermitian part of a matrix which is half the matrix plus its adjoint.
Definition: Data.cpp:2700
friend Data condEvalWorker(escript::Data &mask, escript::Data &trueval, escript::Data &falseval, S sentinel)
Data oneOver() const
Returns 1./ Data object.
Definition: Data.cpp:1228
virtual DataTypes::RealVectorType & getTypedVectorRW(DataTypes::real_t dummy)
These versions use the type system rather than method name to determine return type.
Definition: DataTagged.cpp:1584
DataTypes::real_t LsupWorker() const
Definition: Data.cpp:2399
Data minval() const
Return the minimum value of each data point of this Data object.
Definition: Data.cpp:2575
Data C_TensorBinaryOperation(Data const &arg_0, Data const &arg_1, ES_optype operation)
Compute a tensor operation with two Data objects.
Definition: Data.cpp:6274
bool m_lazy
Definition: Data.h:1745
DataTypes::real_t sup()
Return the maximum value of this Data object.
Definition: Data.cpp:2228
void setTupleForGlobalDataPoint(int id, int proc, boost::python::object)
Set the value of a global data point.
Definition: Data.cpp:1549
const DataTypes::real_t * getSampleDataRO(DataTypes::RealVectorType::size_type sampleNo, DataTypes::real_t dummy=0) const
Definition: DataReady.h:195
int getNumSamples() const
Return the number of samples.
Definition: Data.h:510
boost::python::object __div__(const boost::python::object &right)
wrapper for python divide operation
Definition: Data.cpp:5928
virtual const DataTypes::RealVectorType & getTypedVectorRO(DataTypes::real_t dummy) const =0
Data abs() const
Return the absolute value of each data point of this Data object.
Definition: Data.cpp:2129
Data grad() const
Definition: Data.cpp:1332
Data dp_algorithm(BinaryFunction operation, DataTypes::real_t initial_value) const
Reduce each data-point in this Data object using the given operation. Return a Data object with the s...
Definition: Data.h:2271
void copy(const Data &other)
Make this object a deep copy of "other".
Definition: Data.cpp:759
Wraps an expression tree of other DataObjects. The data will be evaluated when required.
Definition: DataLazy.h:47
std::complex< real_t > cplx_t
complex data type
Definition: DataTypes.h:53
DataTypes::real_t inf_const() const
Definition: Data.cpp:2253
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:2873
void init_from_data_and_fs(const Data &inData, const FunctionSpace &functionspace)
Definition: Data.cpp:433
int getTagNumber(int dpno)
Return the tag number associated with the given data-point.
Definition: Data.cpp:3592
Domain_ptr getDomainPython() const
Return the function space domain. Internal use only! This gets around some python difficulties by cas...
Definition: FunctionSpace.cpp:108
Definition: DataException.h:26
void complicate()
Definition: Data.cpp:5976
Data nonuniforminterp(boost::python::object in, boost::python::object out, bool check_boundaries)
Definition: Data.cpp:5151
virtual const DataTypes::RealVectorType & getTypedVectorRO(DataTypes::real_t dummy) const
Definition: DataExpanded.cpp:1444
Data maxval_nonlazy() const
Definition: Data.cpp:2555
DataTypes::real_t reductionOpVector(const DataTypes::RealVectorType &left, const DataTypes::ShapeType &leftShape, DataTypes::RealVectorType::size_type offset, BinaryFunction operation, DataTypes::real_t initial_value)
Perform the given data point reduction operation on the data point specified by the given offset into...
Definition: DataVectorOps.h:1393
Data condEval(escript::Data &mask, escript::Data &trueval, escript::Data &falseval)
Definition: Data.cpp:5567
Data powO(const boost::python::object &right) const
Return the given power of each data point of this boost python object.
Definition: Data.cpp:3249
Data whereNonNegative() const
Return a Data with a 1 for +ive or 0 values and a 0 for -ive values.
Definition: Data.cpp:1257
Data exp() const
Return the exponential function of each data point of this Data object.
Definition: Data.cpp:2155
Data operator*(const AbstractSystemMatrix &left, const Data &right)
Definition: AbstractSystemMatrix.cpp:59
void exclusiveWrite()
if another object is sharing out member data make a copy to work with instead. This code should only ...
Definition: Data.h:1831
Data whereNegative() const
Return a Data with a 1 for -ive values and a 0 for +ive or 0 values.
Definition: Data.cpp:1246
bool isTagged() const
Return true if this Data is tagged.
Definition: Data.cpp:1107
Definition: DataAbstract.h:61
bool hasNoSamples() const
Return true if this object contains no samples. This is not the same as isEmpty() ...
Definition: Data.h:714
Data conjugate() const
For complex values return the conjugate values. For non-complex data return a copy.
Definition: Data.cpp:1928
const_Domain_ptr getDomain() const
Return the domain.
Definition: Data.h:464
Data interpolateFromTable2D(const WrappedArray &table, DataTypes::real_t Amin, DataTypes::real_t Astep, DataTypes::real_t undef, Data &B, DataTypes::real_t Bmin, DataTypes::real_t Bstep, bool check_boundaries)
Definition: Data.cpp:4812
const ElementType & const_reference
Definition: DataVectorAlt.h:51
void dump(const std::string fileName) const
dumps the object into a netCDF file
Definition: Data.cpp:5336
void setProtection()
switches on update protection
Definition: Data.cpp:1148
const FunctionSpace & getFunctionSpace() const
Return the function space associated with this Data object.
Definition: DataAbstract.h:634
Data besselFirstKind(int order)
Return the Bessel function of the first kind for each data point of this Data object.
Definition: Data.cpp:1821
boost::shared_ptr< DataReady > DataReady_ptr
Definition: DataAbstract.h:56
boost::python::object __rdiv__(const boost::python::object &right)
wrapper for python reverse divide operation
Definition: Data.cpp:5952
bool hasNaN()
returns return true if data contains NaN.
Definition: Data.cpp:2353
DataTypes::RealVectorType::size_type getDataOffset(int sampleNo, int dataPointNo)
Return the offset for the given sample and point within the sample.
Definition: Data.h:670
Data operator/(const Data &left, const Data &right)
Operator/ Takes two Data objects.
Definition: Data.cpp:3297
bool isExpanded() const
Return true if this Data is expanded.
Definition: Data.cpp:1093
Data whereNonPositive() const
Return a Data with a 1 for -ive or 0 values and a 0 for +ive values.
Definition: Data.cpp:1268
Data()
Default constructor. Creates a DataEmpty object.
Definition: Data.cpp:249
void copyWithMask(const Data &other, const Data &mask)
Copy other Data object into this Data object where mask is positive.
Definition: Data.cpp:833
boost::python::object __mul__(const boost::python::object &right)
wrapper for python multiply operation
Definition: Data.cpp:5904
size_t getNumberOfTaggedValues() const
For tagged Data returns the number of tags with values. For non-tagged data will return 0 (even Data ...
Definition: Data.cpp:5778
Data neg() const
Return the negation of each data point of this Data object.
Definition: Data.cpp:2136
const boost::python::tuple getShapeTuple() const
Return the data point shape as a tuple of integers.
Definition: Data.cpp:727
const DataMapType & getTagLookup() const
getTagLookup
Definition: DataTagged.h:741
Data truedivO(const boost::python::object &right)
Newer style division operator for python.
Definition: Data.h:2161
bool isReady() const
Return true if this data is ready.
Definition: Data.cpp:1135
Data symmetric() const
Return the symmetric part of a matrix which is half the matrix plus its transpose.
Definition: Data.cpp:2641
#define ESYS_ASSERT(a, b)
EsysAssert is a MACRO that will throw an exception if the boolean condition specified is false...
Definition: Assert.h:78
Definition: WrappedArray.h:31
DataTypes::RealVectorType & getExpandedVectorReference(DataTypes::real_t dummy=0)
Ensures that the Data is expanded and returns its underlying vector Does not check for exclusive writ...
Definition: Data.cpp:5760
int getNumDataPoints() const
Return the number of data points.
Definition: Data.h:500
Data antisymmetric() const
Return the antisymmetric part of a matrix which is half the matrix minus its transpose.
Definition: Data.cpp:2664
Data sin() const
Return the sin of each data point of this Data object.
Definition: Data.cpp:2010
~Data()
Destructor.
Definition: Data.cpp:604
void unaryOp2(UnaryFunction operation)
Perform the given unary operation on every element of every data point in this Data object...
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:1305
const boost::python::object getValueOfDataPointAsTuple(int dataPointNo)
Return the value of a data point as a python tuple.
Definition: Data.cpp:1502
#define S(_J_, _I_)
Definition: ShapeFunctions.cpp:121
Data interpolateFromTable3D(const WrappedArray &table, DataTypes::real_t Amin, DataTypes::real_t Astep, DataTypes::real_t undef, Data &B, DataTypes::real_t Bmin, DataTypes::real_t Bstep, Data &C, DataTypes::real_t Cmin, DataTypes::real_t Cstep, bool check_boundaries)
Definition: Data.cpp:4968
bool isProtected() const
Returns true, if the data object is protected against update.
Definition: Data.cpp:1154
Data log() const
Return the natural log of each data point of this Data object.
Definition: Data.cpp:2114
Data truedivD(const Data &right)
Newer style division operator for python.
Definition: Data.h:2149
boost::shared_ptr< const AbstractDomain > const_Domain_ptr
Definition: AbstractDomain.h:39
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:1286
DataTypes::real_t Lsup()
Return the maximum absolute value of this Data object.
Definition: Data.cpp:2179
Data & operator-=(const Data &right)
Overloaded operator -=.
Definition: Data.cpp:3129
boost::python::object integrateToTuple_const() const
Calculate the integral over the function space domain as a python tuple.
Definition: Data.cpp:1734
double real_t
type of all real-valued scalars in escript
Definition: DataTypes.h:50
void calc_maxGlobalDataPoint(int &ProcNo, int &DataPointNo) const
Definition: Data.cpp:3012
Data eigenvalues() const
Return the eigenvalues of the symmetric part at each data point of this Data object in increasing val...
Definition: Data.cpp:2846
Data bessel(int order, DataTypes::real_t(*besselfunc)(int, DataTypes::real_t))
Bessel worker function.
Definition: Data.cpp:1835
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:3513
const boost::python::object getValueOfGlobalDataPointAsTuple(int procNo, int dataPointNo)
Return a data point across all processors as a python tuple.
Definition: Data.cpp:1681
Data swapaxes(const int axis0, const int axis1) const
swaps the components axis0 and axis1
Definition: Data.cpp:2584
Domain_ptr getDomainPython() const
Return the domain. TODO: For internal use only. This should be removed.
Definition: Data.h:478
Data tan() const
Return the tan of each data point of this Data object.
Definition: Data.cpp:2024
Data & operator+=(const Data &right)
Overloaded operator +=.
Definition: Data.cpp:3093