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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » fortran 调用pardiso求解 两次结果不一样
141/1返回列表
查看: 1024  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

飞过海_92

金虫 (小有名气)

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

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

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

» 猜你喜欢

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2016-11-08 20:13:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jove1782

木虫 (正式写手)
★
飞过海_92(金币+1): 谢谢参与
厉害

发自小木虫Android客户端
回复此楼
6楼2016-11-08 20:44:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖
简单回复
tan9142楼
2016-11-08 20:31   回复  
飞过海_92(金币+1): 谢谢参与
zhchzhsh20763楼
2016-11-08 20:32   回复  
飞过海_92(金币+1): 谢谢参与
祝福宝。 发自小木虫Android客户端
tzynew4楼
2016-11-08 20:39   回复  
飞过海_92(金币+1): 谢谢参与
5047754905楼
2016-11-08 20:43   回复  
飞过海_92(金币+1): 谢谢参与
发自小木虫Android客户端
dsctg7楼
2016-11-08 21:14   回复  
飞过海_92(金币+1): 谢谢参与
syhorchid10楼
2016-11-10 00:58   回复  
飞过海_92(金币+1): 谢谢参与
JK4511楼
2016-11-10 04:38   回复  
飞过海_92(金币+1): 谢谢参与
yjj198311212楼
2016-11-10 08:21   回复  
飞过海_92(金币+1): 谢谢参与
xhmaohan13楼
2016-11-10 08:25   回复  
相关版块跳转 我要订阅楼主 飞过海_92 的主题更新
141/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭