Re: Octave/C++ matrix Inv() comparison

[ionone]

Also to do it simpler, here are the numbers in a 20x20 matrix in Octave and
in C++

the funny thing is that when i take the inverse in octave, ihave the exact
answer, but if i take the inverse with Lapack i have something totally
different...Not even close as before. Dunno...

R=[ 0.100197 0.100187 0.100155 0.100103 0.100030 0.099936 0.099822 0.099687
0.099531 0.099355 0.099158 0.098942 0.098705 0.098448 0.098171 0.097875
0.097559 0.097224 0.096870 0.096496;

0.100187 0.100197 0.100187 0.100155 0.100103 0.100030 0.099936 0.099822
0.099687 0.099531 0.099355 0.099158 0.098942 0.098705 0.098448 0.098171
0.097875 0.097559 0.097224 0.096870;

0.100155 0.100187 0.100197 0.100187 0.100155 0.100103 0.100030 0.099936
0.099822 0.099687 0.099531 0.099355 0.099158 0.098942 0.098705 0.098448
0.098171 0.097875 0.097559 0.097224;

0.100103 0.100155 0.100187 0.100197 0.100187 0.100155 0.100103 0.100030
0.099936 0.099822 0.099687 0.099531 0.099355 0.099158 0.098942 0.098705
0.098448 0.098171 0.097875 0.097559;

0.100030 0.100103 0.100155 0.100187 0.100197 0.100187 0.100155 0.100103
0.100030 0.099936 0.099822 0.099687 0.099531 0.099355 0.099158 0.098942
0.098705 0.098448 0.098171 0.097875;

0.099936 0.100030 0.100103 0.100155 0.100187 0.100197 0.100187 0.100155
0.100103 0.100030 0.099936 0.099822 0.099687 0.099531 0.099355 0.099158
0.098942 0.098705 0.098448 0.098171;

0.099822 0.099936 0.100030 0.100103 0.100155 0.100187 0.100197 0.100187
0.100155 0.100103 0.100030 0.099936 0.099822 0.099687 0.099531 0.099355
0.099158 0.098942 0.098705 0.098448;

0.099687 0.099822 0.099936 0.100030 0.100103 0.100155 0.100187 0.100197
0.100187 0.100155 0.100103 0.100030 0.099936 0.099822 0.099687 0.099531
0.099355 0.099158 0.098942 0.098705;

0.099531 0.099687 0.099822 0.099936 0.100030 0.100103 0.100155 0.100187
0.100197 0.100187 0.100155 0.100103 0.100030 0.099936 0.099822 0.099687
0.099531 0.099355 0.099158 0.098942;

0.099355 0.099531 0.099687 0.099822 0.099936 0.100030 0.100103 0.100155
0.100187 0.100197 0.100187 0.100155 0.100103 0.100030 0.099936 0.099822
0.099687 0.099531 0.099355 0.099158;

0.099158 0.099355 0.099531 0.099687 0.099822 0.099936 0.100030 0.100103
0.100155 0.100187 0.100197 0.100187 0.100155 0.100103 0.100030 0.099936
0.099822 0.099687 0.099531 0.099355;

0.098942 0.099158 0.099355 0.099531 0.099687 0.099822 0.099936 0.100030
0.100103 0.100155 0.100187 0.100197 0.100187 0.100155 0.100103 0.100030
0.099936 0.099822 0.099687 0.099531;

0.098705 0.098942 0.099158 0.099355 0.099531 0.099687 0.099822 0.099936
0.100030 0.100103 0.100155 0.100187 0.100197 0.100187 0.100155 0.100103
0.100030 0.099936 0.099822 0.099687;

0.098448 0.098705 0.098942 0.099158 0.099355 0.099531 0.099687 0.099822
0.099936 0.100030 0.100103 0.100155 0.100187 0.100197 0.100187 0.100155
0.100103 0.100030 0.099936 0.099822;

0.098171 0.098448 0.098705 0.098942 0.099158 0.099355 0.099531 0.099687
0.099822 0.099936 0.100030 0.100103 0.100155 0.100187 0.100197 0.100187
0.100155 0.100103 0.100030 0.099936;

0.097875 0.098171 0.098448 0.098705 0.098942 0.099158 0.099355 0.099531
0.099687 0.099822 0.099936 0.100030 0.100103 0.100155 0.100187 0.100197
0.100187 0.100155 0.100103 0.100030;

0.097559 0.097875 0.098171 0.098448 0.098705 0.098942 0.099158 0.099355
0.099531 0.099687 0.099822 0.099936 0.100030 0.100103 0.100155 0.100187
0.100197 0.100187 0.100155 0.100103;

0.097224 0.097559 0.097875 0.098171 0.098448 0.098705 0.098942 0.099158
0.099355 0.099531 0.099687 0.099822 0.099936 0.100030 0.100103 0.100155
0.100187 0.100197 0.100187 0.100155;

0.096870 0.097224 0.097559 0.097875 0.098171 0.098448 0.098705 0.098942
0.099158 0.099355 0.099531 0.099687 0.099822 0.099936 0.100030 0.100103
0.100155 0.100187 0.100197 0.100187;

0.096496 0.096870 0.097224 0.097559 0.097875 0.098171 0.098448 0.098705
0.098942 0.099158 0.099355 0.099531 0.099687 0.099822 0.099936 0.100030
0.100103 0.100155 0.100187 0.100197]


and here is the same in C++ (i'll save you some time this way)
double R21[20][20]={{0.100197 ,  0.100187 ,  0.100155  , 0.100103 , 
0.100030 ,  0.099936 ,  0.099822 , 0.099687 ,  0.099531 , 0.099355 , 
0.099158 ,  0.098942 ,  0.098705 ,  0.098448 , 0.098171  , 0.097875 , 
0.097559 ,  0.097224 ,  0.096870 , 0.096496},
   {0.100187 ,  0.100197 ,  0.100187  , 0.100155 ,  0.100103 , 0.100030 , 
0.099936 , 0.099822 ,  0.099687 ,  0.099531 , 0.099355 ,  0.099158 , 
0.098942 ,  0.098705 , 0.098448  , 0.098171 ,  0.097875 ,  0.097559 , 
0.097224 ,  0.096870},
   {0.100155 ,  0.100187 ,  0.100197  , 0.100187 ,  0.100155 , 0.100103 , 
0.100030 , 0.099936 ,  0.099822 ,  0.099687 , 0.099531 ,  0.099355 , 
0.099158 ,  0.098942 , 0.098705  , 0.098448 ,  0.098171 ,  0.097875 , 
0.097559 ,  0.097224},
   {0.100103 ,  0.100155 ,  0.100187  , 0.100197 ,  0.100187 , 0.100155 , 
0.100103 , 0.100030 ,  0.099936 ,  0.099822 , 0.099687 ,  0.099531 , 
0.099355 ,  0.099158 , 0.098942  , 0.098705 ,  0.098448 ,  0.098171 , 
0.097875 ,  0.097559},
   {0.100030 ,  0.100103 ,  0.100155  , 0.100187 ,  0.100197 , 0.100187 , 
0.100155 , 0.100103 ,  0.100030 ,  0.099936 , 0.099822 ,  0.099687 , 
0.099531 ,  0.099355 , 0.099158  , 0.098942 ,  0.098705 ,  0.098448 , 
0.098171 ,  0.097875},
   {0.099936 ,  0.100030 ,  0.100103  , 0.100155 ,  0.100187 , 0.100197 , 
0.100187 , 0.100155 ,  0.100103 ,  0.100030 , 0.099936 ,  0.099822 , 
0.099687 ,  0.099531 , 0.099355  , 0.099158 ,  0.098942 ,  0.098705 , 
0.098448 ,  0.098171},
   {0.099822 ,  0.099936 ,  0.100030  , 0.100103 ,  0.100155 , 0.100187 , 
0.100197 , 0.100187 ,  0.100155 ,  0.100103 , 0.100030 ,  0.099936 , 
0.099822 ,  0.099687 , 0.099531  , 0.099355 ,  0.099158 ,  0.098942 , 
0.098705 ,  0.098448},
   {0.099687 ,  0.099822 ,  0.099936  , 0.100030 ,  0.100103 , 0.100155 , 
0.100187 , 0.100197 ,  0.100187 ,  0.100155 , 0.100103 ,  0.100030 , 
0.099936 ,  0.099822 , 0.099687  , 0.099531 ,  0.099355 ,  0.099158 , 
0.098942 ,  0.098705},
   {0.099531 ,  0.099687 ,  0.099822  , 0.099936 ,  0.100030 , 0.100103 , 
0.100155 , 0.100187 ,  0.100197 ,  0.100187 , 0.100155 ,  0.100103 , 
0.100030 ,  0.099936 , 0.099822  , 0.099687 ,  0.099531 ,  0.099355 , 
0.099158 ,  0.098942},
   {0.099355 ,  0.099531 ,  0.099687  , 0.099822 ,  0.099936 , 0.100030 , 
0.100103 , 0.100155 ,  0.100187 ,  0.100197 , 0.100187 ,  0.100155 , 
0.100103 ,  0.100030 , 0.099936  , 0.099822 ,  0.099687 ,  0.099531 , 
0.099355 ,  0.099158},
   {0.099158 ,  0.099355 ,  0.099531  , 0.099687 ,  0.099822 , 0.099936 , 
0.100030 , 0.100103 ,  0.100155 ,  0.100187 , 0.100197 ,  0.100187 , 
0.100155 ,  0.100103 , 0.100030  , 0.099936 ,  0.099822 ,  0.099687 , 
0.099531 ,  0.099355},
   {0.098942 ,  0.099158 ,  0.099355  , 0.099531 ,  0.099687 , 0.099822 , 
0.099936 , 0.100030 ,  0.100103 ,  0.100155 , 0.100187 ,  0.100197 , 
0.100187 ,  0.100155 , 0.100103  , 0.100030 ,  0.099936 ,  0.099822 , 
0.099687 ,  0.099531},
   {0.098705 ,  0.098942 ,  0.099158  , 0.099355 ,  0.099531 , 0.099687 , 
0.099822 , 0.099936 ,  0.100030 ,  0.100103 , 0.100155 ,  0.100187 , 
0.100197 ,  0.100187 , 0.100155  , 0.100103 ,  0.100030 ,  0.099936 , 
0.099822 ,  0.099687},
   {0.098448 ,  0.098705 ,  0.098942  , 0.099158 ,  0.099355 , 0.099531 , 
0.099687 , 0.099822 ,  0.099936 ,  0.100030 , 0.100103 ,  0.100155 , 
0.100187 ,  0.100197 , 0.100187  , 0.100155 ,  0.100103 ,  0.100030 , 
0.099936 ,  0.099822},
   {0.098171 ,  0.098448 ,  0.098705  , 0.098942 ,  0.099158 , 0.099355 , 
0.099531 , 0.099687 ,  0.099822 ,  0.099936 , 0.100030 ,  0.100103 , 
0.100155 ,  0.100187 , 0.100197  , 0.100187 ,  0.100155 ,  0.100103 , 
0.100030 ,  0.099936},
   {0.097875 ,  0.098171 ,  0.098448  , 0.098705 ,  0.098942 , 0.099158 , 
0.099355 , 0.099531 ,  0.099687 ,  0.099822 , 0.099936 ,  0.100030 , 
0.100103 ,  0.100155 , 0.100187  , 0.100197 ,  0.100187 ,  0.100155 , 
0.100103 ,  0.100030},
   {0.097559 ,  0.097875 ,  0.098171  , 0.098448 ,  0.098705 , 0.098942 , 
0.099158 , 0.099355 ,  0.099531 ,  0.099687 , 0.099822 ,  0.099936 , 
0.100030 ,  0.100103 , 0.100155  , 0.100187 ,  0.100197 ,  0.100187 , 
0.100155 ,  0.100103},
   {0.097224 ,  0.097559 ,  0.097875  , 0.098171 ,  0.098448 , 0.098705 , 
0.098942 , 0.099158 ,  0.099355 ,  0.099531 , 0.099687 ,  0.099822 , 
0.099936 ,  0.100030 , 0.100103  , 0.100155 ,  0.100187 ,  0.100197 , 
0.100187 ,  0.100155},
   {0.096870 ,  0.097224 ,  0.097559  , 0.097875 ,  0.098171 , 0.098448 , 
0.098705 , 0.098942 ,  0.099158 ,  0.099355 , 0.099531 ,  0.099687 , 
0.099822 ,  0.099936 , 0.100030  , 0.100103 ,  0.100155 ,  0.100187 , 
0.100197 ,  0.100187},
   {0.096496 ,  0.096870 ,  0.097224  , 0.097559 ,  0.097875 , 0.098171 , 
0.098448 , 0.098705 ,  0.098942 ,  0.099158 , 0.099355 ,  0.099531 , 
0.099687 ,  0.099822 , 0.099936  , 0.100030 ,  0.100103 ,  0.100155 , 
0.100187 ,  0.100197}}; 



--
View this message in context: 
http://octave.1599824.n4.nabble.com/Octave-C-matrix-Inv-comparison-tp4655291p4655512.html
Sent from the Octave - Maintainers mailing list archive at Nabble.com.