Salome HOME
Boost of expression evaluator DataArrayDouble::applyFunc* + DataArrayDouble::applyFun...
[modules/med.git] / src / MEDCoupling / MEDCouplingMemArray.hxx
1 // Copyright (C) 2007-2014  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 __PARAMEDMEM_MEDCOUPLINGMEMARRAY_HXX__
22 #define __PARAMEDMEM_MEDCOUPLINGMEMARRAY_HXX__
23
24 #include "MEDCoupling.hxx"
25 #include "MEDCouplingTimeLabel.hxx"
26 #include "MEDCouplingRefCountObject.hxx"
27 #include "InterpKernelException.hxx"
28 #include "BBTreePts.txx"
29
30 #include <string>
31 #include <vector>
32 #include <iterator>
33
34 namespace ParaMEDMEM
35 {
36   template<class T>
37   class MEDCouplingPointer
38   {
39   public:
40     MEDCouplingPointer():_internal(0),_external(0) { }
41     void null() { _internal=0; _external=0; }
42     bool isNull() const { return _internal==0 && _external==0; }
43     void setInternal(T *pointer);
44     void setExternal(const T *pointer);
45     const T *getConstPointer() const { if(_internal) return _internal; else return _external; }
46     const T *getConstPointerLoc(std::size_t offset) const { if(_internal) return _internal+offset; else return _external+offset; }
47     T *getPointer() { if(_internal) return _internal; if(_external) throw INTERP_KERNEL::Exception("Trying to write on an external pointer."); else return 0; }
48   private:
49     T *_internal;
50     const T *_external;
51   };
52
53   template<class T>
54   class MemArray
55   {
56   public:
57     typedef void (*Deallocator)(void *,void *);
58   public:
59     MemArray():_nb_of_elem(0),_nb_of_elem_alloc(0),_ownership(false),_dealloc(0),_param_for_deallocator(0) { }
60     MemArray(const MemArray<T>& other);
61     bool isNull() const { return _pointer.isNull(); }
62     const T *getConstPointerLoc(std::size_t offset) const { return _pointer.getConstPointerLoc(offset); }
63     const T *getConstPointer() const { return _pointer.getConstPointer(); }
64     std::size_t getNbOfElem() const { return _nb_of_elem; }
65     std::size_t getNbOfElemAllocated() const { return _nb_of_elem_alloc; }
66     T *getPointer() { return _pointer.getPointer(); }
67     MemArray<T> &operator=(const MemArray<T>& other);
68     T operator[](std::size_t id) const { return _pointer.getConstPointer()[id]; }
69     T& operator[](std::size_t id) { return _pointer.getPointer()[id]; }
70     bool isEqual(const MemArray<T>& other, T prec, std::string& reason) const;
71     void repr(int sl, std::ostream& stream) const;
72     bool reprHeader(int sl, std::ostream& stream) const;
73     void reprZip(int sl, std::ostream& stream) const;
74     void fillWithValue(const T& val);
75     T *fromNoInterlace(int nbOfComp) const;
76     T *toNoInterlace(int nbOfComp) const;
77     void sort(bool asc);
78     void reverse(int nbOfComp);
79     void alloc(std::size_t nbOfElements);
80     void reserve(std::size_t newNbOfElements);
81     void reAlloc(std::size_t newNbOfElements);
82     void useArray(const T *array, bool ownership, DeallocType type, std::size_t nbOfElem);
83     void useExternalArrayWithRWAccess(const T *array, std::size_t nbOfElem);
84     void writeOnPlace(std::size_t id, T element0, const T *others, std::size_t sizeOfOthers);
85     template<class InputIterator>
86     void insertAtTheEnd(InputIterator first, InputIterator last);
87     void pushBack(T elem);
88     T popBack();
89     void pack() const;
90     bool isDeallocatorCalled() const { return _ownership; }
91     Deallocator getDeallocator() const { return _dealloc; }
92     void setSpecificDeallocator(Deallocator dealloc) { _dealloc=dealloc; }
93     void setParameterForDeallocator(void *param) { _param_for_deallocator=param; }
94     void *getParameterForDeallocator() const { return _param_for_deallocator; }
95     void destroy();
96     ~MemArray() { destroy(); }
97   public:
98     static void CPPDeallocator(void *pt, void *param);
99     static void CDeallocator(void *pt, void *param);
100   private:
101     static void DestroyPointer(T *pt, Deallocator dealloc, void *param);
102     static Deallocator BuildFromType(DeallocType type);
103   private:
104     std::size_t _nb_of_elem;
105     std::size_t _nb_of_elem_alloc;
106     bool _ownership;
107     MEDCouplingPointer<T> _pointer;
108     Deallocator _dealloc;
109     void *_param_for_deallocator;
110   };
111
112   class DataArrayInt;
113   class DataArrayByte;
114
115   class DataArray : public RefCountObject, public TimeLabel
116   {
117   public:
118     MEDCOUPLING_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
119     MEDCOUPLING_EXPORT std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
120     MEDCOUPLING_EXPORT void setName(const std::string& name);
121     MEDCOUPLING_EXPORT void copyStringInfoFrom(const DataArray& other);
122     MEDCOUPLING_EXPORT void copyPartOfStringInfoFrom(const DataArray& other, const std::vector<int>& compoIds);
123     MEDCOUPLING_EXPORT void copyPartOfStringInfoFrom2(const std::vector<int>& compoIds, const DataArray& other);
124     MEDCOUPLING_EXPORT bool areInfoEqualsIfNotWhy(const DataArray& other, std::string& reason) const;
125     MEDCOUPLING_EXPORT bool areInfoEquals(const DataArray& other) const;
126     MEDCOUPLING_EXPORT std::string cppRepr(const std::string& varName) const;
127     MEDCOUPLING_EXPORT std::string getName() const { return _name; }
128     MEDCOUPLING_EXPORT const std::vector<std::string> &getInfoOnComponents() const { return _info_on_compo; }
129     MEDCOUPLING_EXPORT std::vector<std::string> &getInfoOnComponents() { return _info_on_compo; }
130     MEDCOUPLING_EXPORT void setInfoOnComponents(const std::vector<std::string>& info);
131     MEDCOUPLING_EXPORT void setInfoAndChangeNbOfCompo(const std::vector<std::string>& info);
132     MEDCOUPLING_EXPORT std::vector<std::string> getVarsOnComponent() const;
133     MEDCOUPLING_EXPORT std::vector<std::string> getUnitsOnComponent() const;
134     MEDCOUPLING_EXPORT std::string getInfoOnComponent(int i) const;
135     MEDCOUPLING_EXPORT std::string getVarOnComponent(int i) const;
136     MEDCOUPLING_EXPORT std::string getUnitOnComponent(int i) const;
137     MEDCOUPLING_EXPORT void setInfoOnComponent(int i, const std::string& info);
138     MEDCOUPLING_EXPORT int getNumberOfComponents() const { return (int)_info_on_compo.size(); }
139     MEDCOUPLING_EXPORT void setPartOfValuesBase3(const DataArray *aBase, const int *bgTuples, const int *endTuples, int bgComp, int endComp, int stepComp, bool strictCompoCompare=true);
140     MEDCOUPLING_EXPORT virtual DataArray *deepCpy() const = 0;
141     MEDCOUPLING_EXPORT virtual bool isAllocated() const = 0;
142     MEDCOUPLING_EXPORT virtual void checkAllocated() const = 0;
143     MEDCOUPLING_EXPORT virtual void desallocate() = 0;
144     MEDCOUPLING_EXPORT virtual int getNumberOfTuples() const = 0;
145     MEDCOUPLING_EXPORT virtual std::size_t getNbOfElems() const = 0;
146     MEDCOUPLING_EXPORT virtual std::size_t getNbOfElemAllocated() const = 0;
147     MEDCOUPLING_EXPORT virtual void alloc(int nbOfTuple, int nbOfCompo=1) = 0;
148     MEDCOUPLING_EXPORT virtual void reAlloc(int newNbOfTuple) = 0;
149     MEDCOUPLING_EXPORT virtual void renumberInPlace(const int *old2New) = 0;
150     MEDCOUPLING_EXPORT virtual void renumberInPlaceR(const int *new2Old) = 0;
151     MEDCOUPLING_EXPORT virtual void setContigPartOfSelectedValues(int tupleIdStart, const DataArray *aBase, const DataArrayInt *tuplesSelec) = 0;
152     MEDCOUPLING_EXPORT virtual void setContigPartOfSelectedValues2(int tupleIdStart, const DataArray *aBase, int bg, int end2, int step) = 0;
153     MEDCOUPLING_EXPORT virtual DataArray *selectByTupleRanges(const std::vector<std::pair<int,int> >& ranges) const = 0;
154     MEDCOUPLING_EXPORT virtual DataArray *keepSelectedComponents(const std::vector<int>& compoIds) const = 0;
155     MEDCOUPLING_EXPORT virtual DataArray *selectByTupleId(const int *new2OldBg, const int *new2OldEnd) const = 0;
156     MEDCOUPLING_EXPORT virtual DataArray *selectByTupleIdSafe(const int *new2OldBg, const int *new2OldEnd) const = 0;
157     MEDCOUPLING_EXPORT virtual DataArray *selectByTupleId2(int bg, int end2, int step) const = 0;
158     MEDCOUPLING_EXPORT virtual void rearrange(int newNbOfCompo) = 0;
159     MEDCOUPLING_EXPORT void checkNbOfTuples(int nbOfTuples, const std::string& msg) const;
160     MEDCOUPLING_EXPORT void checkNbOfComps(int nbOfCompo, const std::string& msg) const;
161     MEDCOUPLING_EXPORT void checkNbOfTuplesAndComp(const DataArray& other, const std::string& msg) const;
162     MEDCOUPLING_EXPORT void checkNbOfTuplesAndComp(int nbOfTuples, int nbOfCompo, const std::string& msg) const;
163     MEDCOUPLING_EXPORT void checkNbOfElems(std::size_t nbOfElems, const std::string& msg) const;
164     MEDCOUPLING_EXPORT static void GetSlice(int start, int stop, int step, int sliceId, int nbOfSlices, int& startSlice, int& stopSlice);
165     MEDCOUPLING_EXPORT static int GetNumberOfItemGivenBES(int begin, int end, int step, const std::string& msg);
166     MEDCOUPLING_EXPORT static int GetNumberOfItemGivenBESRelative(int begin, int end, int step, const std::string& msg);
167     MEDCOUPLING_EXPORT static int GetPosOfItemGivenBESRelativeNoThrow(int value, int begin, int end, int step);
168     MEDCOUPLING_EXPORT static std::string GetVarNameFromInfo(const std::string& info);
169     MEDCOUPLING_EXPORT static std::string GetUnitFromInfo(const std::string& info);
170     MEDCOUPLING_EXPORT static std::string BuildInfoFromVarAndUnit(const std::string& var, const std::string& unit);
171     MEDCOUPLING_EXPORT static DataArray *Aggregate(const std::vector<const DataArray *>& arrs);
172     MEDCOUPLING_EXPORT virtual void reprStream(std::ostream& stream) const = 0;
173     MEDCOUPLING_EXPORT virtual void reprZipStream(std::ostream& stream) const = 0;
174     MEDCOUPLING_EXPORT virtual void reprWithoutNameStream(std::ostream& stream) const;
175     MEDCOUPLING_EXPORT virtual void reprZipWithoutNameStream(std::ostream& stream) const = 0;
176     MEDCOUPLING_EXPORT virtual void reprCppStream(const std::string& varName, std::ostream& stream) const = 0;
177     MEDCOUPLING_EXPORT virtual void reprQuickOverview(std::ostream& stream) const = 0;
178     MEDCOUPLING_EXPORT virtual void reprQuickOverviewData(std::ostream& stream, std::size_t maxNbOfByteInRepr) const = 0;
179   protected:
180     DataArray() { }
181     ~DataArray() { }
182   protected:
183     static void CheckValueInRange(int ref, int value, const std::string& msg);
184     static void CheckValueInRangeEx(int value, int start, int end, const std::string& msg);
185     static void CheckClosingParInRange(int ref, int value, const std::string& msg);
186   protected:
187     std::string _name;
188     std::vector<std::string> _info_on_compo;
189   };
190 }
191
192 #include "MEDCouplingMemArray.txx"
193
194 namespace ParaMEDMEM
195 {
196   class DataArrayInt;
197   class DataArrayDoubleIterator;
198   class DataArrayDouble : public DataArray
199   {
200   public:
201     MEDCOUPLING_EXPORT static DataArrayDouble *New();
202     MEDCOUPLING_EXPORT bool isAllocated() const;
203     MEDCOUPLING_EXPORT void checkAllocated() const;
204     MEDCOUPLING_EXPORT void desallocate();
205     MEDCOUPLING_EXPORT int getNumberOfTuples() const { return _info_on_compo.empty()?0:_mem.getNbOfElem()/getNumberOfComponents(); }
206     MEDCOUPLING_EXPORT std::size_t getNbOfElems() const { return _mem.getNbOfElem(); }
207     MEDCOUPLING_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
208     MEDCOUPLING_EXPORT double doubleValue() const;
209     MEDCOUPLING_EXPORT bool empty() const;
210     MEDCOUPLING_EXPORT DataArrayDouble *deepCpy() const;
211     MEDCOUPLING_EXPORT DataArrayDouble *performCpy(bool deepCpy) const;
212     MEDCOUPLING_EXPORT void cpyFrom(const DataArrayDouble& other);
213     MEDCOUPLING_EXPORT void reserve(std::size_t nbOfElems);
214     MEDCOUPLING_EXPORT void pushBackSilent(double val);
215     MEDCOUPLING_EXPORT void pushBackValsSilent(const double *valsBg, const double *valsEnd);
216     MEDCOUPLING_EXPORT double popBackSilent();
217     MEDCOUPLING_EXPORT void pack() const;
218     MEDCOUPLING_EXPORT std::size_t getNbOfElemAllocated() const { return _mem.getNbOfElemAllocated(); }
219     MEDCOUPLING_EXPORT void alloc(int nbOfTuple, int nbOfCompo=1);
220     MEDCOUPLING_EXPORT void allocIfNecessary(int nbOfTuple, int nbOfCompo);
221     MEDCOUPLING_EXPORT void fillWithZero();
222     MEDCOUPLING_EXPORT void fillWithValue(double val);
223     MEDCOUPLING_EXPORT void iota(double init=0.);
224     MEDCOUPLING_EXPORT bool isUniform(double val, double eps) const;
225     MEDCOUPLING_EXPORT void sort(bool asc=true);
226     MEDCOUPLING_EXPORT void reverse();
227     MEDCOUPLING_EXPORT void checkMonotonic(bool increasing, double eps) const;
228     MEDCOUPLING_EXPORT bool isMonotonic(bool increasing, double eps) const;
229     MEDCOUPLING_EXPORT std::string repr() const;
230     MEDCOUPLING_EXPORT std::string reprZip() const;
231     MEDCOUPLING_EXPORT void writeVTK(std::ostream& ofs, int indent, const std::string& nameInFile, DataArrayByte *byteArr) const;
232     MEDCOUPLING_EXPORT void reprStream(std::ostream& stream) const;
233     MEDCOUPLING_EXPORT void reprZipStream(std::ostream& stream) const;
234     MEDCOUPLING_EXPORT void reprWithoutNameStream(std::ostream& stream) const;
235     MEDCOUPLING_EXPORT void reprZipWithoutNameStream(std::ostream& stream) const;
236     MEDCOUPLING_EXPORT void reprCppStream(const std::string& varName, std::ostream& stream) const;
237     MEDCOUPLING_EXPORT void reprQuickOverview(std::ostream& stream) const;
238     MEDCOUPLING_EXPORT void reprQuickOverviewData(std::ostream& stream, std::size_t maxNbOfByteInRepr) const;
239     MEDCOUPLING_EXPORT bool isEqual(const DataArrayDouble& other, double prec) const;
240     MEDCOUPLING_EXPORT bool isEqualIfNotWhy(const DataArrayDouble& other, double prec, std::string& reason) const;
241     MEDCOUPLING_EXPORT bool isEqualWithoutConsideringStr(const DataArrayDouble& other, double prec) const;
242     MEDCOUPLING_EXPORT void reAlloc(int nbOfTuples);
243     MEDCOUPLING_EXPORT DataArrayInt *convertToIntArr() const;
244     MEDCOUPLING_EXPORT DataArrayDouble *fromNoInterlace() const;
245     MEDCOUPLING_EXPORT DataArrayDouble *toNoInterlace() const;
246     MEDCOUPLING_EXPORT void renumberInPlace(const int *old2New);
247     MEDCOUPLING_EXPORT void renumberInPlaceR(const int *new2Old);
248     MEDCOUPLING_EXPORT DataArrayDouble *renumber(const int *old2New) const;
249     MEDCOUPLING_EXPORT DataArrayDouble *renumberR(const int *new2Old) const;
250     MEDCOUPLING_EXPORT DataArrayDouble *renumberAndReduce(const int *old2New, int newNbOfTuple) const;
251     MEDCOUPLING_EXPORT DataArrayDouble *selectByTupleId(const int *new2OldBg, const int *new2OldEnd) const;
252     MEDCOUPLING_EXPORT DataArrayDouble *selectByTupleIdSafe(const int *new2OldBg, const int *new2OldEnd) const;
253     MEDCOUPLING_EXPORT DataArrayDouble *selectByTupleId2(int bg, int end2, int step) const;
254     MEDCOUPLING_EXPORT DataArray *selectByTupleRanges(const std::vector<std::pair<int,int> >& ranges) const;
255     MEDCOUPLING_EXPORT DataArrayDouble *substr(int tupleIdBg, int tupleIdEnd=-1) const;
256     MEDCOUPLING_EXPORT void rearrange(int newNbOfCompo);
257     MEDCOUPLING_EXPORT void transpose();
258     MEDCOUPLING_EXPORT DataArrayDouble *changeNbOfComponents(int newNbOfComp, double dftValue) const;
259     MEDCOUPLING_EXPORT DataArrayDouble *keepSelectedComponents(const std::vector<int>& compoIds) const;
260     MEDCOUPLING_EXPORT void meldWith(const DataArrayDouble *other);
261     MEDCOUPLING_EXPORT bool areIncludedInMe(const DataArrayDouble *other, double prec, DataArrayInt *&tupleIds) const;
262     MEDCOUPLING_EXPORT void findCommonTuples(double prec, int limitTupleId, DataArrayInt *&comm, DataArrayInt *&commIndex) const;
263     MEDCOUPLING_EXPORT double minimalDistanceTo(const DataArrayDouble *other, int& thisTupleId, int& otherTupleId) const;
264     MEDCOUPLING_EXPORT DataArrayDouble *duplicateEachTupleNTimes(int nbTimes) const;
265     MEDCOUPLING_EXPORT DataArrayDouble *getDifferentValues(double prec, int limitTupleId=-1) const;
266     MEDCOUPLING_EXPORT DataArrayInt *findClosestTupleId(const DataArrayDouble *other) const;
267     MEDCOUPLING_EXPORT DataArrayInt *computeNbOfInteractionsWith(const DataArrayDouble *otherBBoxFrmt, double eps) const;
268     MEDCOUPLING_EXPORT void setSelectedComponents(const DataArrayDouble *a, const std::vector<int>& compoIds);
269     MEDCOUPLING_EXPORT void setPartOfValues1(const DataArrayDouble *a, int bgTuples, int endTuples, int stepTuples, int bgComp, int endComp, int stepComp, bool strictCompoCompare=true);
270     MEDCOUPLING_EXPORT void setPartOfValuesSimple1(double a, int bgTuples, int endTuples, int stepTuples, int bgComp, int endComp, int stepComp);
271     MEDCOUPLING_EXPORT void setPartOfValues2(const DataArrayDouble *a, const int *bgTuples, const int *endTuples, const int *bgComp, const int *endComp, bool strictCompoCompare=true);
272     MEDCOUPLING_EXPORT void setPartOfValuesSimple2(double a, const int *bgTuples, const int *endTuples, const int *bgComp, const int *endComp);
273     MEDCOUPLING_EXPORT void setPartOfValues3(const DataArrayDouble *a, const int *bgTuples, const int *endTuples, int bgComp, int endComp, int stepComp, bool strictCompoCompare=true);
274     MEDCOUPLING_EXPORT void setPartOfValuesSimple3(double a, const int *bgTuples, const int *endTuples, int bgComp, int endComp, int stepComp);
275     MEDCOUPLING_EXPORT void setPartOfValues4(const DataArrayDouble *a, int bgTuples, int endTuples, int stepTuples, const int *bgComp, const int *endComp, bool strictCompoCompare=true);
276     MEDCOUPLING_EXPORT void setPartOfValuesSimple4(double a, int bgTuples, int endTuples, int stepTuples, const int *bgComp, const int *endComp);
277     MEDCOUPLING_EXPORT void setPartOfValuesAdv(const DataArrayDouble *a, const DataArrayInt *tuplesSelec);
278     MEDCOUPLING_EXPORT void setContigPartOfSelectedValues(int tupleIdStart, const DataArray *aBase, const DataArrayInt *tuplesSelec);
279     MEDCOUPLING_EXPORT void setContigPartOfSelectedValues2(int tupleIdStart, const DataArray *aBase, int bg, int end2, int step);
280     MEDCOUPLING_EXPORT void getTuple(int tupleId, double *res) const { std::copy(_mem.getConstPointerLoc(tupleId*_info_on_compo.size()),_mem.getConstPointerLoc((tupleId+1)*_info_on_compo.size()),res); }
281     MEDCOUPLING_EXPORT double getIJ(int tupleId, int compoId) const { return _mem[tupleId*_info_on_compo.size()+compoId]; }
282     MEDCOUPLING_EXPORT double front() const;
283     MEDCOUPLING_EXPORT double back() const;
284     MEDCOUPLING_EXPORT double getIJSafe(int tupleId, int compoId) const;
285     MEDCOUPLING_EXPORT void setIJ(int tupleId, int compoId, double newVal) { _mem[tupleId*_info_on_compo.size()+compoId]=newVal; declareAsNew(); }
286     MEDCOUPLING_EXPORT void setIJSilent(int tupleId, int compoId, double newVal) { _mem[tupleId*_info_on_compo.size()+compoId]=newVal; }
287     MEDCOUPLING_EXPORT double *getPointer() { return _mem.getPointer(); }
288     MEDCOUPLING_EXPORT static void SetArrayIn(DataArrayDouble *newArray, DataArrayDouble* &arrayToSet);
289     MEDCOUPLING_EXPORT const double *getConstPointer() const { return _mem.getConstPointer(); }
290     MEDCOUPLING_EXPORT DataArrayDoubleIterator *iterator();
291     MEDCOUPLING_EXPORT const double *begin() const { return getConstPointer(); }
292     MEDCOUPLING_EXPORT const double *end() const { return getConstPointer()+getNbOfElems(); }
293     MEDCOUPLING_EXPORT void useArray(const double *array, bool ownership, DeallocType type, int nbOfTuple, int nbOfCompo);
294     MEDCOUPLING_EXPORT void useExternalArrayWithRWAccess(const double *array, int nbOfTuple, int nbOfCompo);
295     template<class InputIterator>
296     void insertAtTheEnd(InputIterator first, InputIterator last);
297     MEDCOUPLING_EXPORT void writeOnPlace(std::size_t id, double element0, const double *others, int sizeOfOthers) { _mem.writeOnPlace(id,element0,others,sizeOfOthers); }
298     MEDCOUPLING_EXPORT void checkNoNullValues() const;
299     MEDCOUPLING_EXPORT void getMinMaxPerComponent(double *bounds) const;
300     MEDCOUPLING_EXPORT DataArrayDouble *computeBBoxPerTuple(double epsilon=0.0) const;
301     MEDCOUPLING_EXPORT void computeTupleIdsNearTuples(const DataArrayDouble *other, double eps, DataArrayInt *& c, DataArrayInt *& cI) const;
302     MEDCOUPLING_EXPORT void recenterForMaxPrecision(double eps);
303     MEDCOUPLING_EXPORT double getMaxValue(int& tupleId) const;
304     MEDCOUPLING_EXPORT double getMaxValueInArray() const;
305     MEDCOUPLING_EXPORT double getMinValue(int& tupleId) const;
306     MEDCOUPLING_EXPORT double getMinValueInArray() const;
307     MEDCOUPLING_EXPORT double getMaxValue2(DataArrayInt*& tupleIds) const;
308     MEDCOUPLING_EXPORT double getMinValue2(DataArrayInt*& tupleIds) const;
309     MEDCOUPLING_EXPORT int count(double value, double eps) const;
310     MEDCOUPLING_EXPORT double getAverageValue() const;
311     MEDCOUPLING_EXPORT double norm2() const;
312     MEDCOUPLING_EXPORT double normMax() const;
313     MEDCOUPLING_EXPORT double normMin() const;
314     MEDCOUPLING_EXPORT void accumulate(double *res) const;
315     MEDCOUPLING_EXPORT double accumulate(int compId) const;
316     MEDCOUPLING_EXPORT DataArrayDouble *accumulatePerChunck(const int *bgOfIndex, const int *endOfIndex) const;
317     MEDCOUPLING_EXPORT double distanceToTuple(const double *tupleBg, const double *tupleEnd, int& tupleId) const;
318     MEDCOUPLING_EXPORT DataArrayDouble *fromPolarToCart() const;
319     MEDCOUPLING_EXPORT DataArrayDouble *fromCylToCart() const;
320     MEDCOUPLING_EXPORT DataArrayDouble *fromSpherToCart() const;
321     MEDCOUPLING_EXPORT DataArrayDouble *doublyContractedProduct() const;
322     MEDCOUPLING_EXPORT DataArrayDouble *determinant() const;
323     MEDCOUPLING_EXPORT DataArrayDouble *eigenValues() const;
324     MEDCOUPLING_EXPORT DataArrayDouble *eigenVectors() const;
325     MEDCOUPLING_EXPORT DataArrayDouble *inverse() const;
326     MEDCOUPLING_EXPORT DataArrayDouble *trace() const;
327     MEDCOUPLING_EXPORT DataArrayDouble *deviator() const;
328     MEDCOUPLING_EXPORT DataArrayDouble *magnitude() const;
329     MEDCOUPLING_EXPORT DataArrayDouble *sumPerTuple() const;
330     MEDCOUPLING_EXPORT DataArrayDouble *maxPerTuple() const;
331     MEDCOUPLING_EXPORT DataArrayDouble *maxPerTupleWithCompoId(DataArrayInt* &compoIdOfMaxPerTuple) const;
332     MEDCOUPLING_EXPORT DataArrayDouble *buildEuclidianDistanceDenseMatrix() const;
333     MEDCOUPLING_EXPORT DataArrayDouble *buildEuclidianDistanceDenseMatrixWith(const DataArrayDouble *other) const;
334     MEDCOUPLING_EXPORT void sortPerTuple(bool asc);
335     MEDCOUPLING_EXPORT void abs();
336     MEDCOUPLING_EXPORT DataArrayDouble *computeAbs() const;
337     MEDCOUPLING_EXPORT void applyLin(double a, double b, int compoId);
338     MEDCOUPLING_EXPORT void applyLin(double a, double b);
339     MEDCOUPLING_EXPORT void applyInv(double numerator);
340     MEDCOUPLING_EXPORT void applyPow(double val);
341     MEDCOUPLING_EXPORT void applyRPow(double val);
342     MEDCOUPLING_EXPORT DataArrayDouble *negate() const;
343     MEDCOUPLING_EXPORT DataArrayDouble *applyFunc(int nbOfComp, FunctionToEvaluate func) const;
344     MEDCOUPLING_EXPORT DataArrayDouble *applyFunc(int nbOfComp, const std::string& func, bool isSafe=true) const;
345     MEDCOUPLING_EXPORT DataArrayDouble *applyFunc(const std::string& func, bool isSafe=true) const;
346     MEDCOUPLING_EXPORT void applyFuncOnThis(const std::string& func, bool isSafe=true);
347     MEDCOUPLING_EXPORT DataArrayDouble *applyFunc2(int nbOfComp, const std::string& func, bool isSafe=true) const;
348     MEDCOUPLING_EXPORT DataArrayDouble *applyFunc3(int nbOfComp, const std::vector<std::string>& varsOrder, const std::string& func, bool isSafe=true) const;
349     MEDCOUPLING_EXPORT void applyFuncFast32(const std::string& func);
350     MEDCOUPLING_EXPORT void applyFuncFast64(const std::string& func);
351     MEDCOUPLING_EXPORT DataArrayInt *getIdsInRange(double vmin, double vmax) const;
352     MEDCOUPLING_EXPORT DataArrayInt *getIdsNotInRange(double vmin, double vmax) const;
353     MEDCOUPLING_EXPORT static DataArrayDouble *Aggregate(const DataArrayDouble *a1, const DataArrayDouble *a2);
354     MEDCOUPLING_EXPORT static DataArrayDouble *Aggregate(const std::vector<const DataArrayDouble *>& arr);
355     MEDCOUPLING_EXPORT static DataArrayDouble *Meld(const DataArrayDouble *a1, const DataArrayDouble *a2);
356     MEDCOUPLING_EXPORT static DataArrayDouble *Meld(const std::vector<const DataArrayDouble *>& arr);
357     MEDCOUPLING_EXPORT static DataArrayDouble *Dot(const DataArrayDouble *a1, const DataArrayDouble *a2);
358     MEDCOUPLING_EXPORT static DataArrayDouble *CrossProduct(const DataArrayDouble *a1, const DataArrayDouble *a2);
359     MEDCOUPLING_EXPORT static DataArrayDouble *Max(const DataArrayDouble *a1, const DataArrayDouble *a2);
360     MEDCOUPLING_EXPORT static DataArrayDouble *Min(const DataArrayDouble *a1, const DataArrayDouble *a2);
361     MEDCOUPLING_EXPORT static DataArrayDouble *Add(const DataArrayDouble *a1, const DataArrayDouble *a2);
362     MEDCOUPLING_EXPORT void addEqual(const DataArrayDouble *other);
363     MEDCOUPLING_EXPORT static DataArrayDouble *Substract(const DataArrayDouble *a1, const DataArrayDouble *a2);
364     MEDCOUPLING_EXPORT void substractEqual(const DataArrayDouble *other);
365     MEDCOUPLING_EXPORT static DataArrayDouble *Multiply(const DataArrayDouble *a1, const DataArrayDouble *a2);
366     MEDCOUPLING_EXPORT void multiplyEqual(const DataArrayDouble *other);
367     MEDCOUPLING_EXPORT static DataArrayDouble *Divide(const DataArrayDouble *a1, const DataArrayDouble *a2);
368     MEDCOUPLING_EXPORT void divideEqual(const DataArrayDouble *other);
369     MEDCOUPLING_EXPORT static DataArrayDouble *Pow(const DataArrayDouble *a1, const DataArrayDouble *a2);
370     MEDCOUPLING_EXPORT void powEqual(const DataArrayDouble *other);
371     MEDCOUPLING_EXPORT void updateTime() const { }
372     MEDCOUPLING_EXPORT MemArray<double>& accessToMemArray() { return _mem; }
373     MEDCOUPLING_EXPORT const MemArray<double>& accessToMemArray() const { return _mem; }
374     MEDCOUPLING_EXPORT std::vector<bool> toVectorOfBool(double eps) const;
375   public:
376     MEDCOUPLING_EXPORT void getTinySerializationIntInformation(std::vector<int>& tinyInfo) const;
377     MEDCOUPLING_EXPORT void getTinySerializationStrInformation(std::vector<std::string>& tinyInfo) const;
378     MEDCOUPLING_EXPORT bool resizeForUnserialization(const std::vector<int>& tinyInfoI);
379     MEDCOUPLING_EXPORT void finishUnserialization(const std::vector<int>& tinyInfoI, const std::vector<std::string>& tinyInfoS);
380   public:
381     template<int SPACEDIM>
382     void findCommonTuplesAlg(const double *bbox, int nbNodes, int limitNodeId, double prec, DataArrayInt *c, DataArrayInt *cI) const;
383     template<int SPACEDIM>
384     static void FindClosestTupleIdAlg(const BBTreePts<SPACEDIM,int>& myTree, double dist, const double *pos, int nbOfTuples, const double *thisPt, int thisNbOfTuples, int *res);
385     template<int SPACEDIM>
386     static void FindTupleIdsNearTuplesAlg(const BBTreePts<SPACEDIM,int>& myTree, const double *pos, int nbOfTuples, double eps,
387                                           DataArrayInt *c, DataArrayInt *cI);
388   private:
389     ~DataArrayDouble() { }
390     DataArrayDouble() { }
391   private:
392     MemArray<double> _mem;
393   };
394
395   class DataArrayDoubleTuple;
396
397   class DataArrayDoubleIterator
398   {
399   public:
400     MEDCOUPLING_EXPORT DataArrayDoubleIterator(DataArrayDouble *da);
401     MEDCOUPLING_EXPORT ~DataArrayDoubleIterator();
402     MEDCOUPLING_EXPORT DataArrayDoubleTuple *nextt();
403   private:
404     DataArrayDouble *_da;
405     double *_pt;
406     int _tuple_id;
407     int _nb_comp;
408     int _nb_tuple;
409   };
410
411   class DataArrayDoubleTuple
412   {
413   public:
414     MEDCOUPLING_EXPORT DataArrayDoubleTuple(double *pt, int nbOfComp);
415     MEDCOUPLING_EXPORT std::string repr() const;
416     MEDCOUPLING_EXPORT int getNumberOfCompo() const { return _nb_of_compo; }
417     MEDCOUPLING_EXPORT const double *getConstPointer() const { return  _pt; }
418     MEDCOUPLING_EXPORT double *getPointer() { return _pt; }
419     MEDCOUPLING_EXPORT double doubleValue() const;
420     MEDCOUPLING_EXPORT DataArrayDouble *buildDADouble(int nbOfTuples, int nbOfCompo) const;
421   private:
422     double *_pt;
423     int _nb_of_compo;
424   };
425
426   class DataArrayIntIterator;
427
428   class DataArrayInt : public DataArray
429   {
430   public:
431     MEDCOUPLING_EXPORT static DataArrayInt *New();
432     MEDCOUPLING_EXPORT bool isAllocated() const;
433     MEDCOUPLING_EXPORT void checkAllocated() const;
434     MEDCOUPLING_EXPORT void desallocate();
435     MEDCOUPLING_EXPORT int getNumberOfTuples() const { return _info_on_compo.empty()?0:_mem.getNbOfElem()/getNumberOfComponents(); }
436     MEDCOUPLING_EXPORT std::size_t getNbOfElems() const { return _mem.getNbOfElem(); }
437     MEDCOUPLING_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
438     MEDCOUPLING_EXPORT int intValue() const;
439     MEDCOUPLING_EXPORT int getHashCode() const;
440     MEDCOUPLING_EXPORT bool empty() const;
441     MEDCOUPLING_EXPORT DataArrayInt *deepCpy() const;
442     MEDCOUPLING_EXPORT DataArrayInt *performCpy(bool deepCpy) const;
443     MEDCOUPLING_EXPORT void cpyFrom(const DataArrayInt& other);
444     MEDCOUPLING_EXPORT void reserve(std::size_t nbOfElems);
445     MEDCOUPLING_EXPORT void pushBackSilent(int val);
446     MEDCOUPLING_EXPORT void pushBackValsSilent(const int *valsBg, const int *valsEnd);
447     MEDCOUPLING_EXPORT int popBackSilent();
448     MEDCOUPLING_EXPORT void pack() const;
449     MEDCOUPLING_EXPORT std::size_t getNbOfElemAllocated() const { return _mem.getNbOfElemAllocated(); }
450     MEDCOUPLING_EXPORT void alloc(int nbOfTuple, int nbOfCompo=1);
451     MEDCOUPLING_EXPORT void allocIfNecessary(int nbOfTuple, int nbOfCompo);
452     MEDCOUPLING_EXPORT bool isEqual(const DataArrayInt& other) const;
453     MEDCOUPLING_EXPORT bool isEqualIfNotWhy(const DataArrayInt& other, std::string& reason) const;
454     MEDCOUPLING_EXPORT bool isEqualWithoutConsideringStr(const DataArrayInt& other) const;
455     MEDCOUPLING_EXPORT bool isEqualWithoutConsideringStrAndOrder(const DataArrayInt& other) const;
456     MEDCOUPLING_EXPORT bool isFittingWith(const std::vector<bool>& v) const;
457     MEDCOUPLING_EXPORT DataArrayInt *buildPermutationArr(const DataArrayInt& other) const;
458     MEDCOUPLING_EXPORT DataArrayInt *sumPerTuple() const;
459     MEDCOUPLING_EXPORT void sort(bool asc=true);
460     MEDCOUPLING_EXPORT void reverse();
461     MEDCOUPLING_EXPORT void checkMonotonic(bool increasing) const;
462     MEDCOUPLING_EXPORT bool isMonotonic(bool increasing) const;
463     MEDCOUPLING_EXPORT void checkStrictlyMonotonic(bool increasing) const;
464     MEDCOUPLING_EXPORT bool isStrictlyMonotonic(bool increasing) const;
465     MEDCOUPLING_EXPORT void fillWithZero();
466     MEDCOUPLING_EXPORT void fillWithValue(int val);
467     MEDCOUPLING_EXPORT void iota(int init=0);
468     MEDCOUPLING_EXPORT std::string repr() const;
469     MEDCOUPLING_EXPORT std::string reprZip() const;
470     MEDCOUPLING_EXPORT void writeVTK(std::ostream& ofs, int indent, const std::string& type, const std::string& nameInFile, DataArrayByte *byteArr) const;
471     MEDCOUPLING_EXPORT void reprStream(std::ostream& stream) const;
472     MEDCOUPLING_EXPORT void reprZipStream(std::ostream& stream) const;
473     MEDCOUPLING_EXPORT void reprWithoutNameStream(std::ostream& stream) const;
474     MEDCOUPLING_EXPORT void reprZipWithoutNameStream(std::ostream& stream) const;
475     MEDCOUPLING_EXPORT void reprCppStream(const std::string& varName, std::ostream& stream) const;
476     MEDCOUPLING_EXPORT void reprQuickOverview(std::ostream& stream) const;
477     MEDCOUPLING_EXPORT void reprQuickOverviewData(std::ostream& stream, std::size_t maxNbOfByteInRepr) const;
478     MEDCOUPLING_EXPORT void transformWithIndArr(const int *indArrBg, const int *indArrEnd);
479     MEDCOUPLING_EXPORT DataArrayInt *transformWithIndArrR(const int *indArrBg, const int *indArrEnd) const;
480     MEDCOUPLING_EXPORT void splitByValueRange(const int *arrBg, const int *arrEnd,
481                                               DataArrayInt *& castArr, DataArrayInt *& rankInsideCast, DataArrayInt *& castsPresent) const;
482     MEDCOUPLING_EXPORT DataArrayInt *invertArrayO2N2N2O(int newNbOfElem) const;
483     MEDCOUPLING_EXPORT DataArrayInt *invertArrayN2O2O2N(int oldNbOfElem) const;
484     MEDCOUPLING_EXPORT DataArrayInt *invertArrayO2N2N2OBis(int newNbOfElem) const;
485     MEDCOUPLING_EXPORT void reAlloc(int nbOfTuples);
486     MEDCOUPLING_EXPORT DataArrayDouble *convertToDblArr() const;
487     MEDCOUPLING_EXPORT DataArrayInt *fromNoInterlace() const;
488     MEDCOUPLING_EXPORT DataArrayInt *toNoInterlace() const;
489     MEDCOUPLING_EXPORT void renumberInPlace(const int *old2New);
490     MEDCOUPLING_EXPORT void renumberInPlaceR(const int *new2Old);
491     MEDCOUPLING_EXPORT DataArrayInt *renumber(const int *old2New) const;
492     MEDCOUPLING_EXPORT DataArrayInt *renumberR(const int *new2Old) const;
493     MEDCOUPLING_EXPORT DataArrayInt *renumberAndReduce(const int *old2NewBg, int newNbOfTuple) const;
494     MEDCOUPLING_EXPORT DataArrayInt *selectByTupleId(const int *new2OldBg, const int *new2OldEnd) const;
495     MEDCOUPLING_EXPORT DataArrayInt *selectByTupleIdSafe(const int *new2OldBg, const int *new2OldEnd) const;
496     MEDCOUPLING_EXPORT DataArrayInt *selectByTupleId2(int bg, int end, int step) const;
497     MEDCOUPLING_EXPORT DataArray *selectByTupleRanges(const std::vector<std::pair<int,int> >& ranges) const;
498     MEDCOUPLING_EXPORT DataArrayInt *checkAndPreparePermutation() const;
499     MEDCOUPLING_EXPORT static DataArrayInt *FindPermutationFromFirstToSecond(const DataArrayInt *ids1, const DataArrayInt *ids2);
500     MEDCOUPLING_EXPORT void changeSurjectiveFormat(int targetNb, DataArrayInt *&arr, DataArrayInt *&arrI) const;
501     MEDCOUPLING_EXPORT static DataArrayInt *BuildOld2NewArrayFromSurjectiveFormat2(int nbOfOldTuples, const int *arr, const int *arrIBg, const int *arrIEnd, int &newNbOfTuples);
502     MEDCOUPLING_EXPORT DataArrayInt *buildPermArrPerLevel() const;
503     MEDCOUPLING_EXPORT bool isIdentity() const;
504     MEDCOUPLING_EXPORT bool isUniform(int val) const;
505     MEDCOUPLING_EXPORT DataArrayInt *substr(int tupleIdBg, int tupleIdEnd=-1) const;
506     MEDCOUPLING_EXPORT void rearrange(int newNbOfCompo);
507     MEDCOUPLING_EXPORT void transpose();
508     MEDCOUPLING_EXPORT DataArrayInt *changeNbOfComponents(int newNbOfComp, int dftValue) const;
509     MEDCOUPLING_EXPORT DataArrayInt *keepSelectedComponents(const std::vector<int>& compoIds) const;
510     MEDCOUPLING_EXPORT void meldWith(const DataArrayInt *other);
511     MEDCOUPLING_EXPORT void setSelectedComponents(const DataArrayInt *a, const std::vector<int>& compoIds);
512     MEDCOUPLING_EXPORT void setPartOfValues1(const DataArrayInt *a, int bgTuples, int endTuples, int stepTuples, int bgComp, int endComp, int stepComp, bool strictCompoCompare=true);
513     MEDCOUPLING_EXPORT void setPartOfValuesSimple1(int a, int bgTuples, int endTuples, int stepTuples, int bgComp, int endComp, int stepComp);
514     MEDCOUPLING_EXPORT void setPartOfValues2(const DataArrayInt *a, const int *bgTuples, const int *endTuples, const int *bgComp, const int *endComp, bool strictCompoCompare=true);
515     MEDCOUPLING_EXPORT void setPartOfValuesSimple2(int a, const int *bgTuples, const int *endTuples, const int *bgComp, const int *endComp);
516     MEDCOUPLING_EXPORT void setPartOfValues3(const DataArrayInt *a, const int *bgTuples, const int *endTuples, int bgComp, int endComp, int stepComp, bool strictCompoCompare=true);
517     MEDCOUPLING_EXPORT void setPartOfValuesSimple3(int a, const int *bgTuples, const int *endTuples, int bgComp, int endComp, int stepComp);
518     MEDCOUPLING_EXPORT void setPartOfValues4(const DataArrayInt *a, int bgTuples, int endTuples, int stepTuples, const int *bgComp, const int *endComp, bool strictCompoCompare=true);
519     MEDCOUPLING_EXPORT void setPartOfValuesSimple4(int a, int bgTuples, int endTuples, int stepTuples, const int *bgComp, const int *endComp);
520     MEDCOUPLING_EXPORT void setPartOfValuesAdv(const DataArrayInt *a, const DataArrayInt *tuplesSelec);
521     MEDCOUPLING_EXPORT void setContigPartOfSelectedValues(int tupleIdStart, const DataArray *aBase, const DataArrayInt *tuplesSelec);
522     MEDCOUPLING_EXPORT void setContigPartOfSelectedValues2(int tupleIdStart, const DataArray *aBase, int bg, int end2, int step);
523     MEDCOUPLING_EXPORT void getTuple(int tupleId, int *res) const { std::copy(_mem.getConstPointerLoc(tupleId*_info_on_compo.size()),_mem.getConstPointerLoc((tupleId+1)*_info_on_compo.size()),res); }
524     MEDCOUPLING_EXPORT int getIJ(int tupleId, int compoId) const { return _mem[tupleId*_info_on_compo.size()+compoId]; }
525     MEDCOUPLING_EXPORT int getIJSafe(int tupleId, int compoId) const;
526     MEDCOUPLING_EXPORT int front() const;
527     MEDCOUPLING_EXPORT int back() const;
528     MEDCOUPLING_EXPORT void setIJ(int tupleId, int compoId, int newVal) { _mem[tupleId*_info_on_compo.size()+compoId]=newVal; declareAsNew(); }
529     MEDCOUPLING_EXPORT void setIJSilent(int tupleId, int compoId, int newVal) { _mem[tupleId*_info_on_compo.size()+compoId]=newVal; }
530     MEDCOUPLING_EXPORT int *getPointer() { return _mem.getPointer(); }
531     MEDCOUPLING_EXPORT static void SetArrayIn(DataArrayInt *newArray, DataArrayInt* &arrayToSet);
532     MEDCOUPLING_EXPORT const int *getConstPointer() const { return _mem.getConstPointer(); }
533     MEDCOUPLING_EXPORT DataArrayIntIterator *iterator();
534     MEDCOUPLING_EXPORT const int *begin() const { return getConstPointer(); }
535     MEDCOUPLING_EXPORT const int *end() const { return getConstPointer()+getNbOfElems(); }
536     MEDCOUPLING_EXPORT DataArrayInt *getIdsEqual(int val) const;
537     MEDCOUPLING_EXPORT DataArrayInt *getIdsNotEqual(int val) const;
538     MEDCOUPLING_EXPORT DataArrayInt *getIdsEqualList(const int *valsBg, const int *valsEnd) const;
539     MEDCOUPLING_EXPORT DataArrayInt *getIdsNotEqualList(const int *valsBg, const int *valsEnd) const;
540     MEDCOUPLING_EXPORT DataArrayInt *getIdsEqualTuple(const int *tupleBg, const int *tupleEnd) const;
541     MEDCOUPLING_EXPORT int changeValue(int oldValue, int newValue);
542     MEDCOUPLING_EXPORT int locateTuple(const std::vector<int>& tupl) const;
543     MEDCOUPLING_EXPORT int locateValue(int value) const;
544     MEDCOUPLING_EXPORT int locateValue(const std::vector<int>& vals) const;
545     MEDCOUPLING_EXPORT int search(const std::vector<int>& vals) const;
546     MEDCOUPLING_EXPORT bool presenceOfTuple(const std::vector<int>& tupl) const;
547     MEDCOUPLING_EXPORT bool presenceOfValue(int value) const;
548     MEDCOUPLING_EXPORT bool presenceOfValue(const std::vector<int>& vals) const;
549     MEDCOUPLING_EXPORT int count(int value) const;
550     MEDCOUPLING_EXPORT void accumulate(int *res) const;
551     MEDCOUPLING_EXPORT int accumulate(int compId) const;
552     MEDCOUPLING_EXPORT DataArrayInt *accumulatePerChunck(const int *bgOfIndex, const int *endOfIndex) const;
553     MEDCOUPLING_EXPORT int getMaxValue(int& tupleId) const;
554     MEDCOUPLING_EXPORT int getMaxValueInArray() const;
555     MEDCOUPLING_EXPORT int getMinValue(int& tupleId) const;
556     MEDCOUPLING_EXPORT int getMinValueInArray() const;
557     MEDCOUPLING_EXPORT void abs();
558     MEDCOUPLING_EXPORT DataArrayInt *computeAbs() const;
559     MEDCOUPLING_EXPORT void applyLin(int a, int b, int compoId);
560     MEDCOUPLING_EXPORT void applyLin(int a, int b);
561     MEDCOUPLING_EXPORT void applyInv(int numerator);
562     MEDCOUPLING_EXPORT DataArrayInt *negate() const;
563     MEDCOUPLING_EXPORT void applyDivideBy(int val);
564     MEDCOUPLING_EXPORT void applyModulus(int val);
565     MEDCOUPLING_EXPORT void applyRModulus(int val);
566     MEDCOUPLING_EXPORT void applyPow(int val);
567     MEDCOUPLING_EXPORT void applyRPow(int val);
568     MEDCOUPLING_EXPORT DataArrayInt *getIdsInRange(int vmin, int vmax) const;
569     MEDCOUPLING_EXPORT DataArrayInt *getIdsNotInRange(int vmin, int vmax) const;
570     MEDCOUPLING_EXPORT bool checkAllIdsInRange(int vmin, int vmax) const;
571     MEDCOUPLING_EXPORT static DataArrayInt *Aggregate(const DataArrayInt *a1, const DataArrayInt *a2, int offsetA2);
572     MEDCOUPLING_EXPORT static DataArrayInt *Aggregate(const std::vector<const DataArrayInt *>& arr);
573     MEDCOUPLING_EXPORT static DataArrayInt *AggregateIndexes(const std::vector<const DataArrayInt *>& arrs);
574     MEDCOUPLING_EXPORT static DataArrayInt *Meld(const DataArrayInt *a1, const DataArrayInt *a2);
575     MEDCOUPLING_EXPORT static DataArrayInt *Meld(const std::vector<const DataArrayInt *>& arr);
576     MEDCOUPLING_EXPORT static DataArrayInt *MakePartition(const std::vector<const DataArrayInt *>& groups, int newNb, std::vector< std::vector<int> >& fidsOfGroups);
577     MEDCOUPLING_EXPORT static DataArrayInt *BuildUnion(const std::vector<const DataArrayInt *>& arr);
578     MEDCOUPLING_EXPORT static DataArrayInt *BuildIntersection(const std::vector<const DataArrayInt *>& arr);
579     MEDCOUPLING_EXPORT static DataArrayInt *BuildListOfSwitchedOn(const std::vector<bool>& v);
580     MEDCOUPLING_EXPORT static DataArrayInt *BuildListOfSwitchedOff(const std::vector<bool>& v);
581     MEDCOUPLING_EXPORT static void PutIntoToSkylineFrmt(const std::vector< std::vector<int> >& v, DataArrayInt *& data, DataArrayInt *& dataIndex);
582     MEDCOUPLING_EXPORT DataArrayInt *buildComplement(int nbOfElement) const;
583     MEDCOUPLING_EXPORT DataArrayInt *buildSubstraction(const DataArrayInt *other) const;
584     MEDCOUPLING_EXPORT DataArrayInt *buildSubstractionOptimized(const DataArrayInt *other) const;
585     MEDCOUPLING_EXPORT DataArrayInt *buildUnion(const DataArrayInt *other) const;
586     MEDCOUPLING_EXPORT DataArrayInt *buildIntersection(const DataArrayInt *other) const;
587     MEDCOUPLING_EXPORT DataArrayInt *buildUnique() const;
588     MEDCOUPLING_EXPORT DataArrayInt *deltaShiftIndex() const;
589     MEDCOUPLING_EXPORT void computeOffsets();
590     MEDCOUPLING_EXPORT void computeOffsets2();
591     MEDCOUPLING_EXPORT void searchRangesInListOfIds(const DataArrayInt *listOfIds, DataArrayInt *& rangeIdsFetched, DataArrayInt *& idsInInputListThatFetch) const;
592     MEDCOUPLING_EXPORT DataArrayInt *buildExplicitArrByRanges(const DataArrayInt *offsets) const;
593     MEDCOUPLING_EXPORT DataArrayInt *buildExplicitArrOfSliceOnScaledArr(int begin, int stop, int step) const;
594     MEDCOUPLING_EXPORT DataArrayInt *findRangeIdForEachTuple(const DataArrayInt *ranges) const;
595     MEDCOUPLING_EXPORT DataArrayInt *findIdInRangeForEachTuple(const DataArrayInt *ranges) const;
596     MEDCOUPLING_EXPORT DataArrayInt *duplicateEachTupleNTimes(int nbTimes) const;
597     MEDCOUPLING_EXPORT DataArrayInt *getDifferentValues() const;
598     MEDCOUPLING_EXPORT std::vector<DataArrayInt *> partitionByDifferentValues(std::vector<int>& differentIds) const;
599     MEDCOUPLING_EXPORT std::vector< std::pair<int,int> > splitInBalancedSlices(int nbOfSlices) const;
600     MEDCOUPLING_EXPORT void useArray(const int *array, bool ownership, DeallocType type, int nbOfTuple, int nbOfCompo);
601     MEDCOUPLING_EXPORT void useExternalArrayWithRWAccess(const int *array, int nbOfTuple, int nbOfCompo);
602     template<class InputIterator>
603     void insertAtTheEnd(InputIterator first, InputIterator last);
604     MEDCOUPLING_EXPORT void writeOnPlace(std::size_t id, int element0, const int *others, int sizeOfOthers) { _mem.writeOnPlace(id,element0,others,sizeOfOthers); }
605     MEDCOUPLING_EXPORT static DataArrayInt *Add(const DataArrayInt *a1, const DataArrayInt *a2);
606     MEDCOUPLING_EXPORT void addEqual(const DataArrayInt *other);
607     MEDCOUPLING_EXPORT static DataArrayInt *Substract(const DataArrayInt *a1, const DataArrayInt *a2);
608     MEDCOUPLING_EXPORT void substractEqual(const DataArrayInt *other);
609     MEDCOUPLING_EXPORT static DataArrayInt *Multiply(const DataArrayInt *a1, const DataArrayInt *a2);
610     MEDCOUPLING_EXPORT void multiplyEqual(const DataArrayInt *other);
611     MEDCOUPLING_EXPORT static DataArrayInt *Divide(const DataArrayInt *a1, const DataArrayInt *a2);
612     MEDCOUPLING_EXPORT void divideEqual(const DataArrayInt *other);
613     MEDCOUPLING_EXPORT static DataArrayInt *Modulus(const DataArrayInt *a1, const DataArrayInt *a2);
614     MEDCOUPLING_EXPORT void modulusEqual(const DataArrayInt *other);
615     MEDCOUPLING_EXPORT static DataArrayInt *Pow(const DataArrayInt *a1, const DataArrayInt *a2);
616     MEDCOUPLING_EXPORT void powEqual(const DataArrayInt *other);
617     MEDCOUPLING_EXPORT void updateTime() const { }
618     MEDCOUPLING_EXPORT MemArray<int>& accessToMemArray() { return _mem; }
619     MEDCOUPLING_EXPORT const MemArray<int>& accessToMemArray() const { return _mem; }
620   public:
621     MEDCOUPLING_EXPORT static int *CheckAndPreparePermutation(const int *start, const int *end);
622     MEDCOUPLING_EXPORT static DataArrayInt *Range(int begin, int end, int step);
623   public:
624     MEDCOUPLING_EXPORT void getTinySerializationIntInformation(std::vector<int>& tinyInfo) const;
625     MEDCOUPLING_EXPORT void getTinySerializationStrInformation(std::vector<std::string>& tinyInfo) const;
626     MEDCOUPLING_EXPORT bool resizeForUnserialization(const std::vector<int>& tinyInfoI);
627     MEDCOUPLING_EXPORT void finishUnserialization(const std::vector<int>& tinyInfoI, const std::vector<std::string>& tinyInfoS);
628   private:
629     ~DataArrayInt() { }
630     DataArrayInt() { }
631   private:
632     MemArray<int> _mem;
633   };
634
635   class DataArrayIntTuple;
636
637   class DataArrayIntIterator
638   {
639   public:
640     MEDCOUPLING_EXPORT DataArrayIntIterator(DataArrayInt *da);
641     MEDCOUPLING_EXPORT ~DataArrayIntIterator();
642     MEDCOUPLING_EXPORT DataArrayIntTuple *nextt();
643   private:
644     DataArrayInt *_da;
645     int *_pt;
646     int _tuple_id;
647     int _nb_comp;
648     int _nb_tuple;
649   };
650
651   class DataArrayIntTuple
652   {
653   public:
654     MEDCOUPLING_EXPORT DataArrayIntTuple(int *pt, int nbOfComp);
655     MEDCOUPLING_EXPORT std::string repr() const;
656     MEDCOUPLING_EXPORT int getNumberOfCompo() const { return _nb_of_compo; }
657     MEDCOUPLING_EXPORT const int *getConstPointer() const { return  _pt; }
658     MEDCOUPLING_EXPORT int *getPointer() { return _pt; }
659     MEDCOUPLING_EXPORT int intValue() const;
660     MEDCOUPLING_EXPORT DataArrayInt *buildDAInt(int nbOfTuples, int nbOfCompo) const;
661   private:
662     int *_pt;
663     int _nb_of_compo;
664   };
665
666   class DataArrayChar : public DataArray
667   {
668   public:
669     MEDCOUPLING_EXPORT virtual DataArrayChar *buildEmptySpecializedDAChar() const = 0;
670     MEDCOUPLING_EXPORT bool isAllocated() const;
671     MEDCOUPLING_EXPORT void checkAllocated() const;
672     MEDCOUPLING_EXPORT void desallocate();
673     MEDCOUPLING_EXPORT int getNumberOfTuples() const { return _info_on_compo.empty()?0:_mem.getNbOfElem()/getNumberOfComponents(); }
674     MEDCOUPLING_EXPORT std::size_t getNbOfElems() const { return _mem.getNbOfElem(); }
675     MEDCOUPLING_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
676     MEDCOUPLING_EXPORT int getHashCode() const;
677     MEDCOUPLING_EXPORT bool empty() const;
678     MEDCOUPLING_EXPORT void cpyFrom(const DataArrayChar& other);
679     MEDCOUPLING_EXPORT void reserve(std::size_t nbOfElems);
680     MEDCOUPLING_EXPORT void pushBackSilent(char val);
681     MEDCOUPLING_EXPORT void pushBackValsSilent(const char *valsBg, const char *valsEnd);
682     MEDCOUPLING_EXPORT char popBackSilent();
683     MEDCOUPLING_EXPORT void pack() const;
684     MEDCOUPLING_EXPORT std::size_t getNbOfElemAllocated() const { return _mem.getNbOfElemAllocated(); }
685     MEDCOUPLING_EXPORT void alloc(int nbOfTuple, int nbOfCompo=1);
686     MEDCOUPLING_EXPORT void allocIfNecessary(int nbOfTuple, int nbOfCompo);
687     MEDCOUPLING_EXPORT bool isEqual(const DataArrayChar& other) const;
688     MEDCOUPLING_EXPORT virtual bool isEqualIfNotWhy(const DataArrayChar& other, std::string& reason) const;
689     MEDCOUPLING_EXPORT bool isEqualWithoutConsideringStr(const DataArrayChar& other) const;
690     MEDCOUPLING_EXPORT void reverse();
691     MEDCOUPLING_EXPORT void fillWithZero();
692     MEDCOUPLING_EXPORT void fillWithValue(char val);
693     MEDCOUPLING_EXPORT std::string repr() const;
694     MEDCOUPLING_EXPORT std::string reprZip() const;
695     MEDCOUPLING_EXPORT void reAlloc(int nbOfTuples);
696     MEDCOUPLING_EXPORT DataArrayInt *convertToIntArr() const;
697     MEDCOUPLING_EXPORT void renumberInPlace(const int *old2New);
698     MEDCOUPLING_EXPORT void renumberInPlaceR(const int *new2Old);
699     MEDCOUPLING_EXPORT DataArrayChar *renumber(const int *old2New) const;
700     MEDCOUPLING_EXPORT DataArrayChar *renumberR(const int *new2Old) const;
701     MEDCOUPLING_EXPORT DataArrayChar *renumberAndReduce(const int *old2NewBg, int newNbOfTuple) const;
702     MEDCOUPLING_EXPORT DataArrayChar *selectByTupleId(const int *new2OldBg, const int *new2OldEnd) const;
703     MEDCOUPLING_EXPORT DataArrayChar *selectByTupleIdSafe(const int *new2OldBg, const int *new2OldEnd) const;
704     MEDCOUPLING_EXPORT DataArrayChar *selectByTupleId2(int bg, int end, int step) const;
705     MEDCOUPLING_EXPORT bool isUniform(char val) const;
706     MEDCOUPLING_EXPORT void rearrange(int newNbOfCompo);
707     MEDCOUPLING_EXPORT DataArrayChar *substr(int tupleIdBg, int tupleIdEnd=-1) const;
708     MEDCOUPLING_EXPORT DataArrayChar *changeNbOfComponents(int newNbOfComp, char dftValue) const;
709     MEDCOUPLING_EXPORT DataArrayChar *keepSelectedComponents(const std::vector<int>& compoIds) const;
710     MEDCOUPLING_EXPORT void meldWith(const DataArrayChar *other);
711     MEDCOUPLING_EXPORT void setPartOfValues1(const DataArrayChar *a, int bgTuples, int endTuples, int stepTuples, int bgComp, int endComp, int stepComp, bool strictCompoCompare=true);
712     MEDCOUPLING_EXPORT void setPartOfValuesSimple1(char a, int bgTuples, int endTuples, int stepTuples, int bgComp, int endComp, int stepComp);
713     MEDCOUPLING_EXPORT void setPartOfValues2(const DataArrayChar *a, const int *bgTuples, const int *endTuples, const int *bgComp, const int *endComp, bool strictCompoCompare=true);
714     MEDCOUPLING_EXPORT void setPartOfValuesSimple2(char a, const int *bgTuples, const int *endTuples, const int *bgComp, const int *endComp);
715     MEDCOUPLING_EXPORT void setPartOfValues3(const DataArrayChar *a, const int *bgTuples, const int *endTuples, int bgComp, int endComp, int stepComp, bool strictCompoCompare=true);
716     MEDCOUPLING_EXPORT void setPartOfValuesSimple3(char a, const int *bgTuples, const int *endTuples, int bgComp, int endComp, int stepComp);
717     MEDCOUPLING_EXPORT void setPartOfValues4(const DataArrayChar *a, int bgTuples, int endTuples, int stepTuples, const int *bgComp, const int *endComp, bool strictCompoCompare=true);
718     MEDCOUPLING_EXPORT void setPartOfValuesSimple4(char a, int bgTuples, int endTuples, int stepTuples, const int *bgComp, const int *endComp);
719     MEDCOUPLING_EXPORT void setPartOfValuesAdv(const DataArrayChar *a, const DataArrayChar *tuplesSelec);
720     MEDCOUPLING_EXPORT void setContigPartOfSelectedValues(int tupleIdStart, const DataArray *aBase, const DataArrayInt *tuplesSelec);
721     MEDCOUPLING_EXPORT void setContigPartOfSelectedValues2(int tupleIdStart, const DataArray *aBase, int bg, int end2, int step);
722     MEDCOUPLING_EXPORT DataArray *selectByTupleRanges(const std::vector<std::pair<int,int> >& ranges) const;
723     MEDCOUPLING_EXPORT void getTuple(int tupleId, char *res) const { std::copy(_mem.getConstPointerLoc(tupleId*_info_on_compo.size()),_mem.getConstPointerLoc((tupleId+1)*_info_on_compo.size()),res); }
724     MEDCOUPLING_EXPORT char getIJ(int tupleId, int compoId) const { return _mem[tupleId*_info_on_compo.size()+compoId]; }
725     MEDCOUPLING_EXPORT char getIJSafe(int tupleId, int compoId) const;
726     MEDCOUPLING_EXPORT char front() const;
727     MEDCOUPLING_EXPORT char back() const;
728     MEDCOUPLING_EXPORT void setIJ(int tupleId, int compoId, char newVal) { _mem[tupleId*_info_on_compo.size()+compoId]=newVal; declareAsNew(); }
729     MEDCOUPLING_EXPORT void setIJSilent(int tupleId, int compoId, char newVal) { _mem[tupleId*_info_on_compo.size()+compoId]=newVal; }
730     MEDCOUPLING_EXPORT char *getPointer() { return _mem.getPointer(); }
731     MEDCOUPLING_EXPORT const char *getConstPointer() const { return _mem.getConstPointer(); }
732     MEDCOUPLING_EXPORT const char *begin() const { return getConstPointer(); }
733     MEDCOUPLING_EXPORT const char *end() const { return getConstPointer()+getNbOfElems(); }
734     MEDCOUPLING_EXPORT DataArrayInt *getIdsEqual(char val) const;
735     MEDCOUPLING_EXPORT DataArrayInt *getIdsNotEqual(char val) const;
736     MEDCOUPLING_EXPORT int search(const std::vector<char>& vals) const;
737     MEDCOUPLING_EXPORT int locateTuple(const std::vector<char>& tupl) const;
738     MEDCOUPLING_EXPORT int locateValue(char value) const;
739     MEDCOUPLING_EXPORT int locateValue(const std::vector<char>& vals) const;
740     MEDCOUPLING_EXPORT bool presenceOfTuple(const std::vector<char>& tupl) const;
741     MEDCOUPLING_EXPORT bool presenceOfValue(char value) const;
742     MEDCOUPLING_EXPORT bool presenceOfValue(const std::vector<char>& vals) const;
743     MEDCOUPLING_EXPORT char getMaxValue(int& tupleId) const;
744     MEDCOUPLING_EXPORT char getMaxValueInArray() const;
745     MEDCOUPLING_EXPORT char getMinValue(int& tupleId) const;
746     MEDCOUPLING_EXPORT char getMinValueInArray() const;
747     MEDCOUPLING_EXPORT DataArrayInt *getIdsInRange(char vmin, char vmax) const;
748     MEDCOUPLING_EXPORT static DataArrayChar *Aggregate(const DataArrayChar *a1, const DataArrayChar *a2);
749     MEDCOUPLING_EXPORT static DataArrayChar *Aggregate(const std::vector<const DataArrayChar *>& arr);
750     MEDCOUPLING_EXPORT static DataArrayChar *Meld(const DataArrayChar *a1, const DataArrayChar *a2);
751     MEDCOUPLING_EXPORT static DataArrayChar *Meld(const std::vector<const DataArrayChar *>& arr);
752     MEDCOUPLING_EXPORT void useArray(const char *array, bool ownership, DeallocType type, int nbOfTuple, int nbOfCompo);
753     template<class InputIterator>
754     void insertAtTheEnd(InputIterator first, InputIterator last);
755     MEDCOUPLING_EXPORT void useExternalArrayWithRWAccess(const char *array, int nbOfTuple, int nbOfCompo);
756     MEDCOUPLING_EXPORT void updateTime() const { }
757     MEDCOUPLING_EXPORT MemArray<char>& accessToMemArray() { return _mem; }
758     MEDCOUPLING_EXPORT const MemArray<char>& accessToMemArray() const { return _mem; }
759   public:
760     //MEDCOUPLING_EXPORT void getTinySerializationIntInformation(std::vector<int>& tinyInfo) const;
761     //MEDCOUPLING_EXPORT void getTinySerializationStrInformation(std::vector<std::string>& tinyInfo) const;
762     //MEDCOUPLING_EXPORT bool resizeForUnserialization(const std::vector<int>& tinyInfoI);
763     //MEDCOUPLING_EXPORT void finishUnserialization(const std::vector<int>& tinyInfoI, const std::vector<std::string>& tinyInfoS);
764   protected:
765     DataArrayChar() { }
766   protected:
767     MemArray<char> _mem;
768   };
769
770   class DataArrayByteIterator;
771
772   class DataArrayByte : public DataArrayChar
773   {
774   public:
775     MEDCOUPLING_EXPORT static DataArrayByte *New();
776     MEDCOUPLING_EXPORT DataArrayChar *buildEmptySpecializedDAChar() const;
777     MEDCOUPLING_EXPORT DataArrayByteIterator *iterator();
778     MEDCOUPLING_EXPORT DataArrayByte *deepCpy() const;
779     MEDCOUPLING_EXPORT DataArrayByte *performCpy(bool deepCpy) const;
780     MEDCOUPLING_EXPORT char byteValue() const;
781     MEDCOUPLING_EXPORT void reprStream(std::ostream& stream) const;
782     MEDCOUPLING_EXPORT void reprZipStream(std::ostream& stream) const;
783     MEDCOUPLING_EXPORT void reprWithoutNameStream(std::ostream& stream) const;
784     MEDCOUPLING_EXPORT void reprZipWithoutNameStream(std::ostream& stream) const;
785     MEDCOUPLING_EXPORT void reprCppStream(const std::string& varName, std::ostream& stream) const;
786     MEDCOUPLING_EXPORT void reprQuickOverview(std::ostream& stream) const;
787     MEDCOUPLING_EXPORT void reprQuickOverviewData(std::ostream& stream, std::size_t maxNbOfByteInRepr) const;
788     MEDCOUPLING_EXPORT bool isEqualIfNotWhy(const DataArrayChar& other, std::string& reason) const;
789     MEDCOUPLING_EXPORT std::vector<bool> toVectorOfBool() const;
790   private:
791     ~DataArrayByte() { }
792     DataArrayByte() { }
793   };
794
795   class DataArrayByteTuple;
796
797   class DataArrayByteIterator
798   {
799   public:
800     MEDCOUPLING_EXPORT DataArrayByteIterator(DataArrayByte *da);
801     MEDCOUPLING_EXPORT ~DataArrayByteIterator();
802     MEDCOUPLING_EXPORT DataArrayByteTuple *nextt();
803   private:
804     DataArrayByte *_da;
805     char *_pt;
806     int _tuple_id;
807     int _nb_comp;
808     int _nb_tuple;
809   };
810
811   class DataArrayByteTuple
812   {
813   public:
814     MEDCOUPLING_EXPORT DataArrayByteTuple(char *pt, int nbOfComp);
815     MEDCOUPLING_EXPORT std::string repr() const;
816     MEDCOUPLING_EXPORT int getNumberOfCompo() const { return _nb_of_compo; }
817     MEDCOUPLING_EXPORT const char *getConstPointer() const { return  _pt; }
818     MEDCOUPLING_EXPORT char *getPointer() { return _pt; }
819     MEDCOUPLING_EXPORT char byteValue() const;
820     MEDCOUPLING_EXPORT DataArrayByte *buildDAByte(int nbOfTuples, int nbOfCompo) const;
821   private:
822     char *_pt;
823     int _nb_of_compo;
824   };
825
826   class DataArrayAsciiCharIterator;
827
828   class DataArrayAsciiChar : public DataArrayChar
829   {
830   public:
831     MEDCOUPLING_EXPORT static DataArrayAsciiChar *New();
832     MEDCOUPLING_EXPORT static DataArrayAsciiChar *New(const std::string& st);
833     MEDCOUPLING_EXPORT static DataArrayAsciiChar *New(const std::vector<std::string>& vst, char defaultChar);
834     MEDCOUPLING_EXPORT DataArrayChar *buildEmptySpecializedDAChar() const;
835     MEDCOUPLING_EXPORT DataArrayAsciiCharIterator *iterator();
836     MEDCOUPLING_EXPORT DataArrayAsciiChar *deepCpy() const;
837     MEDCOUPLING_EXPORT DataArrayAsciiChar *performCpy(bool deepCpy) const;
838     MEDCOUPLING_EXPORT char asciiCharValue() const;
839     MEDCOUPLING_EXPORT void reprStream(std::ostream& stream) const;
840     MEDCOUPLING_EXPORT void reprZipStream(std::ostream& stream) const;
841     MEDCOUPLING_EXPORT void reprWithoutNameStream(std::ostream& stream) const;
842     MEDCOUPLING_EXPORT void reprZipWithoutNameStream(std::ostream& stream) const;
843     MEDCOUPLING_EXPORT void reprCppStream(const std::string& varName, std::ostream& stream) const;
844     MEDCOUPLING_EXPORT void reprQuickOverview(std::ostream& stream) const;
845     MEDCOUPLING_EXPORT void reprQuickOverviewData(std::ostream& stream, std::size_t maxNbOfByteInRepr) const;
846     MEDCOUPLING_EXPORT bool isEqualIfNotWhy(const DataArrayChar& other, std::string& reason) const;
847   private:
848     ~DataArrayAsciiChar() { }
849     DataArrayAsciiChar() { }
850     DataArrayAsciiChar(const std::string& st);
851     DataArrayAsciiChar(const std::vector<std::string>& vst, char defaultChar);
852   };
853
854   class DataArrayAsciiCharTuple;
855
856   class DataArrayAsciiCharIterator
857   {
858   public:
859     MEDCOUPLING_EXPORT DataArrayAsciiCharIterator(DataArrayAsciiChar *da);
860     MEDCOUPLING_EXPORT ~DataArrayAsciiCharIterator();
861     MEDCOUPLING_EXPORT DataArrayAsciiCharTuple *nextt();
862   private:
863     DataArrayAsciiChar *_da;
864     char *_pt;
865     int _tuple_id;
866     int _nb_comp;
867     int _nb_tuple;
868   };
869
870   class DataArrayAsciiCharTuple
871   {
872   public:
873     MEDCOUPLING_EXPORT DataArrayAsciiCharTuple(char *pt, int nbOfComp);
874     MEDCOUPLING_EXPORT std::string repr() const;
875     MEDCOUPLING_EXPORT int getNumberOfCompo() const { return _nb_of_compo; }
876     MEDCOUPLING_EXPORT const char *getConstPointer() const { return  _pt; }
877     MEDCOUPLING_EXPORT char *getPointer() { return _pt; }
878     MEDCOUPLING_EXPORT char asciiCharValue() const;
879     MEDCOUPLING_EXPORT DataArrayAsciiChar *buildDAAsciiChar(int nbOfTuples, int nbOfCompo) const;
880   private:
881     char *_pt;
882     int _nb_of_compo;
883   };
884
885   template<class InputIterator>
886   void DataArrayDouble::insertAtTheEnd(InputIterator first, InputIterator last)
887   {
888     int nbCompo(getNumberOfComponents());
889     if(nbCompo==1)
890       _mem.insertAtTheEnd(first,last);
891     else if(nbCompo==0)
892       {
893         _info_on_compo.resize(1);
894         _mem.insertAtTheEnd(first,last);
895       }
896     else
897       throw INTERP_KERNEL::Exception("DataArrayDouble::insertAtTheEnd : not available for DataArrayDouble with number of components different than 1 !");
898   }
899
900   template<class InputIterator>
901   void DataArrayInt::insertAtTheEnd(InputIterator first, InputIterator last)
902   {
903     int nbCompo(getNumberOfComponents());
904     if(nbCompo==1)
905       _mem.insertAtTheEnd(first,last);
906     else if(nbCompo==0)
907       {
908         _info_on_compo.resize(1);
909         _mem.insertAtTheEnd(first,last);
910       }
911     else
912       throw INTERP_KERNEL::Exception("DataArrayInt::insertAtTheEnd : not available for DataArrayInt with number of components different than 1 !");
913   }
914
915   template<class InputIterator>
916   void DataArrayChar::insertAtTheEnd(InputIterator first, InputIterator last)
917   {
918     int nbCompo(getNumberOfComponents());
919     if(nbCompo==1)
920       _mem.insertAtTheEnd(first,last);
921     else if(nbCompo==0)
922       {
923         _info_on_compo.resize(1);
924         _mem.insertAtTheEnd(first,last);
925       }
926     else
927       throw INTERP_KERNEL::Exception("DataArrayChar::insertAtTheEnd : not available for DataArrayChar with number of components different than 1 !");
928   }
929 }
930
931 #endif