#include <string>
#include <MEDCouplingFieldDouble.hxx>
+#include <MEDCouplingCMesh.hxx>
+#include <MEDFileField1TS.hxx>
+#include "MEDFileMesh.hxx"
using namespace std;
using namespace MEDCoupling;
CPPUNIT_ASSERT_EQUAL( 1, concF1.getNumberOfComponents() );
CPPUNIT_ASSERT_EQUAL( 0.0, concF1.getTime() );
CPPUNIT_ASSERT_EQUAL( 24, concF1.getNumberOfElements() );
+
+ Mesh M2(0.,1.,2,0.,1.,2,1);
+ Field concF2("CONCENTRATION",CELLS,M2) ;
+ for (int j=0;j<concF2.getNumberOfComponents();j++)
+ for (int i=0;i<concF2.getNumberOfElements();i++)
+ concF2(i,j)=i+j;
+
+ for (int j=0;j<concF2.getNumberOfComponents();j++)
+ for (int i=0;i<concF2.getNumberOfElements();i++)
+ CPPUNIT_ASSERT_EQUAL( double(i+j), concF2(i,j) );
+ CPPUNIT_ASSERT_EQUAL( 1, concF2.getNumberOfComponents() );
+ CPPUNIT_ASSERT_EQUAL( 0.0, concF2.getTime() );
+ CPPUNIT_ASSERT_EQUAL( 8, concF2.getNumberOfElements() );
+ CPPUNIT_ASSERT(concF2.meshNotDeleted());
+ concF2.writeMED("FieldConcF2");//This saves the mesh and the values of iteration 0 at time t=0
+// concF2.deleteMEDCouplingUMesh();//medcouplingmesh is no longer needed as the mesh was already saved in the previous line
+ concF2.setTime(0.5,1);//Increase the time to 0.5 and the iteration to 1
+ for (int j=0;j<concF2.getNumberOfComponents();j++)
+ for (int i=0;i<concF2.getNumberOfElements();i++)
+ concF2(i,j)=i*j;
+ concF2.writeMED("FieldConcF2", true);//This saves only the values of iteration 1 at time t=0.5. The previous values are not deleted
+
+ Mesh M3(0.0,1.0,2,0.,1.,2,0.,1.,2);
+ Field concF3("CONCENTRATION",FACES,M3) ;
+ for (int j=0;j<concF3.getNumberOfComponents();j++)
+ for (int i=0;i<concF3.getNumberOfElements();i++)
+ concF3(i,j)=i+j;
+
+ for (int j=0;j<concF3.getNumberOfComponents();j++)
+ for (int i=0;i<concF3.getNumberOfElements();i++)
+ CPPUNIT_ASSERT_EQUAL( double(i+j), concF3(i,j) );
+ CPPUNIT_ASSERT_EQUAL( 1, concF3.getNumberOfComponents() );
+ CPPUNIT_ASSERT_EQUAL( 0.0, concF3.getTime() );
+ CPPUNIT_ASSERT_EQUAL( 36, concF3.getNumberOfElements() );
+ CPPUNIT_ASSERT(concF3.meshNotDeleted());
+
+ //Load the Field CONCENTRATION in the file fileNameMED
+ Field concF4(fileNameMED,CELLS,"CONCENTRATION",0,0);
+ CPPUNIT_ASSERT_EQUAL( 2, concF4.getNumberOfComponents() );
+ CPPUNIT_ASSERT_EQUAL( 1.2, concF4.getTime() );
+ CPPUNIT_ASSERT_EQUAL( 50, concF4.getNumberOfElements() );
+ for (int j=0;j<concF4.getNumberOfComponents();j++)
+ for (int i=0;i<concF4.getNumberOfElements();i++)
+ CPPUNIT_ASSERT_EQUAL( double(i+j), concF4(i,j) );
+ CPPUNIT_ASSERT(concF4.meshNotDeleted());
+ concF4.writeMED("FieldConcF4");//This saves the mesh and the values of iteration 0 at time t=0
+// concF4.deleteMEDCouplingUMesh();//medcouplingmesh is no longer needed as the mesh was already saved in the previous line
+ concF4.setTime(0.5,1);//Increase the time to 0.5 and the iteration to 1
+ for (int j=0;j<concF4.getNumberOfComponents();j++)
+ for (int i=0;i<concF4.getNumberOfElements();i++)
+ concF4(i,j)=i*j;
+ concF4.writeMED("FieldConcF4", false);//This saves only the values of iteration 1 at time t=0.5. The previous values are not deleted
+
+ //Create a constant field on the mesh fileNameMEDn
+ Field concF5(fileNameMEDn,NODES,std::vector<double> (3,1),"CONSTANT_Field");
+ CPPUNIT_ASSERT_EQUAL( 3, concF5.getNumberOfComponents() );
+ CPPUNIT_ASSERT_EQUAL( 0., concF5.getTime() );
+ CPPUNIT_ASSERT_EQUAL( 66, concF5.getNumberOfElements() );
+ for (int j=0;j<concF5.getNumberOfComponents();j++)
+ for (int i=0;i<concF5.getNumberOfElements();i++)
+ CPPUNIT_ASSERT_EQUAL( 1., concF5(i,j) );
+ CPPUNIT_ASSERT(concF5.meshNotDeleted());
+ (concF5.getMesh()).writeMED("FieldConcF5");//This saves only the mesh
+ cout<<"Mesh name : " << concF5.getMesh().getName()<<endl;
+ cout<<"Field name : " << concF5.getName()<<endl;
+ concF5.setTime(0.5,1);//Increase the time to 0.5 and the iteration to 1
+ concF5.writeMED("FieldConcF5",false);//This saves the mesh and the values of iteration 0 at time t=0
+ //concF5.deleteMEDCouplingUMesh();//medcouplingmesh is no longer needed as the mesh was already saved in the previous line
+ for (int j=0;j<concF5.getNumberOfComponents();j++)
+ for (int i=0;i<concF5.getNumberOfElements();i++)
+ concF5(i,j)=i*j;
+ //concF5.writeMED("FieldConcF5", false);//This saves only the values of iteration 1 at time t=0.5. The previous values are not deleted
+
+
+ /* 2D image mesh */
+ //int _spaceDim=2;
+ //double *originPtr = new double[_spaceDim];
+ //double *dxyzPtr = new double[_spaceDim];
+ //mcIdType *nodeStrctPtr = new mcIdType[_spaceDim];
+
+ //originPtr[0]=0;
+ //originPtr[1]=0;
+ //nodeStrctPtr[0]=3;
+ //nodeStrctPtr[1]=3;
+ //dxyzPtr[0]=1;
+ //dxyzPtr[1]=1;
+
+ //MEDCouplingIMesh * _mesh=MEDCouplingIMesh::New("test",
+ //_spaceDim,
+ //nodeStrctPtr,
+ //nodeStrctPtr+_spaceDim,
+ //originPtr,
+ //originPtr+_spaceDim,
+ //dxyzPtr,
+ //dxyzPtr+_spaceDim);
+ //MEDCouplingUMesh * m1 = _mesh->buildUnstructured();
+ //m1->setName("mesh");
+
+ //MEDCouplingFieldDouble * f = MEDCouplingFieldDouble::New(ON_CELLS, ONE_TIME);
+ //f->setMesh(m1);
+ //f->setName("F");
+ //*f=0;
+ //f->setTime(0.0,0,0);
+
+ //MEDFileField1TS * ff;
+ //ff->setFieldNoProfileSBT(f);
+
+ /* 1D image mesh */
+ //int _spaceDim=1;
+ //double *originPtr = new double[_spaceDim];
+ //double *dxyzPtr = new double[_spaceDim];
+ //mcIdType *nodeStrctPtr = new mcIdType[_spaceDim];
+
+ //originPtr[0]=0;
+ //nodeStrctPtr[0]=3;
+ //dxyzPtr[0]=1;
+
+ //MEDCouplingIMesh * _mesh=MEDCouplingIMesh::New("test",
+ //_spaceDim,
+ //nodeStrctPtr,
+ //nodeStrctPtr+_spaceDim,
+ //originPtr,
+ //originPtr+_spaceDim,
+ //dxyzPtr,
+ //dxyzPtr+_spaceDim);
+ //MEDCouplingUMesh * m1 = _mesh->buildUnstructured();
+ //m1->setName("mesh");
+
+ //MEDCouplingFieldDouble * f = MEDCouplingFieldDouble::New(ON_CELLS, ONE_TIME);
+ //f->setMesh(m1);
+ //f->setName("F");
+ //*f=0;
+ //f->setTime(0.0,0,0);
+
+ //MEDFileField1TS * ff;
+ //ff->setFieldNoProfileSBT(f);
+
+ /* 2D cartesian mesh */
+ //Dataarray
+ double XCoords[3]={0.,1.,2.};
+ double YCoords[3]={0.,1.,2.};
+ MEDCoupling::DataArrayDouble *arrX=MEDCoupling::DataArrayDouble::New();
+ arrX->alloc(3,1);
+ std::copy(XCoords,XCoords+3,arrX->getPointer());
+ arrX->setInfoOnComponent(0,"X [m]");
+ MEDCoupling::DataArrayDouble *arrY=MEDCoupling::DataArrayDouble::New();
+ arrY->alloc(3,1);
+ std::copy(YCoords,YCoords+3,arrY->getPointer());
+ arrY->setInfoOnComponent(0,"Y [m]");
+ //Mesh
+ MEDCoupling::MEDCouplingCMesh *mesh=MEDCoupling::MEDCouplingCMesh::New("My2D_CMesh");
+ mesh->setCoords(arrX,arrY);
+ arrX->decrRef();
+ arrY->decrRef();
+ MEDCouplingUMesh * m1 = mesh->buildUnstructured();
+ m1->setName("mesh");
+ //Field
+ MEDCouplingFieldDouble * f = MEDCouplingFieldDouble::New(ON_CELLS, ONE_TIME);
+ f->setMesh(m1);
+ f->setName("F");
+ *f=0;
+ f->setTime(0.0,0,0);
+ //MEDFileField1TS
+ MEDFileField1TS * ff;
+ //ff->setFieldNoProfileSBT(f);
}
