CPPLapack
|
//============================================================================= /*! dsysv_check */ void dsysv_check_vector() { cout << "############ check dsysv vector ############" << endl; srand(time(NULL)); int N(3); //// make dsymatrix A //// CPPL::dsymatrix A(N); for(int i=0; i<A.n; i++){ for(int j=0; j<A.n; j++){ A(i,j) =double( rand() /(RAND_MAX/10) ); }} //// make dcovector y //// CPPL::dcovector y(N); for(int i=0; i<y.l; i++){ y(i) =double( rand() /(RAND_MAX/10) ); } //// make A_original and y_original //// CPPL::dsymatrix A_original(A); CPPL::dcovector y_original(y); cout << "A_original=\n" << A_original << endl; cout << "y_original=\n" << y_original << endl; //// solve Ax=y //// A.dsysv(y); //// print A, y and A_original*y //// cout << "A=\n" << A << endl; cout << "y=\n" << y << endl; cout << "A_original*y=\n" << A_original*y << endl; } //============================================================================= void dsysv_check_matrix() { cout << "############ check dsysv matrix ############" << endl; srand(time(NULL)); int N(3); //// make dsymatrix A //// CPPL::dsymatrix A(N); for(int i=0; i<A.n; i++){ for(int j=0; j<A.n; j++){ A(i,j) =double( rand() /(RAND_MAX/10) ); }} //// make dsymatrix Y //// CPPL::dgematrix Y(N,N); for(int i=0; i<Y.m; i++){ for(int j=0; j<Y.n; j++){ Y(i,j) =double( rand() /(RAND_MAX/10) ); }} //// make A_original and Y_original //// CPPL::dsymatrix A_original(A); CPPL::dgematrix Y_original(Y); cout << "A_original=\n" << A_original << endl; cout << "Y_original=\n" << Y_original << endl; //// solve AY=B //// A.dsysv(Y); //// print A, Y and A_original*Y //// cout << "A=\n" << A << endl; cout << "Y=\n" << Y << endl; cout << "A_original*Y=\n" << A_original*Y << endl; }
//============================================================================= void dsyev_check_value() { cout << "############ check dsyev value ############" << endl; srand(time(NULL)); int N(3); //// make dsymatrix A //// CPPL::dsymatrix A(N); for(int i=0; i<N; i++){ for(int j=0; j<N; j++){ A(i,j) =double( rand() /(RAND_MAX/10) ); }} //// make wr wi vr //// vector<double> w; //// make A_original //// CPPL::dsymatrix A_original(A); //// dsyev //// A.dsyev(w); //A.dsyev(w,1); //// print //// cout << "A=\n" << A << endl; cout << "A_original=\n" << A_original << endl; for(int i=0; i<N; i++){ cout << "#### " << i << "th eigen ####" << endl; cout << "w=" << w[i] <<endl; } } //============================================================================= void dsyev_check_right() { cout << "############ check dsyev right ############" << endl; srand(time(NULL)); int N(3); //// make dsymatrix A //// CPPL::dsymatrix A(N); for(int i=0; i<A.n; i++){ for(int j=0; j<A.n; j++){ A(i,j) =double( rand() /(RAND_MAX/10) ); }} //// make wr wi vr //// vector<double> w; vector<CPPL::dcovector> v; //// make A_original //// CPPL::dsymatrix A_original(A); //// dsyev //// A.dsyev(w,v); //// print //// cout << "A_original=\n" << A_original << endl; for(int i=0; i<A.n; i++){ cout << "#### " << i << "th eigen ####" << endl; cout << "w=" << w[i] <<endl; cout << "v=\n" << v[i] <<endl; cout << "[A]*{x} -lambda*{x} = (Should be zeros)\n" << A_original*v[i] -w[i]*v[i] << endl; } } //============================================================================= void dsyev_check_left() { cout << "############ check dsyev left ############" << endl; srand(time(NULL)); int N(3); //// make dsymatrix A //// CPPL::dsymatrix A(N); for(int i=0; i<A.n; i++){ for(int j=0; j<A.n; j++){ A(i,j) =double( rand() /(RAND_MAX/10) ); }} //// make wr wi vl //// vector<double> w; vector<CPPL::drovector> v; //// make A_original //// CPPL::dsymatrix A_original(A); //// dsyev //// A.dsyev(w, v); //// print //// cout << "A_original=\n" << A_original << endl; for(int i=0; i<A.n; i++){ cout << "#### " << i << "th eigen ####" << endl; cout << "w = " << w[i] << endl; cout << "v = " << v[i] << endl; cout << "{x}*[A] -{x}*lambda = (Should be zeros)\n" << v[i]*A_original -v[i]*w[i] << endl; } }
//====================================================================[include] #include <iostream> #include <cstdlib> #include <ctime> #include "cpplapack.h" using namespace std; #include "dsysv_check.hpp" #include "dsyev_check.hpp" #include "dsygv_check.hpp" //=======================================================================[main] /*! main */ int main(int argc, char** argv) { //////// dsysv //////// dsysv_check_vector(); dsysv_check_matrix(); //////// dsyev //////// dsyev_check_value(); dsyev_check_right(); dsyev_check_left(); //////// dsygv //////// dsygv_check_value(); dsygv_check_both(); return 0; }