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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » Fortran 求解复数矩阵SVD时，用到cgesvd，结果与matlab不同，为什么？
162/2返回列表上一页12
查看: 3325  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

somomo91

专家顾问 (职业作家)
引用回帖:
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

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

刚刚想编译一下FORTRAN版的MEX，没想到他们只支持 Intel Visual Fortran，
再下一个 Intel Visual Fortran，太麻烦了，
不试了，
赞一下 回复此楼
12楼2013-06-27 18:04:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

长路漫漫

新虫 (初入文坛)
引用回帖:
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

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

墙里墙外各不同
回复此楼

» 本帖已获得的红花（最新10朵）

长路漫漫
14楼2013-06-28 17:20:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

长路漫漫

新虫 (初入文坛)
送红花一朵
引用回帖:
14楼: Originally posted by somomo91 at 2013-06-28 17:20:39
墙里墙外各不同...


哦～～
谢谢你，我明白了～你真是个很nice的人哇～
我要送你一朵花，本来比较拮据的，可是还是要送你
表示感谢
赞一下 回复此楼
15楼2013-06-28 17:51:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

清风点兵

新虫 (小有名气)
请教您 http://muchong.com/t-11731506-1
赞一下 回复此楼
16楼2017-10-11 22:47:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 长路漫漫 的主题更新
162/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:43 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 11/550 2026-08-23 02:55 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 00:52 by OEbVnUOu01ol
[基金申请] 2026国自然函评费到账 +11 羊腰板 2026-08-21 11/550 2026-08-22 22:03 by lsbin3733
[论文投稿] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 3/150 2026-08-22 22:03 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 3/150 2026-08-22 22:03 by OEbVnUOu01ol
[教师之家] 跳槽后在研项目怎么办？ +3 简单化xn 2026-08-22 7/350 2026-08-22 19:46 by 简单化xn
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今天系统多次维护，明天很可能放榜！ +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
信息提示
关闭
请填处理意见
关闭