Changeset 2557 in Sophya for trunk/SophyaProg/Tests/tsttminv.cc

Timestamp:

Jul 21, 2004, 6:23:00 PM (21 years ago)

Author:

cmv

Message:

cosmetique... cmv 21/07/04

File:

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

trunk/SophyaProg/Tests/tsttminv.cc

@@ -21 +21 @@
 
 // --- Choix de ce que doit faire Lapack
-#ifdef USE_LAPACK
 #define ALSO_LAPACK_INV
 #define ALSO_LAPACK_INV_SYM
@@ -27 +26 @@
 #define ALSO_LAPACK_EV
 #define ALSO_LAPACK_EV_SYM
-#endif
 
 //////////////////////////////////////////////////
@@ -161 +159 @@
         <<"-- Only GausPiv --"<<endl
         <<"-t 0/1/2 : data scaling 0=no, 1=global (def), 2=row-by-row"<<endl
-        <<"-d : also compute determinant"<<endl;
+        <<"-d : also compute determinant (ne pas utiliser si N est grand)"<<endl;
     return(-1);
   }
@@ -257 +255 @@
 
 
-#ifdef ALSO_LAPACK_INV
-////////////////////////////////////
-///////// Test avec Lapack /////////
-////////////////////////////////////
+//////////////////////////////////////////////
+///////// Test Inversion avec Lapack /////////
+//////////////////////////////////////////////
+#if defined(USE_LAPACK) && defined(ALSO_LAPACK_INV)
 {
 cout<<"\n=========================================="<<endl;
@@ -282 +280 @@
 
 
-#ifdef ALSO_LAPACK_INV_SYM
-////////////////////////////////////////
-///////// Test avec Lapack sym /////////
-////////////////////////////////////////
+//////////////////////////////////////////////////
+///////// Test Inversion avec Lapack sym /////////
+//////////////////////////////////////////////////
+#if defined(USE_LAPACK) && defined(ALSO_LAPACK_INV_SYM)
 {
 cout<<"\n=========================================="<<endl;
@@ -307 +305 @@
 
 
-#ifdef ALSO_LAPACK_INV_LSS
 ////////////////////////////////////////////////
 ///////// Test avec Lapack LeastSquare /////////
 ////////////////////////////////////////////////
+#if defined(USE_LAPACK) && defined(ALSO_LAPACK_INV_LSS)
 {
 cout<<"\n=========================================="<<endl;
@@ -332 +330 @@
 
 
-#ifdef ALSO_LAPACK_EV
 ////////////////////////////////////////////////
 ///////// Test avec Lapack pour EV /////////
 ////////////////////////////////////////////////
+#if defined(USE_LAPACK) && defined(ALSO_LAPACK_EV)
 {
 cout<<"\n=========================================="<<endl;
@@ -347 +345 @@
 PrtTim("--- End of LapackEigen decompositon ---");
 if(nprline>0) {cout<<Eval; cout<<endl; cout<<Evec; cout<<endl;}
-#ifndef COMPLEX
 //-- Find the complex conjugate pairs
+#if !defined(COMPLEX)
 TVector< uint_2 > Evalconj(N); Evalconj = 0;
 int_4 nconj=0;
@@ -369 +367 @@
   for(int i=0;i<N;i++) {
     complex<TYPER> Evec_il;
-    #ifdef COMPLEX
+    #if defined(COMPLEX)
     Evec_il = Evec(i,l);
     #else
@@ -378 +376 @@
     for(int j=0;j<N;j++) {
       complex<TYPER> Evec_jl;
-      #ifdef COMPLEX
+      #if defined(COMPLEX)
       Evec_jl = Evec(j,l);
       #else
@@ -398 +396 @@
 
 
-#ifdef ALSO_LAPACK_EV_SYM
 ////////////////////////////////////////////////
 ///////// Test avec Lapack sym pour EV /////////
 ////////////////////////////////////////////////
+#if defined(USE_LAPACK) && defined(ALSO_LAPACK_EV_SYM)
 {
 cout<<"\n=========================================="<<endl;
@@ -428 +426 @@
 
 
-#ifdef USE_GAUSPIV
-////////////////////////////////////
-///////// Test avec GausPiv /////////
-////////////////////////////////////
+///////////////////////////////////////////////
+///////// Test Inversion avec GausPiv /////////
+///////////////////////////////////////////////
+#if defined(USE_GAUSPIV)
 {
 cout<<"\n==========================================\n"