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
999 const boost::python::tuple
1009 const boost::python::tuple
1062 trace(
int axis_offset)
const;
1090 const boost::python::tuple
1099 swapaxes(
const int axis0,
const int axis1)
const;
1313 powO(
const boost::python::object& right)
const;
1324 rpowO(
const boost::python::object& left)
const;
1386 boost::python::object
__add__(
const boost::python::object& right);
1393 boost::python::object
__sub__(
const boost::python::object& right);
1399 boost::python::object
__rsub__(
const boost::python::object& right);
1405 boost::python::object
__mul__(
const boost::python::object& right);
1411 boost::python::object
__div__(
const boost::python::object& right);
1417 boost::python::object
__rdiv__(
const boost::python::object& right);
1448 getItem(
const boost::python::object& key)
const;
1462 setItemD(
const boost::python::object& key,
1466 setItemO(
const boost::python::object& key,
1467 const boost::python::object& value);
1476 template <
class UnaryFunction>
1608 template <
class BinaryOp>
1625 boost::python::object
1662 template <
class BinaryFunction>
1675 template <
class BinaryFunction>
1746 ESYS_ASSERT(dr!=0,
"error casting to DataReady.");
1754 ESYS_ASSERT(dr!=0,
"error casting to DataReady.");
1766 ESYS_ASSERT(dr.get()!=0,
"error casting to DataReady.");
1774 ESYS_ASSERT(dr.get()!=0,
"error casting to DataReady.");
1794 #ifdef SLOWSHARECHECK 1795 return m_data->isShared();
1797 return !m_data.unique();
1806 if (omp_in_parallel())
1808 throw DataException(
"Please do not call forceResolve() in a parallel region.");
1822 if (omp_in_parallel())
1824 throw DataException(
"Programming error. Please do not run exclusiveWrite() in multi-threaded sections.");
1834 m_data->exclusivewritecalled=
true;
1845 std::ostringstream oss;
1846 oss <<
"Programming error. ExclusiveWrite required - please call requireWrite() isLazy=" <<
isLazy() <<
" isShared()=" <<
isShared();
1864 #ifdef IKNOWWHATIMDOING 1868 friend Data randomData(
const boost::python::tuple& shape,
const FunctionSpace& what,
long seed,
const boost::python::tuple& filter);
1873 #ifdef IKNOWWHATIMDOING 1888 long seed,
const boost::python::tuple& filter);
1914 throw DataException(
"Error, attempt to acquire RW access to lazy data. Please call requireWrite() first.");
1917 if (!
getReady()->exclusivewritecalled)
1919 throw DataException(
"Error, call to Data::getSampleDataRW without a preceeding call to requireWrite/exclusiveWrite.");
1931 throw DataException(
"Error, attempt to acquire RW access to lazy data. Please call requireWrite() first.");
1934 if (!
getReady()->exclusivewritecalled)
1936 throw DataException(
"Error, call to Data::getSampleDataRW without a preceeding call to requireWrite/exclusiveWrite.");
1952 return &((*res)[offset]);
1964 throw DataException(
"Programming error: complex lazy objects are not supported.");
1976 throw DataException(
"Programmer error - getDataRO must not be called on Lazy Data.");
1994 throw DataException(
"Programmer error - getDataRO must not be called on Lazy Data.");
2138 return *
this / right;
2176 template <
class BinaryFunction>
2183 ESYS_ASSERT(leftC!=0,
"Programming error - casting to DataExpanded.");
2191 const auto& vec=data.
getTypedVectorRO(
typename BinaryFunction::first_argument_type(0));
2195 #pragma omp parallel private(local_current_value) 2197 local_current_value=initial_value;
2198 #pragma omp for private(i,j) schedule(static) 2199 for (i=0;i<numSamples;i++) {
2200 for (j=0;j<numDPPSample;j++) {
2205 #pragma omp critical 2206 global_current_value=operation(global_current_value,local_current_value);
2208 return global_current_value;
2211 ESYS_ASSERT(leftC!=0,
"Programming error - casting to DataTagged.");
2216 const auto& vec=data.
getTypedVectorRO(
typename BinaryFunction::first_argument_type(0));
2220 for (std::list<int>::const_iterator i=used.begin();i!=used.end();++i)
2229 DataTagged::DataMapType::const_iterator it=lookup.find(tag);
2230 if (it!=lookup.end())
2236 return current_value;
2239 ESYS_ASSERT(leftC!=0,
"Programming error - casting to DataConstant.");
2242 throw DataException(
"Error - Operations (algorithm) not permitted on instances of DataEmpty.");
2244 throw DataException(
"Error - Operations not permitted on instances of DataLazy.");
2246 throw DataException(
"Error - Data encapsulates an unknown type.");
2258 template <
class BinaryFunction>
2264 throw DataException(
"Error - Operations (dp_algorithm) not permitted on instances of DataEmpty.");
2270 ESYS_ASSERT(dataE!=0,
"Programming error - casting data to DataExpanded.");
2271 ESYS_ASSERT(resultE!=0,
"Programming error - casting result to DataExpanded.");
2285 #pragma omp parallel for private(i,j) schedule(static) 2286 for (i=0;i<numSamples;i++) {
2287 for (j=0;j<numDPPSample;j++) {
2298 ESYS_ASSERT(dataT!=0,
"Programming error - casting data to DataTagged.");
2307 for (DataTagged::DataMapType::const_iterator i=lookup.begin(); i!=lookup.end(); i++) {
2317 return Data(resultT);
2323 ESYS_ASSERT(dataC!=0,
"Programming error - casting data to DataConstant.");
2324 ESYS_ASSERT(resultC!=0,
"Programming error - casting result to DataConstant.");
2333 throw DataException(
"Error - Operations not permitted on instances of DataLazy.");
2335 throw DataException(
"Error - Data encapsulates an unknown type.");
2360 #endif // __ESCRIPT_DATA_H__ Definition: FunctionSpace.h:34
std::ostream & operator<<(std::ostream &o, const Data &data)
Output operator.
Definition: Data.cpp:3427
void typeMatchLeft(Data &right) const
Convert the data type of the RHS to match this.
Definition: Data.cpp:3273
DataAbstract_ptr m_data
Definition: Data.h:1738
DataTypes::RealVectorType::size_type getLength() const
Return the number of doubles stored for this Data.
Definition: Data.cpp:1241
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:1252
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:1999
void setItemO(const boost::python::object &key, const boost::python::object &value)
Definition: Data.cpp:3234
Data antihermitian() const
Return the anti-hermitian part of a matrix which is half the matrix minus its hermitian.
Definition: Data.cpp:2564
Data interpolateFromTable1D(const WrappedArray &table, DataTypes::real_t Amin, DataTypes::real_t Astep, DataTypes::real_t undef, bool check_boundaries)
Definition: Data.cpp:4505
boost::python::object __add__(const boost::python::object &right)
wrapper for python add operation
Definition: Data.cpp:5615
const_DataReady_ptr getReadyPtr() const
Definition: Data.h:1771
void checkExclusiveWrite()
checks if caller can have exclusive write to the object
Definition: Data.h:1841
boost::python::object __rsub__(const boost::python::object &right)
wrapper for python reverse subtract operation
Definition: Data.cpp:5663
std::string toString() const
Write the data as a string. For large amounts of data, a summary is printed.
Definition: Data.cpp:4384
void setToZero()
set all values to zero
Definition: Data.cpp:742
DataTypes::real_t rpow(DataTypes::real_t x, DataTypes::real_t y)
Definition: Data.h:2010
ElementType & reference
Definition: DataVectorAlt.h:50
DataTypes::real_t inf()
Return the minimum value of this Data object.
Definition: Data.cpp:2104
void replaceNaN(DataTypes::real_t value)
replaces all NaN values with value
Definition: Data.cpp:2201
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:713
Data & operator*=(const Data &right)
Overloaded operator *=.
Definition: Data.cpp:2981
void tag()
If possible convert this Data to DataTagged. This will only allow Constant data to be converted to ta...
Definition: Data.cpp:1087
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:557
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:1227
boost::python::object integrateToTuple()
Calculate the integral over the function space domain as a python tuple.
Definition: Data.cpp:1637
Data getSlice(const DataTypes::RegionType ®ion) const
Return a Data object containing the specified slice of this Data object.
Definition: Data.cpp:3226
DataTypes::real_t supWorker() const
Definition: Data.cpp:2293
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:1959
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:1234
Definition: AbstractContinuousDomain.cpp:22
Data imag() const
Definition: Data.cpp:1826
Data atan() const
Return the atan of each data point of this Data object.
Definition: Data.cpp:1884
Data log10() const
Return the log to base 10 of each data point of this Data object.
Definition: Data.cpp:1945
Data pos() const
Return the identity of each data point of this Data object. Simply returns this object unmodified...
Definition: Data.cpp:1980
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:3242
Data operator-(const Data &left, const Data &right)
Operator- Takes two Data objects.
Definition: Data.cpp:3103
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:5187
DataTypes::real_t sup_const() const
Definition: Data.cpp:2051
Data besselSecondKind(int order)
Return the Bessel function of the second kind for each data point of this Data object.
Definition: Data.cpp:1689
void print(void)
print the data values to stdout. Used for debugging
Definition: Data.cpp:5127
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:4477
void typeMatchRight(const Data &right)
Convert the data type of this to match the RHS.
Definition: Data.cpp:3293
void setSlice(const Data &value, const DataTypes::RegionType ®ion)
Copy the specified slice from the given value into this Data object.
Definition: Data.cpp:3258
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:1433
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:4453
bool isLazy() const
Return true if this Data is lazy.
Definition: Data.cpp:1029
Data C_TensorUnaryOperation(Data const &arg_0, escript::ES_optype operation, DataTypes::real_t tol=0)
Definition: Data.cpp:5769
DataAbstract_ptr borrowDataPtr(void) const
Definition: Data.cpp:4369
Data real() const
Definition: Data.cpp:1807
virtual const DataTypes::RealVectorType & getTypedVectorRO(DataTypes::real_t dummy) const
Definition: DataConstant.cpp:602
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:3321
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:1384
Data sinh() const
Return the sinh of each data point of this Data object.
Definition: Data.cpp:1891
Data atanh() const
Return the atanh of each data point of this Data object.
Definition: Data.cpp:1938
DataReady_ptr getReadyPtr()
Definition: Data.h:1763
void TensorSelfUpdateBinaryOperation(const Data &right, escript::ES_optype operation)
Definition: Data.cpp:6220
DataTypes::RealVectorType::reference getDataByTagRW(int tag, DataTypes::RealVectorType::size_type i, DataTypes::real_t dummy=0)
getDataByTag
Definition: DataTagged.cpp:875
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:3067
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:4496
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:3378
Data nonuniformslope(boost::python::object in, boost::python::object out, bool check_boundaries)
Definition: Data.cpp:5048
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:5639
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:4487
Data acos() const
Return the acos of each data point of this Data object.
Definition: Data.cpp:1876
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:1200
boost::shared_ptr< DataAbstract > DataAbstract_ptr
Definition: DataAbstract.h:51
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:3212
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:586
Data maxval() const
Return the maximum value of each data point of this Data object.
Definition: Data.cpp:2403
void setValueOfDataPointC(int dataPointNo, const DataTypes::cplx_t)
Definition: Data.cpp:1550
Data cos() const
Return the cos of each data point of this Data object.
Definition: Data.cpp:1855
const DataTypes::real_t * getDataRO(DataTypes::real_t dummy=0) const
Return a pointer to the beginning of the underlying data.
Definition: Data.h:1972
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:3434
bool m_protected
Definition: Data.h:1732
Data cosh() const
Return the cosh of each data point of this Data object.
Definition: Data.cpp:1898
void setValueOfDataPoint(int dataPointNo, const DataTypes::real_t)
sets the values of a data-point on this process
Definition: Data.cpp:1529
Data interpolate(const FunctionSpace &functionspace) const
Interpolates this onto the given functionspace and returns the result as a Data object.
Definition: Data.cpp:1194
Data matrixInverse() const
return inverse of matricies.
Definition: Data.cpp:3039
Data & operator/=(const Data &right)
Overloaded operator /=.
Definition: Data.cpp:3008
DataReady * getReady()
Definition: Data.h:1751
Data copySelf() const
Return a pointer to a deep copy of this object.
Definition: Data.cpp:697
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:1113
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:1015
Data tanh() const
Return the tanh of each data point of this Data object.
Definition: Data.cpp:1905
void calc_minGlobalDataPoint(int &ProcNo, int &DataPointNo) const
Definition: Data.cpp:2750
void operandCheck(const Data &right) const
Check *this and the right operand are compatible. Throws an exception if they aren't.
Definition: Data.h:1652
DataAbstract * borrowData(void) const
return the object produced by the factory, which is a DataConstant or DataExpanded TODO Ownership of ...
Definition: Data.cpp:4362
Data erf() const
Return the error function erf of each data point of this Data object.
Definition: Data.cpp:1913
const DataTypes::ShapeType & getShape() const
Return the shape information for the point data.
Definition: DataAbstract.h:631
int getNumSamples() const
Return the number of samples.
Definition: DataAbstract.h:612
int getNumDPPSample() const
Return the number of data points per sample.
Definition: DataAbstract.h:601
DataTypes::real_t infWorker() const
Definition: Data.cpp:2336
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:481
Data rtruedivO(const boost::python::object &left)
Newer style division operator for python.
Definition: Data.h:2161
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:2828
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:1945
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:3081
void replaceNaNPython(boost::python::object obj)
replaces all NaN values with value
Definition: Data.cpp:2221
Data trace(int axis_offset) const
Return the trace of a matrix.
Definition: Data.cpp:2604
DataTypes::RealVectorType::size_type getOffsetForTag(int tag) const
getOffsetForTag
Definition: DataTagged.cpp:853
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:2736
Data minval_nonlazy() const
Definition: Data.cpp:2381
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:5163
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:724
Data transpose(int axis_offset) const
Transpose each data point of this Data object around the given axis.
Definition: Data.cpp:2659
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:5175
DataTypes::RealVectorType::reference getDataAtOffsetRW(DataTypes::RealVectorType::size_type i, DataTypes::real_t dummy=0)
Definition: Data.cpp:4422
void expand()
Whatever the current Data type make this into a DataExpanded.
Definition: Data.cpp:1063
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:2179
Data operator+(const Data &left, const Data &right)
Operator+ Takes two Data objects.
Definition: Data.cpp:3093
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:1180
DataTypes::RealVectorType::reference getDataPointRW(int sampleNo, int dataPointNo)
Return a reference into the DataVector which points to the specified data point.
Definition: Data.cpp:4469
Data asinh() const
Return the asinh of each data point of this Data object.
Definition: Data.cpp:1924
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:595
ES_optype
Definition: ES_optype.h:26
Data asin() const
Return the asin of each data point of this Data object.
Definition: Data.cpp:1869
bool isConstant() const
Return true if this Data is constant.
Definition: Data.cpp:1022
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:1486
bool isComplex() const
True if components of this data are stored as complex.
Definition: Data.cpp:1043
bool isShared() const
Definition: Data.h:1792
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:1532
DataReady_ptr borrowReadyPtr(void) const
Definition: Data.cpp:4376
Data acosh() const
Return the acosh of each data point of this Data object.
Definition: Data.cpp:1931
Data wherePositive() const
Return a Data with a 1 for +ive values and a 0 for 0 or -ive values.
Definition: Data.cpp:1136
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:1001
Data & operator=(const Data &other)
Definition: Data.cpp:2946
void forceResolve()
Definition: Data.h:1801
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:1122
DataTypes::real_t Lsup_const() const
Definition: Data.cpp:2006
int getNumDataPointsPerSample() const
Return the number of data points per sample.
Definition: Data.h:521
const DataReady * getReady() const
Definition: Data.h:1743
DataTypes::RealVectorType::size_type getDefaultOffset() const
Returns the offset in the structure which stores the default value.
Definition: DataTagged.h:693
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:1910
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:2537
Data oneOver() const
Returns 1./ Data object.
Definition: Data.cpp:1129
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:1377
DataTypes::real_t LsupWorker() const
Definition: Data.cpp:2236
Data minval() const
Return the minimum value of each data point of this Data object.
Definition: Data.cpp:2412
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:6044
bool m_lazy
Definition: Data.h:1733
DataTypes::real_t sup()
Return the maximum value of this Data object.
Definition: Data.cpp:2065
boost::python::object integrateWorker() const
Definition: Data.cpp:1648
void setTupleForGlobalDataPoint(int id, int proc, boost::python::object)
Set the value of a global data point.
Definition: Data.cpp:1441
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:5712
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:1966
Data grad() const
Definition: Data.cpp:1223
friend Data condEval(escript::Data &mask, escript::Data &trueval, escript::Data &falseval)
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:2261
void copy(const Data &other)
Make this object a deep copy of "other".
Definition: Data.cpp:704
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:2090
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:2710
void init_from_data_and_fs(const Data &inData, const FunctionSpace &functionspace)
Definition: Data.cpp:433
DataTypes::RealVectorType::const_reference getDataAtOffsetRO(DataTypes::RealVectorType::size_type i, DataTypes::real_t dummy=0)
Return a pointer to the beginning of the datapoint at the specified offset. TODO Eventually these sho...
Definition: Data.cpp:4430
int getTagNumber(int dpno)
Return the tag number associated with the given data-point.
Definition: Data.cpp:3417
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:5760
Data nonuniforminterp(boost::python::object in, boost::python::object out, bool check_boundaries)
Definition: Data.cpp:4972
virtual const DataTypes::RealVectorType & getTypedVectorRO(DataTypes::real_t dummy) const
Definition: DataExpanded.cpp:1233
Data maxval_nonlazy() const
Definition: Data.cpp:2392
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:1361
Data powO(const boost::python::object &right) const
Return the given power of each data point of this boost python object.
Definition: Data.cpp:3074
Data whereNonNegative() const
Return a Data with a 1 for +ive or 0 values and a 0 for -ive values.
Definition: Data.cpp:1158
Data exp() const
Return the exponential function of each data point of this Data object.
Definition: Data.cpp:1992
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:1819
Data whereNegative() const
Return a Data with a 1 for -ive values and a 0 for +ive or 0 values.
Definition: Data.cpp:1147
bool isTagged() const
Return true if this Data is tagged.
Definition: Data.cpp:1008
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:1787
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:4633
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:5141
void setProtection()
switches on update protection
Definition: Data.cpp:1049
const FunctionSpace & getFunctionSpace() const
Return the function space associated with this Data object.
Definition: DataAbstract.h:624
Data besselFirstKind(int order)
Return the Bessel function of the first kind for each data point of this Data object.
Definition: Data.cpp:1682
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:5736
bool hasNaN()
returns return true if data contains NaN.
Definition: Data.cpp:2190
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:3122
bool isExpanded() const
Return true if this Data is expanded.
Definition: Data.cpp:994
Data whereNonPositive() const
Return a Data with a 1 for -ive or 0 values and a 0 for +ive values.
Definition: Data.cpp:1169
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:764
boost::python::object __mul__(const boost::python::object &right)
wrapper for python multiply operation
Definition: Data.cpp:5688
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:5562
Data neg() const
Return the negation of each data point of this Data object.
Definition: Data.cpp:1973
const boost::python::tuple getShapeTuple() const
Return the data point shape as a tuple of integers.
Definition: Data.cpp:672
const DataMapType & getTagLookup() const
getTagLookup
Definition: DataTagged.h:731
Data truedivO(const boost::python::object &right)
Newer style division operator for python.
Definition: Data.h:2148
bool isReady() const
Return true if this data is ready.
Definition: Data.cpp:1036
Data symmetric() const
Return the symmetric part of a matrix which is half the matrix plus its transpose.
Definition: Data.cpp:2478
#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:5544
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:2501
Data sin() const
Return the sin of each data point of this Data object.
Definition: Data.cpp:1848
~Data()
Destructor.
Definition: Data.cpp:549
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:1206
const boost::python::object getValueOfDataPointAsTuple(int dataPointNo)
Return the value of a data point as a python tuple.
Definition: Data.cpp:1394
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:4789
bool isProtected() const
Returns true, if the data object is protected against update.
Definition: Data.cpp:1055
Data log() const
Return the natural log of each data point of this Data object.
Definition: Data.cpp:1952
Data truedivD(const Data &right)
Newer style division operator for python.
Definition: Data.h:2136
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:1187
DataTypes::real_t Lsup()
Return the maximum absolute value of this Data object.
Definition: Data.cpp:2016
Data & operator-=(const Data &right)
Overloaded operator -=.
Definition: Data.cpp:2954
boost::python::object integrateToTuple_const() const
Calculate the integral over the function space domain as a python tuple.
Definition: Data.cpp:1626
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:2838
Data eigenvalues() const
Return the eigenvalues of the symmetric part at each data point of this Data object in increasing val...
Definition: Data.cpp:2683
Data bessel(int order, DataTypes::real_t(*besselfunc)(int, DataTypes::real_t))
Bessel worker function.
Definition: Data.cpp:1696
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:3338
const boost::python::object getValueOfGlobalDataPointAsTuple(int procNo, int dataPointNo)
Return a data point across all processors as a python tuple.
Definition: Data.cpp:1573
Data swapaxes(const int axis0, const int axis1) const
swaps the components axis0 and axis1
Definition: Data.cpp:2421
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:1862
Data & operator+=(const Data &right)
Overloaded operator +=.
Definition: Data.cpp:2918