Changeset 2584 in Sophya for trunk/Eval

Timestamp:

Jul 30, 2004, 1:04:18 PM (21 years ago)

Author:

ansari

Message:

Ajout test simple boucle (fop0) ds cpupower.c - Reza 30/07/2004

File:

• trunk/Eval/Speed/cpupower.c (modified) (6 diffs)

trunk/Eval/Speed/cpupower.c

@@ -17 +17 @@
 #define  T_Type  double 
 
-static int SZ;                    /* Taille de tableau */
-static double N_OP;               /* Nb operations */
-static  int OPE;                  /* Choix Operation fop_1/2/3/4    */
-static int ckprt;                 /* > 0 : Print check apres calcul */
+static int SZ=20000;                    /* Taille de tableau */
+static double N_OP=0;               /* Nb operations */
+static  int OPE=0;                  /* Choix Operation fop_1/2/3/4    */
+static int ckprt=0;                 /* > 0 : Print check apres calcul */
 
 /*  Tableaux X,Y,Z   */
@@ -26 +26 @@
 static T_Type * y;
 static T_Type * z;
+
+void fop_0()
+{
+  T_Type s=0;
+  int i,k;
+  printf("--- fop_0: Simple Loop: z[k] = x[k]*y[i] --- \n");
+  for(k=0; k<SZ; k++) z[k] = x[k]*y[k];
+  N_OP += SZ;
+  return;
+}
 
 void fop_1()
@@ -101 +111 @@
 {
   int maxnprt = 1000;
+  int nloop = 1000;
   int i,nprt;
   double mflops;
@@ -109 +120 @@
   if (narg < 2) {
     printf("--- Programme cpupower: (Puissance de calcul) ---- \n");
-    printf("  Usage cpupower Op=1/2/3 [Size=20000] [CkPrt=0] \n");
+    printf("  Usage cpupower Op=0/1/2/3/4/5 [Size=20000] [NLoop=1000] [CkPrt=0]\n");
+    printf("  Op=0 Simple Loop: z[k] = x[k]*y[k]  effectue NLoop fois\n");
+    printf("    Size=10^6 par defaut, NLoop=100 \n");
     printf("  Op=1 Double Loop: z[k] = Somme_i(x[k]*y[i])  \n");
     printf("    Op=1 *,+=2 op / tour de boucle \n");
@@ -120 +133 @@
     printf("  Op=5 Double Loop: z[k] = Somme_i(x[k]*sin(y[i])+log(fabs(x[i])+0.1)+y[k]) \n");
     printf("    Op=5   ~ 100 op float /  tour de boucle \n");
+    printf("    Test 4,5 Size=2000 par defaut  \n");
     return 1;
   } 
   OPE = atoi(arg[1]);
   SZ = 20000;
+  if (OPE == 0) SZ = 1000000;
+  if (OPE >= 4) SZ = 2000;
   if (narg > 2)  SZ = atoi(arg[2]);
+  nloop = 100;
+  if (narg > 3)  nloop = atoi(arg[3]);
   ckprt = 0;
-  if (narg > 3)  ckprt = atoi(arg[3]);
+  if (narg > 3)  ckprt = atoi(arg[4]);
 
   printf("::::::: cpupower: OPE=%d SZ= %d ::::::: \n", OPE,SZ);
@@ -136 +154 @@
     x[i] = (T_Type)((random()%10000)-5000);  
     y[i] = (T_Type)((random()%10000)-5000);  
-    y[i] = (T_Type)(0);
+    z[i] = (T_Type)(0);
   }
   PrtTim("--Fin malloc+init x,y,z ");
+  N_OP = 0;
   if (OPE == 5)  fop_5();
   else if (OPE == 4)  fop_4();
   else if (OPE == 3)  fop_3();
   else if (OPE == 2)  fop_2();
-  else fop_1();
+  else if (OPE == 1)  fop_1();
+  else for(i=0; i<nloop; i++) fop_0();
   PrtTim("---Fin OpeDoubleBoucle ");
   mflops = N_OP/ GetPartialCPUTime()*1.e-6;