From d60fe0cde471ced5d2a9b7d95bfb3812e521bf8c Mon Sep 17 00:00:00 2001
From: ageay <ageay>
Date: Tue, 25 Aug 2009 07:00:26 +0000
Subject: [PATCH] Add one test.

---
 src/ParaMEDMEM/Test/Makefile.am               |   1 +
 src/ParaMEDMEM/Test/ParaMEDMEMTest.hxx        |   2 +
 .../Test/ParaMEDMEMTest_Gauthier1.cxx         | 298 ++++++++++++++++++
 3 files changed, 301 insertions(+)
 create mode 100644 src/ParaMEDMEM/Test/ParaMEDMEMTest_Gauthier1.cxx

diff --git a/src/ParaMEDMEM/Test/Makefile.am b/src/ParaMEDMEM/Test/Makefile.am
index 71a2e0132..ba5d75adc 100644
--- a/src/ParaMEDMEM/Test/Makefile.am
+++ b/src/ParaMEDMEM/Test/Makefile.am
@@ -37,6 +37,7 @@ dist_libParaMEDMEMTest_la_SOURCES = \
 	ParaMEDMEMTest_StructuredCoincidentDEC.cxx \
 	ParaMEDMEMTest_MEDLoader.cxx \
 	ParaMEDMEMTest_ICocoTrio.cxx \
+	ParaMEDMEMTest_Gauthier1.cxx \
 	MPIAccessDECTest.cxx \
 	test_AllToAllDEC.cxx \
 	test_AllToAllvDEC.cxx \
diff --git a/src/ParaMEDMEM/Test/ParaMEDMEMTest.hxx b/src/ParaMEDMEM/Test/ParaMEDMEMTest.hxx
index 4d94f7f86..9e33a62b6 100644
--- a/src/ParaMEDMEM/Test/ParaMEDMEMTest.hxx
+++ b/src/ParaMEDMEM/Test/ParaMEDMEMTest.hxx
@@ -61,6 +61,7 @@ class ParaMEDMEMTest : public CppUnit::TestFixture
   CPPUNIT_TEST(testStructuredCoincidentDEC);
   CPPUNIT_TEST(testStructuredCoincidentDEC);
   CPPUNIT_TEST(testICocoTrio1);
+  CPPUNIT_TEST(testGauthier1);
   CPPUNIT_TEST(testMEDLoaderRead1);
   CPPUNIT_TEST(testMEDLoaderPolygonRead);
   CPPUNIT_TEST(testMEDLoaderPolyhedronRead);
@@ -105,6 +106,7 @@ public:
   void testAsynchronousFastSourceIntersectionDEC_2D();
   //
   void testICocoTrio1();
+  void testGauthier1();
   //
   void testMEDLoaderRead1();
   void testMEDLoaderPolygonRead();
diff --git a/src/ParaMEDMEM/Test/ParaMEDMEMTest_Gauthier1.cxx b/src/ParaMEDMEM/Test/ParaMEDMEMTest_Gauthier1.cxx
new file mode 100644
index 000000000..b104583af
--- /dev/null
+++ b/src/ParaMEDMEM/Test/ParaMEDMEMTest_Gauthier1.cxx
@@ -0,0 +1,298 @@
+#include "ParaMEDMEMTest.hxx"
+#include <cppunit/TestAssert.h>
+
+#include <string>
+#include "CommInterface.hxx"
+#include "ProcessorGroup.hxx"
+#include "MPIProcessorGroup.hxx"
+#include "DEC.hxx"
+#include "IntersectionDEC.hxx"
+#include <set>
+#include <time.h>
+#include "ICoCoTrioField.hxx"
+#include <iostream>
+#include <assert.h>
+#include <math.h>
+
+using namespace std;
+using namespace ParaMEDMEM;
+using namespace ICoCo;
+
+void afficheGauthier1( const TrioField&   field, const double *vals, int lgth)
+{
+  CPPUNIT_ASSERT_EQUAL(lgth,field._nb_elems);
+  for (int ele=0;ele<field._nb_elems;ele++)
+    CPPUNIT_ASSERT_DOUBLES_EQUAL(vals[ele],field._field[ele],1e-12);
+}
+
+void remplit_coordGauthier1(double* coords)
+{
+  double angle,epaisseur;
+  angle=0*45*(asin(1)/90);
+  epaisseur=1e-0;
+  coords[0*3+0]=0.;
+  coords[0*3+1]=0.;
+  coords[0*3+2]=0.;
+  
+  coords[1*3+0]=cos(angle);
+  coords[1*3+1]=0.;
+  coords[1*3+2]=sin(angle);
+  
+    
+  coords[2*3+0]=-sin(angle);
+  coords[2*3+1]=0.;
+  coords[2*3+2]=cos(angle);
+  
+  for (int d=0;d<3;d++)
+    coords[3*3+d]=coords[1*3+d]+ coords[2*3+d];
+  
+  for (int i=4;i<8;i++)
+    {
+      for (int d=0;d<3;d++)
+        coords[i*3+d]=coords[(i-4)*3+d];
+      coords[i*3+1]+=epaisseur;
+    }
+
+}
+
+void init_quadGauthier1(TrioField& champ_quad,int is_master)
+{
+  
+  champ_quad.setName("champ_quad");
+  champ_quad._space_dim=3;
+  champ_quad._mesh_dim=2;
+  champ_quad._nodes_per_elem=4;
+  champ_quad._itnumber=0;
+  champ_quad._time1=0;
+  champ_quad._time2=1;
+  champ_quad._nb_field_components=1;
+
+  if (is_master)
+    {
+      champ_quad._nbnodes=8;
+      champ_quad._nb_elems=2;
+      
+      champ_quad._coords=new double[champ_quad._nbnodes*champ_quad._space_dim];
+      //memcpy(afield._coords,sommets.addr(),champ_quad._nbnodes*champ_quad._space_dim*sizeof(double));
+      
+      remplit_coordGauthier1(champ_quad._coords);
+  
+  
+      champ_quad._connectivity=new int[champ_quad._nb_elems*champ_quad._nodes_per_elem];
+      champ_quad._connectivity[0*champ_quad._nodes_per_elem+0]=0;
+      champ_quad._connectivity[0*champ_quad._nodes_per_elem+1]=1;
+      champ_quad._connectivity[0*champ_quad._nodes_per_elem+2]=3;
+      champ_quad._connectivity[0*champ_quad._nodes_per_elem+3]=2;
+      champ_quad._connectivity[1*champ_quad._nodes_per_elem+0]=4;
+      champ_quad._connectivity[1*champ_quad._nodes_per_elem+1]=5;
+      champ_quad._connectivity[1*champ_quad._nodes_per_elem+2]=7;
+      champ_quad._connectivity[1*champ_quad._nodes_per_elem+3]=6;
+      
+    }
+  else
+    {
+      champ_quad._nbnodes=0;
+      champ_quad._nb_elems=0;
+      champ_quad._coords=new double[champ_quad._nbnodes*champ_quad._space_dim];
+    
+    }
+  champ_quad._has_field_ownership=false;
+  champ_quad._field=0;
+  //champ_quad._field=new double[champ_quad._nb_elems];
+  //  assert(champ_quad._nb_field_components==1);
+}
+void init_triangleGauthier1(TrioField& champ_triangle,int is_master)
+{
+   
+  champ_triangle.setName("champ_triangle");
+  champ_triangle._space_dim=3;
+  champ_triangle._mesh_dim=2;
+  champ_triangle._nodes_per_elem=3;
+  champ_triangle._itnumber=0;
+  champ_triangle._time1=0;
+  champ_triangle._time2=1;
+  champ_triangle._nb_field_components=1;
+
+  if (is_master)
+    {
+      champ_triangle._nb_elems=4;
+      champ_triangle._nbnodes=8;
+    
+      champ_triangle._coords=new double[champ_triangle._nbnodes*champ_triangle._space_dim];
+      //memcpy(afield._coords,sommets.addr(),champ_triangle._nbnodes*champ_triangle._space_dim*sizeof(double));
+      remplit_coordGauthier1(champ_triangle._coords);
+      
+      champ_triangle._connectivity=new int[champ_triangle._nb_elems*champ_triangle._nodes_per_elem];
+      champ_triangle._connectivity[0*champ_triangle._nodes_per_elem+0]=0;
+      champ_triangle._connectivity[0*champ_triangle._nodes_per_elem+1]=1;
+      champ_triangle._connectivity[0*champ_triangle._nodes_per_elem+2]=2;
+      champ_triangle._connectivity[1*champ_triangle._nodes_per_elem+0]=1;
+      champ_triangle._connectivity[1*champ_triangle._nodes_per_elem+1]=2;
+      champ_triangle._connectivity[1*champ_triangle._nodes_per_elem+2]=3;
+      
+      champ_triangle._connectivity[2*champ_triangle._nodes_per_elem+0]=4;
+      champ_triangle._connectivity[2*champ_triangle._nodes_per_elem+1]=5;
+      champ_triangle._connectivity[2*champ_triangle._nodes_per_elem+2]=7;
+      champ_triangle._connectivity[3*champ_triangle._nodes_per_elem+0]=4;
+      champ_triangle._connectivity[3*champ_triangle._nodes_per_elem+1]=6;
+      champ_triangle._connectivity[3*champ_triangle._nodes_per_elem+2]=7;
+    }
+  else
+    {
+      champ_triangle._nb_elems=0;
+      champ_triangle._nbnodes=0;
+      champ_triangle._coords=new double[champ_triangle._nbnodes*champ_triangle._space_dim];
+    
+    }
+  champ_triangle._has_field_ownership=false;
+  // champ_triangle._field=new double[champ_triangle._nb_elems];
+  champ_triangle._field=0;
+  
+}
+
+
+void ParaMEDMEMTest::testGauthier1()
+{
+  int num_cas=0;
+  int rank, size;
+  MPI_Comm_rank(MPI_COMM_WORLD,&rank);
+  MPI_Comm_size(MPI_COMM_WORLD,&size);
+  
+  int is_master=0;
+
+  CommInterface comm;
+  set<int> emetteur_ids;
+  set<int> recepteur_ids;
+  emetteur_ids.insert(0);
+  if(size!=4)
+    return;
+  recepteur_ids.insert(1);
+  if (size >2) 
+    recepteur_ids.insert(2);
+  if (size >2) 
+    emetteur_ids.insert(3);
+  if ((rank==0)||(rank==1)) 
+    is_master=1;
+  
+  MPIProcessorGroup recepteur_group(comm,recepteur_ids);
+  MPIProcessorGroup emetteur_group(comm,emetteur_ids);
+
+
+  string cas;
+  if (recepteur_group.containsMyRank())
+    {
+      cas="recepteur";
+      //freopen("recpeteur.out","w",stdout);
+      //freopen("recepteur.err","w",stderr);
+      
+    }
+  else
+    {
+      cas="emetteur";
+      // freopen("emetteur.out","w",stdout);
+      //freopen("emetteur.err","w",stderr);
+    }
+  double expected[8][4]={
+    {1.,1.,1.,1.},
+    {40., 40., 1., 1.},
+    {1.,1.,1e200,1e200},
+    {40.,1.,1e200,1e200},
+    {1.,1.,1.,1.},
+    {40.,1.,1.,1.},
+    {1.,1.,1e200,1e200},
+    {20.5,1.,1e200,1e200}
+  };
+
+  int expectedLgth[8]={4,4,2,2,4,4,2,2};
+  
+  for (int send=0;send<2;send++)
+    for (int rec=0;rec<2;rec++)
+      {
+        IntersectionDEC dec_emetteur(emetteur_group, recepteur_group);
+        dec_emetteur.setOrientation(2);
+        TrioField champ_emetteur, champ_recepteur;
+   
+        if (send==0)
+          init_quadGauthier1(champ_emetteur,is_master);
+        else
+          init_triangleGauthier1(champ_emetteur,is_master);
+        if (rec==0)
+          init_triangleGauthier1(champ_recepteur,is_master);
+        else
+          init_quadGauthier1(champ_recepteur,is_master);
+  
+        if (cas=="emetteur") 
+          {
+            champ_emetteur._field=new double[champ_emetteur._nb_elems];
+            for (int ele=0;ele<champ_emetteur._nb_elems;ele++)
+              champ_emetteur._field[ele]=1;
+      
+            champ_emetteur._has_field_ownership=true;
+          }
+  
+  
+        MPI_Barrier(MPI_COMM_WORLD);
+
+        clock_t clock0= clock ();
+        int compti=0;
+
+        bool init=true; // first time step ??
+        bool stop=false;
+        //boucle sur les pas de quads
+        while (!stop) {
+  
+          compti++;
+          clock_t clocki= clock ();
+          //cout << compti << " CLOCK " << (clocki-clock0)*1.e-6 << endl; 
+          for (int non_unif=0;non_unif<2;non_unif++)
+            {
+              // if (champ_recepteur._field)
+              //   delete [] champ_recepteur._field;
+              champ_recepteur._field=0;
+              // champ_recepteur._has_field_ownership=false;
+  
+
+  
+              if (cas=="emetteur") 
+                {
+                  if (non_unif)
+                    if(rank!=3)
+                      champ_emetteur._field[0]=40;
+                }
+              bool ok=false; // Is the time interval successfully solved ?
+    
+              // Loop on the time interval tries
+              if(1) {
+      
+
+                if (cas=="emetteur")
+                  dec_emetteur.attachLocalField((ICoCo::Field*) &champ_emetteur);
+                else
+                  dec_emetteur.attachLocalField((ICoCo::Field*) &champ_recepteur);
+
+
+                if(init) dec_emetteur.synchronize();
+                init=false;
+
+                if (cas=="emetteur") {
+                  //    affiche(champ_emetteur);
+                  dec_emetteur.sendData();
+                }
+                else if (cas=="recepteur")
+                  {
+                    dec_emetteur.recvData();
+                    if (is_master)
+                      afficheGauthier1(champ_recepteur,expected[num_cas],expectedLgth[num_cas]);
+                  }
+                else
+                  throw 0;
+                MPI_Barrier(MPI_COMM_WORLD);
+              }
+              stop=true;
+              num_cas++;
+            }
+          // destruction des champs, des DEC, et des tableaux associés
+        }
+      }
+}
+
-- 
2.39.2