MEDMEM::SUPPORT* support;
ParaMEDMEM::ParaMESH* paramesh;
ParaMEDMEM::ParaFIELD* parafield;
+ ParaMEDMEM::ParaSUPPORT* parasupport;
+ double * value;
+ ICoCo::Field* icocofield ;
string data_dir = getenv("MED_ROOT_DIR");
string tmp_dir = getenv("TMP");
paramesh=new ParaMESH (*mesh,*source_group,"source mesh");
- ParaMEDMEM::ParaSUPPORT* parasupport=new UnstructuredParaSUPPORT( support,*source_group);
+// ParaMEDMEM::ParaSUPPORT* parasupport=new UnstructuredParaSUPPORT( support,*source_group);
+ parasupport=new UnstructuredParaSUPPORT( support,*source_group);
ParaMEDMEM::ComponentTopology comptopo;
parafield = new ParaFIELD(parasupport, comptopo);
int nb_local=support->getNumberOfElements(MED_EN::MED_ALL_ELEMENTS);
- double * value= new double[nb_local];
+// double * value= new double[nb_local];
+ value= new double[nb_local];
for(int ielem=0; ielem<nb_local;ielem++)
value[ielem]=1.0;
parafield->getField()->setValue(value);
- ICoCo::Field* icocofield=new ICoCo::MEDField(paramesh,parafield);
+// ICoCo::Field* icocofield=new ICoCo::MEDField(paramesh,parafield);
+ icocofield=new ICoCo::MEDField(paramesh,parafield);
dec.attachLocalField(icocofield);
}
support=new MEDMEM::SUPPORT(mesh,"all elements",MED_EN::MED_CELL);
paramesh=new ParaMESH (*mesh,*target_group,"target mesh");
- ParaMEDMEM::ParaSUPPORT* parasupport=new UnstructuredParaSUPPORT(support,*target_group);
+// ParaMEDMEM::ParaSUPPORT* parasupport=new UnstructuredParaSUPPORT(support,*target_group);
+ parasupport=new UnstructuredParaSUPPORT(support,*target_group);
ParaMEDMEM::ComponentTopology comptopo;
parafield = new ParaFIELD(parasupport, comptopo);
int nb_local=support->getNumberOfElements(MED_EN::MED_ALL_ELEMENTS);
- double * value= new double[nb_local];
+// double * value= new double[nb_local];
+ value= new double[nb_local];
for(int ielem=0; ielem<nb_local;ielem++)
value[ielem]=0.0;
parafield->getField()->setValue(value);
- ICoCo::Field* icocofield=new ICoCo::MEDField(paramesh,parafield);
+// ICoCo::Field* icocofield=new ICoCo::MEDField(paramesh,parafield);
+ icocofield=new ICoCo::MEDField(paramesh,parafield);
dec.attachLocalField(icocofield);
}
delete source_group;
delete target_group;
delete self_group;
-
+ delete mesh;
+ delete support;
+ delete paramesh;
+ delete parafield;
+ delete parasupport;
+ delete [] value;
+ delete icocofield;
MPI_Barrier(MPI_COMM_WORLD);
cout << "end of IntersectionDEC_2D test"<<endl;