home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD2.mdf / c / library / dos / diverses / leda / prog / basic / matrix1.c < prev    next >
Encoding:
C/C++ Source or Header  |  1991-11-15  |  931 b   |  57 lines
  1. #include <LEDA/vector.h>
  2. #include <LEDA/matrix.h>
  3.  
  4. main()
  5.  
  6. { 
  7.   int d = read_int("dimension = ");
  8.  
  9.   matrix       A(d,d);
  10.   matrix       I(d,d);
  11.   double       det;
  12.  
  13.  
  14.   init_random();
  15.  
  16.   for(int i=0;i<d;i++)
  17.     for(int j=0;j<d;j++)
  18.       A(i,j) = double(random(-1000,1000))/100;
  19.  
  20.   float T = used_time();
  21.  
  22.  
  23.   cout << "A.inv():  ";
  24.   cout.flush();
  25.   I =  A.inv();
  26.   cout << string("%5.2f sec\n",used_time(T));
  27.  
  28.   cout << "A.det():  ";
  29.   cout.flush();
  30.   det = A.det();
  31.   cout << string("%5.2f sec \n",used_time(T));
  32.   newline;
  33.  
  34. /*
  35.   cout << "A*A.inv() = \n" << A*I << "\n";
  36.   newline;
  37. */
  38.  
  39.   while (Yes("A.solve(v) ? "))
  40.   {
  41.     vector v(d), x(d), y(d);
  42.  
  43.     for(int i=0;i<d;i++) v[i] = double(random(-1000,1000))/100;
  44.     
  45.     used_time(T);
  46.     x = A.solve(v);
  47.     cout << string("time for solve:  %5.2f \n",used_time(T));
  48.     newline;
  49.  
  50.     cout << "v   = " << v   << "\n";
  51.     cout << "A*x = " << A*x << "\n";
  52.     newline;
  53.   }
  54.  
  55.  
  56. }
  57.