求强大的C/C++矩阵运算类

[yangshucheng1986]

最近在处理数据时需要进行大矩阵的求逆运算，碰到了下面的两个矩阵求逆问题.
第一个问题：
有 A = {
1.70484441204335100000 , 0.89959011009135648000 , 0.03036127331075771900 ,
0.89959011009135648000 , 2.12292572413248370000 , 0.19604287211824964000 ,
0.03036127331075771900 , 0.19604287211824964000 , 0.02020286524503244600 }

利用我目前使用的矩阵运算类，求逆，得到逆矩阵：
AI = {
-1828419220961452.50 5014677334459694.00 -45913220746131040.00
5014677334459694.00 -13753404296155084.00 125922974768663220.00
-45913220746131040.00 125922974768663220.00 -1152921504606847000.00 }

然后将两个矩阵相乘
E = A*AI = {
0.500000000000 1.000000000000 -0.125000000000
0.000000000000 4.000000000000 0.000000000000
0.000000000000 0.000000000000 0.000000000000 }

结果显示求逆精度明显不够！

第二个问题：
B = {
2.08260837917549590000 , 0.37352310083165263000 , -0.73572925695810321000 ,
0.37352310083165263000 ,0.42182944496501185000 , 0.36553608709804863000 ,
-0.73572925695810321000 ,0.36553608709804863000 , 0.95741142753192932000 }
这个是求逆失败！

感觉是求逆函数的运算精度不够，求解决方法和强大矩阵运算类！！！！！！！！！！！！！！！！！！