| 查看: 581 | 回复: 2 | ||
| 本帖产生 1 个 博学EPI ，点击这里进行查看 | ||
[求助]
fortran编程
|
||
| 求大神帮忙解答一下：fortran中如何将一个二维数组求行列式，然后将该行列式的值作为另一个数组的元素 |
回复此楼» 猜你喜欢
为什么网上很多人说本周 12号出结果
已经有4人回复
-
基金中了
已经有15人回复
-
帮忙看看fileCode
已经有8人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有3人回复
-
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
静等基金结果
已经有16人回复
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
100