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

Timestamp:

Jul 23, 2004, 12:54:27 PM (21 years ago)

Author:

cmv

Message:

add test of lapack SVD decomp by Divide and Conquer (cmv 23/07/04)

File:

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

trunk/SophyaProg/Tests/tsttminv.cc

@@ -27 +27 @@
 #define ALSO_LAPACK_EV_SYM
 #define ALSO_LAPACK_SVD
+#define ALSO_LAPACK_SVD_DC
 
 //////////////////////////////////////////////////
@@ -338 +339 @@
 cout<<"------------ Eigen decompositon LapackEigen "<<endl;
 A=Ainput;
-TMatrix< TYPE > Evec(N,N); Evec = (TYPE) 0;
-TVector< complex<r_8> > Eval(N); Eval = complex<r_8>(0,0);
+TMatrix< TYPE > Evec;
+TVector< complex<r_8> > Eval;
 //-- Decompositon
 int_4 info = LapackEigen(A,Eval,Evec,true);
@@ -434 +435 @@
 cout<<"------------ SVD decompositon LapackSVD "<<endl;
 A=Ainput;
-TVector< TYPE > S(N,N); TMatrix< TYPE > U(N,N); TMatrix< TYPE > VT(N,N); 
+TVector< TYPE > S; TMatrix< TYPE > U; TMatrix< TYPE > VT; 
 //-- Decompositon
 int_4 info = LapackSVD(A,S,U,VT);
@@ -453 +454 @@
 Check_Mat_Zero(AmUSVt);
 PrtTim("--- End of LapackSVD Test ---");
+}
+#endif
+
+
+///////////////////////////////////////////////////
+///////////// Test avec Lapack SVD_DC /////////////
+///////////////////////////////////////////////////
+#if defined(USE_LAPACK) && defined(ALSO_LAPACK_SVD_DC)
+{
+cout<<"\n=========================================="<<endl;
+cout<<"------------ SVD decompositon LapackSVD_DC "<<endl;
+A=Ainput;
+TVector< TYPE > S; TMatrix< TYPE > U; TMatrix< TYPE > VT; 
+//-- Decompositon
+int_4 info = LapackSVD_DC(A,S,U,VT);
+cout<<"info="<<info<<endl;
+PrtTim("--- End of LapackSVD_DC decompositon ---");
+if(nprline>0) {cout<<S; cout<<endl;}
+double smax = ABS_VAL(S(0)), smin = ABS_VAL(S(N-1));
+cout<<" Smin = |"<<S(N-1)<<"| = "<<smin<<", "
+    <<" Smax = |"<<S(0)<<"| = "<<smax<<", "
+    <<" --> Smin/Smax = "<<smin/smax<<endl;
+//-- A = U*S*VT
+cout<<"Compute A = U*S*VT"<<endl;
+TMatrix< TYPE > AmUSVt(N,N); AmUSVt = U;
+for(int i=0;i<N;i++) for(int j=0;j<N;j++) AmUSVt(i,j) *= S(j);
+AmUSVt *= VT; AmUSVt -= Ainput;
+if(nprline>0) {cout<<AmUSVt; cout<<endl;} 
+//-- Check
+Check_Mat_Zero(AmUSVt);
+PrtTim("--- End of LapackSVD_DC Test ---");
 }
 #endif