数值计算——矩阵求逆误差测试

N= 500

float  no pivoting: 

23:19:13:~/..er/inverse$ ./main 500
*****************************************************

// float
no pivoting:
no_pivot_inverse:0.088006 seconds
no_pivot_mult:0.148009 seconds
A_original*Ainv:
       0.941271     -0.0576429    -0.00558718      -0.105509      0.0770019      0.0143333
     -0.0730097       0.926172    -0.00830387      -0.129601       0.088998      0.0141948
      -0.113874      -0.111164       0.990121      -0.201627       0.139809      0.0223503
      -0.138538      -0.135908      -0.012533        0.75218        0.17064      0.0292929
      0.0119933       0.007765    -0.00265124      0.0239916       0.985565    -0.00296608
      -0.053791     -0.0581362      -0.010759      -0.102463      0.0724007        1.01128
no_pivot_sub_norm:0 seconds
(no pivot)N=500:norm(A*Ainv-I):44.36532593
*****************************************************
*****************************************************

//float
pivoting:
pivot_inverse:0.096006 seconds
pivot_Mult:0.144009 seconds
A_original*Ainv:
       0.999998    6.85453e-06   -2.08616e-06    1.81794e-06   -2.87592e-06     1.0252e-05
    2.32272e-06       0.999986   -2.63564e-06    4.11719e-05    1.54972e-05   -3.00407e-05
    1.13025e-05   -4.70281e-05       0.999982    0.000106543    6.32703e-05   -0.000158131
   -1.03563e-06    -1.3411e-05   -3.29316e