24小时热门版块排行榜    

查看: 1169  |  回复: 13
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

飞过海_92

金虫 (小有名气)


[交流] fortran 调用pardiso求解 两次结果不一样

如题,第一次:全部实数稀疏矩阵运算,mtype=11 ,实数非对称矩阵
第二次:将实数矩阵变成了复数矩阵,(实部等于原来的实部,虚部全部为零)复数稀疏矩阵运算,mtype=13 ,复数非对称矩阵

这个按理应该是没有影响的吧,
但解出来结果不一样,虽然复数的解出来虚部还是全部为0
回复此楼

» 猜你喜欢

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

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

jove1782

木虫 (正式写手)



飞过海_92(金币+1): 谢谢参与
6楼2016-11-08 20:44:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

wz_work

木虫 (小有名气)



飞过海_92(金币+1): 谢谢参与
你只给问题和自己的结论,看不到代码也没法说话啊。。。
8楼2016-11-09 10:00:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

飞过海_92

金虫 (小有名气)


引用回帖:
8楼: Originally posted by wz_work at 2016-11-09 10:00:53
你只给问题和自己的结论,看不到代码也没法说话啊。。。

找到原因了,复数定义complex(kind(1d0)),实数定义real*8,
将实数复制给复数实部是精度变了
9楼2016-11-09 15:51:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wz_work

木虫 (小有名气)


引用回帖:
9楼: Originally posted by 飞过海_92 at 2016-11-09 15:51:29
找到原因了,复数定义complex(kind(1d0)),实数定义real*8,
将实数复制给复数实部是精度变了...

找到问题就好咯
14楼2016-11-10 09:24:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
5047754905楼
2016-11-08 20:43   回复  
飞过海_92(金币+1): 谢谢参与
发自小木虫Android客户端
2016-11-08 20:32   回复  
飞过海_92(金币+1): 谢谢参与
祝福宝。 发自小木虫Android客户端
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见