-// Copyright (C) 2007-2020 CEA/DEN, EDF R&D
+// Copyright (C) 2007-2022 CEA/DEN, EDF R&D
//
// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
#include "ParaMESH.hxx"
#include "ParaFIELD.hxx"
#include "ComponentTopology.hxx"
-#include "ICoCoMEDField.hxx"
#include "ParaMEDLoader.hxx"
+#include "ICoCoMEDDoubleField.hxx"
#include "MEDLoader.hxx"
+#include "MEDCouplingUMesh.hxx"
#include "TestInterpKernelUtils.hxx"
MEDCoupling::InterpKernelDEC dec (*source_group,*target_group);
- MEDCoupling::MEDCouplingUMesh* mesh;
- MEDCoupling::ParaMESH* paramesh;
- MEDCoupling::ParaFIELD* parafield;
- ICoCo::MEDField* icocofield ;
+ MEDCoupling::MEDCouplingUMesh* mesh = nullptr;
+ MEDCoupling::ParaMESH* paramesh = nullptr;
+ MEDCoupling::ParaFIELD* parafield = nullptr;
+ ICoCo::MEDDoubleField* icocofield = nullptr;
string filename_xml1 = "square1_split";
string filename_xml2 = "square2_split";
value[ielem]=1.0;
// ICoCo::Field* icocofield=new ICoCo::MEDField(paramesh,parafield);
- icocofield=new ICoCo::MEDField(parafield->getField());
+ icocofield=new ICoCo::MEDDoubleField(parafield->getField());
dec.setMethod(srcMeth);
dec.attachLocalField(icocofield);
}
for(int ielem=0; ielem<nb_local;ielem++)
value[ielem]=0.0;
// ICoCo::Field* icocofield=new ICoCo::MEDField(paramesh,parafield);
- icocofield=new ICoCo::MEDField(parafield->getField());
+ icocofield=new ICoCo::MEDDoubleField(parafield->getField());
dec.setMethod(targetMeth);
dec.attachLocalField(icocofield);
}
MEDCoupling::InterpKernelDEC dec (*source_group,*target_group);
- MEDCoupling::MEDCouplingUMesh* mesh;
- MEDCoupling::MEDCouplingFieldDouble* mcfield;
+ MEDCoupling::MEDCouplingUMesh* mesh = nullptr;
+ MEDCoupling::MEDCouplingFieldDouble* mcfield = nullptr;
string filename_xml1 = "square1_split";
string filename_xml2 = "square2_split";
MEDCoupling::InterpKernelDEC dec (*source_group,*target_group);
- MEDCoupling::MEDCouplingUMesh* mesh;
- MEDCoupling::ParaMESH* paramesh;
- MEDCoupling::ParaFIELD* parafield;
- ICoCo::MEDField* icocofield ;
+ MEDCoupling::MEDCouplingUMesh* mesh = nullptr;
+ MEDCoupling::ParaMESH* paramesh = nullptr;
+ MEDCoupling::ParaFIELD* parafield = nullptr;
+ ICoCo::MEDDoubleField* icocofield = nullptr;
char * tmp_dir_c = getenv("TMP");
string tmp_dir;
value[ielem]=1.0;
// ICoCo::Field* icocofield=new ICoCo::MEDField(paramesh,parafield);
- icocofield=new ICoCo::MEDField(parafield->getField());
+ icocofield=new ICoCo::MEDDoubleField(parafield->getField());
dec.setMethod(srcMeth);
dec.attachLocalField(icocofield);
}
for(int ielem=0; ielem<nb_local;ielem++)
value[ielem]=0.0;
// ICoCo::Field* icocofield=new ICoCo::MEDField(paramesh,parafield);
- icocofield=new ICoCo::MEDField(parafield->getField());
+ icocofield=new ICoCo::MEDDoubleField(parafield->getField());
dec.setMethod(targetMeth);
dec.attachLocalField(icocofield);
}
mesh->setCoords(myCoords);
myCoords->decrRef();
paramesh=new ParaMESH(mesh,*target_group,targetMeshName);
- DataArrayIdType *da=DataArrayIdType::New();
- const mcIdType globalNumberingP2[5]={0,1,2,3,4};
- da->useArray(globalNumberingP2,false,DeallocType::CPP_DEALLOC,5,1);
+ std::vector<mcIdType> globalNumberingP2 = {0,1,2,3,4};
+ MCAuto<DataArrayIdType> da=DataArrayIdType::New(); da->alloc(5,1);
+ std::copy(globalNumberingP2.begin(), globalNumberingP2.end(), da->rwBegin());
paramesh->setNodeGlobal(da);
- da->decrRef();
}
if(rank==3)
{
mesh->setCoords(myCoords);
myCoords->decrRef();
paramesh=new ParaMESH(mesh,*target_group,targetMeshName);
- DataArrayIdType *da=DataArrayIdType::New();
- const mcIdType globalNumberingP3[3]={4,2,5};
- da->useArray(globalNumberingP3,false,DeallocType::CPP_DEALLOC,3,1);
+ std::vector<mcIdType> globalNumberingP3 = {4,2,5};
+ MCAuto<DataArrayIdType> da=DataArrayIdType::New(); da->alloc(3,1);
+ std::copy(globalNumberingP3.begin(), globalNumberingP3.end(), da->rwBegin());
paramesh->setNodeGlobal(da);
- da->decrRef();
}
if(rank==4)
{
mesh->setCoords(myCoords);
myCoords->decrRef();
paramesh=new ParaMESH(mesh,*target_group,targetMeshName);
- DataArrayIdType *da=DataArrayIdType::New();
- const mcIdType globalNumberingP4[6]={3,6,7,4,8,5};
- da->useArray(globalNumberingP4,false,DeallocType::CPP_DEALLOC,6,1);
+ std::vector<mcIdType> globalNumberingP4 = {3,6,7,4,8,5};
+ MCAuto<DataArrayIdType> da=DataArrayIdType::New(); da->alloc(6,1);
+ std::copy(globalNumberingP4.begin(), globalNumberingP4.end(), da->rwBegin());
paramesh->setNodeGlobal(da);
- da->decrRef();
}
MEDCoupling::ComponentTopology comptopo;
parafieldP0 = new ParaFIELD(ON_CELLS,NO_TIME,paramesh, comptopo);
MEDCoupling::InterpKernelDEC dec (*source_group,*target_group);
- MEDCoupling::MEDCouplingUMesh* mesh;
- MEDCoupling::ParaMESH* paramesh;
- MEDCoupling::ParaFIELD* parafield;
-
- ICoCo::MEDField* icocofield ;
+ MEDCoupling::MEDCouplingUMesh* mesh = nullptr;
+ MEDCoupling::ParaMESH* paramesh = nullptr;
+ MEDCoupling::ParaFIELD* parafield = nullptr;
+ ICoCo::MEDDoubleField* icocofield = nullptr;
char * tmp_dir_c = getenv("TMP");
string tmp_dir;
value[ielem]=0.0;
// ICoCo::Field* icocofield=new ICoCo::MEDField(paramesh,parafield);
- icocofield=new ICoCo::MEDField(parafield->getField());
+ icocofield=new ICoCo::MEDDoubleField(parafield->getField());
dec.attachLocalField(icocofield);
for(int ielem=0; ielem<nb_local;ielem++)
value[ielem]=0.0;
// ICoCo::Field* icocofield=new ICoCo::MEDField(paramesh,parafield);
- icocofield=new ICoCo::MEDField(parafield->getField());
+ icocofield=new ICoCo::MEDDoubleField(parafield->getField());
dec.attachLocalField(icocofield);
}