1 // Copyright (C) 2007-2016 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, or (at your option) any later version.
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 __MEDCOUPLINGTYPEMAPS_I__
22 #define __MEDCOUPLINGTYPEMAPS_I__
24 #include "MEDCouplingDataArrayTypemaps.i"
25 #include "MEDCouplingUMesh.hxx"
26 #include "MEDCouplingCMesh.hxx"
27 #include "MEDCouplingIMesh.hxx"
28 #include "MEDCouplingCurveLinearMesh.hxx"
29 #include "MEDCouplingMappedExtrudedMesh.hxx"
30 #include "MEDCoupling1GTUMesh.hxx"
31 #include "MEDCouplingFieldDiscretization.hxx"
32 #include "MEDCouplingMultiFields.hxx"
33 #include "MEDCouplingPartDefinition.hxx"
34 #include "MEDCouplingCartesianAMRMesh.hxx"
36 static PyObject *convertMesh(MEDCoupling::MEDCouplingMesh *mesh, int owner)
44 if(dynamic_cast<MEDCoupling::MEDCouplingUMesh *>(mesh))
45 ret=SWIG_NewPointerObj((void*)mesh,SWIGTYPE_p_MEDCoupling__MEDCouplingUMesh,owner);
46 if(dynamic_cast<MEDCoupling::MEDCoupling1SGTUMesh *>(mesh))
47 ret=SWIG_NewPointerObj((void*)mesh,SWIGTYPE_p_MEDCoupling__MEDCoupling1SGTUMesh,owner);
48 if(dynamic_cast<MEDCoupling::MEDCoupling1DGTUMesh *>(mesh))
49 ret=SWIG_NewPointerObj((void*)mesh,SWIGTYPE_p_MEDCoupling__MEDCoupling1DGTUMesh,owner);
50 if(dynamic_cast<MEDCoupling::MEDCouplingMappedExtrudedMesh *>(mesh))
51 ret=SWIG_NewPointerObj((void*)mesh,SWIGTYPE_p_MEDCoupling__MEDCouplingMappedExtrudedMesh,owner);
52 if(dynamic_cast<MEDCoupling::MEDCouplingCMesh *>(mesh))
53 ret=SWIG_NewPointerObj((void*)mesh,SWIGTYPE_p_MEDCoupling__MEDCouplingCMesh,owner);
54 if(dynamic_cast<MEDCoupling::MEDCouplingCurveLinearMesh *>(mesh))
55 ret=SWIG_NewPointerObj((void*)mesh,SWIGTYPE_p_MEDCoupling__MEDCouplingCurveLinearMesh,owner);
56 if(dynamic_cast<MEDCoupling::MEDCouplingIMesh *>(mesh))
57 ret=SWIG_NewPointerObj((void*)mesh,SWIGTYPE_p_MEDCoupling__MEDCouplingIMesh,owner);
59 throw INTERP_KERNEL::Exception("Not recognized type of mesh on downcast !");
63 static PyObject *convertFieldDiscretization(MEDCoupling::MEDCouplingFieldDiscretization *fd, int owner)
71 if(dynamic_cast<MEDCoupling::MEDCouplingFieldDiscretizationP0 *>(fd))
72 ret=SWIG_NewPointerObj(reinterpret_cast<void*>(fd),SWIGTYPE_p_MEDCoupling__MEDCouplingFieldDiscretizationP0,owner);
73 if(dynamic_cast<MEDCoupling::MEDCouplingFieldDiscretizationP1 *>(fd))
74 ret=SWIG_NewPointerObj(reinterpret_cast<void*>(fd),SWIGTYPE_p_MEDCoupling__MEDCouplingFieldDiscretizationP1,owner);
75 if(dynamic_cast<MEDCoupling::MEDCouplingFieldDiscretizationGauss *>(fd))
76 ret=SWIG_NewPointerObj(reinterpret_cast<void*>(fd),SWIGTYPE_p_MEDCoupling__MEDCouplingFieldDiscretizationGauss,owner);
77 if(dynamic_cast<MEDCoupling::MEDCouplingFieldDiscretizationGaussNE *>(fd))
78 ret=SWIG_NewPointerObj(reinterpret_cast<void*>(fd),SWIGTYPE_p_MEDCoupling__MEDCouplingFieldDiscretizationGaussNE,owner);
79 if(dynamic_cast<MEDCoupling::MEDCouplingFieldDiscretizationKriging *>(fd))
80 ret=SWIG_NewPointerObj(reinterpret_cast<void*>(fd),SWIGTYPE_p_MEDCoupling__MEDCouplingFieldDiscretizationKriging,owner);
82 throw INTERP_KERNEL::Exception("Not recognized type of field discretization on downcast !");
86 static PyObject *convertField(MEDCoupling::MEDCouplingField *f, int owner)
94 if(dynamic_cast<MEDCoupling::MEDCouplingFieldDouble *>(f))
95 ret=SWIG_NewPointerObj(reinterpret_cast<void*>(f),SWIGTYPE_p_MEDCoupling__MEDCouplingFieldDouble,owner);
96 if(dynamic_cast<MEDCoupling::MEDCouplingFieldInt *>(f))
97 ret=SWIG_NewPointerObj(reinterpret_cast<void*>(f),SWIGTYPE_p_MEDCoupling__MEDCouplingFieldInt,owner);
98 if(dynamic_cast<MEDCoupling::MEDCouplingFieldFloat *>(f))
99 ret=SWIG_NewPointerObj(reinterpret_cast<void*>(f),SWIGTYPE_p_MEDCoupling__MEDCouplingFieldFloat,owner);
101 throw INTERP_KERNEL::Exception("Not recognized type of field on downcast !");
105 static PyObject* convertMultiFields(MEDCoupling::MEDCouplingMultiFields *mfs, int owner)
113 if(dynamic_cast<MEDCoupling::MEDCouplingFieldOverTime *>(mfs))
114 ret=SWIG_NewPointerObj((void*)mfs,SWIGTYPE_p_MEDCoupling__MEDCouplingFieldOverTime,owner);
116 ret=SWIG_NewPointerObj((void*)mfs,SWIGTYPE_p_MEDCoupling__MEDCouplingMultiFields,owner);
120 static PyObject *convertCartesianAMRMesh(MEDCoupling::MEDCouplingCartesianAMRMeshGen *mesh, int owner)
127 if(dynamic_cast<MEDCoupling::MEDCouplingCartesianAMRMeshSub *>(mesh))
129 return SWIG_NewPointerObj(reinterpret_cast<void*>(mesh),SWIGTYPE_p_MEDCoupling__MEDCouplingCartesianAMRMeshSub,owner);
131 if(dynamic_cast<MEDCoupling::MEDCouplingCartesianAMRMesh *>(mesh))
133 return SWIG_NewPointerObj(reinterpret_cast<void*>(mesh),SWIGTYPE_p_MEDCoupling__MEDCouplingCartesianAMRMesh,owner);
135 throw INTERP_KERNEL::Exception("convertCartesianAMRMesh wrap : unrecognized type of cartesian AMR mesh !");
138 static PyObject *convertDataForGodFather(MEDCoupling::MEDCouplingDataForGodFather *data, int owner)
145 if(dynamic_cast<MEDCoupling::MEDCouplingAMRAttribute *>(data))
147 return SWIG_NewPointerObj(reinterpret_cast<void*>(data),SWIGTYPE_p_MEDCoupling__MEDCouplingAMRAttribute,owner);
149 throw INTERP_KERNEL::Exception("convertDataForGodFather wrap : unrecognized data type for AMR !");
152 static PyObject *convertCartesianAMRPatch(MEDCoupling::MEDCouplingCartesianAMRPatchGen *patch, int owner) throw(INTERP_KERNEL::Exception)
159 if(dynamic_cast<MEDCoupling::MEDCouplingCartesianAMRPatchGF *>(patch))
161 return SWIG_NewPointerObj(reinterpret_cast<void*>(patch),SWIGTYPE_p_MEDCoupling__MEDCouplingCartesianAMRPatchGF,owner);
163 if(dynamic_cast<MEDCoupling::MEDCouplingCartesianAMRPatch *>(patch))
165 return SWIG_NewPointerObj(reinterpret_cast<void*>(patch),SWIGTYPE_p_MEDCoupling__MEDCouplingCartesianAMRPatch,owner);
167 throw INTERP_KERNEL::Exception("convertCartesianAMRPatch wrap : unrecognized type of cartesian AMR patch !");
170 static MEDCoupling::MEDCouplingFieldDouble *MEDCoupling_MEDCouplingFieldDouble___add__Impl(MEDCoupling::MEDCouplingFieldDouble *self, PyObject *obj) throw(INTERP_KERNEL::Exception)
172 const char msg[]="Unexpected situation in MEDCouplingFieldDouble.__add__ ! Expecting a not null MEDCouplingFieldDouble or DataArrayDouble or DataArrayDoubleTuple instance, or a list of double, or a double.";
173 const char msg2[]="in MEDCouplingFieldDouble.__add__ : self field has no Array of values set !";
176 if(SWIG_IsOK(SWIG_ConvertPtr(obj,&argp,SWIGTYPE_p_MEDCoupling__MEDCouplingFieldDouble,0|0)))
178 MEDCoupling::MEDCouplingFieldDouble *other=reinterpret_cast< MEDCoupling::MEDCouplingFieldDouble * >(argp);
180 return (*self)+(*other);
182 throw INTERP_KERNEL::Exception(msg);
186 MEDCoupling::DataArrayDouble *a;
187 MEDCoupling::DataArrayDoubleTuple *aa;
188 std::vector<double> bb;
190 convertDoubleStarLikePyObjToCpp_2(obj,sw,val,a,aa,bb);
195 if(!self->getArray())
196 throw INTERP_KERNEL::Exception(msg2);
197 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> ret=self->getArray()->deepCopy();
198 ret->applyLin(1.,val);
199 MEDCoupling::MCAuto<MEDCoupling::MEDCouplingFieldDouble> ret2=self->clone(false);
205 if(!self->getArray())
206 throw INTERP_KERNEL::Exception(msg2);
207 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> ret=MEDCoupling::DataArrayDouble::Add(self->getArray(),a);
208 MEDCoupling::MCAuto<MEDCoupling::MEDCouplingFieldDouble> ret2=self->clone(false);
214 if(!self->getArray())
215 throw INTERP_KERNEL::Exception(msg2);
216 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> aaa=aa->buildDADouble(1,self->getNumberOfComponents());
217 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> ret=MEDCoupling::DataArrayDouble::Add(self->getArray(),aaa);
218 MEDCoupling::MCAuto<MEDCoupling::MEDCouplingFieldDouble> ret2=self->clone(false);
224 if(!self->getArray())
225 throw INTERP_KERNEL::Exception(msg2);
226 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> aaa=MEDCoupling::DataArrayDouble::New(); aaa->useArray(&bb[0],false,MEDCoupling::CPP_DEALLOC,1,(int)bb.size());
227 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> ret=MEDCoupling::DataArrayDouble::Add(self->getArray(),aaa);
228 MEDCoupling::MCAuto<MEDCoupling::MEDCouplingFieldDouble> ret2=self->clone(false);
233 { throw INTERP_KERNEL::Exception(msg); }
237 static MEDCoupling::MEDCouplingFieldDouble *MEDCoupling_MEDCouplingFieldDouble___radd__Impl(MEDCoupling::MEDCouplingFieldDouble *self, PyObject *obj) throw(INTERP_KERNEL::Exception)
239 return MEDCoupling_MEDCouplingFieldDouble___add__Impl(self,obj);
242 static MEDCoupling::MEDCouplingFieldDouble *MEDCoupling_MEDCouplingFieldDouble___rsub__Impl(MEDCoupling::MEDCouplingFieldDouble *self, PyObject *obj) throw(INTERP_KERNEL::Exception)
244 const char msg[]="Unexpected situation in MEDCouplingFieldDouble.__rsub__ ! Expecting a not null MEDCouplingFieldDouble or DataArrayDouble or DataArrayDoubleTuple instance, or a list of double, or a double.";
245 const char msg2[]="in MEDCouplingFieldDouble.__rsub__ : self field has no Array of values set !";
248 if(SWIG_IsOK(SWIG_ConvertPtr(obj,&argp,SWIGTYPE_p_MEDCoupling__MEDCouplingFieldDouble,0|0)))
250 MEDCoupling::MEDCouplingFieldDouble *other=reinterpret_cast< MEDCoupling::MEDCouplingFieldDouble * >(argp);
252 return (*other)-(*self);
254 throw INTERP_KERNEL::Exception(msg);
258 MEDCoupling::DataArrayDouble *a;
259 MEDCoupling::DataArrayDoubleTuple *aa;
260 std::vector<double> bb;
262 convertDoubleStarLikePyObjToCpp_2(obj,sw,val,a,aa,bb);
267 if(!self->getArray())
268 throw INTERP_KERNEL::Exception(msg2);
269 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> ret=self->getArray()->deepCopy();
270 ret->applyLin(-1.,val);
271 MEDCoupling::MCAuto<MEDCoupling::MEDCouplingFieldDouble> ret2=self->clone(false);
277 if(!self->getArray())
278 throw INTERP_KERNEL::Exception(msg2);
279 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> ret=MEDCoupling::DataArrayDouble::Substract(a,self->getArray());
280 MEDCoupling::MCAuto<MEDCoupling::MEDCouplingFieldDouble> ret2=self->clone(false);
286 if(!self->getArray())
287 throw INTERP_KERNEL::Exception(msg2);
288 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> aaa=aa->buildDADouble(1,self->getNumberOfComponents());
289 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> ret=MEDCoupling::DataArrayDouble::Substract(aaa,self->getArray());
290 MEDCoupling::MCAuto<MEDCoupling::MEDCouplingFieldDouble> ret2=self->clone(false);
296 if(!self->getArray())
297 throw INTERP_KERNEL::Exception(msg2);
298 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> aaa=MEDCoupling::DataArrayDouble::New(); aaa->useArray(&bb[0],false,MEDCoupling::CPP_DEALLOC,1,(int)bb.size());
299 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> ret=MEDCoupling::DataArrayDouble::Substract(aaa,self->getArray());
300 MEDCoupling::MCAuto<MEDCoupling::MEDCouplingFieldDouble> ret2=self->clone(false);
305 { throw INTERP_KERNEL::Exception(msg); }
309 static MEDCoupling::MEDCouplingFieldDouble *MEDCoupling_MEDCouplingFieldDouble___mul__Impl(MEDCoupling::MEDCouplingFieldDouble *self, PyObject *obj) throw(INTERP_KERNEL::Exception)
311 const char msg[]="Unexpected situation in MEDCouplingFieldDouble.__mul__ ! Expecting a not null MEDCouplingFieldDouble or DataArrayDouble or DataArrayDoubleTuple instance, or a list of double, or a double.";
312 const char msg2[]="in MEDCouplingFieldDouble.__mul__ : self field has no Array of values set !";
315 if(SWIG_IsOK(SWIG_ConvertPtr(obj,&argp,SWIGTYPE_p_MEDCoupling__MEDCouplingFieldDouble,0|0)))
317 MEDCoupling::MEDCouplingFieldDouble *other=reinterpret_cast< MEDCoupling::MEDCouplingFieldDouble * >(argp);
319 return (*self)*(*other);
321 throw INTERP_KERNEL::Exception(msg);
325 MEDCoupling::DataArrayDouble *a;
326 MEDCoupling::DataArrayDoubleTuple *aa;
327 std::vector<double> bb;
329 convertDoubleStarLikePyObjToCpp_2(obj,sw,val,a,aa,bb);
334 if(!self->getArray())
335 throw INTERP_KERNEL::Exception(msg2);
336 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> ret=self->getArray()->deepCopy();
337 ret->applyLin(val,0.);
338 MEDCoupling::MCAuto<MEDCoupling::MEDCouplingFieldDouble> ret2=self->clone(false);
344 if(!self->getArray())
345 throw INTERP_KERNEL::Exception(msg2);
346 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> ret=MEDCoupling::DataArrayDouble::Multiply(self->getArray(),a);
347 MEDCoupling::MCAuto<MEDCoupling::MEDCouplingFieldDouble> ret2=self->clone(false);
353 if(!self->getArray())
354 throw INTERP_KERNEL::Exception(msg2);
355 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> aaa=aa->buildDADouble(1,self->getNumberOfComponents());
356 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> ret=MEDCoupling::DataArrayDouble::Multiply(self->getArray(),aaa);
357 MEDCoupling::MCAuto<MEDCoupling::MEDCouplingFieldDouble> ret2=self->clone(false);
363 if(!self->getArray())
364 throw INTERP_KERNEL::Exception(msg2);
365 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> aaa=MEDCoupling::DataArrayDouble::New(); aaa->useArray(&bb[0],false,MEDCoupling::CPP_DEALLOC,1,(int)bb.size());
366 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> ret=MEDCoupling::DataArrayDouble::Multiply(self->getArray(),aaa);
367 MEDCoupling::MCAuto<MEDCoupling::MEDCouplingFieldDouble> ret2=self->clone(false);
372 { throw INTERP_KERNEL::Exception(msg); }
376 MEDCoupling::MEDCouplingFieldDouble *MEDCoupling_MEDCouplingFieldDouble___rmul__Impl(MEDCoupling::MEDCouplingFieldDouble *self, PyObject *obj) throw(INTERP_KERNEL::Exception)
378 return MEDCoupling_MEDCouplingFieldDouble___mul__Impl(self,obj);
381 MEDCoupling::MEDCouplingFieldDouble *MEDCoupling_MEDCouplingFieldDouble___rdiv__Impl(MEDCoupling::MEDCouplingFieldDouble *self, PyObject *obj) throw(INTERP_KERNEL::Exception)
383 const char msg[]="Unexpected situation in MEDCouplingFieldDouble.__rdiv__ ! Expecting a not null MEDCouplingFieldDouble or DataArrayDouble or DataArrayDoubleTuple instance, or a list of double, or a double.";
384 const char msg2[]="in MEDCouplingFieldDouble.__div__ : self field has no Array of values set !";
387 if(SWIG_IsOK(SWIG_ConvertPtr(obj,&argp,SWIGTYPE_p_MEDCoupling__MEDCouplingFieldDouble,0|0)))
389 MEDCoupling::MEDCouplingFieldDouble *other=reinterpret_cast< MEDCoupling::MEDCouplingFieldDouble * >(argp);
391 return (*other)/(*self);
393 throw INTERP_KERNEL::Exception(msg);
397 MEDCoupling::DataArrayDouble *a;
398 MEDCoupling::DataArrayDoubleTuple *aa;
399 std::vector<double> bb;
401 convertDoubleStarLikePyObjToCpp_2(obj,sw,val,a,aa,bb);
406 if(!self->getArray())
407 throw INTERP_KERNEL::Exception(msg2);
408 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> ret=self->getArray()->deepCopy();
410 MEDCoupling::MCAuto<MEDCoupling::MEDCouplingFieldDouble> ret2=self->clone(false);
416 if(!self->getArray())
417 throw INTERP_KERNEL::Exception(msg2);
418 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> ret=MEDCoupling::DataArrayDouble::Divide(a,self->getArray());
419 MEDCoupling::MCAuto<MEDCoupling::MEDCouplingFieldDouble> ret2=self->clone(false);
425 if(!self->getArray())
426 throw INTERP_KERNEL::Exception(msg2);
427 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> aaa=aa->buildDADouble(1,self->getNumberOfComponents());
428 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> ret=MEDCoupling::DataArrayDouble::Divide(aaa,self->getArray());
429 MEDCoupling::MCAuto<MEDCoupling::MEDCouplingFieldDouble> ret2=self->clone(false);
435 if(!self->getArray())
436 throw INTERP_KERNEL::Exception(msg2);
437 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> aaa=MEDCoupling::DataArrayDouble::New(); aaa->useArray(&bb[0],false,MEDCoupling::CPP_DEALLOC,1,(int)bb.size());
438 MEDCoupling::MCAuto<MEDCoupling::DataArrayDouble> ret=MEDCoupling::DataArrayDouble::Divide(aaa,self->getArray());
439 MEDCoupling::MCAuto<MEDCoupling::MEDCouplingFieldDouble> ret2=self->clone(false);
444 { throw INTERP_KERNEL::Exception(msg); }
449 typename MEDCoupling::Traits<T>::FieldType *fieldT_buildSubPart(const MEDCoupling::MEDCouplingFieldT<T> *self, PyObject *li)
453 std::vector<int> multiVal;
454 std::pair<int, std::pair<int,int> > slic;
455 MEDCoupling::DataArrayInt *daIntTyypp=0;
456 const MEDCoupling::MEDCouplingMesh *mesh=self->getMesh();
458 throw INTERP_KERNEL::Exception("MEDCouplingFieldDouble::buildSubPart : field lies on a null mesh !");
459 int nbc=mesh->getNumberOfCells();
460 convertIntStarOrSliceLikePyObjToCpp(li,nbc,sw,singleVal,multiVal,slic,daIntTyypp);
467 std::ostringstream oss;
468 oss << "Requesting for cell id " << singleVal << " having only " << nbc << " cells !";
469 throw INTERP_KERNEL::Exception(oss.str().c_str());
472 return self->buildSubPart(&singleVal,&singleVal+1);
477 int tmp=nbc+singleVal;
478 return self->buildSubPart(&tmp,&tmp+1);
482 std::ostringstream oss;
483 oss << "Requesting for cell id " << singleVal << " having only " << nbc << " cells !";
484 throw INTERP_KERNEL::Exception(oss.str().c_str());
490 return self->buildSubPart(&multiVal[0],&multiVal[0]+multiVal.size());
494 return self->buildSubPartRange(slic.first,slic.second.first,slic.second.second);
499 throw INTERP_KERNEL::Exception("MEDCouplingFieldDouble::buildSubPart : null instance has been given in input !");
500 daIntTyypp->checkAllocated();
501 return self->buildSubPart(daIntTyypp->begin(),daIntTyypp->end());
504 throw INTERP_KERNEL::Exception("MEDCouplingFieldDouble::buildSubPart : unrecognized type in input ! Possibilities are : int, list or tuple of int DataArrayInt instance !");
509 typename MEDCoupling::Traits<T>::FieldType *fieldT__getitem__(const MEDCoupling::MEDCouplingFieldT<T> *self, PyObject *li)
511 const char msg[]="MEDCouplingFieldDouble::__getitem__ : invalid call Available API are : \n-myField[dataArrayInt]\n-myField[slice]\n-myField[pythonListOfCellIds]\n-myField[integer]\n-myField[dataArrayInt,1]\n-myField[slice,1]\n-myField[pythonListOfCellIds,1]\n-myField[integer,1]\n";
512 if(PyTuple_Check(li))
514 Py_ssize_t sz=PyTuple_Size(li);
516 throw INTERP_KERNEL::Exception(msg);
517 PyObject *elt0=PyTuple_GetItem(li,0),*elt1=PyTuple_GetItem(li,1);
520 std::vector<int> multiVal;
521 std::pair<int, std::pair<int,int> > slic;
522 MEDCoupling::DataArrayInt *daIntTyypp=0;
523 if(!self->getArray())
524 throw INTERP_KERNEL::Exception("MEDCouplingFieldDouble::__getitem__ : no array set on field to deduce number of components !");
526 { convertIntStarOrSliceLikePyObjToCpp(elt1,self->getArray()->getNumberOfComponents(),sw,singleVal,multiVal,slic,daIntTyypp); }
527 catch(INTERP_KERNEL::Exception& e)
528 { std::ostringstream oss; oss << "MEDCouplingFieldDouble::__getitem__ : invalid type in 2nd parameter (compo) !" << e.what(); throw INTERP_KERNEL::Exception(oss.str().c_str()); }
529 typename MEDCoupling::MCAuto< typename MEDCoupling::Traits<T>::FieldType > ret0(fieldT_buildSubPart<T>(self,elt0));
530 typename MEDCoupling::Traits<T>::ArrayType *ret0Arr=ret0->getArray();
532 throw INTERP_KERNEL::Exception("MEDCouplingFieldDouble::__getitem__ : no array exists to apply restriction on component on it !");
537 std::vector<int> v2(1,singleVal);
538 MEDCoupling::MCAuto< typename MEDCoupling::Traits<T>::ArrayType > aarr(ret0Arr->keepSelectedComponents(v2));
539 ret0->setArray(aarr);
544 MEDCoupling::MCAuto< typename MEDCoupling::Traits<T>::ArrayType > aarr(ret0Arr->keepSelectedComponents(multiVal));
545 ret0->setArray(aarr);
550 int nbOfComp(MEDCoupling::DataArray::GetNumberOfItemGivenBESRelative(slic.first,slic.second.first,slic.second.second,"MEDCouplingFieldDouble::__getitem__ : invalid range in 2nd parameter (components) !"));
551 std::vector<int> v2(nbOfComp);
552 for(int i=0;i<nbOfComp;i++)
553 v2[i]=slic.first+i*slic.second.second;
554 MEDCoupling::MCAuto< typename MEDCoupling::Traits<T>::ArrayType > aarr(ret0Arr->keepSelectedComponents(v2));
555 ret0->setArray(aarr);
559 throw INTERP_KERNEL::Exception(msg);
563 return fieldT_buildSubPart<T>(self,li);
566 template<class FIELDT>
567 PyObject *field_getTinySerializationInformation(const FIELDT *self)
569 std::vector<double> a0;
571 std::vector<std::string> a2;
572 self->getTinySerializationDbleInformation(a0);
573 self->getTinySerializationIntInformation(a1);
574 self->getTinySerializationStrInformation(a2);
576 PyObject *ret(PyTuple_New(3));
577 PyTuple_SetItem(ret,0,convertDblArrToPyList2(a0));
578 PyTuple_SetItem(ret,1,convertIntArrToPyList2(a1));
580 PyObject *ret2(PyList_New(sz));
582 for(int i=0;i<sz;i++)
583 PyList_SetItem(ret2,i,PyString_FromString(a2[i].c_str()));
585 PyTuple_SetItem(ret,2,ret2);
590 PyObject *field_serialize(const typename MEDCoupling::Traits<T>::FieldType *self)
592 MEDCoupling::DataArrayInt *ret0(0);
593 std::vector<typename MEDCoupling::Traits<T>::ArrayType *> ret1;
594 self->serialize(ret0,ret1);
597 std::size_t sz(ret1.size());
598 PyObject *ret(PyTuple_New(2));
599 PyTuple_SetItem(ret,0,SWIG_NewPointerObj(SWIG_as_voidptr(ret0),SWIGTYPE_p_MEDCoupling__DataArrayInt, SWIG_POINTER_OWN | 0 ));
600 PyObject *ret1Py(PyList_New(sz));
601 for(std::size_t i=0;i<sz;i++)
605 PyList_SetItem(ret1Py,i,convertArray(ret1[i],SWIG_POINTER_OWN | 0));
607 PyTuple_SetItem(ret,1,ret1Py);
611 template<class FIELDT>
612 PyObject *field__getnewargs__(FIELDT *self)
614 self->checkConsistencyLight();
615 PyObject *ret(PyTuple_New(1));
616 PyObject *ret0(PyDict_New());
618 PyObject *a(PyInt_FromLong(0)),*b(PyInt_FromLong(self->getTypeOfField())),*c(PyInt_FromLong(self->getTimeDiscretization()));
619 PyObject *d(PyTuple_New(2)); PyTuple_SetItem(d,0,b); PyTuple_SetItem(d,1,c);
620 PyDict_SetItem(ret0,a,d);
621 Py_DECREF(a); Py_DECREF(d);
623 PyTuple_SetItem(ret,0,ret0);
627 template<class FIELDT>
628 PyObject *field__getstate__(const FIELDT *self, PyObject *(*tinyserial)(const FIELDT *), PyObject *(*bigserial)(const FIELDT *))
630 self->checkConsistencyLight();
631 PyObject *ret0(tinyserial(self));
632 PyObject *ret1(bigserial(self));
633 const MEDCoupling::MEDCouplingMesh *mesh(self->getMesh());
636 PyObject *ret(PyTuple_New(3));
637 PyTuple_SetItem(ret,0,ret0);
638 PyTuple_SetItem(ret,1,ret1);
639 PyTuple_SetItem(ret,2,convertMesh(const_cast<MEDCoupling::MEDCouplingMesh *>(mesh),SWIG_POINTER_OWN | 0 ));
644 void field__setstate__(typename MEDCoupling::Traits<T>::FieldType *self, PyObject *inp)
646 static const char MSG[]="MEDCouplingFieldDouble.__setstate__ : expected input is a tuple of size 3 !";
647 if(!PyTuple_Check(inp))
648 throw INTERP_KERNEL::Exception(MSG);
649 int sz(PyTuple_Size(inp));
651 throw INTERP_KERNEL::Exception(MSG);
653 PyObject *elt2(PyTuple_GetItem(inp,2));
655 int status(SWIG_ConvertPtr(elt2,&argp,SWIGTYPE_p_MEDCoupling__MEDCouplingMesh,0|0));
656 if(!SWIG_IsOK(status))
657 throw INTERP_KERNEL::Exception(MSG);
658 self->setMesh(reinterpret_cast< const MEDCoupling::MEDCouplingMesh * >(argp));
660 PyObject *elt0(PyTuple_GetItem(inp,0));
661 PyObject *elt1(PyTuple_GetItem(inp,1));
662 std::vector<double> a0;
664 std::vector<std::string> a2;
665 MEDCoupling::DataArrayInt *b0(0);
666 std::vector<typename MEDCoupling::Traits<T>::ArrayType *>b1;
668 if(!PyTuple_Check(elt0) && PyTuple_Size(elt0)!=3)
669 throw INTERP_KERNEL::Exception(MSG);
670 PyObject *a0py(PyTuple_GetItem(elt0,0)),*a1py(PyTuple_GetItem(elt0,1)),*a2py(PyTuple_GetItem(elt0,2));
672 fillArrayWithPyListDbl3(a0py,tmp,a0);
673 convertPyToNewIntArr3(a1py,a1);
674 fillStringVector(a2py,a2);
677 if(!PyTuple_Check(elt1) && PyTuple_Size(elt1)!=2)
678 throw INTERP_KERNEL::Exception(MSG);
679 PyObject *b0py(PyTuple_GetItem(elt1,0)),*b1py(PyTuple_GetItem(elt1,1));
681 int status(SWIG_ConvertPtr(b0py,&argp,SWIGTYPE_p_MEDCoupling__DataArrayInt,0|0));
682 if(!SWIG_IsOK(status))
683 throw INTERP_KERNEL::Exception(MSG);
684 b0=reinterpret_cast<MEDCoupling::DataArrayInt *>(argp);
685 convertFromPyObjVectorOfObj<typename MEDCoupling::Traits<T>::ArrayType *>(b1py,SWIGTITraits<T>::TI,MEDCoupling::Traits<T>::ArrayTypeName,b1);
687 self->checkForUnserialization(a1,b0,b1);
688 // useless here to call resizeForUnserialization because arrays are well resized.
689 self->finishUnserialization(a1,a0,a2);