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

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

*****************************************************