src/MEDCoupling/MEDCouplingFieldDouble.hxx

   1 // Copyright (C) 2007-2012  CEA/DEN, EDF R&D
   2 //
   3 // This library is free software; you can redistribute it and/or
   4 // modify it under the terms of the GNU Lesser General Public
   5 // License as published by the Free Software Foundation; either
   6 // version 2.1 of the License.
   7 //
   8 // This library is distributed in the hope that it will be useful,
   9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 // Lesser General Public License for more details.
  12 //
  13 // You should have received a copy of the GNU Lesser General Public
  14 // License along with this library; if not, write to the Free Software
  15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  16 //
  17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  18 //
  19 // Author : Anthony Geay (CEA/DEN)
  20
  21 #ifndef __PARAMEDMEM_MEDCOUPLINGFIELDDOUBLE_HXX__
  22 #define __PARAMEDMEM_MEDCOUPLINGFIELDDOUBLE_HXX__
  23
  24 #include "MEDCoupling.hxx"
  25 #include "MEDCouplingField.hxx"
  26 #include "MEDCouplingTimeDiscretization.hxx"
  27 #include "MEDCouplingMemArray.hxx"
  28
  29 namespace ParaMEDMEM
  30 {
  31   class MEDCouplingFieldTemplate;
  32
  33   class MEDCOUPLING_EXPORT MEDCouplingFieldDouble : public MEDCouplingField
  34   {
  35   public:
  36     static MEDCouplingFieldDouble *New(TypeOfField type, TypeOfTimeDiscretization td=NO_TIME);
  37     static MEDCouplingFieldDouble *New(const MEDCouplingFieldTemplate *ft, TypeOfTimeDiscretization td=NO_TIME);
  38     void setTimeUnit(const char *unit);
  39     const char *getTimeUnit() const;
  40     void copyTinyStringsFrom(const MEDCouplingFieldDouble *other) throw(INTERP_KERNEL::Exception);
  41     void copyTinyAttrFrom(const MEDCouplingFieldDouble *other) throw(INTERP_KERNEL::Exception);
  42     std::string simpleRepr() const;
  43     std::string advancedRepr() const;
  44     bool isEqualIfNotWhy(const MEDCouplingField *other, double meshPrec, double valsPrec, std::string& reason) const throw(INTERP_KERNEL::Exception);
  45     bool isEqualWithoutConsideringStr(const MEDCouplingField *other, double meshPrec, double valsPrec) const;
  46     bool areCompatibleForMerge(const MEDCouplingField *other) const;
  47     bool areStrictlyCompatible(const MEDCouplingField *other) const;
  48     bool areCompatibleForMul(const MEDCouplingField *other) const;
  49     bool areCompatibleForDiv(const MEDCouplingField *other) const;
  50     bool areCompatibleForMeld(const MEDCouplingFieldDouble *other) const;
  51     void renumberCells(const int *old2NewBg, bool check=true) throw(INTERP_KERNEL::Exception);
  52     void renumberCellsWithoutMesh(const int *old2NewBg, bool check=true) throw(INTERP_KERNEL::Exception);
  53     void renumberNodes(const int *old2NewBg) throw(INTERP_KERNEL::Exception);
  54     void renumberNodesWithoutMesh(const int *old2NewBg, int newNbOfNodes, double eps=1e-15) throw(INTERP_KERNEL::Exception);
  55     DataArrayInt *getIdsInRange(double vmin, double vmax) const throw(INTERP_KERNEL::Exception);
  56     MEDCouplingFieldDouble *buildSubPart(const DataArrayInt *part) const throw(INTERP_KERNEL::Exception);
  57     MEDCouplingFieldDouble *buildSubPart(const int *partBg, const int *partEnd) const throw(INTERP_KERNEL::Exception);
  58     MEDCouplingFieldDouble *deepCpy() const;
  59     MEDCouplingFieldDouble *clone(bool recDeepCpy) const;
  60     MEDCouplingFieldDouble *cloneWithMesh(bool recDeepCpy) const;
  61     MEDCouplingFieldDouble *buildNewTimeReprFromThis(TypeOfTimeDiscretization td, bool deepCopy) const;
  62     TypeOfTimeDiscretization getTimeDiscretization() const;
  63     void checkCoherency() const throw(INTERP_KERNEL::Exception);
  64     void setNature(NatureOfField nat) throw(INTERP_KERNEL::Exception);
  65     void setTimeTolerance(double val) { _time_discr->setTimeTolerance(val); }
  66     double getTimeTolerance() const { return _time_discr->getTimeTolerance(); }
  67     void setIteration(int it) throw(INTERP_KERNEL::Exception) { _time_discr->setIteration(it); }
  68     void setEndIteration(int it) throw(INTERP_KERNEL::Exception) { _time_discr->setEndIteration(it); }
  69     void setOrder(int order) throw(INTERP_KERNEL::Exception) { _time_discr->setOrder(order); }
  70     void setEndOrder(int order) throw(INTERP_KERNEL::Exception) { _time_discr->setEndOrder(order); }
  71     void setTimeValue(double val) throw(INTERP_KERNEL::Exception) { _time_discr->setTimeValue(val); }
  72     void setEndTimeValue(double val) throw(INTERP_KERNEL::Exception) { _time_discr->setEndTimeValue(val); }
  73     void setTime(double val, int iteration, int order) { _time_discr->setTime(val,iteration,order); }
  74     void setStartTime(double val, int iteration, int order) { _time_discr->setStartTime(val,iteration,order); }
  75     void setEndTime(double val, int iteration, int order) { _time_discr->setEndTime(val,iteration,order); }
  76     double getTime(int& iteration, int& order) const { return _time_discr->getTime(iteration,order); }
  77     double getStartTime(int& iteration, int& order) const { return _time_discr->getStartTime(iteration,order); }
  78     double getEndTime(int& iteration, int& order) const { return _time_discr->getEndTime(iteration,order); }
  79     double getIJ(int tupleId, int compoId) const { return getArray()->getIJ(tupleId,compoId); }
  80     double getIJK(int cellId, int nodeIdInCell, int compoId) const;
  81     void setArray(DataArrayDouble *array);
  82     void setEndArray(DataArrayDouble *array);
  83     void setArrays(const std::vector<DataArrayDouble *>& arrs) throw(INTERP_KERNEL::Exception);
  84     const DataArrayDouble *getArray() const { return _time_discr->getArray(); }
  85     DataArrayDouble *getArray() { return _time_discr->getArray(); }
  86     const DataArrayDouble *getEndArray() const { return _time_discr->getEndArray(); }
  87     DataArrayDouble *getEndArray() { return _time_discr->getEndArray(); }
  88     std::vector<DataArrayDouble *> getArrays() const { std::vector<DataArrayDouble *> ret; _time_discr->getArrays(ret); return ret; }
  89     double accumulate(int compId) const;
  90     void accumulate(double *res) const;
  91     double getMaxValue() const throw(INTERP_KERNEL::Exception);
  92     double getMaxValue2(DataArrayInt*& tupleIds) const throw(INTERP_KERNEL::Exception);
  93     double getMinValue() const throw(INTERP_KERNEL::Exception);
  94     double getMinValue2(DataArrayInt*& tupleIds) const throw(INTERP_KERNEL::Exception);
  95     double getAverageValue() const throw(INTERP_KERNEL::Exception);
  96     double norm2() const throw(INTERP_KERNEL::Exception);
  97     double normMax() const throw(INTERP_KERNEL::Exception);
  98     double getWeightedAverageValue() const throw(INTERP_KERNEL::Exception);
  99     double normL1(int compId) const throw(INTERP_KERNEL::Exception);
 100     void normL1(double *res) const throw(INTERP_KERNEL::Exception);
 101     double normL2(int compId) const throw(INTERP_KERNEL::Exception);
 102     void normL2(double *res) const throw(INTERP_KERNEL::Exception);
 103     double integral(int compId, bool isWAbs) const throw(INTERP_KERNEL::Exception);
 104     void integral(bool isWAbs, double *res) const throw(INTERP_KERNEL::Exception);
 105     void getValueOnPos(int i, int j, int k, double *res) const throw(INTERP_KERNEL::Exception);
 106     void getValueOn(const double *spaceLoc, double *res) const throw(INTERP_KERNEL::Exception);
 107     void getValueOn(const double *spaceLoc, double time, double *res) const throw(INTERP_KERNEL::Exception);
 108     DataArrayDouble *getValueOnMulti(const double *spaceLoc, int nbOfPoints) const throw(INTERP_KERNEL::Exception);
 109     //! \b temporary
 110     void applyLin(double a, double b, int compoId);
 111     MEDCouplingFieldDouble &operator=(double value) throw(INTERP_KERNEL::Exception);
 112     void fillFromAnalytic(int nbOfComp, FunctionToEvaluate func) throw(INTERP_KERNEL::Exception);
 113     void fillFromAnalytic(int nbOfComp, const char *func) throw(INTERP_KERNEL::Exception);
 114     void fillFromAnalytic2(int nbOfComp, const char *func) throw(INTERP_KERNEL::Exception);
 115     void fillFromAnalytic3(int nbOfComp, const std::vector<std::string>& varsOrder, const char *func) throw(INTERP_KERNEL::Exception);
 116     void applyFunc(int nbOfComp, FunctionToEvaluate func);
 117     void applyFunc(int nbOfComp, double val);
 118     void applyFunc(int nbOfComp, const char *func) throw(INTERP_KERNEL::Exception);
 119     void applyFunc2(int nbOfComp, const char *func) throw(INTERP_KERNEL::Exception);
 120     void applyFunc3(int nbOfComp, const std::vector<std::string>& varsOrder, const char *func) throw(INTERP_KERNEL::Exception);
 121     void applyFunc(const char *func) throw(INTERP_KERNEL::Exception);
 122     void applyFuncFast32(const char *func) throw(INTERP_KERNEL::Exception);
 123     void applyFuncFast64(const char *func) throw(INTERP_KERNEL::Exception);
 124     int getNumberOfComponents() const throw(INTERP_KERNEL::Exception);
 125     int getNumberOfTuples() const throw(INTERP_KERNEL::Exception);
 126     int getNumberOfValues() const throw(INTERP_KERNEL::Exception);
 127     void updateTime() const;
 128     //
 129     void getTinySerializationIntInformation(std::vector<int>& tinyInfo) const;
 130     void getTinySerializationDbleInformation(std::vector<double>& tinyInfo) const;
 131     void getTinySerializationStrInformation(std::vector<std::string>& tinyInfo) const;
 132     void resizeForUnserialization(const std::vector<int>& tinyInfoI, DataArrayInt *&dataInt, std::vector<DataArrayDouble *>& arrays);
 133     void finishUnserialization(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD, const std::vector<std::string>& tinyInfoS);
 134     void serialize(DataArrayInt *&dataInt, std::vector<DataArrayDouble *>& arrays) const;
 135     //
 136     void changeUnderlyingMesh(const MEDCouplingMesh *other, int levOfCheck, double prec) throw(INTERP_KERNEL::Exception);
 137     void substractInPlaceDM(const MEDCouplingFieldDouble *f, int levOfCheck, double prec) throw(INTERP_KERNEL::Exception);
 138     bool mergeNodes(double eps, double epsOnVals=1e-15) throw(INTERP_KERNEL::Exception);
 139     bool mergeNodes2(double eps, double epsOnVals=1e-15) throw(INTERP_KERNEL::Exception);
 140     bool zipCoords(double epsOnVals=1e-15) throw(INTERP_KERNEL::Exception);
 141     bool zipConnectivity(int compType, double epsOnVals=1e-15) throw(INTERP_KERNEL::Exception);
 142     MEDCouplingFieldDouble *extractSlice3D(const double *origin, const double *vec, double eps) const throw(INTERP_KERNEL::Exception);
 143     bool simplexize(int policy) throw(INTERP_KERNEL::Exception);
 144     MEDCouplingFieldDouble *doublyContractedProduct() const throw(INTERP_KERNEL::Exception);
 145     MEDCouplingFieldDouble *determinant() const throw(INTERP_KERNEL::Exception);
 146     MEDCouplingFieldDouble *eigenValues() const throw(INTERP_KERNEL::Exception);
 147     MEDCouplingFieldDouble *eigenVectors() const throw(INTERP_KERNEL::Exception);
 148     MEDCouplingFieldDouble *inverse() const throw(INTERP_KERNEL::Exception);
 149     MEDCouplingFieldDouble *trace() const throw(INTERP_KERNEL::Exception);
 150     MEDCouplingFieldDouble *deviator() const throw(INTERP_KERNEL::Exception);
 151     MEDCouplingFieldDouble *magnitude() const throw(INTERP_KERNEL::Exception);
 152     MEDCouplingFieldDouble *maxPerTuple() const throw(INTERP_KERNEL::Exception);
 153     void changeNbOfComponents(int newNbOfComp, double dftValue=0.) throw(INTERP_KERNEL::Exception);
 154     MEDCouplingFieldDouble *keepSelectedComponents(const std::vector<int>& compoIds) const throw(INTERP_KERNEL::Exception);
 155     void setSelectedComponents(const MEDCouplingFieldDouble *f, const std::vector<int>& compoIds) throw(INTERP_KERNEL::Exception);
 156     void sortPerTuple(bool asc) throw(INTERP_KERNEL::Exception);
 157     static MEDCouplingFieldDouble *MergeFields(const MEDCouplingFieldDouble *f1, const MEDCouplingFieldDouble *f2) throw(INTERP_KERNEL::Exception);
 158     static MEDCouplingFieldDouble *MergeFields(const std::vector<const MEDCouplingFieldDouble *>& a) throw(INTERP_KERNEL::Exception);
 159     static MEDCouplingFieldDouble *MeldFields(const MEDCouplingFieldDouble *f1, const MEDCouplingFieldDouble *f2) throw(INTERP_KERNEL::Exception);
 160     static MEDCouplingFieldDouble *DotFields(const MEDCouplingFieldDouble *f1, const MEDCouplingFieldDouble *f2) throw(INTERP_KERNEL::Exception);
 161     MEDCouplingFieldDouble *dot(const MEDCouplingFieldDouble& other) const throw(INTERP_KERNEL::Exception) { return DotFields(this,&other); }
 162     static MEDCouplingFieldDouble *CrossProductFields(const MEDCouplingFieldDouble *f1, const MEDCouplingFieldDouble *f2) throw(INTERP_KERNEL::Exception);
 163     MEDCouplingFieldDouble *crossProduct(const MEDCouplingFieldDouble& other) const throw(INTERP_KERNEL::Exception) { return CrossProductFields(this,&other); }
 164     static MEDCouplingFieldDouble *MaxFields(const MEDCouplingFieldDouble *f1, const MEDCouplingFieldDouble *f2) throw(INTERP_KERNEL::Exception);
 165     MEDCouplingFieldDouble *max(const MEDCouplingFieldDouble& other) const throw(INTERP_KERNEL::Exception) { return MaxFields(this,&other); }
 166     static MEDCouplingFieldDouble *MinFields(const MEDCouplingFieldDouble *f1, const MEDCouplingFieldDouble *f2) throw(INTERP_KERNEL::Exception);
 167     MEDCouplingFieldDouble *min(const MEDCouplingFieldDouble& other) const throw(INTERP_KERNEL::Exception) { return MinFields(this,&other); }
 168     MEDCouplingFieldDouble *operator+(const MEDCouplingFieldDouble& other) const throw(INTERP_KERNEL::Exception) { return AddFields(this,&other); }
 169     const MEDCouplingFieldDouble &operator+=(const MEDCouplingFieldDouble& other) throw(INTERP_KERNEL::Exception);
 170     static MEDCouplingFieldDouble *AddFields(const MEDCouplingFieldDouble *f1, const MEDCouplingFieldDouble *f2) throw(INTERP_KERNEL::Exception);
 171     MEDCouplingFieldDouble *operator-(const MEDCouplingFieldDouble& other) const throw(INTERP_KERNEL::Exception) { return SubstractFields(this,&other); }
 172     const MEDCouplingFieldDouble &operator-=(const MEDCouplingFieldDouble& other) throw(INTERP_KERNEL::Exception);
 173     static MEDCouplingFieldDouble *SubstractFields(const MEDCouplingFieldDouble *f1, const MEDCouplingFieldDouble *f2) throw(INTERP_KERNEL::Exception);
 174     MEDCouplingFieldDouble *operator*(const MEDCouplingFieldDouble& other) const throw(INTERP_KERNEL::Exception) { return MultiplyFields(this,&other); }
 175     const MEDCouplingFieldDouble &operator*=(const MEDCouplingFieldDouble& other) throw(INTERP_KERNEL::Exception);
 176     static MEDCouplingFieldDouble *MultiplyFields(const MEDCouplingFieldDouble *f1, const MEDCouplingFieldDouble *f2) throw(INTERP_KERNEL::Exception);
 177     MEDCouplingFieldDouble *operator/(const MEDCouplingFieldDouble& other) const throw(INTERP_KERNEL::Exception) { return DivideFields(this,&other); }
 178     const MEDCouplingFieldDouble &operator/=(const MEDCouplingFieldDouble& other) throw(INTERP_KERNEL::Exception);
 179     static MEDCouplingFieldDouble *DivideFields(const MEDCouplingFieldDouble *f1, const MEDCouplingFieldDouble *f2) throw(INTERP_KERNEL::Exception);
 180     static void WriteVTK(const char *fileName, const std::vector<const MEDCouplingFieldDouble *>& fs) throw(INTERP_KERNEL::Exception);
 181   public:
 182     const MEDCouplingTimeDiscretization *getTimeDiscretizationUnderGround() const { return _time_discr; }
 183     MEDCouplingTimeDiscretization *getTimeDiscretizationUnderGround() { return _time_discr; }
 184   private:
 185     MEDCouplingFieldDouble(TypeOfField type, TypeOfTimeDiscretization td);
 186     MEDCouplingFieldDouble(const MEDCouplingFieldTemplate *ft, TypeOfTimeDiscretization td);
 187     MEDCouplingFieldDouble(const MEDCouplingFieldDouble& other, bool deepCopy);
 188     MEDCouplingFieldDouble(NatureOfField n, MEDCouplingTimeDiscretization *td, MEDCouplingFieldDiscretization *type);
 189     ~MEDCouplingFieldDouble();
 190   private:
 191     MEDCouplingTimeDiscretization *_time_discr;
 192   };
 193 }
 194
 195 #endif