Menu

[r151]: / trunk / test / vc_sample / main.cpp  Maximize  Restore  History

Download this file

28 lines (23 with data), 772 Bytes

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
#include <cpplapack.h>
#include <ctime>
int main(int /*argc*/, char** /*argv*/)
{
//////// declare A and y ////////
const size_t size(3);
CPPL::dgematrix A(size,size); //3x3 matrix
CPPL::dcovector y(size); //3 column vector
//////// set components of A and y at random ////////
srand(unsigned(time(NULL)));
for(size_t i=0; i<size; i++){
for(size_t j=0; j<size; j++){
A(i,j) =double(rand())/double(RAND_MAX);
}
y(i) =double(rand())/double(RAND_MAX);
}
std::cout << "A=\n" << A << std::endl; //print A
std::cout << "y=\n" << y << std::endl; //print y
//////// solve A*x=y ////////
A.dgesv(y); //call dgesv of LAPACK (y becomes x)
std::cout << "x=\n" << y << std::endl; //print x
return 0;
}
Want the latest updates on software, tech news, and AI?
Get latest updates about software, tech news, and AI from SourceForge directly in your inbox once a month.