CPPLapack
|
//====================================================================[include] #include <iostream> #include <cstdlib> #include <ctime> #include "cpplapack.h" using namespace std; //=======================================================================[main] /*! main */ int main(int argc, char** argv) { srand(time(NULL)); int L(4); //// make dcovector x //// CPPL::dcovector x(L); for(int i=0; i<x.l; i++){ x(i) =double( rand() /(RAND_MAX/10) ); } //// print x in two ways //// cout << "x =\n" << x << endl; for(int i=0; i<x.l; i++){ cout << "x(" << i << ") =" << x(i) << endl; } //// make dcovector y //// CPPL::dcovector y(x); //// print y in two ways //// cout << "y =\n" << y << endl; for(int i=0; i<y.l; i++){ cout << "y(" << i << ") =" << y(i) << endl; } //// print x+y //// cout << "x+y=\n" << x+y << endl; //// write/read //// x.write( "tmp.txt" ); CPPL::dcovector z; z.read( "tmp.txt" ); cout << "z-x =\n" << z-x << "<-Should be zero." << endl; //// const //// CPPL::dcovector X(CPPL::dcovector(3).set(0,1).set(1,2).set(2,3)); cout << "X=\n" << X << endl; return 0; }
//====================================================================[include] #include <iostream> #include <cstdlib> #include <ctime> #include "cpplapack.h" using namespace std; //=======================================================================[main] /*! main */ int main(int argc, char** argv) { srand(time(NULL)); int L(4); //// make drovector x //// CPPL::drovector x(L); for(int i=0; i<x.l; i++){ x(i) =double( rand() /(RAND_MAX/10) ); } //// print x in two ways //// cout << "x =\n" << x << endl; for(int i=0; i<x.l; i++){ cout << "x(" << i << ") =" << x(i) << endl; } //// make drovector y //// CPPL::drovector y(x); //// print y in two ways //// cout << "y =\n" << y << endl; for(int i=0; i<y.l; i++){ cout << "y(" << i << ") =" << y(i) << endl; } //// print x+y //// cout << "x+y=\n" << x+y << endl; //// write/read //// x.write( "tmp.txt" ); CPPL::drovector z; z.read( "tmp.txt" ); cout << "z-x =\n" << z-x << "<-Should be zero." << endl; return 0; }