关于矩阵计算

shshsh_0510 · 发表于 2008-4-15 14:52:15

您需要登录才可以下载或查看，没有账号？欢迎注册

×

前些日子有个问题需要计算一个矩阵的逆。但由于矩阵条件数很大，很病态。用maple算，它最高支持8字节浮点数。当时想编个任意精度的，但不通此道，不知各位高手有什么看法。能否推荐个现成的库或相关资料？谢了先。

无心人 · 发表于 2008-4-15 14:54:47

米有哦
现有的库似乎都精度不高

除非你自己用GMP写一个库

无心人 · 发表于 2008-4-15 14:57:08

对了
规模多大？

shshsh_0510 · 发表于 2008-4-15 15:13:43

100*100

无心人 · 发表于 2008-4-15 15:20:36

倒倒倒

是100行100列否？
最笨的办法也不过是100万次运算吧

mathe · 发表于 2008-4-15 16:26:38

呵呵，随便找一个矩阵程序程序，然后将double全局替换成mpf_class,然后加上#include <gmpxx.h>,链接gmp编译一把看看

mathe · 发表于 2008-4-15 16:28:05

关于矩阵运算库的源代码，可以查找LaPack, 这个是比较通用的库（Linux下面直接有，不过都是Fortran接口）

shshsh_0510 · 发表于 2008-4-15 16:43:52

谢谢math。gmp没用过，是可以任意精度的吗？回去学习一下

无心人 · 发表于 2008-4-15 16:45:44

不是任意精度的
但是很高精度的

mathe · 发表于 2008-4-16 10:05:44

呵呵，gmp其实也支持任意精度运算（有理数运算），那要使用mpq_class.不过具体到本题计算机内存是否受得了就不能保证了

账号		自动登录	找回密码
密码			欢迎注册

[讨论] 关于矩阵计算