| 查看: 494 | 回复: 2 | ||
| 本帖产生 1 个 博学EPI ，点击这里进行查看 | ||
[求助]
fortran编程
|
||
| 求大神帮忙解答一下：fortran中如何将一个二维数组求行列式，然后将该行列式的值作为另一个数组的元素 |
回复此楼» 猜你喜欢
垃圾破二本职称评审标准
已经有19人回复
-
职称评审没过，求安慰
已经有53人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有5人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
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