回 帖 发 新 帖 刷新版面

主题:矩阵求逆的子程序里数据溢出,求高手指点,急急急!!!!!!!!!!!

我编写了一个阻尼最小二乘法模拟波速的程序,迭代了四次显示数据溢出了,经查看发现数据溢出发生在矩阵求逆的子程序处,被求逆的矩阵是36*36的矩阵,矩阵中的大多数元素都是-11次方左右的数量级,请高手帮忙分析下这种溢出问题该怎么解决。急急急

回复列表 (共3个回复)

沙发

这么小得量级 需要real*16了, 用ivf吧,在32位的机子上也可以,另 矩阵求逆可以调用imsl库,应该 比自己写的要好

板凳

是的,自己也可以编成的。

3 楼

虽然 -11次方接近real*4的指数限(-31次方),但个人觉得可能还是你的程序有问题。

建议:一个算法子程序编写好以后,一定要先算几个标准例子。

我来回复

您尚未登录,请登录后再回复。点此登录或注册