现在的位置: 首页 > 综合 > 正文

数值计算——矩阵求逆误差测试(n=2000)

2013年07月20日 ⁄ 综合 ⁄ 共 3187字 ⁄ 字号 评论关闭

23:37:01:~/..er/inverse$ ./main 2000
*****************************************************
no pivoting:
no_pivot_inverse:10.924682 seconds
no_pivot_mult:67.480217 seconds
A_original*Ainv:
        1.00568     -0.0088769     -0.0140972    -0.00617585     0.00524703     -0.0247897
     -0.0119802        1.00609     0.00365259    -0.00677496      0.0130823     -0.0080104
     -0.0128222    -0.00213862       0.979183      -0.008476      0.0305491     -0.0298309
     -0.0236932     -0.0021551     -0.0295723       0.993011      0.0449296     -0.0380806
      0.0172102     0.00938947      0.0163724      0.0109373       0.989297    -0.00549839
     -0.0101167    -0.00248698     -0.0186263    -0.00551816      0.0225024       0.973027
no_pivot_sub_norm:0.004001 seconds
(no pivot)N=2000:norm(A*Ainv-I):39.28850555
*****************************************************
*****************************************************
pivoting:
pivot_inverse:10.984687 seconds
pivot_Mult:67.696231 seconds
A_original*Ainv:
        1.00001   -1.89058e-05   -4.31342e-06   -2.10684e-05   -1.20411e-05     -1.413e-05
   -0.000140131        1.00017   -1.95652e-05    8.83937e-05    1.54078e-05    3.45707e-05
    2.19345e-05   -1.51992e-05        1.00001   -2.58684e-05    -9.6187e-06   -2.58088e-05
    0.000131607   -0.000198007    2.01538e-05       0.999871    -3.9719e-05   -7.89165e-05
    0.000182509   -0.000258565    1.26064e-05   -0.000145733       0.999962   -0.000104785
    2.27094e-05   -3.88026e-05    8.60542e-06   -2.16961e-05   -9.98378e-06       0.999961
pivot_sub_norm:0.004 seconds
(pivot)N=2000:norm(A*Ainv-I):0.1356286108
*****************************************************
------------------------------------------------------------------------------
------------------------------------------------------------------------------
------------------------------------------------------------------------------
------------------------------------------------------------------------------
*****************************************************
no pivoting(double):
no_pivot_inverse:21.437339 seconds
no_pivot_mult:78.980937 seconds
A_original*Ainv:
              1   -6.34344e-11    4.55167e-11    2.89133e-11   -3.87063e-11    1.32871e-10
    7.98329e-11              1    5.56323e-11    4.26117e-11   -4.24686e-11      1.475e-10
    8.67711e-11   -8.34241e-11              1     4.4937e-11   -4.64317e-11    1.69494e-10
    1.20053e-10   -1.15211e-10    8.11088e-11              1     -5.084e-11    2.12459e-10
    1.38203e-10   -1.36513e-10    1.00567e-10    6.94243e-11              1    2.27682e-10
    9.78832e-11   -1.07153e-10    7.18054e-11    4.29636e-11   -4.56121e-11              1
no_pivot_sub_norm:0.004 seconds
(no pivot)N=2000:norm(A*Ainv-I):1.746190785e-07
*****************************************************
*****************************************************
pivoting(double):
pivot_inverse:21.685356 seconds
pivot_Mult:78.932932 seconds
A_original*Ainv:
              1    -5.7343e-14    4.91759e-14   -4.10193e-14    4.32154e-14    4.10019e-14
   -1.77636e-14              1    1.50713e-14    4.52416e-14    2.54519e-14   -1.43219e-14
   -1.08247e-14    4.58522e-14              1    2.47302e-14    3.08711e-14   -1.33227e-14
    8.35096e-14    -3.1232e-14    6.28577e-14              1    5.05715e-14    5.20053e-14
   -9.88654e-14    1.42636e-13   -4.99323e-14      8.607e-14              1   -8.16153e-14
   -3.18148e-14    8.51368e-14    7.66921e-15    2.38004e-14    1.50444e-14              1
pivot_sub_norm:0.004 seconds
(pivot)N=2000:norm(A*Ainv-I):2.97245601e-10
*****************************************************

抱歉!评论已关闭.