Salome HOME
2b88a67c543127ae1cec29046149b686c25d3ffe
[tools/medcoupling.git] / src / MEDCoupling / MEDCouplingTimeDiscretization.hxx
1 // Copyright (C) 2007-2016  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, or (at your option) any later version.
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 __MEDCOUPLINGTIMEDISCRETIZATION_HXX__
22 #define __MEDCOUPLINGTIMEDISCRETIZATION_HXX__
23
24 #include "MEDCoupling.hxx"
25 #include "MEDCouplingTimeLabel.hxx"
26 #include "MEDCouplingRefCountObject.hxx"
27 #include "InterpKernelException.hxx"
28 #include "MEDCouplingTraits.hxx"
29
30 #include <vector>
31
32 namespace MEDCoupling
33 {
34   class MEDCouplingMesh;
35   class DataArrayDouble;
36   class TimeLabel;
37
38   template<class T>
39   class MEDCouplingTimeDiscretizationTemplate : public TimeLabel, public BigMemoryObject
40   {
41   public:
42     MEDCOUPLING_EXPORT void updateTime() const;
43     MEDCOUPLING_EXPORT virtual void setArray(typename Traits<T>::ArrayType *array, TimeLabel *owner);
44     MEDCOUPLING_EXPORT typename Traits<T>::ArrayType *getArray() { return _array; }
45     MEDCOUPLING_EXPORT const typename Traits<T>::ArrayType *getArray() const { return _array; }
46     MEDCOUPLING_EXPORT void setTimeUnit(const std::string& unit) { _time_unit=unit; }
47     MEDCOUPLING_EXPORT std::string getTimeUnit() const { return _time_unit; }
48     MEDCOUPLING_EXPORT void setTimeTolerance(double val) { _time_tolerance=val; }
49     MEDCOUPLING_EXPORT double getTimeTolerance() const { return _time_tolerance; }
50     MEDCOUPLING_EXPORT double getTime(int& iteration, int& order) const { return getStartTime(iteration,order); }
51     MEDCOUPLING_EXPORT void setTime(double time, int iteration, int order) { setStartTime(time,iteration,order); }
52     MEDCOUPLING_EXPORT void setIteration(int it) { setStartIteration(it); }
53     MEDCOUPLING_EXPORT void setOrder(int order) { setStartOrder(order); }
54     MEDCOUPLING_EXPORT void setTimeValue(double val) { setStartTimeValue(val); }
55     MEDCOUPLING_EXPORT virtual double getStartTime(int& iteration, int& order) const = 0;
56     MEDCOUPLING_EXPORT virtual void setStartIteration(int it) = 0;
57     MEDCOUPLING_EXPORT virtual void setStartOrder(int order) = 0;
58     MEDCOUPLING_EXPORT virtual void setStartTime(double time, int iteration, int order) = 0;
59     MEDCOUPLING_EXPORT virtual void setStartTimeValue(double time) = 0;
60     MEDCOUPLING_EXPORT virtual std::string getStringRepr() const = 0;
61     MEDCOUPLING_EXPORT virtual TypeOfTimeDiscretization getEnum() const = 0;
62     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretizationTemplate<T> *performCopyOrIncrRef(bool deepCopy) const = 0;
63     //
64     MEDCOUPLING_EXPORT virtual void copyTinyAttrFrom(const MEDCouplingTimeDiscretizationTemplate<T>& other);
65     MEDCOUPLING_EXPORT virtual void copyTinyStringsFrom(const MEDCouplingTimeDiscretizationTemplate<T>& other);
66     MEDCOUPLING_EXPORT virtual void checkConsistencyLight() const;
67     MEDCOUPLING_EXPORT virtual std::size_t getHeapMemorySizeWithoutChildren() const;
68     MEDCOUPLING_EXPORT virtual std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
69     //
70     MEDCOUPLING_EXPORT virtual bool areStrictlyCompatible(const MEDCouplingTimeDiscretizationTemplate<T> *other, std::string& reason) const;
71     MEDCOUPLING_EXPORT virtual bool areCompatible(const MEDCouplingTimeDiscretizationTemplate<T> *other) const;
72     MEDCOUPLING_EXPORT virtual bool areStrictlyCompatibleForMul(const MEDCouplingTimeDiscretizationTemplate<T> *other) const;
73     MEDCOUPLING_EXPORT virtual bool areStrictlyCompatibleForDiv(const MEDCouplingTimeDiscretizationTemplate<T> *other) const;
74   protected:
75     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretizationTemplate();
76     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretizationTemplate(const MEDCouplingTimeDiscretizationTemplate<T>& other, bool deepCopy);
77     MEDCOUPLING_EXPORT virtual ~MEDCouplingTimeDiscretizationTemplate();
78   protected:
79     std::string _time_unit;
80     double _time_tolerance;
81     typename Traits<T>::ArrayType *_array;
82   protected:
83     static const double TIME_TOLERANCE_DFT;
84   };
85
86   class MEDCouplingTimeKeeper
87   {
88   public:
89     MEDCOUPLING_EXPORT MEDCouplingTimeKeeper():_time(0.),_iteration(-1),_order(-1) { }
90     MEDCOUPLING_EXPORT double getAllInfo(int& iteration, int& order) const { iteration=_iteration; order=_order; return _time; }
91     MEDCOUPLING_EXPORT void setAllInfo(double time, int iteration, int order) { _time=time; _iteration=iteration; _order=order; }
92     MEDCOUPLING_EXPORT int getIteration() const { return _iteration; }
93     MEDCOUPLING_EXPORT void setIteration(int it) { _iteration=it; }
94     MEDCOUPLING_EXPORT int getOrder() const { return _order; }
95     MEDCOUPLING_EXPORT void setOrder(int order) { _order=order; }
96     MEDCOUPLING_EXPORT double getTimeValue() const { return _time; }
97     MEDCOUPLING_EXPORT void setTimeValue(double time) { _time=time; }
98     MEDCOUPLING_EXPORT bool isEqualIfNotWhy(const MEDCouplingTimeKeeper& other, double prec, std::string& reason) const;
99     MEDCOUPLING_EXPORT bool isEqual(const MEDCouplingTimeKeeper& other, double prec) const;
100     MEDCOUPLING_EXPORT void copyFrom(const MEDCouplingTimeKeeper& other);
101     MEDCOUPLING_EXPORT void checkTimePresence(double time, double eps) const;
102   private:
103     double _time;
104     int _iteration;
105     int _order;
106   };
107
108   class MEDCouplingTimeDiscretization : public MEDCouplingTimeDiscretizationTemplate<double>
109   {
110   protected:
111     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization() { }
112     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization(const MEDCouplingTimeDiscretization& other, bool deepCopy):MEDCouplingTimeDiscretizationTemplate<double>(other,deepCopy) { }
113   public:
114     MEDCOUPLING_EXPORT static MEDCouplingTimeDiscretization *New(TypeOfTimeDiscretization type);
115     MEDCOUPLING_EXPORT virtual bool areCompatibleForMeld(const MEDCouplingTimeDiscretization *other) const;
116     MEDCOUPLING_EXPORT virtual bool isEqualIfNotWhy(const MEDCouplingTimeDiscretization *other, double prec, std::string& reason) const;
117     MEDCOUPLING_EXPORT virtual bool isEqual(const MEDCouplingTimeDiscretization *other, double prec) const;
118     MEDCOUPLING_EXPORT virtual bool isEqualWithoutConsideringStr(const MEDCouplingTimeDiscretization *other, double prec) const;
119     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *buildNewTimeReprFromThis(TypeOfTimeDiscretization type, bool deepCopy) const;
120     MEDCOUPLING_EXPORT virtual void synchronizeTimeWith(const MEDCouplingMesh *mesh) = 0;
121     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *aggregate(const MEDCouplingTimeDiscretization *other) const = 0;
122     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *aggregate(const std::vector<const MEDCouplingTimeDiscretization *>& other) const = 0;
123     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *meld(const MEDCouplingTimeDiscretization *other) const = 0;
124     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *dot(const MEDCouplingTimeDiscretization *other) const = 0;
125     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *crossProduct(const MEDCouplingTimeDiscretization *other) const = 0;
126     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *max(const MEDCouplingTimeDiscretization *other) const = 0;
127     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *min(const MEDCouplingTimeDiscretization *other) const = 0;
128     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *add(const MEDCouplingTimeDiscretization *other) const = 0;
129     MEDCOUPLING_EXPORT virtual void addEqual(const MEDCouplingTimeDiscretization *other) = 0;
130     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *substract(const MEDCouplingTimeDiscretization *other) const = 0;
131     MEDCOUPLING_EXPORT virtual void substractEqual(const MEDCouplingTimeDiscretization *other) = 0;
132     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *multiply(const MEDCouplingTimeDiscretization *other) const = 0;
133     MEDCOUPLING_EXPORT virtual void multiplyEqual(const MEDCouplingTimeDiscretization *other) = 0;
134     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *divide(const MEDCouplingTimeDiscretization *other) const = 0;
135     MEDCOUPLING_EXPORT virtual void divideEqual(const MEDCouplingTimeDiscretization *other) = 0;
136     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *pow(const MEDCouplingTimeDiscretization *other) const = 0;
137     MEDCOUPLING_EXPORT virtual void powEqual(const MEDCouplingTimeDiscretization *other) = 0;
138     MEDCOUPLING_EXPORT virtual void getTinySerializationIntInformation(std::vector<int>& tinyInfo) const;
139     MEDCOUPLING_EXPORT virtual void getTinySerializationDbleInformation(std::vector<double>& tinyInfo) const;
140     MEDCOUPLING_EXPORT virtual void getTinySerializationStrInformation(std::vector<std::string>& tinyInfo) const;
141     MEDCOUPLING_EXPORT virtual void resizeForUnserialization(const std::vector<int>& tinyInfoI, std::vector<DataArrayDouble *>& arrays);
142     MEDCOUPLING_EXPORT virtual void checkForUnserialization(const std::vector<int>& tinyInfoI, const std::vector<DataArrayDouble *>& arrays);
143     MEDCOUPLING_EXPORT virtual void finishUnserialization(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD, const std::vector<std::string>& tinyInfoS);
144     MEDCOUPLING_EXPORT virtual void getTinySerializationIntInformation2(std::vector<int>& tinyInfo) const = 0;
145     MEDCOUPLING_EXPORT virtual void getTinySerializationDbleInformation2(std::vector<double>& tinyInfo) const = 0;
146     MEDCOUPLING_EXPORT virtual void finishUnserialization2(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD) = 0;
147     MEDCOUPLING_EXPORT virtual void checkNoTimePresence() const = 0;
148     MEDCOUPLING_EXPORT virtual void checkTimePresence(double time) const = 0;
149     MEDCOUPLING_EXPORT virtual void setEndArray(DataArrayDouble *array, TimeLabel *owner);
150     MEDCOUPLING_EXPORT virtual void setArrays(const std::vector<DataArrayDouble *>& arrays, TimeLabel *owner);
151     MEDCOUPLING_EXPORT virtual const DataArrayDouble *getEndArray() const;
152     MEDCOUPLING_EXPORT virtual DataArrayDouble *getEndArray();
153     MEDCOUPLING_EXPORT virtual std::vector< const DataArrayDouble *> getArraysForTime(double time) const = 0;
154     MEDCOUPLING_EXPORT virtual void getValueForTime(double time, const std::vector<double>& vals, double *res) const = 0; 
155     MEDCOUPLING_EXPORT virtual void getArrays(std::vector<DataArrayDouble *>& arrays) const;
156     MEDCOUPLING_EXPORT virtual bool isBefore(const MEDCouplingTimeDiscretization *other) const;
157     MEDCOUPLING_EXPORT virtual bool isStrictlyBefore(const MEDCouplingTimeDiscretization *other) const;
158     MEDCOUPLING_EXPORT virtual double getEndTime(int& iteration, int& order) const = 0;
159     MEDCOUPLING_EXPORT virtual void setEndIteration(int it) = 0;
160     MEDCOUPLING_EXPORT virtual void setEndOrder(int order) = 0;
161     MEDCOUPLING_EXPORT virtual void setEndTimeValue(double time) = 0;
162     MEDCOUPLING_EXPORT virtual void setEndTime(double time, int iteration, int order) = 0;
163     MEDCOUPLING_EXPORT virtual void getValueOnTime(int eltId, double time, double *value) const = 0;
164     MEDCOUPLING_EXPORT virtual void getValueOnDiscTime(int eltId, int iteration, int order, double *value) const = 0;
165     //
166     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *doublyContractedProduct() const;
167     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *determinant() const;
168     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *eigenValues() const;
169     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *eigenVectors() const;
170     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *inverse() const;
171     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *trace() const;
172     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *deviator() const;
173     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *magnitude() const;
174     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *negate() const;
175     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *maxPerTuple() const;
176     MEDCOUPLING_EXPORT virtual MEDCouplingTimeDiscretization *keepSelectedComponents(const std::vector<int>& compoIds) const;
177     MEDCOUPLING_EXPORT virtual void setSelectedComponents(const MEDCouplingTimeDiscretization *other, const std::vector<int>& compoIds);
178     MEDCOUPLING_EXPORT virtual void changeNbOfComponents(int newNbOfComp, double dftValue);
179     MEDCOUPLING_EXPORT virtual void sortPerTuple(bool asc);
180     MEDCOUPLING_EXPORT virtual void setUniformValue(int nbOfTuple, int nbOfCompo, double value);
181     MEDCOUPLING_EXPORT virtual void setOrCreateUniformValueOnAllComponents(int nbOfTuple, double value);
182     MEDCOUPLING_EXPORT virtual void applyLin(double a, double b, int compoId);
183     MEDCOUPLING_EXPORT virtual void applyLin(double a, double b);
184     MEDCOUPLING_EXPORT virtual void applyFunc(int nbOfComp, FunctionToEvaluate func);
185     MEDCOUPLING_EXPORT virtual void applyFunc(int nbOfComp, const std::string& func);
186     MEDCOUPLING_EXPORT virtual void applyFuncCompo(int nbOfComp, const std::string& func);
187     MEDCOUPLING_EXPORT virtual void applyFuncNamedCompo(int nbOfComp, const std::vector<std::string>& varsOrder, const std::string& func);
188     MEDCOUPLING_EXPORT virtual void applyFunc(const std::string& func);
189     MEDCOUPLING_EXPORT virtual void applyFuncFast32(const std::string& func);
190     MEDCOUPLING_EXPORT virtual void applyFuncFast64(const std::string& func);
191     MEDCOUPLING_EXPORT virtual void fillFromAnalytic(const DataArrayDouble *loc, int nbOfComp, FunctionToEvaluate func);
192     MEDCOUPLING_EXPORT virtual void fillFromAnalytic(const DataArrayDouble *loc, int nbOfComp, const std::string& func);
193     MEDCOUPLING_EXPORT virtual void fillFromAnalyticCompo(const DataArrayDouble *loc, int nbOfComp, const std::string& func);
194     MEDCOUPLING_EXPORT virtual void fillFromAnalyticNamedCompo(const DataArrayDouble *loc, int nbOfComp, const std::vector<std::string>& varsOrder, const std::string& func);
195   };
196   
197   class MEDCouplingTimeDiscretizationInt : public MEDCouplingTimeDiscretizationTemplate<int>
198   {
199   public:
200     MEDCouplingTimeDiscretizationInt();
201     MEDCouplingTimeDiscretizationInt(const MEDCouplingTimeDiscretizationInt& other, bool deepCopy);
202     static MEDCouplingTimeDiscretizationInt *New(TypeOfTimeDiscretization type);
203     double getStartTime(int& iteration, int& order) const { return _tk.getAllInfo(iteration,order); }
204     void setStartIteration(int it) { _tk.setIteration(it); }
205     void setStartOrder(int order) { _tk.setOrder(order); }
206     void setStartTimeValue(double time) { _tk.setTimeValue(time); } 
207     void setStartTime(double time, int iteration, int order) { _tk.setAllInfo(time,iteration,order); }
208     std::string getStringRepr() const;
209     TypeOfTimeDiscretization getEnum() const { return DISCRETIZATION; }
210     MEDCouplingTimeDiscretizationInt *performCopyOrIncrRef(bool deepCopy) const;
211   private:
212     static const TypeOfTimeDiscretization DISCRETIZATION=ONE_TIME;
213     MEDCOUPLING_EXPORT static const char REPR[];
214   protected:
215     MEDCouplingTimeKeeper _tk;
216   };
217   
218   class MEDCouplingNoTimeLabel : public MEDCouplingTimeDiscretization
219   {
220   public:
221     MEDCOUPLING_EXPORT MEDCouplingNoTimeLabel();
222     MEDCOUPLING_EXPORT MEDCouplingNoTimeLabel(const MEDCouplingTimeDiscretization& other, bool deepCopy);
223     MEDCOUPLING_EXPORT std::string getStringRepr() const;
224     MEDCOUPLING_EXPORT TypeOfTimeDiscretization getEnum() const { return DISCRETIZATION; }
225     MEDCOUPLING_EXPORT void synchronizeTimeWith(const MEDCouplingMesh *mesh);
226     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *aggregate(const MEDCouplingTimeDiscretization *other) const;
227     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *aggregate(const std::vector<const MEDCouplingTimeDiscretization *>& other) const;
228     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *meld(const MEDCouplingTimeDiscretization *other) const;
229     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *dot(const MEDCouplingTimeDiscretization *other) const;
230     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *crossProduct(const MEDCouplingTimeDiscretization *other) const;
231     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *max(const MEDCouplingTimeDiscretization *other) const;
232     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *min(const MEDCouplingTimeDiscretization *other) const;
233     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *add(const MEDCouplingTimeDiscretization *other) const;
234     MEDCOUPLING_EXPORT void addEqual(const MEDCouplingTimeDiscretization *other);
235     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *substract(const MEDCouplingTimeDiscretization *other) const;
236     MEDCOUPLING_EXPORT void substractEqual(const MEDCouplingTimeDiscretization *other);
237     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *multiply(const MEDCouplingTimeDiscretization *other) const;
238     MEDCOUPLING_EXPORT void multiplyEqual(const MEDCouplingTimeDiscretization *other);
239     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *divide(const MEDCouplingTimeDiscretization *other) const;
240     MEDCOUPLING_EXPORT void divideEqual(const MEDCouplingTimeDiscretization *other);
241     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *pow(const MEDCouplingTimeDiscretization *other) const;
242     MEDCOUPLING_EXPORT void powEqual(const MEDCouplingTimeDiscretization *other);
243     MEDCOUPLING_EXPORT bool isEqualIfNotWhy(const MEDCouplingTimeDiscretization *other, double prec, std::string& reason) const;
244     MEDCOUPLING_EXPORT bool isEqualWithoutConsideringStr(const MEDCouplingTimeDiscretization *other, double prec) const;
245     MEDCOUPLING_EXPORT bool areCompatible(const MEDCouplingTimeDiscretizationTemplate<double> *other) const;
246     MEDCOUPLING_EXPORT bool areStrictlyCompatible(const MEDCouplingTimeDiscretizationTemplate<double> *other, std::string& reason) const;
247     MEDCOUPLING_EXPORT bool areStrictlyCompatibleForMul(const MEDCouplingTimeDiscretizationTemplate<double> *other) const;
248     MEDCOUPLING_EXPORT bool areStrictlyCompatibleForDiv(const MEDCouplingTimeDiscretizationTemplate<double> *other) const;
249     MEDCOUPLING_EXPORT bool areCompatibleForMeld(const MEDCouplingTimeDiscretization *other) const;
250     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *performCopyOrIncrRef(bool deepCopy) const;
251     MEDCOUPLING_EXPORT void checkNoTimePresence() const { }
252     MEDCOUPLING_EXPORT void checkTimePresence(double time) const;
253     MEDCOUPLING_EXPORT std::vector< const DataArrayDouble *> getArraysForTime(double time) const;
254     MEDCOUPLING_EXPORT void getValueForTime(double time, const std::vector<double>& vals, double *res) const;
255     MEDCOUPLING_EXPORT bool isBefore(const MEDCouplingTimeDiscretization *other) const;
256     MEDCOUPLING_EXPORT bool isStrictlyBefore(const MEDCouplingTimeDiscretization *other) const;
257     MEDCOUPLING_EXPORT double getStartTime(int& iteration, int& order) const;
258     MEDCOUPLING_EXPORT double getEndTime(int& iteration, int& order) const;
259     MEDCOUPLING_EXPORT void setStartIteration(int it);
260     MEDCOUPLING_EXPORT void setEndIteration(int it);
261     MEDCOUPLING_EXPORT void setStartOrder(int order);
262     MEDCOUPLING_EXPORT void setEndOrder(int order);
263     MEDCOUPLING_EXPORT void setStartTimeValue(double time);
264     MEDCOUPLING_EXPORT void setEndTimeValue(double time);
265     MEDCOUPLING_EXPORT void setStartTime(double time, int iteration, int order);
266     MEDCOUPLING_EXPORT void setEndTime(double time, int iteration, int order);
267     MEDCOUPLING_EXPORT void getValueOnTime(int eltId, double time, double *value) const;
268     MEDCOUPLING_EXPORT void getValueOnDiscTime(int eltId, int iteration, int order, double *value) const;
269     MEDCOUPLING_EXPORT void getTinySerializationIntInformation2(std::vector<int>& tinyInfo) const;
270     MEDCOUPLING_EXPORT void getTinySerializationDbleInformation2(std::vector<double>& tinyInfo) const;
271     MEDCOUPLING_EXPORT void finishUnserialization2(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD);
272   public:
273     static const TypeOfTimeDiscretization DISCRETIZATION=NO_TIME;
274     MEDCOUPLING_EXPORT static const char REPR[];
275   private:
276     static const char EXCEPTION_MSG[];
277   };
278
279   class MEDCouplingWithTimeStep : public MEDCouplingTimeDiscretization
280   {
281   protected:
282     MEDCOUPLING_EXPORT MEDCouplingWithTimeStep(const MEDCouplingWithTimeStep& other, bool deepCopy);
283   public:
284     MEDCOUPLING_EXPORT MEDCouplingWithTimeStep();
285     MEDCOUPLING_EXPORT std::string getStringRepr() const;
286     MEDCOUPLING_EXPORT void copyTinyAttrFrom(const MEDCouplingTimeDiscretizationTemplate<double>& other);
287     MEDCOUPLING_EXPORT TypeOfTimeDiscretization getEnum() const { return DISCRETIZATION; }
288     MEDCOUPLING_EXPORT void synchronizeTimeWith(const MEDCouplingMesh *mesh);
289     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *aggregate(const MEDCouplingTimeDiscretization *other) const;
290     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *aggregate(const std::vector<const MEDCouplingTimeDiscretization *>& other) const;
291     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *meld(const MEDCouplingTimeDiscretization *other) const;
292     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *dot(const MEDCouplingTimeDiscretization *other) const;
293     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *crossProduct(const MEDCouplingTimeDiscretization *other) const;
294     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *max(const MEDCouplingTimeDiscretization *other) const;
295     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *min(const MEDCouplingTimeDiscretization *other) const;
296     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *add(const MEDCouplingTimeDiscretization *other) const;
297     MEDCOUPLING_EXPORT void addEqual(const MEDCouplingTimeDiscretization *other);
298     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *substract(const MEDCouplingTimeDiscretization *other) const;
299     MEDCOUPLING_EXPORT void substractEqual(const MEDCouplingTimeDiscretization *other);
300     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *multiply(const MEDCouplingTimeDiscretization *other) const;
301     MEDCOUPLING_EXPORT void multiplyEqual(const MEDCouplingTimeDiscretization *other);
302     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *divide(const MEDCouplingTimeDiscretization *other) const;
303     MEDCOUPLING_EXPORT void divideEqual(const MEDCouplingTimeDiscretization *other);
304     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *pow(const MEDCouplingTimeDiscretization *other) const;
305     MEDCOUPLING_EXPORT void powEqual(const MEDCouplingTimeDiscretization *other);
306     MEDCOUPLING_EXPORT bool isEqualIfNotWhy(const MEDCouplingTimeDiscretization *other, double prec, std::string& reason) const;
307     MEDCOUPLING_EXPORT bool isEqualWithoutConsideringStr(const MEDCouplingTimeDiscretization *other, double prec) const;
308     MEDCOUPLING_EXPORT bool areCompatible(const MEDCouplingTimeDiscretizationTemplate<double> *other) const;
309     MEDCOUPLING_EXPORT bool areStrictlyCompatible(const MEDCouplingTimeDiscretizationTemplate<double> *other, std::string& reason) const;
310     MEDCOUPLING_EXPORT bool areStrictlyCompatibleForMul(const MEDCouplingTimeDiscretizationTemplate<double> *other) const;
311     MEDCOUPLING_EXPORT bool areStrictlyCompatibleForDiv(const MEDCouplingTimeDiscretizationTemplate<double> *other) const;
312     MEDCOUPLING_EXPORT bool areCompatibleForMeld(const MEDCouplingTimeDiscretization *other) const;
313     MEDCOUPLING_EXPORT void getTinySerializationIntInformation(std::vector<int>& tinyInfo) const;
314     MEDCOUPLING_EXPORT void getTinySerializationDbleInformation(std::vector<double>& tinyInfo) const;
315     MEDCOUPLING_EXPORT void finishUnserialization(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD, const std::vector<std::string>& tinyInfoS);
316     MEDCOUPLING_EXPORT void getTinySerializationIntInformation2(std::vector<int>& tinyInfo) const;
317     MEDCOUPLING_EXPORT void getTinySerializationDbleInformation2(std::vector<double>& tinyInfo) const;
318     MEDCOUPLING_EXPORT void finishUnserialization2(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD);
319     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *performCopyOrIncrRef(bool deepCopy) const;
320     MEDCOUPLING_EXPORT void checkNoTimePresence() const;
321     MEDCOUPLING_EXPORT void checkTimePresence(double time) const;
322     MEDCOUPLING_EXPORT void setStartTime(double time, int iteration, int order) { _tk.setAllInfo(time,iteration,order); }
323     MEDCOUPLING_EXPORT void setEndTime(double time, int iteration, int order) { _tk.setAllInfo(time,iteration,order); }
324     MEDCOUPLING_EXPORT double getStartTime(int& iteration, int& order) const { return _tk.getAllInfo(iteration,order); }
325     MEDCOUPLING_EXPORT double getEndTime(int& iteration, int& order) const { return _tk.getAllInfo(iteration,order); }
326     MEDCOUPLING_EXPORT void setStartIteration(int it) { _tk.setIteration(it); }
327     MEDCOUPLING_EXPORT void setEndIteration(int it) { _tk.setIteration(it); }
328     MEDCOUPLING_EXPORT void setStartOrder(int order) { _tk.setOrder(order); }
329     MEDCOUPLING_EXPORT void setEndOrder(int order) { _tk.setOrder(order); }
330     MEDCOUPLING_EXPORT void setStartTimeValue(double time) { _tk.setTimeValue(time); }
331     MEDCOUPLING_EXPORT void setEndTimeValue(double time) { _tk.setTimeValue(time); }
332     MEDCOUPLING_EXPORT std::vector< const DataArrayDouble *> getArraysForTime(double time) const;
333     MEDCOUPLING_EXPORT void getValueForTime(double time, const std::vector<double>& vals, double *res) const;
334     MEDCOUPLING_EXPORT void getValueOnTime(int eltId, double time, double *value) const;
335     MEDCOUPLING_EXPORT void getValueOnDiscTime(int eltId, int iteration, int order, double *value) const;
336   public:
337     static const TypeOfTimeDiscretization DISCRETIZATION=ONE_TIME;
338     MEDCOUPLING_EXPORT static const char REPR[];
339   private:
340     static const char EXCEPTION_MSG[];
341   protected:
342     MEDCouplingTimeKeeper _tk;
343   };
344
345   class MEDCouplingConstOnTimeInterval : public MEDCouplingTimeDiscretization
346   {
347   protected:
348     MEDCOUPLING_EXPORT MEDCouplingConstOnTimeInterval(const MEDCouplingConstOnTimeInterval& other, bool deepCopy);
349   public:
350     MEDCOUPLING_EXPORT MEDCouplingConstOnTimeInterval();
351     MEDCOUPLING_EXPORT void copyTinyAttrFrom(const MEDCouplingTimeDiscretizationTemplate<double>& other);
352     MEDCOUPLING_EXPORT void getTinySerializationIntInformation(std::vector<int>& tinyInfo) const;
353     MEDCOUPLING_EXPORT void getTinySerializationDbleInformation(std::vector<double>& tinyInfo) const;
354     MEDCOUPLING_EXPORT void finishUnserialization(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD, const std::vector<std::string>& tinyInfoS);
355     MEDCOUPLING_EXPORT void getTinySerializationIntInformation2(std::vector<int>& tinyInfo) const;
356     MEDCOUPLING_EXPORT void getTinySerializationDbleInformation2(std::vector<double>& tinyInfo) const;
357     MEDCOUPLING_EXPORT void finishUnserialization2(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD);
358     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *performCopyOrIncrRef(bool deepCopy) const;
359     MEDCOUPLING_EXPORT bool areCompatible(const MEDCouplingTimeDiscretizationTemplate<double> *other) const;
360     MEDCOUPLING_EXPORT bool areStrictlyCompatible(const MEDCouplingTimeDiscretizationTemplate<double> *other, std::string& reason) const;
361     MEDCOUPLING_EXPORT bool areStrictlyCompatibleForMul(const MEDCouplingTimeDiscretizationTemplate<double> *other) const;
362     MEDCOUPLING_EXPORT bool areStrictlyCompatibleForDiv(const MEDCouplingTimeDiscretizationTemplate<double> *other) const;
363     MEDCOUPLING_EXPORT bool areCompatibleForMeld(const MEDCouplingTimeDiscretization *other) const;
364     MEDCOUPLING_EXPORT bool isEqualIfNotWhy(const MEDCouplingTimeDiscretization *other, double prec, std::string& reason) const;
365     MEDCOUPLING_EXPORT bool isEqualWithoutConsideringStr(const MEDCouplingTimeDiscretization *other, double prec) const;
366     MEDCOUPLING_EXPORT std::vector< const DataArrayDouble *> getArraysForTime(double time) const;
367     MEDCOUPLING_EXPORT void getValueForTime(double time, const std::vector<double>& vals, double *res) const;
368     MEDCOUPLING_EXPORT void getValueOnTime(int eltId, double time, double *value) const;
369     MEDCOUPLING_EXPORT void getValueOnDiscTime(int eltId, int iteration, int order, double *value) const;
370     MEDCOUPLING_EXPORT TypeOfTimeDiscretization getEnum() const { return DISCRETIZATION; }
371     MEDCOUPLING_EXPORT void synchronizeTimeWith(const MEDCouplingMesh *mesh);
372     MEDCOUPLING_EXPORT std::string getStringRepr() const;
373     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *aggregate(const MEDCouplingTimeDiscretization *other) const;
374     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *aggregate(const std::vector<const MEDCouplingTimeDiscretization *>& other) const;
375     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *meld(const MEDCouplingTimeDiscretization *other) const;
376     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *dot(const MEDCouplingTimeDiscretization *other) const;
377     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *crossProduct(const MEDCouplingTimeDiscretization *other) const;
378     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *max(const MEDCouplingTimeDiscretization *other) const;
379     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *min(const MEDCouplingTimeDiscretization *other) const;
380     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *add(const MEDCouplingTimeDiscretization *other) const;
381     MEDCOUPLING_EXPORT void addEqual(const MEDCouplingTimeDiscretization *other);
382     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *substract(const MEDCouplingTimeDiscretization *other) const;
383     MEDCOUPLING_EXPORT void substractEqual(const MEDCouplingTimeDiscretization *other);
384     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *multiply(const MEDCouplingTimeDiscretization *other) const;
385     MEDCOUPLING_EXPORT void multiplyEqual(const MEDCouplingTimeDiscretization *other);
386     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *divide(const MEDCouplingTimeDiscretization *other) const;
387     MEDCOUPLING_EXPORT void divideEqual(const MEDCouplingTimeDiscretization *other);
388     MEDCouplingTimeDiscretization *pow(const MEDCouplingTimeDiscretization *other) const;
389     MEDCOUPLING_EXPORT void powEqual(const MEDCouplingTimeDiscretization *other);
390     MEDCOUPLING_EXPORT void setStartTime(double time, int iteration, int order) { _start.setAllInfo(time,iteration,order); }
391     MEDCOUPLING_EXPORT void setEndTime(double time, int iteration, int order) { _end.setAllInfo(time,iteration,order); }
392     MEDCOUPLING_EXPORT double getStartTime(int& iteration, int& order) const { return _start.getAllInfo(iteration,order); }
393     MEDCOUPLING_EXPORT double getEndTime(int& iteration, int& order) const { return _end.getAllInfo(iteration,order); }
394     MEDCOUPLING_EXPORT void setStartIteration(int it) { _start.setIteration(it); }
395     MEDCOUPLING_EXPORT void setEndIteration(int it) { _end.setIteration(it); }
396     MEDCOUPLING_EXPORT void setStartOrder(int order) { _start.setOrder(order); }
397     MEDCOUPLING_EXPORT void setEndOrder(int order) { _end.setOrder(order); }
398     MEDCOUPLING_EXPORT void setStartTimeValue(double time) { _start.setTimeValue(time); }
399     MEDCOUPLING_EXPORT void setEndTimeValue(double time) { _end.setTimeValue(time); }
400     MEDCOUPLING_EXPORT void checkNoTimePresence() const;
401     MEDCOUPLING_EXPORT void checkTimePresence(double time) const;
402   public:
403     static const TypeOfTimeDiscretization DISCRETIZATION=CONST_ON_TIME_INTERVAL;
404     MEDCOUPLING_EXPORT static const char REPR[];
405   private:
406     static const char EXCEPTION_MSG[];
407   protected:
408     MEDCouplingTimeKeeper _start;
409     MEDCouplingTimeKeeper _end;
410   };
411
412   class MEDCouplingTwoTimeSteps : public MEDCouplingTimeDiscretization
413   {
414   protected:
415     MEDCOUPLING_EXPORT MEDCouplingTwoTimeSteps(const MEDCouplingTwoTimeSteps& other, bool deepCopy);
416     MEDCOUPLING_EXPORT MEDCouplingTwoTimeSteps();
417     MEDCOUPLING_EXPORT ~MEDCouplingTwoTimeSteps();
418   public:
419     MEDCOUPLING_EXPORT void updateTime() const;
420     MEDCOUPLING_EXPORT void synchronizeTimeWith(const MEDCouplingMesh *mesh);
421     MEDCOUPLING_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
422     MEDCOUPLING_EXPORT std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
423     MEDCOUPLING_EXPORT void copyTinyAttrFrom(const MEDCouplingTimeDiscretizationTemplate<double>& other);
424     MEDCOUPLING_EXPORT void copyTinyStringsFrom(const MEDCouplingTimeDiscretizationTemplate<double>& other);
425     MEDCOUPLING_EXPORT const DataArrayDouble *getEndArray() const;
426     MEDCOUPLING_EXPORT DataArrayDouble *getEndArray();
427     MEDCOUPLING_EXPORT void checkConsistencyLight() const;
428     MEDCOUPLING_EXPORT bool isEqualIfNotWhy(const MEDCouplingTimeDiscretization *other, double prec, std::string& reason) const;
429     MEDCOUPLING_EXPORT bool isEqualWithoutConsideringStr(const MEDCouplingTimeDiscretization *other, double prec) const;
430     MEDCOUPLING_EXPORT void checkNoTimePresence() const;
431     MEDCOUPLING_EXPORT void checkTimePresence(double time) const;
432     MEDCOUPLING_EXPORT void getArrays(std::vector<DataArrayDouble *>& arrays) const;
433     MEDCOUPLING_EXPORT void setEndArray(DataArrayDouble *array, TimeLabel *owner);
434     MEDCOUPLING_EXPORT void setStartTime(double time, int iteration, int order) { _start_time=time; _start_iteration=iteration; _start_order=order; }
435     MEDCOUPLING_EXPORT void setEndTime(double time, int iteration, int order) { _end_time=time; _end_iteration=iteration; _end_order=order; }
436     MEDCOUPLING_EXPORT double getStartTime(int& iteration, int& order) const { iteration=_start_iteration; order=_start_order; return _start_time; }
437     MEDCOUPLING_EXPORT double getEndTime(int& iteration, int& order) const { iteration=_end_iteration; order=_end_order; return _end_time; }
438     MEDCOUPLING_EXPORT void setStartIteration(int it) { _start_iteration=it; }
439     MEDCOUPLING_EXPORT void setEndIteration(int it) { _end_iteration=it; }
440     MEDCOUPLING_EXPORT void setStartOrder(int order) { _start_order=order; }
441     MEDCOUPLING_EXPORT void setEndOrder(int order) { _end_order=order; }
442     MEDCOUPLING_EXPORT void setStartTimeValue(double time) { _start_time=time; }
443     MEDCOUPLING_EXPORT void setEndTimeValue(double time) { _end_time=time; }
444     MEDCOUPLING_EXPORT void getTinySerializationIntInformation(std::vector<int>& tinyInfo) const;
445     MEDCOUPLING_EXPORT void getTinySerializationDbleInformation(std::vector<double>& tinyInfo) const;
446     MEDCOUPLING_EXPORT void getTinySerializationStrInformation(std::vector<std::string>& tinyInfo) const;
447     MEDCOUPLING_EXPORT void resizeForUnserialization(const std::vector<int>& tinyInfoI, std::vector<DataArrayDouble *>& arrays);
448     MEDCOUPLING_EXPORT void checkForUnserialization(const std::vector<int>& tinyInfoI, const std::vector<DataArrayDouble *>& arrays);
449     MEDCOUPLING_EXPORT void finishUnserialization(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD, const std::vector<std::string>& tinyInfoS);
450     MEDCOUPLING_EXPORT void getTinySerializationIntInformation2(std::vector<int>& tinyInfo) const;
451     MEDCOUPLING_EXPORT void getTinySerializationDbleInformation2(std::vector<double>& tinyInfo) const;
452     MEDCOUPLING_EXPORT void finishUnserialization2(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD);
453     MEDCOUPLING_EXPORT std::vector< const DataArrayDouble *> getArraysForTime(double time) const;
454     MEDCOUPLING_EXPORT void setArrays(const std::vector<DataArrayDouble *>& arrays, TimeLabel *owner);
455   protected:
456     static const char EXCEPTION_MSG[];
457   protected:
458     double _start_time;
459     double _end_time;
460     int _start_iteration;
461     int _end_iteration;
462     int _start_order;
463     int _end_order;
464     DataArrayDouble *_end_array;
465   };
466
467   class MEDCouplingLinearTime : public MEDCouplingTwoTimeSteps
468   {
469   protected:
470     MEDCOUPLING_EXPORT MEDCouplingLinearTime(const MEDCouplingLinearTime& other, bool deepCopy);
471   public:
472     MEDCOUPLING_EXPORT MEDCouplingLinearTime();
473     MEDCOUPLING_EXPORT std::string getStringRepr() const;
474     MEDCOUPLING_EXPORT TypeOfTimeDiscretization getEnum() const { return DISCRETIZATION; }
475     MEDCOUPLING_EXPORT void checkConsistencyLight() const;
476     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *performCopyOrIncrRef(bool deepCopy) const;
477     MEDCOUPLING_EXPORT bool areCompatible(const MEDCouplingTimeDiscretizationTemplate<double> *other) const;
478     MEDCOUPLING_EXPORT bool areStrictlyCompatible(const MEDCouplingTimeDiscretizationTemplate<double> *other, std::string& reason) const;
479     MEDCOUPLING_EXPORT bool areStrictlyCompatibleForMul(const MEDCouplingTimeDiscretizationTemplate<double> *other) const;
480     MEDCOUPLING_EXPORT bool areStrictlyCompatibleForDiv(const MEDCouplingTimeDiscretizationTemplate<double> *other) const;
481     MEDCOUPLING_EXPORT bool areCompatibleForMeld(const MEDCouplingTimeDiscretization *other) const;
482     MEDCOUPLING_EXPORT void getValueForTime(double time, const std::vector<double>& vals, double *res) const;
483     MEDCOUPLING_EXPORT void getValueOnTime(int eltId, double time, double *value) const;
484     MEDCOUPLING_EXPORT void getValueOnDiscTime(int eltId, int iteration, int order, double *value) const;
485     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *aggregate(const MEDCouplingTimeDiscretization *other) const;
486     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *aggregate(const std::vector<const MEDCouplingTimeDiscretization *>& other) const;
487     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *meld(const MEDCouplingTimeDiscretization *other) const;
488     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *dot(const MEDCouplingTimeDiscretization *other) const;
489     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *crossProduct(const MEDCouplingTimeDiscretization *other) const;
490     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *max(const MEDCouplingTimeDiscretization *other) const;
491     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *min(const MEDCouplingTimeDiscretization *other) const;
492     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *add(const MEDCouplingTimeDiscretization *other) const;
493     MEDCOUPLING_EXPORT void addEqual(const MEDCouplingTimeDiscretization *other);
494     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *substract(const MEDCouplingTimeDiscretization *other) const;
495     MEDCOUPLING_EXPORT void substractEqual(const MEDCouplingTimeDiscretization *other);
496     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *multiply(const MEDCouplingTimeDiscretization *other) const;
497     MEDCOUPLING_EXPORT void multiplyEqual(const MEDCouplingTimeDiscretization *other);
498     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *divide(const MEDCouplingTimeDiscretization *other) const;
499     MEDCOUPLING_EXPORT void divideEqual(const MEDCouplingTimeDiscretization *other);
500     MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization *pow(const MEDCouplingTimeDiscretization *other) const;
501     MEDCOUPLING_EXPORT void powEqual(const MEDCouplingTimeDiscretization *other);
502   public:
503     static const TypeOfTimeDiscretization DISCRETIZATION=LINEAR_TIME;
504     MEDCOUPLING_EXPORT static const char REPR[];
505   };
506 }
507
508 #endif