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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 逆矩阵运算
北京石油化工学院2026年研究生招生接收调剂公告
212/3返回列表上一页123下一页
查看: 2882  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 11408 321分求调剂 +3 huchun12138 2026-03-30 4/200 2026-04-01 22:48 by guanxin1001
[考研] 266求调剂 +7 阳阳哇塞 2026-04-01 7/350 2026-04-01 22:27 by barlinike
[考研] 261求调剂 +3 明仔· 2026-04-01 3/150 2026-04-01 20:52 by cq2548
[考研] 353求调剂 +4 拉钩不许变 2026-04-01 4/200 2026-04-01 18:10 by 记事本2026
[考研] 311求调剂 +11 蓝月亮亮 2026-03-30 11/550 2026-04-01 16:33 by 七度不信任
[考研] 一志愿郑大085600，310分求调剂 +6 李潇可 2026-03-26 6/300 2026-04-01 14:44 by chenqifeng666
[考研] 求调剂 +5 零八# 2026-03-27 5/250 2026-04-01 14:40 by yulian1987
[考研] 283求调剂 +9 A child 2026-03-28 9/450 2026-04-01 14:20 by Jaylen.
[考研] 一志愿同济大学323分（080500）求调剂 +4 yikeniu 2026-04-01 4/200 2026-04-01 14:06 by asdfzly
[考研] 086000生物与医药298调剂求助 +4 元元青青 2026-03-31 6/300 2026-04-01 11:13 by syh9288
[考研] 一志愿：西北大学，英一数一408-284分求调剂 +7 12.27 2026-03-27 7/350 2026-03-31 21:59 by lbsjt
[考研] 354求调剂 +3 lxb598 2026-03-31 4/200 2026-03-31 13:42 by sophie2180
[考研] 085600，专业课化工原理，320分求调剂 +6 大馋小子 2026-03-29 6/300 2026-03-31 10:03 by 氯化亚硝酰
[考研] 调剂 +4 GK72 2026-03-30 4/200 2026-03-30 20:32 by dick_runner
[考研] 求调剂 +10 家佳佳佳佳佳 2026-03-29 10/500 2026-03-30 18:34 by 544594351
[考研] 297求调剂 +17 田洪有 2026-03-26 18/900 2026-03-30 18:32 by nothing投稿中
[考研] 一志愿南开大学0710生物学359求调剂 +5 兔兔兔111223314 2026-03-29 7/350 2026-03-30 18:29 by 兔兔兔111223314
[考研] 本科新能源科学与工程，一志愿华理能动285求调剂 +3 AZMK 2026-03-27 5/250 2026-03-28 16:19 by xxxsssccc
[考研] 304求调剂 +6 曼殊2266 2026-03-27 6/300 2026-03-28 14:10 by 唐沐儿
[考研] 调剂 +4 柚柚yoyo 2026-03-26 4/200 2026-03-26 20:43 by fmesaito
信息提示
关闭
请填处理意见
关闭