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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 逆矩阵运算
212/3返回列表上一页123下一页
查看: 2976  |  回复: 20
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

kimileegdut

捐助贵宾 (小有名气)
送红花一朵
引用回帖:
8楼: Originally posted by pippi6 at 2015-05-13 21:27:12
我试了一下，gaussr_inverse 应该工作。 逆矩阵和原矩阵的乘积与单位矩阵的最大差值为  2.442490654175344E-015。


subroutine  InverseMatrix
  implicit none
  integer, parameter :: nn=100,m=20
  inte ...

还有一个问题想请教您，就是我的主程序里面的数据都是单精度，但是子程序里面的计算涉及到双精度，我在调用子程序前需要把数据先变为双精度吗？然后计算出来的逆矩阵，在放回到主程序被应用的话，是不是要再变回单精度？这样会有问题吗？
赞一下 回复此楼
11楼2015-05-14 11:45:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobiao007

木虫 (职业作家)

中国特色

【答案】应助回帖

当然应该用老外编写的经过几十年考验的Lapack软件包来处理矩阵计算问题
赞一下 回复此楼
我同意叔本华的观点，人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活，逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
12楼2015-05-14 13:46:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kimileegdut

捐助贵宾 (小有名气)
引用回帖:
12楼: Originally posted by baobiao007 at 2015-05-14 13:46:14
当然应该用老外编写的经过几十年考验的Lapack软件包来处理矩阵计算问题

导师喜欢fortran,没办法。。。
赞一下 回复此楼
13楼2015-05-14 14:49:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobiao007

木虫 (职业作家)

中国特色
引用回帖:
13楼: Originally posted by kimileegdut at 2015-05-14 14:49:38
导师喜欢fortran,没办法。。。...

lapack就是fortran写的
赞一下 回复此楼
我同意叔本华的观点，人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活，逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
14楼2015-05-14 15:30:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询
引用回帖:
10楼: Originally posted by kimileegdut at 2015-05-14 10:47:47
你好，我还有个问题想请教一下您，就是上面计算b=a*c的时候，为什么不直接用矩阵相乘语句b=matmul(a,c),而使用了循环？...

无所谓，应该也可以的
赞一下 回复此楼
15楼2015-05-14 21:58:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询
引用回帖:
9楼: Originally posted by kimileegdut at 2015-05-14 09:48:38
好的，很感谢你啊！可能是我调用子程序的时候形参和实参没有对应上！gaussr_inverse 的形参表各个变量是代表什么？...

! How to use
!  real*8 :: c(n_dim,n_dim),e(n_dim,n_dim)
!  integer :: n=4
!     n           size of the matrix to be inversed
!     n_dim       array size
!     e           unit matrix
!     c           matrix to be inversed
!     call gaussr_inverse(c,n,n_dim,e,n,n_dim)
赞一下 回复此楼
16楼2015-05-14 21:59:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kimileegdut

捐助贵宾 (小有名气)
引用回帖:
15楼: Originally posted by pippi6 at 2015-05-14 21:58:42
无所谓，应该也可以的...

ok!还有一个问题想请教您，就是我的主程序里面的数据都是单精度，但是子程序里面的计算涉及到双精度，我在调用子程序前需要把数据先变为双精度吗？然后计算出来的逆矩阵，在放回到主程序被应用的话，是不是要再变回单精度？这样会有问题吗？
赞一下 回复此楼
17楼2015-05-14 23:08:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询
引用回帖:
17楼: Originally posted by kimileegdut at 2015-05-14 23:08:17
ok!还有一个问题想请教您，就是我的主程序里面的数据都是单精度，但是子程序里面的计算涉及到双精度，我在调用子程序前需要把数据先变为双精度吗？然后计算出来的逆矩阵，在放回到主程序被应用的话，是不是要再变回 ...

其实，像你自己所意识到的一样，就是形参要对应上，包括数据类型。如果子程序要求双精度，那么在主程序里调用时也需要提供双精度变量。双精度和单精度在主程序里可以互相赋值没问题，只是从双精度到单精度损失精度而已。如果你不在乎，就没有关系。
赞一下 回复此楼
18楼2015-05-15 06:58:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kimileegdut

捐助贵宾 (小有名气)
引用回帖:
8楼: Originally posted by pippi6 at 2015-05-13 21:27:12
我试了一下，gaussr_inverse 应该工作。 逆矩阵和原矩阵的乘积与单位矩阵的最大差值为  2.442490654175344E-015。


subroutine  InverseMatrix
  implicit none
  integer, parameter :: nn=100,m=20
  inte ...

你好，我还有一个不理解的地方，就是在调用子程序计算逆矩阵前，为什么要先把原来的矩阵元素值赋给b?如果没有b=a，这一句，调用子程序时实参用a，计算出来的结果有误，这什么原因导致的？
赞一下 回复此楼
19楼2015-05-15 10:54:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询
引用回帖:
19楼: Originally posted by kimileegdut at 2015-05-15 10:54:50
你好，我还有一个不理解的地方，就是在调用子程序计算逆矩阵前，为什么要先把原来的矩阵元素值赋给b?如果没有b=a，这一句，调用子程序时实参用a，计算出来的结果有误，这什么原因导致的？...

因为原矩阵a会被破坏
赞一下 回复此楼
20楼2015-05-15 11:44:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kimileegdut 的主题更新
212/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今年审到国自然15份，谈谈感受 +16 国自然国社科中 2026-05-17 16/800 2026-05-18 14:58 by gy116024
[文学芳草园] 半夜喝咖啡 +3 myrtle 2026-05-15 5/250 2026-05-18 01:03 by 小沈2018
[找工作] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 ky2p12rrjj 2026-05-15 4/200 2026-05-17 19:47 by Equinoxhua
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 xx7gd5zq4e 2026-05-15 6/300 2026-05-17 19:36 by Equinoxhua
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 cjf4bx70cj 2026-05-14 7/350 2026-05-17 18:49 by Equinoxhua
[公派出国] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 5/250 2026-05-17 18:39 by Equinoxhua
[考博] 找博士生导师 +4 小代想上岸 2026-05-15 5/250 2026-05-17 17:44 by 脱颖而出
[基金申请] 这年头没有找到涵评专家，还有中面上的可能吗 +13 dd921ww 2026-05-12 16/800 2026-05-17 12:38 by 于轩
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 08:11 by 11n4dfd8yn
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 7/350 2026-05-17 06:55 by 11n4dfd8yn
[博后之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 k37jurhrau 2026-05-16 4/200 2026-05-17 01:35 by ue3ir18jc3
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 3/150 2026-05-17 00:45 by ue3ir18jc3
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:25 by ue3ir18jc3
[高分子] 本人最近太闲了，谁有问题可以提，每天会统一回复 +9 一切都是空工 2026-05-12 20/1000 2026-05-16 19:52 by Equinoxhua
[有机交流] 求助2，4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-16 13:57 by vcdazktkjx
[文学芳草园] 风把牡丹吹跑了 +5 myrtle 2026-05-12 9/450 2026-05-15 15:27 by myrtle
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
信息提示
关闭
请填处理意见
关闭