Changeset 3570 in Sophya for trunk/SophyaProg

Timestamp:

Feb 7, 2009, 10:30:29 PM (17 years ago)

Author:

ansari

Message:

Ajouts nouveaux tests avec Rc sur operations TArray, Reza 07/02/2009

File:

• trunk/SophyaProg/Tests/arrt.cc (modified) (1 diff)

trunk/SophyaProg/Tests/arrt.cc

@@ -525 +525 @@
     rc += 4;
   }
-
+  
+  if (prtlev > 0)  
+    cout << " CheckArry: Starting Test-4 ArrayOp <float> " << endl;  nerr = 0;
+  NX = 7;
+  NY = 5;
+  NZ = 3;  
+  TArray<float> FA(NX,NY,NZ), FB, FM1(NX,NY,NZ), FM2(NX,NY,NZ);
+  for(int iz=0; iz<NZ; iz++)  
+    FM1(Range::all(), Range::all(), Range(iz, iz)) = M_PI*(iz+0.75);
+  for(int iy=0; iy<NY; iy++)  
+    FM2(Range::all(), Range(iy, iy), Range::all()) = (iy+2.22)*78.654;
+
+  FA = RandomSequence(RandomSequence::Gaussian, 7., 2.7);
+  FB = FA;  
+
+  float cst[3] = {6., 4., 8.2};
+
+  for(int ix=0; ix<NX; ix++)
+    for(int iy=0; iy<NY; iy++) 
+      for(int iz=0; iz<NZ; iz++)
+        FB(ix,iy,iz) = FB(ix,iy,iz)*FM1(ix,iy,iz)+cst[0]*FM2(ix,iy,iz);         
+
+  nerr = 0;
+  TArray<float> FD = ((FA&&FM1) + (cst[0]*FM2))/cst[1];
+  TArray<float> FE = (FD*cst[1] - FB);
+  float fmin, fmax;
+  FE.MinMax(fmin, fmax);
+  
+//  cout << FM1 << FM2;
+//  cout << FB << FD ; 
+//  cout << " ---- Resultat FE: " << FE;
+  if ((fabs(fmin)>1.e-5) || (fabs(fmax)>1e-5)) nerr++;
+  if (nerr == 0)
+    cout << "(4) CheckTArray <int>  ArrayOp <float> OK " << endl; 
+  else {
+    cout << "(4) CheckTArray <int>  ArrayOp <float>  ERROR  Rc += 8 !!! fmin=" 
+         << fmin << " fmax=" << fmax << endl; 
+    rc += 8;
+  }
+
+  if (prtlev > 0)  
+    cout << " CheckArry: Starting Test-5 PPF_Array IO " << endl;  
+  nerr = 0;
+  {
+    POutPersist po("arrt_ck.ppf");
+    TArray<float> FAS = FA(Range::all(), Range::all(), Range(0,0));
+    po << FA << FAS;    
+  }
+  {
+    PInPersist pin("arrt_ck.ppf");
+    TArray<float> RFA, RFAS;
+    pin >> RFA >> RFAS;
+    TArray<float> DA = RFA - FA;
+    TArray<float> DAS = RFAS - FA(Range::all(), Range::all(), Range(0,0));
+    DA.MinMax(fmin, fmax);
+    if ((fabs(fmin)>1.e-8) || (fabs(fmax)>1e-8)) nerr++;
+    DAS.MinMax(fmin, fmax);
+    if ((fabs(fmin)>1.e-8) || (fabs(fmax)>1e-8)) nerr++;
+  }
+  if (nerr == 0)
+    cout << "(5) CheckTArray <float>  PPF_Array IO OK " << endl; 
+  else {
+    cout << "(5) CheckTArray <float>  PPF_Array IO  ERROR  Rc += 16 !!!  nerr=" 
+         << nerr << endl; 
+    rc += 16;
+  }
+
+  if (prtlev > 0)  
+    cout << " CheckArry: Starting Test-6 Matrix inverse/multiply " << endl;  
+  nerr = 0;
+ 
+  sa_size_t MSZ = 32;
+  Matrix A(MSZ,MSZ), B(MSZ, MSZ);
+  A = 0.;
+  for(sa_size_t cc=0; cc<MSZ; cc+=2) {
+        sa_size_t ca = cc;
+        sa_size_t ra = MSZ-ca-2;
+        A(ra,ca) = 1.;   B(ca, ra) = 0.5;
+        A(ra+1,ca) = 1.;   B(ca, ra+1) = 0.5;
+        A(ra,ca+1) = 1.;   B(ca+1, ra) = 0.5;
+        A(ra+1,ca+1) = -1.;   B(ca+1, ra+1) = -0.5;
+  }
+//  cout << A << B;
+  Matrix C = A*B;
+  Matrix Identite(MSZ,MSZ);
+  Identite = IdentityMatrix(1., MSZ);
+  Matrix D = C-Identite;
+//  cout << "---- C=A*B=" << C ;
+  Matrix Ainv = Inverse(A);
+  Matrix E =B-Ainv;
+  double dmin, dmax;
+  D.MinMax(dmin, dmax);
+  if ((fabs(dmin)>1.e-7) || (fabs(dmax)>1e-7)) nerr++;
+  E.MinMax(dmin, dmax);
+  if ((fabs(dmin)>1.e-7) || (fabs(dmax)>1e-7)) nerr++;
+  if (nerr == 0)
+    cout << "(6) CheckTArray <double>  Matrix inverse/multiply OK " << endl; 
+  else {
+    cout << "(6) CheckTArray <double>  Matrix inverse/multiply  ERROR  Rc += 32 !!!  nerr=" 
+         << nerr << endl; 
+    rc += 32;
+  }
+ 
   return rc;
 }