| 查看: 520 | 回复: 2 | ||
| 本帖产生 1 个 博学EPI ，点击这里进行查看 | ||
[求助]
fortran编程
|
||
| 求大神帮忙解答一下：fortran中如何将一个二维数组求行列式，然后将该行列式的值作为另一个数组的元素 |
回复此楼» 猜你喜欢
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有7人回复
-
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有7人回复
-
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有10人回复
-
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）
已经有9人回复
-
有院领导为了换新车，用横向课题经费买了俩车
已经有10人回复
-
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性
已经有6人回复
-
酰胺脱乙酰基
已经有13人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有10人回复
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