1 // Copyright (C) 2007-2013 CEA/DEN, EDF R&D
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.
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.
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
17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
19 // Author : Anthony Geay (CEA/DEN)
21 #ifndef __PARAMEDMEM_MEDCOUPLINGTIMEDISCRETIZATION_HXX__
22 #define __PARAMEDMEM_MEDCOUPLINGTIMEDISCRETIZATION_HXX__
24 #include "MEDCoupling.hxx"
25 #include "MEDCouplingTimeLabel.hxx"
26 #include "MEDCouplingRefCountObject.hxx"
27 #include "InterpKernelException.hxx"
33 class MEDCouplingMesh;
34 class DataArrayDouble;
37 class MEDCOUPLING_EXPORT MEDCouplingTimeDiscretization : public TimeLabel
40 MEDCouplingTimeDiscretization();
41 MEDCouplingTimeDiscretization(const MEDCouplingTimeDiscretization& other, bool deepCpy);
43 void updateTime() const;
44 virtual std::size_t getHeapMemorySize() const;
45 static MEDCouplingTimeDiscretization *New(TypeOfTimeDiscretization type);
46 void setTimeUnit(const char *unit) { _time_unit=unit; }
47 const char *getTimeUnit() const { return _time_unit.c_str(); }
48 virtual void copyTinyAttrFrom(const MEDCouplingTimeDiscretization& other) throw(INTERP_KERNEL::Exception);
49 virtual void copyTinyStringsFrom(const MEDCouplingTimeDiscretization& other);
50 virtual void checkCoherency() const throw(INTERP_KERNEL::Exception);
51 virtual bool areCompatible(const MEDCouplingTimeDiscretization *other) const;
52 virtual bool areStrictlyCompatible(const MEDCouplingTimeDiscretization *other, std::string& reason) const;
53 virtual bool areStrictlyCompatibleForMul(const MEDCouplingTimeDiscretization *other) const;
54 virtual bool areStrictlyCompatibleForDiv(const MEDCouplingTimeDiscretization *other) const;
55 virtual bool areCompatibleForMeld(const MEDCouplingTimeDiscretization *other) const;
56 virtual bool isEqualIfNotWhy(const MEDCouplingTimeDiscretization *other, double prec, std::string& reason) const;
57 virtual bool isEqual(const MEDCouplingTimeDiscretization *other, double prec) const;
58 virtual bool isEqualWithoutConsideringStr(const MEDCouplingTimeDiscretization *other, double prec) const;
59 virtual MEDCouplingTimeDiscretization *buildNewTimeReprFromThis(TypeOfTimeDiscretization type, bool deepCpy) const;
60 virtual std::string getStringRepr() const = 0;
61 virtual TypeOfTimeDiscretization getEnum() const = 0;
62 virtual void synchronizeTimeWith(const MEDCouplingMesh *mesh) throw(INTERP_KERNEL::Exception) = 0;
63 virtual MEDCouplingTimeDiscretization *aggregate(const MEDCouplingTimeDiscretization *other) const = 0;
64 virtual MEDCouplingTimeDiscretization *aggregate(const std::vector<const MEDCouplingTimeDiscretization *>& other) const = 0;
65 virtual MEDCouplingTimeDiscretization *meld(const MEDCouplingTimeDiscretization *other) const = 0;
66 virtual MEDCouplingTimeDiscretization *dot(const MEDCouplingTimeDiscretization *other) const = 0;
67 virtual MEDCouplingTimeDiscretization *crossProduct(const MEDCouplingTimeDiscretization *other) const = 0;
68 virtual MEDCouplingTimeDiscretization *max(const MEDCouplingTimeDiscretization *other) const = 0;
69 virtual MEDCouplingTimeDiscretization *min(const MEDCouplingTimeDiscretization *other) const = 0;
70 virtual MEDCouplingTimeDiscretization *add(const MEDCouplingTimeDiscretization *other) const = 0;
71 virtual void addEqual(const MEDCouplingTimeDiscretization *other) = 0;
72 virtual MEDCouplingTimeDiscretization *substract(const MEDCouplingTimeDiscretization *other) const = 0;
73 virtual void substractEqual(const MEDCouplingTimeDiscretization *other) = 0;
74 virtual MEDCouplingTimeDiscretization *multiply(const MEDCouplingTimeDiscretization *other) const = 0;
75 virtual void multiplyEqual(const MEDCouplingTimeDiscretization *other) = 0;
76 virtual MEDCouplingTimeDiscretization *divide(const MEDCouplingTimeDiscretization *other) const = 0;
77 virtual void divideEqual(const MEDCouplingTimeDiscretization *other) = 0;
78 virtual void getTinySerializationIntInformation(std::vector<int>& tinyInfo) const;
79 virtual void getTinySerializationDbleInformation(std::vector<double>& tinyInfo) const;
80 virtual void getTinySerializationStrInformation(std::vector<std::string>& tinyInfo) const;
81 virtual void resizeForUnserialization(const std::vector<int>& tinyInfoI, std::vector<DataArrayDouble *>& arrays);
82 virtual void finishUnserialization(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD, const std::vector<std::string>& tinyInfoS);
83 virtual void getTinySerializationIntInformation2(std::vector<int>& tinyInfo) const = 0;
84 virtual void getTinySerializationDbleInformation2(std::vector<double>& tinyInfo) const = 0;
85 virtual void finishUnserialization2(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD) = 0;
86 virtual MEDCouplingTimeDiscretization *performCpy(bool deepCpy) const = 0;
87 void setTimeTolerance(double val) { _time_tolerance=val; }
88 double getTimeTolerance() const { return _time_tolerance; }
89 virtual void checkNoTimePresence() const throw(INTERP_KERNEL::Exception) = 0;
90 virtual void checkTimePresence(double time) const throw(INTERP_KERNEL::Exception) = 0;
91 virtual void setArray(DataArrayDouble *array, TimeLabel *owner);
92 virtual void setEndArray(DataArrayDouble *array, TimeLabel *owner);
93 virtual void setArrays(const std::vector<DataArrayDouble *>& arrays, TimeLabel *owner) throw(INTERP_KERNEL::Exception);
94 DataArrayDouble *getArray() { return _array; }
95 const DataArrayDouble *getArray() const { return _array; }
96 virtual const DataArrayDouble *getEndArray() const;
97 virtual DataArrayDouble *getEndArray();
98 virtual std::vector< const DataArrayDouble *> getArraysForTime(double time) const throw(INTERP_KERNEL::Exception) = 0;
99 virtual void getValueForTime(double time, const std::vector<double>& vals, double *res) const = 0;
100 virtual void getArrays(std::vector<DataArrayDouble *>& arrays) const;
101 virtual bool isBefore(const MEDCouplingTimeDiscretization *other) const throw(INTERP_KERNEL::Exception);
102 virtual bool isStrictlyBefore(const MEDCouplingTimeDiscretization *other) const throw(INTERP_KERNEL::Exception);
103 double getTime(int& iteration, int& order) const throw(INTERP_KERNEL::Exception) { return getStartTime(iteration,order); }
104 virtual double getStartTime(int& iteration, int& order) const throw(INTERP_KERNEL::Exception) = 0;
105 virtual double getEndTime(int& iteration, int& order) const throw(INTERP_KERNEL::Exception) = 0;
106 void setTime(double time, int iteration, int order) throw(INTERP_KERNEL::Exception) { setStartTime(time,iteration,order); }
107 void setIteration(int it) throw(INTERP_KERNEL::Exception) { setStartIteration(it); }
108 void setOrder(int order) throw(INTERP_KERNEL::Exception) { setStartOrder(order); }
109 void setTimeValue(double val) throw(INTERP_KERNEL::Exception) { setStartTimeValue(val); }
110 virtual void setStartIteration(int it) throw(INTERP_KERNEL::Exception) = 0;
111 virtual void setEndIteration(int it) throw(INTERP_KERNEL::Exception) = 0;
112 virtual void setStartOrder(int order) throw(INTERP_KERNEL::Exception) = 0;
113 virtual void setEndOrder(int order) throw(INTERP_KERNEL::Exception) = 0;
114 virtual void setStartTimeValue(double time) throw(INTERP_KERNEL::Exception) = 0;
115 virtual void setEndTimeValue(double time) throw(INTERP_KERNEL::Exception) = 0;
116 virtual void setStartTime(double time, int iteration, int order) throw(INTERP_KERNEL::Exception) = 0;
117 virtual void setEndTime(double time, int iteration, int order) throw(INTERP_KERNEL::Exception) = 0;
118 virtual void getValueOnTime(int eltId, double time, double *value) const throw(INTERP_KERNEL::Exception) = 0;
119 virtual void getValueOnDiscTime(int eltId, int iteration, int order, double *value) const throw(INTERP_KERNEL::Exception) = 0;
121 virtual MEDCouplingTimeDiscretization *doublyContractedProduct() const throw(INTERP_KERNEL::Exception);
122 virtual MEDCouplingTimeDiscretization *determinant() const throw(INTERP_KERNEL::Exception);
123 virtual MEDCouplingTimeDiscretization *eigenValues() const throw(INTERP_KERNEL::Exception);
124 virtual MEDCouplingTimeDiscretization *eigenVectors() const throw(INTERP_KERNEL::Exception);
125 virtual MEDCouplingTimeDiscretization *inverse() const throw(INTERP_KERNEL::Exception);
126 virtual MEDCouplingTimeDiscretization *trace() const throw(INTERP_KERNEL::Exception);
127 virtual MEDCouplingTimeDiscretization *deviator() const throw(INTERP_KERNEL::Exception);
128 virtual MEDCouplingTimeDiscretization *magnitude() const throw(INTERP_KERNEL::Exception);
129 virtual MEDCouplingTimeDiscretization *maxPerTuple() const throw(INTERP_KERNEL::Exception);
130 virtual MEDCouplingTimeDiscretization *keepSelectedComponents(const std::vector<int>& compoIds) const throw(INTERP_KERNEL::Exception);
131 virtual void setSelectedComponents(const MEDCouplingTimeDiscretization *other, const std::vector<int>& compoIds) throw(INTERP_KERNEL::Exception);
132 virtual void changeNbOfComponents(int newNbOfComp, double dftValue) throw(INTERP_KERNEL::Exception);
133 virtual void sortPerTuple(bool asc) throw(INTERP_KERNEL::Exception);
134 virtual void setUniformValue(int nbOfTuple, int nbOfCompo, double value);
135 virtual void applyLin(double a, double b, int compoId);
136 virtual void applyFunc(int nbOfComp, FunctionToEvaluate func);
137 virtual void applyFunc(int nbOfComp, const char *func);
138 virtual void applyFunc2(int nbOfComp, const char *func);
139 virtual void applyFunc3(int nbOfComp, const std::vector<std::string>& varsOrder, const char *func);
140 virtual void applyFunc(const char *func);
141 virtual void applyFuncFast32(const char *func);
142 virtual void applyFuncFast64(const char *func);
143 virtual void fillFromAnalytic(const DataArrayDouble *loc, int nbOfComp, FunctionToEvaluate func) throw(INTERP_KERNEL::Exception);
144 virtual void fillFromAnalytic(const DataArrayDouble *loc, int nbOfComp, const char *func) throw(INTERP_KERNEL::Exception);
145 virtual void fillFromAnalytic2(const DataArrayDouble *loc, int nbOfComp, const char *func) throw(INTERP_KERNEL::Exception);
146 virtual void fillFromAnalytic3(const DataArrayDouble *loc, int nbOfComp, const std::vector<std::string>& varsOrder, const char *func) throw(INTERP_KERNEL::Exception);
148 virtual ~MEDCouplingTimeDiscretization();
150 std::string _time_unit;
151 double _time_tolerance;
152 DataArrayDouble *_array;
154 static const double TIME_TOLERANCE_DFT;
157 class MEDCOUPLING_EXPORT MEDCouplingNoTimeLabel : public MEDCouplingTimeDiscretization
160 MEDCouplingNoTimeLabel();
161 MEDCouplingNoTimeLabel(const MEDCouplingTimeDiscretization& other, bool deepCpy);
162 std::string getStringRepr() const;
163 TypeOfTimeDiscretization getEnum() const { return DISCRETIZATION; }
164 void synchronizeTimeWith(const MEDCouplingMesh *mesh) throw(INTERP_KERNEL::Exception);
165 MEDCouplingTimeDiscretization *aggregate(const MEDCouplingTimeDiscretization *other) const;
166 MEDCouplingTimeDiscretization *aggregate(const std::vector<const MEDCouplingTimeDiscretization *>& other) const;
167 MEDCouplingTimeDiscretization *meld(const MEDCouplingTimeDiscretization *other) const;
168 MEDCouplingTimeDiscretization *dot(const MEDCouplingTimeDiscretization *other) const;
169 MEDCouplingTimeDiscretization *crossProduct(const MEDCouplingTimeDiscretization *other) const;
170 MEDCouplingTimeDiscretization *max(const MEDCouplingTimeDiscretization *other) const;
171 MEDCouplingTimeDiscretization *min(const MEDCouplingTimeDiscretization *other) const;
172 MEDCouplingTimeDiscretization *add(const MEDCouplingTimeDiscretization *other) const;
173 void addEqual(const MEDCouplingTimeDiscretization *other);
174 MEDCouplingTimeDiscretization *substract(const MEDCouplingTimeDiscretization *other) const;
175 void substractEqual(const MEDCouplingTimeDiscretization *other);
176 MEDCouplingTimeDiscretization *multiply(const MEDCouplingTimeDiscretization *other) const;
177 void multiplyEqual(const MEDCouplingTimeDiscretization *other);
178 MEDCouplingTimeDiscretization *divide(const MEDCouplingTimeDiscretization *other) const;
179 void divideEqual(const MEDCouplingTimeDiscretization *other);
180 bool isEqualIfNotWhy(const MEDCouplingTimeDiscretization *other, double prec, std::string& reason) const;
181 bool isEqualWithoutConsideringStr(const MEDCouplingTimeDiscretization *other, double prec) const;
182 bool areCompatible(const MEDCouplingTimeDiscretization *other) const;
183 bool areStrictlyCompatible(const MEDCouplingTimeDiscretization *other, std::string& reason) const;
184 bool areStrictlyCompatibleForMul(const MEDCouplingTimeDiscretization *other) const;
185 bool areStrictlyCompatibleForDiv(const MEDCouplingTimeDiscretization *other) const;
186 bool areCompatibleForMeld(const MEDCouplingTimeDiscretization *other) const;
187 MEDCouplingTimeDiscretization *performCpy(bool deepCpy) const;
188 void checkNoTimePresence() const throw(INTERP_KERNEL::Exception) { }
189 void checkTimePresence(double time) const throw(INTERP_KERNEL::Exception);
190 std::vector< const DataArrayDouble *> getArraysForTime(double time) const throw(INTERP_KERNEL::Exception);
191 void getValueForTime(double time, const std::vector<double>& vals, double *res) const;
192 bool isBefore(const MEDCouplingTimeDiscretization *other) const throw(INTERP_KERNEL::Exception);
193 bool isStrictlyBefore(const MEDCouplingTimeDiscretization *other) const throw(INTERP_KERNEL::Exception);
194 double getStartTime(int& iteration, int& order) const throw(INTERP_KERNEL::Exception);
195 double getEndTime(int& iteration, int& order) const throw(INTERP_KERNEL::Exception);
196 void setStartIteration(int it) throw(INTERP_KERNEL::Exception);
197 void setEndIteration(int it) throw(INTERP_KERNEL::Exception);
198 void setStartOrder(int order) throw(INTERP_KERNEL::Exception);
199 void setEndOrder(int order) throw(INTERP_KERNEL::Exception);
200 void setStartTimeValue(double time) throw(INTERP_KERNEL::Exception);
201 void setEndTimeValue(double time) throw(INTERP_KERNEL::Exception);
202 void setStartTime(double time, int iteration, int order) throw(INTERP_KERNEL::Exception);
203 void setEndTime(double time, int iteration, int order) throw(INTERP_KERNEL::Exception);
204 void getValueOnTime(int eltId, double time, double *value) const throw(INTERP_KERNEL::Exception);
205 void getValueOnDiscTime(int eltId, int iteration, int order, double *value) const throw(INTERP_KERNEL::Exception);
206 void getTinySerializationIntInformation2(std::vector<int>& tinyInfo) const;
207 void getTinySerializationDbleInformation2(std::vector<double>& tinyInfo) const;
208 void finishUnserialization2(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD);
210 static const TypeOfTimeDiscretization DISCRETIZATION=NO_TIME;
211 static const char REPR[];
213 static const char EXCEPTION_MSG[];
216 class MEDCOUPLING_EXPORT MEDCouplingWithTimeStep : public MEDCouplingTimeDiscretization
219 MEDCouplingWithTimeStep(const MEDCouplingWithTimeStep& other, bool deepCpy);
221 MEDCouplingWithTimeStep();
222 std::string getStringRepr() const;
223 void copyTinyAttrFrom(const MEDCouplingTimeDiscretization& other) throw(INTERP_KERNEL::Exception);
224 TypeOfTimeDiscretization getEnum() const { return DISCRETIZATION; }
225 void synchronizeTimeWith(const MEDCouplingMesh *mesh) throw(INTERP_KERNEL::Exception);
226 MEDCouplingTimeDiscretization *aggregate(const MEDCouplingTimeDiscretization *other) const;
227 MEDCouplingTimeDiscretization *aggregate(const std::vector<const MEDCouplingTimeDiscretization *>& other) const;
228 MEDCouplingTimeDiscretization *meld(const MEDCouplingTimeDiscretization *other) const;
229 MEDCouplingTimeDiscretization *dot(const MEDCouplingTimeDiscretization *other) const;
230 MEDCouplingTimeDiscretization *crossProduct(const MEDCouplingTimeDiscretization *other) const;
231 MEDCouplingTimeDiscretization *max(const MEDCouplingTimeDiscretization *other) const;
232 MEDCouplingTimeDiscretization *min(const MEDCouplingTimeDiscretization *other) const;
233 MEDCouplingTimeDiscretization *add(const MEDCouplingTimeDiscretization *other) const;
234 void addEqual(const MEDCouplingTimeDiscretization *other);
235 MEDCouplingTimeDiscretization *substract(const MEDCouplingTimeDiscretization *other) const;
236 void substractEqual(const MEDCouplingTimeDiscretization *other);
237 MEDCouplingTimeDiscretization *multiply(const MEDCouplingTimeDiscretization *other) const;
238 void multiplyEqual(const MEDCouplingTimeDiscretization *other);
239 MEDCouplingTimeDiscretization *divide(const MEDCouplingTimeDiscretization *other) const;
240 void divideEqual(const MEDCouplingTimeDiscretization *other);
241 bool isEqualIfNotWhy(const MEDCouplingTimeDiscretization *other, double prec, std::string& reason) const;
242 bool isEqualWithoutConsideringStr(const MEDCouplingTimeDiscretization *other, double prec) const;
243 bool areCompatible(const MEDCouplingTimeDiscretization *other) const;
244 bool areStrictlyCompatible(const MEDCouplingTimeDiscretization *other, std::string& reason) const;
245 bool areStrictlyCompatibleForMul(const MEDCouplingTimeDiscretization *other) const;
246 bool areStrictlyCompatibleForDiv(const MEDCouplingTimeDiscretization *other) const;
247 bool areCompatibleForMeld(const MEDCouplingTimeDiscretization *other) const;
248 void getTinySerializationIntInformation(std::vector<int>& tinyInfo) const;
249 void getTinySerializationDbleInformation(std::vector<double>& tinyInfo) const;
250 void finishUnserialization(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD, const std::vector<std::string>& tinyInfoS);
251 void getTinySerializationIntInformation2(std::vector<int>& tinyInfo) const;
252 void getTinySerializationDbleInformation2(std::vector<double>& tinyInfo) const;
253 void finishUnserialization2(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD);
254 MEDCouplingTimeDiscretization *performCpy(bool deepCpy) const;
255 void checkNoTimePresence() const throw(INTERP_KERNEL::Exception);
256 void checkTimePresence(double time) const throw(INTERP_KERNEL::Exception);
257 void setStartTime(double time, int iteration, int order) throw(INTERP_KERNEL::Exception) { _time=time; _iteration=iteration; _order=order; }
258 void setEndTime(double time, int iteration, int order) throw(INTERP_KERNEL::Exception) { _time=time; _iteration=iteration; _order=order; }
259 double getStartTime(int& iteration, int& order) const throw(INTERP_KERNEL::Exception) { iteration=_iteration; order=_order; return _time; }
260 double getEndTime(int& iteration, int& order) const throw(INTERP_KERNEL::Exception) { iteration=_iteration; order=_order; return _time; }
261 void setStartIteration(int it) throw(INTERP_KERNEL::Exception) { _iteration=it; }
262 void setEndIteration(int it) throw(INTERP_KERNEL::Exception) { _iteration=it; }
263 void setStartOrder(int order) throw(INTERP_KERNEL::Exception) { _order=order; }
264 void setEndOrder(int order) throw(INTERP_KERNEL::Exception) { _order=order; }
265 void setStartTimeValue(double time) throw(INTERP_KERNEL::Exception) { _time=time; }
266 void setEndTimeValue(double time) throw(INTERP_KERNEL::Exception) { _time=time; }
267 std::vector< const DataArrayDouble *> getArraysForTime(double time) const throw(INTERP_KERNEL::Exception);
268 void getValueForTime(double time, const std::vector<double>& vals, double *res) const;
269 void getValueOnTime(int eltId, double time, double *value) const throw(INTERP_KERNEL::Exception);
270 void getValueOnDiscTime(int eltId, int iteration, int order, double *value) const throw(INTERP_KERNEL::Exception);
272 static const TypeOfTimeDiscretization DISCRETIZATION=ONE_TIME;
273 static const char REPR[];
275 static const char EXCEPTION_MSG[];
282 class MEDCOUPLING_EXPORT MEDCouplingConstOnTimeInterval : public MEDCouplingTimeDiscretization
285 MEDCouplingConstOnTimeInterval(const MEDCouplingConstOnTimeInterval& other, bool deepCpy);
287 MEDCouplingConstOnTimeInterval();
288 void copyTinyAttrFrom(const MEDCouplingTimeDiscretization& other) throw(INTERP_KERNEL::Exception);
289 void getTinySerializationIntInformation(std::vector<int>& tinyInfo) const;
290 void getTinySerializationDbleInformation(std::vector<double>& tinyInfo) const;
291 void finishUnserialization(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD, const std::vector<std::string>& tinyInfoS);
292 void getTinySerializationIntInformation2(std::vector<int>& tinyInfo) const;
293 void getTinySerializationDbleInformation2(std::vector<double>& tinyInfo) const;
294 void finishUnserialization2(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD);
295 MEDCouplingTimeDiscretization *performCpy(bool deepCpy) const;
296 bool areCompatible(const MEDCouplingTimeDiscretization *other) const;
297 bool areStrictlyCompatible(const MEDCouplingTimeDiscretization *other, std::string& reason) const;
298 bool areStrictlyCompatibleForMul(const MEDCouplingTimeDiscretization *other) const;
299 bool areStrictlyCompatibleForDiv(const MEDCouplingTimeDiscretization *other) const;
300 bool areCompatibleForMeld(const MEDCouplingTimeDiscretization *other) const;
301 bool isEqualIfNotWhy(const MEDCouplingTimeDiscretization *other, double prec, std::string& reason) const;
302 bool isEqualWithoutConsideringStr(const MEDCouplingTimeDiscretization *other, double prec) const;
303 std::vector< const DataArrayDouble *> getArraysForTime(double time) const throw(INTERP_KERNEL::Exception);
304 void getValueForTime(double time, const std::vector<double>& vals, double *res) const;
305 void getValueOnTime(int eltId, double time, double *value) const throw(INTERP_KERNEL::Exception);
306 void getValueOnDiscTime(int eltId, int iteration, int order, double *value) const throw(INTERP_KERNEL::Exception);
307 TypeOfTimeDiscretization getEnum() const { return DISCRETIZATION; }
308 void synchronizeTimeWith(const MEDCouplingMesh *mesh) throw(INTERP_KERNEL::Exception);
309 std::string getStringRepr() const;
310 MEDCouplingTimeDiscretization *aggregate(const MEDCouplingTimeDiscretization *other) const;
311 MEDCouplingTimeDiscretization *aggregate(const std::vector<const MEDCouplingTimeDiscretization *>& other) const;
312 MEDCouplingTimeDiscretization *meld(const MEDCouplingTimeDiscretization *other) const;
313 MEDCouplingTimeDiscretization *dot(const MEDCouplingTimeDiscretization *other) const;
314 MEDCouplingTimeDiscretization *crossProduct(const MEDCouplingTimeDiscretization *other) const;
315 MEDCouplingTimeDiscretization *max(const MEDCouplingTimeDiscretization *other) const;
316 MEDCouplingTimeDiscretization *min(const MEDCouplingTimeDiscretization *other) const;
317 MEDCouplingTimeDiscretization *add(const MEDCouplingTimeDiscretization *other) const;
318 void addEqual(const MEDCouplingTimeDiscretization *other);
319 MEDCouplingTimeDiscretization *substract(const MEDCouplingTimeDiscretization *other) const;
320 void substractEqual(const MEDCouplingTimeDiscretization *other);
321 MEDCouplingTimeDiscretization *multiply(const MEDCouplingTimeDiscretization *other) const;
322 void multiplyEqual(const MEDCouplingTimeDiscretization *other);
323 MEDCouplingTimeDiscretization *divide(const MEDCouplingTimeDiscretization *other) const;
324 void divideEqual(const MEDCouplingTimeDiscretization *other);
325 void setStartTime(double time, int iteration, int order) throw(INTERP_KERNEL::Exception) { _start_time=time; _start_iteration=iteration; _start_order=order; }
326 void setEndTime(double time, int iteration, int order) throw(INTERP_KERNEL::Exception) { _end_time=time; _end_iteration=iteration; _end_order=order; }
327 double getStartTime(int& iteration, int& order) const throw(INTERP_KERNEL::Exception) { iteration=_start_iteration; order=_start_order; return _start_time; }
328 double getEndTime(int& iteration, int& order) const throw(INTERP_KERNEL::Exception) { iteration=_end_iteration; order=_end_order; return _end_time; }
329 void setStartIteration(int it) throw(INTERP_KERNEL::Exception) { _start_iteration=it; }
330 void setEndIteration(int it) throw(INTERP_KERNEL::Exception) { _end_iteration=it; }
331 void setStartOrder(int order) throw(INTERP_KERNEL::Exception) { _start_order=order; }
332 void setEndOrder(int order) throw(INTERP_KERNEL::Exception) { _end_order=order; }
333 void setStartTimeValue(double time) throw(INTERP_KERNEL::Exception) { _start_time=time; }
334 void setEndTimeValue(double time) throw(INTERP_KERNEL::Exception) { _end_time=time; }
335 void checkNoTimePresence() const throw(INTERP_KERNEL::Exception);
336 void checkTimePresence(double time) const throw(INTERP_KERNEL::Exception);
338 static const TypeOfTimeDiscretization DISCRETIZATION=CONST_ON_TIME_INTERVAL;
339 static const char REPR[];
341 static const char EXCEPTION_MSG[];
345 int _start_iteration;
351 class MEDCOUPLING_EXPORT MEDCouplingTwoTimeSteps : public MEDCouplingTimeDiscretization
354 MEDCouplingTwoTimeSteps(const MEDCouplingTwoTimeSteps& other, bool deepCpy);
355 MEDCouplingTwoTimeSteps();
356 ~MEDCouplingTwoTimeSteps();
358 void updateTime() const;
359 void synchronizeTimeWith(const MEDCouplingMesh *mesh) throw(INTERP_KERNEL::Exception);
360 std::size_t getHeapMemorySize() const;
361 void copyTinyAttrFrom(const MEDCouplingTimeDiscretization& other) throw(INTERP_KERNEL::Exception);
362 void copyTinyStringsFrom(const MEDCouplingTimeDiscretization& other);
363 const DataArrayDouble *getEndArray() const;
364 DataArrayDouble *getEndArray();
365 void checkCoherency() const throw(INTERP_KERNEL::Exception);
366 bool isEqualIfNotWhy(const MEDCouplingTimeDiscretization *other, double prec, std::string& reason) const;
367 bool isEqualWithoutConsideringStr(const MEDCouplingTimeDiscretization *other, double prec) const;
368 void checkNoTimePresence() const throw(INTERP_KERNEL::Exception);
369 void checkTimePresence(double time) const throw(INTERP_KERNEL::Exception);
370 void getArrays(std::vector<DataArrayDouble *>& arrays) const;
371 void setEndArray(DataArrayDouble *array, TimeLabel *owner);
372 void setStartTime(double time, int iteration, int order) throw(INTERP_KERNEL::Exception) { _start_time=time; _start_iteration=iteration; _start_order=order; }
373 void setEndTime(double time, int iteration, int order) throw(INTERP_KERNEL::Exception) { _end_time=time; _end_iteration=iteration; _end_order=order; }
374 double getStartTime(int& iteration, int& order) const throw(INTERP_KERNEL::Exception) { iteration=_start_iteration; order=_start_order; return _start_time; }
375 double getEndTime(int& iteration, int& order) const throw(INTERP_KERNEL::Exception) { iteration=_end_iteration; order=_end_order; return _end_time; }
376 void setStartIteration(int it) throw(INTERP_KERNEL::Exception) { _start_iteration=it; }
377 void setEndIteration(int it) throw(INTERP_KERNEL::Exception) { _end_iteration=it; }
378 void setStartOrder(int order) throw(INTERP_KERNEL::Exception) { _start_order=order; }
379 void setEndOrder(int order) throw(INTERP_KERNEL::Exception) { _end_order=order; }
380 void setStartTimeValue(double time) throw(INTERP_KERNEL::Exception) { _start_time=time; }
381 void setEndTimeValue(double time) throw(INTERP_KERNEL::Exception) { _end_time=time; }
382 void getTinySerializationIntInformation(std::vector<int>& tinyInfo) const;
383 void getTinySerializationDbleInformation(std::vector<double>& tinyInfo) const;
384 void getTinySerializationStrInformation(std::vector<std::string>& tinyInfo) const;
385 void resizeForUnserialization(const std::vector<int>& tinyInfoI, std::vector<DataArrayDouble *>& arrays);
386 void finishUnserialization(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD, const std::vector<std::string>& tinyInfoS);
387 void getTinySerializationIntInformation2(std::vector<int>& tinyInfo) const;
388 void getTinySerializationDbleInformation2(std::vector<double>& tinyInfo) const;
389 void finishUnserialization2(const std::vector<int>& tinyInfoI, const std::vector<double>& tinyInfoD);
390 std::vector< const DataArrayDouble *> getArraysForTime(double time) const throw(INTERP_KERNEL::Exception);
391 void setArrays(const std::vector<DataArrayDouble *>& arrays, TimeLabel *owner) throw(INTERP_KERNEL::Exception);
393 static const char EXCEPTION_MSG[];
397 int _start_iteration;
401 DataArrayDouble *_end_array;
404 class MEDCOUPLING_EXPORT MEDCouplingLinearTime : public MEDCouplingTwoTimeSteps
407 MEDCouplingLinearTime(const MEDCouplingLinearTime& other, bool deepCpy);
409 MEDCouplingLinearTime();
410 std::string getStringRepr() const;
411 TypeOfTimeDiscretization getEnum() const { return DISCRETIZATION; }
412 void checkCoherency() const throw(INTERP_KERNEL::Exception);
413 MEDCouplingTimeDiscretization *performCpy(bool deepCpy) const;
414 bool areCompatible(const MEDCouplingTimeDiscretization *other) const;
415 bool areStrictlyCompatible(const MEDCouplingTimeDiscretization *other, std::string& reason) const;
416 bool areStrictlyCompatibleForMul(const MEDCouplingTimeDiscretization *other) const;
417 bool areStrictlyCompatibleForDiv(const MEDCouplingTimeDiscretization *other) const;
418 bool areCompatibleForMeld(const MEDCouplingTimeDiscretization *other) const;
419 void getValueForTime(double time, const std::vector<double>& vals, double *res) const;
420 void getValueOnTime(int eltId, double time, double *value) const throw(INTERP_KERNEL::Exception);
421 void getValueOnDiscTime(int eltId, int iteration, int order, double *value) const throw(INTERP_KERNEL::Exception);
422 MEDCouplingTimeDiscretization *aggregate(const MEDCouplingTimeDiscretization *other) const;
423 MEDCouplingTimeDiscretization *aggregate(const std::vector<const MEDCouplingTimeDiscretization *>& other) const;
424 MEDCouplingTimeDiscretization *meld(const MEDCouplingTimeDiscretization *other) const;
425 MEDCouplingTimeDiscretization *dot(const MEDCouplingTimeDiscretization *other) const;
426 MEDCouplingTimeDiscretization *crossProduct(const MEDCouplingTimeDiscretization *other) const;
427 MEDCouplingTimeDiscretization *max(const MEDCouplingTimeDiscretization *other) const;
428 MEDCouplingTimeDiscretization *min(const MEDCouplingTimeDiscretization *other) const;
429 MEDCouplingTimeDiscretization *add(const MEDCouplingTimeDiscretization *other) const;
430 void addEqual(const MEDCouplingTimeDiscretization *other);
431 MEDCouplingTimeDiscretization *substract(const MEDCouplingTimeDiscretization *other) const;
432 void substractEqual(const MEDCouplingTimeDiscretization *other);
433 MEDCouplingTimeDiscretization *multiply(const MEDCouplingTimeDiscretization *other) const;
434 void multiplyEqual(const MEDCouplingTimeDiscretization *other);
435 MEDCouplingTimeDiscretization *divide(const MEDCouplingTimeDiscretization *other) const;
436 void divideEqual(const MEDCouplingTimeDiscretization *other);
438 static const TypeOfTimeDiscretization DISCRETIZATION=LINEAR_TIME;
439 static const char REPR[];