Changeset 2534 in Sophya for trunk/Eval

Timestamp:

May 7, 2004, 1:43:53 PM (21 years ago)

Author:

ansari

Message:

Programme cpupower : passage en double et ajout OPE=3 - Reza 7 Mai 2004

Location:

trunk/Eval/Speed

Files:

• cpupower.c (modified) (5 diffs)
• tompCXX.cc (modified) (1 diff)
• vectorCXX.cc (modified) (1 diff)

trunk/Eval/Speed/cpupower.c

@@ -15 +15 @@
 
 /*  Choix de type d'operations  float double int ... */
-#define  T_Type  float 
+#define  T_Type  double 
 
 static int SZ;                    /* Taille de tableau */
@@ -53 +53 @@
 }
 
+void fop_3() 
+{
+  T_Type s;
+  int i,j,k;
+  printf("--- fop_3: Double Loop: z[k] = Somme_i(x[k]*y[i]+x[i]*y[k]-0.85*(y[k]+x[k])) --- \n");
+  for(k=0; k<SZ; k++) {
+    for(i=0; i<SZ; i++)  s += x[k]*y[i]+x[i]*y[k]-0.85*(y[k]+x[k]);
+    z[k] = s;
+    N_OP += 7*SZ;
+  }
+  return;
+}
+
 /* Fonctions de timing (TCPU) - voir en fin de fichier */
 void InitTim(void);                                                             
@@ -70 +83 @@
   if (narg < 2) {
     printf("--- Programme cpupower: (Puissance de calcul) ---- \n");
-    printf("  Usage cpupower Op=1/2 [Size=20000] [CkPrt=0] \n");
-    return;
+    printf("  Usage cpupower Op=1/2/3 [Size=20000] [CkPrt=0] \n");
+    return 1;
   } 
   OPE = atoi(arg[1]);
@@ -90 +103 @@
   }
   PrtTim("--Fin malloc+init x,y,z ");
-  if (OPE == 2)  fop_2();
+  if (OPE == 3)  fop_3();
+  else if (OPE == 2)  fop_2();
   else fop_1();
   PrtTim("---Fin OpeDoubleBoucle ");
   mflops = N_OP/ GetPartialCPUTime()*1.e-6;
-  printf("-> Nb Operations~= %g  MFLOPS~= \n",N_OP,mflops);
+  printf("-> Nb Operations= %g  MFLOPS= %g \n",N_OP,mflops);
   if (ckprt > 0) {
     printf(" CheckPrint - ckprt= %d maxnprt= %d \n", ckprt, maxnprt);
@@ -107 +121 @@
 
   PrtTim("----Fin cpupower");
-  printf(":::::: FIN cpupower N_OP~= %g MFLOPS~= %g ::::::: \n", N_OP,mflops);
+  printf(":::::: FIN cpupower N_OP= %g MFLOPS= %g ::::::: \n", N_OP,mflops);
 
   free(x);
   free(y);
   free(z);
+  return 0;
 }
 

trunk/Eval/Speed/tompCXX.cc

@@ -346 +346 @@
 #pragma omp  single 
   if (k==0) printf("MultB_OMP3() -pragma2- NumThr= %d \n", omp_get_num_threads());
-#pragma omp for private(i) schedule(static) 
+  /* #pragma omp for private(i) schedule(static) */
 #endif
     for(i=0; i<sz; i++) 

trunk/Eval/Speed/vectorCXX.cc

@@ -16 +16 @@
 public:
   MyException(const char * msg) { _msg = msg; }
-//  ~MyException() { }
+  // virtual  ~MyException() { }
   string Msg() { return(_msg); }
 private: