From 6b5612304d161e1de3d76f71fd6b5d503d511f9b Mon Sep 17 00:00:00 2001
From: secher <secher>
Date: Wed, 18 Mar 2009 15:42:09 +0000
Subject: [PATCH] add test for perf measures

---
 src/ParaMEDMEM/Test/Makefile.am   |   4 +-
 src/ParaMEDMEM/Test/test_perf.cxx | 316 ++++++++++++++++++++++++++++++
 2 files changed, 319 insertions(+), 1 deletion(-)
 create mode 100644 src/ParaMEDMEM/Test/test_perf.cxx

diff --git a/src/ParaMEDMEM/Test/Makefile.am b/src/ParaMEDMEM/Test/Makefile.am
index c10555ebe..b00a9498c 100644
--- a/src/ParaMEDMEM/Test/Makefile.am
+++ b/src/ParaMEDMEM/Test/Makefile.am
@@ -72,11 +72,13 @@ libParaMEDMEMTest_la_LDFLAGS= ../MEDLoader/libparamedmemmedloader.la @CPPUNIT_LI
 # Executables targets
 bin_PROGRAMS = TestParaMEDMEM  \
 							 TestMPIAccessDEC \
-		  			   TestMPIAccess
+		  			   TestMPIAccess \
+		test_perf
 
 dist_TestParaMEDMEM_SOURCES = TestParaMEDMEM.cxx
 dist_TestMPIAccessDEC_SOURCES = TestMPIAccessDEC.cxx
 dist_TestMPIAccess_SOURCES = TestMPIAccess.cxx
+dist_test_perf_SOURCES = test_perf.cxx
 
 LDADD= $(MED2_LIBS) $(libMEDMEMTest_la_LDFLAGS) -lm $(MPI_LIBS) \
 	../libparamedmem.la libParaMEDMEMTest.la \
diff --git a/src/ParaMEDMEM/Test/test_perf.cxx b/src/ParaMEDMEM/Test/test_perf.cxx
new file mode 100644
index 000000000..6720d1173
--- /dev/null
+++ b/src/ParaMEDMEM/Test/test_perf.cxx
@@ -0,0 +1,316 @@
+#include <time.h>
+#include <sys/times.h>
+#include <sys/time.h>
+#include "ParaMEDMEMTest.hxx"
+#include <cppunit/TestAssert.h>
+
+#include "CommInterface.hxx"
+#include "ProcessorGroup.hxx"
+#include "MPIProcessorGroup.hxx"
+#include "Topology.hxx"
+#include "DEC.hxx"
+#include "MxN_Mapping.hxx"
+#include "IntersectionDEC.hxx"
+#include "ParaMESH.hxx"
+#include "ParaFIELD.hxx"
+#include "ICoCoMEDField.hxx"
+#include "MEDLoader.hxx"
+ 
+#include <string>
+
+// use this define to enable lines, execution of which leads to Segmentation Fault
+#define ENABLE_FAULTS
+
+// use this define to enable CPPUNIT asserts and fails, showing bugs
+#define ENABLE_FORCED_FAILURES
+
+#ifndef CLK_TCK 
+#include <unistd.h>
+#define CLK_TCK sysconf(_SC_CLK_TCK);
+#endif 
+
+using namespace std;
+using namespace ParaMEDMEM;
+ 
+void testIntersectionDEC_2D(const string& filename1, const string& meshname1,
+			    const string& filename2, const string& meshname2,
+			    int nproc_source, double epsilon, bool tri, bool all);
+void get_time( float *telps, float *tuser, float *tsys, float *tcpu );
+
+int main(int argc, char *argv[])
+{
+  string filename1, filename2;
+  string meshname1, meshname2;
+  int nproc_source=1, rank;
+  double epsilon=1.e-6;
+  int count=0;
+  bool tri=false;
+  bool all=false;
+
+  MPI_Init(&argc,&argv);
+
+  for(int i=1;i<argc;i++){
+    if( strcmp(argv[i],"-f1") == 0 ){
+      filename1 = argv[++i];
+      count++;
+    }
+    else if( strcmp(argv[i],"-f2") == 0 ){
+      filename2 = argv[++i];
+      count++;
+    }
+    else if( strcmp(argv[i],"-m1") == 0 ){
+      meshname1 = argv[++i];
+      count++;
+    }
+    else if( strcmp(argv[i],"-m2") == 0 ){
+      meshname2 = argv[++i];
+      count++;
+    }
+    else if( strcmp(argv[i],"-ns") == 0 ){
+      nproc_source = atoi(argv[++i]);
+    }
+    else if( strcmp(argv[i],"-eps") == 0 ){
+      epsilon = atof(argv[++i]);
+    }
+    else if( strcmp(argv[i],"-tri") == 0 ){
+      tri = true;
+    }
+    else if( strcmp(argv[i],"-all") == 0 ){
+      all = true;
+    }
+  }
+
+  if( count != 4 ){
+    cout << "usage test_perf -f1 filename1 -m1 meshname1 -f2 filename2 -m2 meshname2 (-ns nproc_source -eps epsilon -tri -all)" << endl;
+    exit(0);
+  }
+
+  MPI_Comm_rank(MPI_COMM_WORLD,&rank);
+  testIntersectionDEC_2D(filename1,meshname1,filename2,meshname2,nproc_source,epsilon,tri,all);
+
+  MPI_Finalize();
+}
+
+void testIntersectionDEC_2D(const string& filename_xml1, const string& meshname1,
+			    const string& filename_xml2, const string& meshname2,
+			    int nproc_source, double epsilon, bool tri, bool all)
+{
+  float tcpu, tcpu_u, tcpu_s, telps;
+  int size;
+  int rank;
+  MPI_Comm_size(MPI_COMM_WORLD,&size);
+  MPI_Comm_rank(MPI_COMM_WORLD,&rank);
+ 
+  set<int> self_procs;
+  set<int> procs_source;
+  set<int> procs_target;
+  
+  for (int i=0; i<nproc_source; i++)
+    procs_source.insert(i);
+  for (int i=nproc_source; i<size; i++)
+    procs_target.insert(i);
+  self_procs.insert(rank);
+  
+  ParaMEDMEM::CommInterface interface;
+    
+  ParaMEDMEM::ProcessorGroup* self_group = new ParaMEDMEM::MPIProcessorGroup(interface,self_procs);
+  ParaMEDMEM::ProcessorGroup* target_group = new ParaMEDMEM::MPIProcessorGroup(interface,procs_target);
+  ParaMEDMEM::ProcessorGroup* source_group = new ParaMEDMEM::MPIProcessorGroup(interface,procs_source);
+  
+  //loading the geometry for the source group
+
+  ParaMEDMEM::IntersectionDEC dec (*source_group,*target_group);
+  if(tri)
+    dec.setIntersectionType(INTERP_KERNEL::Triangulation);
+  else
+    dec.setIntersectionType(INTERP_KERNEL::Convex);
+
+  ParaMEDMEM::MEDCouplingUMesh* mesh;
+  ParaMEDMEM::ParaMESH* paramesh;
+  ParaMEDMEM::ParaFIELD* parafield;
+  ICoCo::Field* icocofield ;
+  
+  // To remove tmp files from disk
+  ParaMEDMEMTest_TmpFilesRemover aRemover;
+  
+  MPI_Barrier(MPI_COMM_WORLD);
+  if (source_group->containsMyRank()){
+    string master = filename_xml1;
+      
+    ostringstream strstream;
+    if( nproc_source == 1 )
+      strstream <<master<<".med";
+    else
+      strstream <<master<<rank+1<<".med";
+
+    ostringstream meshname ;
+    if( nproc_source == 1 )
+      meshname<< meshname1;
+    else
+      meshname<< meshname1<<"_"<< rank+1;
+      
+    get_time( &telps, &tcpu_u, &tcpu_s, &tcpu );
+    mesh=MEDLoader::ReadUMeshFromFile(strstream.str().c_str(),meshname.str().c_str(),0);
+    get_time( &telps, &tcpu_u, &tcpu_s, &tcpu );
+    if( rank == 0 )
+      cout << "IO : Telapse = " << telps << " TuserCPU = " << tcpu_u << " TsysCPU = " << tcpu_s << " TCPU = " << tcpu << endl;
+    mesh->incrRef();
+    
+    paramesh=new ParaMESH (mesh,*source_group,"source mesh");
+    
+    ParaMEDMEM::ComponentTopology comptopo;
+    parafield = new ParaFIELD(ON_CELLS, paramesh, comptopo);
+
+    int nb_local=mesh->getNumberOfCells();
+    double *value=parafield->getField()->getArray()->getPointer();
+    for(int ielem=0; ielem<nb_local;ielem++)
+      value[ielem]=1.0;
+    
+    icocofield=new ICoCo::MEDField(paramesh,parafield);
+     
+    dec.attachLocalField(icocofield);
+  }
+  
+  //loading the geometry for the target group
+  if (target_group->containsMyRank()){
+    string master= filename_xml2;
+    ostringstream strstream;
+    if( (size-nproc_source) == 1 )
+      strstream << master<<".med";
+    else
+      strstream << master<<(rank-nproc_source+1)<<".med";
+    ostringstream meshname ;
+    if( (size-nproc_source) == 1 )
+      meshname<< meshname2;
+    else
+      meshname<< meshname2<<"_"<<rank-nproc_source+1;
+      
+    get_time( &telps, &tcpu_u, &tcpu_s, &tcpu );
+    mesh = MEDLoader::ReadUMeshFromFile(strstream.str().c_str(),meshname.str().c_str(),0);
+    get_time( &telps, &tcpu_u, &tcpu_s, &tcpu );
+    mesh->incrRef();
+
+    paramesh=new ParaMESH (mesh,*target_group,"target mesh");
+    ParaMEDMEM::ComponentTopology comptopo;
+    parafield = new ParaFIELD(ON_CELLS,paramesh, comptopo);
+		
+    int nb_local=mesh->getNumberOfCells();
+    double *value=parafield->getField()->getArray()->getPointer();
+    for(int ielem=0; ielem<nb_local;ielem++)
+      value[ielem]=0.0;
+    icocofield=new ICoCo::MEDField(paramesh,parafield);
+      
+    dec.attachLocalField(icocofield);
+  }
+    
+  
+  //attaching a DEC to the source group 
+  double field_before_int;
+  double field_after_int;
+  
+  if (source_group->containsMyRank()){ 
+    field_before_int = parafield->getVolumeIntegral(0);
+    get_time( &telps, &tcpu_u, &tcpu_s, &tcpu );
+    dec.synchronize();
+    get_time( &telps, &tcpu_u, &tcpu_s, &tcpu );
+    if( rank == 0 )
+      cout << "SYNCHRONIZE : Telapse = " << telps << " TuserCPU = " << tcpu_u << " TsysCPU = " << tcpu_s << " TCPU = " << tcpu << endl;
+    cout<<"DEC usage"<<endl;
+    dec.setForcedRenormalization(false);
+    if(all)
+      dec.setAllToAllMethod(PointToPoint);
+
+    get_time( &telps, &tcpu_u, &tcpu_s, &tcpu );
+    dec.sendData();
+    
+    get_time( &telps, &tcpu_u, &tcpu_s, &tcpu );
+    if( rank == 0 )
+      cout << "SEND DATA : Telapse = " << telps << " TuserCPU = " << tcpu_u << " TsysCPU = " << tcpu_s << " TCPU = " << tcpu << endl;
+    dec.recvData();
+     
+    field_after_int = parafield->getVolumeIntegral(0);
+//    CPPUNIT_ASSERT_DOUBLES_EQUAL(field_before_int, field_after_int, epsilon);
+      
+  }
+  
+  //attaching a DEC to the target group
+  if (target_group->containsMyRank()){
+    get_time( &telps, &tcpu_u, &tcpu_s, &tcpu );
+    dec.synchronize();
+    get_time( &telps, &tcpu_u, &tcpu_s, &tcpu );
+    dec.setForcedRenormalization(false);
+    if(all)
+      dec.setAllToAllMethod(PointToPoint);
+
+    get_time( &telps, &tcpu_u, &tcpu_s, &tcpu );
+    dec.recvData();
+    get_time( &telps, &tcpu_u, &tcpu_s, &tcpu );
+    dec.sendData();
+  }
+  
+  get_time( &telps, &tcpu_u, &tcpu_s, &tcpu );
+  if( rank == 0 )
+    cout << "RECV DATA : Telapse = " << telps << " TuserCPU = " << tcpu_u << " TsysCPU = " << tcpu_s << " TCPU = " << tcpu << endl;
+
+  delete source_group;
+  delete target_group;
+  delete self_group;
+  delete paramesh;
+  delete parafield;
+  mesh->decrRef() ;
+  delete icocofield;
+
+  MPI_Barrier(MPI_COMM_WORLD);
+  cout << "end of IntersectionDEC_2D test"<<endl;
+}
+
+void get_time( float *telps, float *tuser, float *tsys, float *tcpu )
+{
+
+  /* Variables declaration */
+  static time_t zsec = 0;
+  static long zusec = 0;
+  time_t nsec;
+  long nusec;
+  static clock_t zclock = 0;
+  clock_t nclock;
+  static clock_t zuser = 0;
+  static clock_t zsys = 0;
+  clock_t nuser, nsys;
+
+  struct timeval tp;
+  struct timezone tzp;
+  struct tms local;
+
+  MPI_Barrier(MPI_COMM_WORLD);
+
+  /* Elapsed time reading */
+
+  gettimeofday(&tp,&tzp);
+  nsec = tp.tv_sec;
+  nusec = tp.tv_usec;
+  *telps = (float)(nsec-zsec) + (float)(nusec-zusec)/(float)CLOCKS_PER_SEC;
+  
+  zsec = nsec;
+  zusec = nusec;
+
+  /* User and system CPU time reading */
+
+  times(&local);
+  nuser = local.tms_utime;
+  nsys = local.tms_stime;
+  *tuser = (float)(nuser-zuser) / (float)CLK_TCK;
+  *tsys = (float)(nsys-zsys) / (float)CLK_TCK;
+
+  zuser = nuser;
+  zsys = nsys;
+
+  /* CPU time reading */
+
+  nclock = clock();
+  *tcpu = (float)(nclock-zclock) / (float)CLOCKS_PER_SEC;
+  zclock = nclock;
+
+}
+
+
-- 
2.39.2