Changeset 995 in Sophya for trunk/SophyaProg/Tests

Timestamp:

May 2, 2000, 5:48:24 PM (25 years ago)

Author:

ansari

Message:

add LaPack comparison 2/5/00 cmv

File:

• trunk/SophyaProg/Tests/tsttminv.cc (modified) (7 diffs)

trunk/SophyaProg/Tests/tsttminv.cc

@@ -1 +1 @@
 // Test de l'inversion de matrice (cmv 14/04/00)
+// Pb numeriques: dapax22 et daplxa115
+//                tsttminv -v 55,100,130,1000000000 -l 0 -a 1 
+
+// if defined COMPLEX , if not REAL
+//#define COMPLEX
+// if defined 32 bits precision, if not 64 bits
+// #define PRECIS32
+#define ALSO_LAPACK
+
 #include "machdefs.h"
 #include <iostream.h>
@@ -12 +21 @@
 #include "srandgen.h"
 
-// if defined COMPLEX , if not REAL
-#define COMPLEX
-// if defined 32 bits precision, if not 64 bits
-// #define PRECIS32
+#ifdef ALSO_LAPACK
+#include "intflapack.h"
+#endif
+
 
 #if defined(COMPLEX)
@@ -79 +88 @@
 //-- Initialization arrays
 SophyaInit();
+#ifdef ALSO_LAPACK
+  BaseArray::SetDefaultMemoryMapping(BaseArray::FortranMemoryMapping);
+#endif
 
 //-- Definition arrays
@@ -85 +97 @@
 TMatrix< TYPE > AiA(N,N);
 TMatrix< TYPE > B(N,N);
-TMatrix< TYPE > C(N,N);
+TMatrix< TYPE > Adum(N,N);
 TVector< int >  Vind(N*N);
 A.Show();
@@ -94 +106 @@
 AiA  = (TYPE) 0;
 B    = (TYPE) 0;
-C    = (TYPE) 0;
+Adum = (TYPE) 0;
 BaseArray::SetMaxPrint(nprline*N,0);
 
@@ -122 +134 @@
 #endif
 cout<<"Nombre de valeurs BIG reelles = "<<nbr<<", imaginaires = "<<nbi<<endl;
+Adum = A;
 
 //-- Print matrice A
@@ -154 +167 @@
 cout<<"Ecart maximum par rapport a 0 hors diagonale = "<<vmax<<endl;
 
+////////////////////////////////////
+///////// Test avec Lapack /////////
+////////////////////////////////////
+#ifdef ALSO_LAPACK
+
+cout<<endl<<"LAPACK Cross-Check"<<endl;
+InvA = IdentityMatrix(1.,N);
+LapackLinSolve(Adum,InvA);
+AiA = A * InvA;
+vmax=-1.;
+for(i=0;i<N;i++) {
+  double absv = ABS_VAL( 1. - AiA(i,i) );
+  if( absv > vmax ) vmax = absv;
+}
+cout<<"Ecart maximum par rapport a 1 sur diagonale = "<<vmax<<endl;
+vmax = -1.;
+for(i=0;i<N;i++) for(j=0;j<N;j++) {
+  if( i == j ) continue;
+  double absv = ABS_VAL( AiA(i,j) );
+  if( absv > vmax ) vmax = absv;
+}
+cout<<"Ecart maximum par rapport a 0 hors diagonale = "<<vmax<<endl;
+
+#endif
+
 exit(0);
 }