GPLIB++: testgeev.cpp Source File

00001 #include "VecMat.h"
00002 #include "UniformRNG.h"
00003 #include <iostream>
00004 #include <fstream>
00005 
00006 
00007 namespace ublas = boost::numeric::ublas;
00008 using namespace std;
00009 int main()
00010 {
00011         UniformRNG Random;
00012         const int msize = 3;
00013         gplib::cvec s(msize);
00014         gplib::cmat vl(msize,msize),vr(msize,msize), in(msize,msize);
00015                         
00016         ifstream infile ("in");
00017         infile >> in;
00018         cout << "In: " << in << endl;
00019         boost::numeric::bindings::lapack::geev(in,s,&vl,&vr,boost::numeric::bindings::lapack::optimal_workspace());
00020                         //boost::numeric::bindings::lapack::geev(in,s,NULL,NULL,boost::numeric::bindings::lapack::optimal_workspace());
00021         cout << "In after: " << in << endl;
00022         cout << "s: " << s << endl;
00023         cout << "vl: " << vl << endl;
00024         cout << "vr: " << vr << endl;
00025 }