24小时热门版块排行榜    

查看: 1684  |  回复: 1

yixijiumeng

银虫 (初入文坛)

[求助] 调用LAPACK中的DSYEVD进行对角化的误差问题。 已有1人参与

菜鸟一枚,并不是专门搞计算的。只是在解一个矩阵时,直接调用了DSYEVD将这个矩阵对角化,求出了本征值与本征矢。不过在和老板聊的时候,他提到了关于这个结果的误差问题。我只知道DSYEVD使用了divide and conquer 算法,所以想问问各位大虾,这是怎样的一种算法,这个DSYEVD子程序是如何控制误差的,它的精度能达到多少??
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fish.yfyh

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
yixijiumeng: 金币+5, 有帮助 2014-05-27 14:26:37
其实,LAPACK库这些子程序是很成熟的,不要太担心精度问题。
MATLAB在求解矩阵特征值的时候,其实也调用了LAPACK库 的。

当然你自己找一个例子,调用DSYEVD,把结果和MATLAB的结果一比就知道差距了。

我曾经测试过多次,至少小数点后10位是一样的。
2楼2014-05-26 01:16:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yixijiumeng 的主题更新
信息提示
请填处理意见