24小时热门版块排行榜    

查看: 1583  |  回复: 10

goaliman

木虫 (著名写手)

[求助] 关于Fortran中的乘法和求逆比较

对于矩阵乘法,可以用A=matmul(B,C),也可以用A=B.x.C,这两个方法的计算结果是否一致?

对于矩阵求逆,可以用A=.i.B,也可以用ZGETRF,这两个方法的计算结果是否一致?

以上A,B和C均为复数矩阵,谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

goaliman

木虫 (著名写手)

引用回帖:
6楼: Originally posted by sunkingxu at 2012-12-17 15:11:21
我用得太多了.x.,.i.这类的,方便简明。对于一般的复数矩阵没有什么区别。

恩,但是我昨天做了一个比较,跟lapack里面的getrf和getri,还有差别的,后者更准确
7楼2012-12-17 16:39:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
Fortran 现在这么高级了啊,.x. 和 .i. 根本就没有见过呢……
2楼2012-12-16 12:20:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
goaliman: 金币+10, 2 2012-12-17 16:39:39
你自己测试一下不就好了吗?不过,依我的直觉,你的 .x. 和 .i. 应该是某个编译器的扩展,不具备通用性,能不用就不用吧……
3楼2012-12-16 12:21:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goaliman

木虫 (著名写手)

.x.和.i.是imsl里面的命令,比较简明
4楼2012-12-16 12:33:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主

楼主自己找几个极端的矩阵例子试试不就知道了?
5楼2012-12-16 13:50:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunkingxu

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
goaliman: 金币+10 2012-12-17 16:41:02
我用得太多了.x.,.i.这类的,方便简明。对于一般的复数矩阵没有什么区别。
6楼2012-12-17 15:11:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goaliman

木虫 (著名写手)

引用回帖:
5楼: Originally posted by nono2009 at 2012-12-16 13:50:09
楼主自己找几个极端的矩阵例子试试不就知道了?

你把剩下的10分转过去吧,我没权限给你转
8楼2012-12-17 16:41:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主

引用回帖:
8楼: Originally posted by goaliman at 2012-12-17 16:41:53
你把剩下的10分转过去吧,我没权限给你转...

谢谢。可我也转不了。
9楼2012-12-17 16:45:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goaliman

木虫 (著名写手)

引用回帖:
9楼: Originally posted by nono2009 at 2012-12-17 16:45:51
谢谢。可我也转不了。...

你回帖应助,我就可以给你评分了
10楼2012-12-17 17:18:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 goaliman 的主题更新
信息提示
请填处理意见