Fortran 求解复数矩阵SVD时，用到cgesvd，结果与matlab不同，为什么？ - 程序语言 - Fortran

版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

somomo91

专家顾问 (职业作家)

专家经验: +3
应助: 144 (高中生)
金币: 19353
红花: 115
帖子: 3854
在线: 1053.2小时
虫号: 2399301
注册: 2013-04-03
专业: 动力学与控制
管辖: 计算模拟

引用回帖:

10楼: Originally posted by 长路漫漫 at 2013-06-27 17:55:49
啊？？有点懵。。。我希望？呵呵，没有希望，因为得到的SVD结果不和matlab近似，所以后来的运算，结果就很不对，所以，如果按这样说，我就是希望matlab对的
呵呵...

正常来说，MATLAB应该是对的
事实上，在MATLAB官方论坛里，曾经提到过这个问题： FORTRAN 下LAPACK SVD 结果和 MATLAB 不一样
有人提到，其中是在FORTRAN下参数 lwork 太小，修改之后，两个结果相同了
MATLAB的预处理，会尽量优化这些参数

Qianyong Chen, I think the segfault was because lwork was too small; doing as follows doesn't seem to crash matlab. The fortranError as well as fortranError_SVD remain both numerically zero.
http://www.mathworks.fr/matlabcentral/fileexchange/16777-lapack

赞一下

回复此楼

11楼2013-06-27 18:00:25

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

somomo91

专家顾问 (职业作家)

专家经验: +3
应助: 144 (高中生)
金币: 19353
红花: 115
帖子: 3854
在线: 1053.2小时
虫号: 2399301
注册: 2013-04-03
专业: 动力学与控制
管辖: 计算模拟

引用回帖:

刚刚想编译一下FORTRAN版的MEX，没想到他们只支持 Intel Visual Fortran，
再下一个 Intel Visual Fortran，太麻烦了，
不试了，

赞一下

回复此楼

12楼2013-06-27 18:04:47

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

长路漫漫

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 67
红花: 2
帖子: 19
在线: 13.8小时
虫号: 2167661
注册: 2012-12-05
专业: 电磁场与波

引用回帖:

12楼: Originally posted by somomo91 at 2013-06-27 18:04:47
刚刚想编译一下FORTRAN版的MEX，没想到他们只支持 Intel Visual Fortran，
再下一个 Intel Visual Fortran，太麻烦了，
不试了，...

哦，还是谢谢你了，只是想问一下，如果你遇到类似的问题，你搜网页查找原因的过程是怎么样？
因为你发的链接，我就没有从Google 搜索到。。。
只是想学习解决问题的方法～
谢谢～

赞一下

回复此楼

13楼2013-06-28 09:34:39

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

somomo91

专家顾问 (职业作家)

专家经验: +3
应助: 144 (高中生)
金币: 19353
红花: 115
帖子: 3854
在线: 1053.2小时
虫号: 2399301
注册: 2013-04-03
专业: 动力学与控制
管辖: 计算模拟

引用回帖:

13楼: Originally posted by 长路漫漫 at 2013-06-28 09:34:39
哦，还是谢谢你了，只是想问一下，如果你遇到类似的问题，你搜网页查找原因的过程是怎么样？
因为你发的链接，我就没有从Google 搜索到。。。
只是想学习解决问题的方法～
谢谢～...

墙里墙外各不同

回复此楼