| 查看: 620 | 回复: 2 | |||
| 本帖产生 1 个 博学EPI ，点击这里进行查看 | |||
[求助]
fortran编程
|
|||
| 求大神帮忙解答一下：fortran中如何将一个二维数组求行列式，然后将该行列式的值作为另一个数组的元素 |
回复此楼» 猜你喜欢
现代”学阀”该如何界定
已经有12人回复
-
国社科系统bug了，是不是要放榜了？
已经有10人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有3人回复
-
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有11人回复
-
申博发邮件
已经有11人回复
-
我的奶奶
已经有3人回复
2楼2016-10-31 23:46:00
|
do i=1,nele,1 PP(1:3)=(/nodx(nod(i,1)), nody(nod(i,1)),nodz(nod(i,1)) /) QQ(1:3)=(/nodx(nod(i,2)), nody(nod(i,2)),nodz(nod(i,2)) /) RR(1:3)=(/nodx(nod(i,3)), nody(nod(i,3)),nodz(nod(i,3)) /) SS(1:3)=(/nodx(nod(i,4)), nody(nod(i,4)),nodz(nod(i,4)) /) a1(1:3,1:3)=(/QQ(1),QQ(2),QQ(3),RR(1),RR(2),RR(3),SS(1),SS(2),SS(3)/) b1(1:3,1:3)=-(/1.0,QQ(2),QQ(3),1.0,RR(2),RR(3),1.0,SS(2),SS(3)/) c1(1:3,1:3)=(/1.0,QQ(1),QQ(3),1.0,RR(1),RR(3),1.0,SS(1),SS(3)/) d1(1:3,1:3)=-(/1.0,QQ(1),QQ(2),1.0,RR(1),RR(2),1.0,SS(1),SS(2)/) aae1(i)=det(a1(1:3,1:3)) bbe1(i)=det(b1(1:3,1:3)) cce1(i)=det(c1(1:3,1:3)) dde1(i)=det(d1(1:3,1:3)) end do （这是我程序中的一部分，数组都以正确定义，但运行时a1,b1,c1,d1部分总会出现 Error: The shapes of the array expressions do not conform. 是什么原因啊？） |
3楼2016-11-01 09:15:53