[求助]如何编写?-----写出各三角形三个顶点标号.[Fortran中文网]

楼主

[求助]如何编写?-----写出各三角形三个顶点标号.

此主题相关图片如下：

图形如上所示.
      结果应该是:1 7 8
8 2 1
8 9 2
2 9 3
3 9 10
10 4 3
10 11 4
4 11 5
5 11 12
12 6  5
刚学fortran ,希望各位高手指点，谢谢！！！

第 2 楼

用户已被锁定

第 3 楼

PROGRAM IJM
COMMON /SM/ II(20),JJ(20),KK(20)
NE=12

CALL B1(1,10,1,7)

OPEN (8,FILE='IJM.DAT',FORM='FORMATTED',STATUS='NEW')
WRITE (8,90) NE
WRITE(8,90)(II(I),I=1,NE)
WRITE(8,90)(JJ(I),I=1,NE)
WRITE(8,90)(KK(I),I=1,NE)
90 FORMAT(1X,10I6)
END

SUBROUTINE B1(I1,I2,L1,L2)
COMMON /SM/ II(20),JJ(20),KK(20)
L11=L1
L22=L2
DO 4 I=I1,I2,2
II(I)=L11
JJ(I)=L22
KK(I)=L11+1
II(I+1)=L11+1
JJ(I+1)=L22
KK(I+1)=L22+1
L11=L11+1
4 L22=L22+1
END

这个倒是可以,但不是三个一对显示的.还需自己找对应关系

第 4 楼

鄙人现在正学有限元，网格划分时也遇见这个问题，要把各单元对应的节点表示出来，还要逆时针的，好麻烦啊，有没有高手指教一下

第 5 楼

用支付宝给pingjianming付款或购买其商品，支付宝交易免手续费、安全、快捷！

说一下我的看法,如果有了上面的图形,也就是说已经知道了每个三角形单元的三个节点号了,现在只是要进行排序而已. 不知道我的理解是否正确,如果是.
下一步,设置两个矢量,分别为 xy1，xy2，xy1（1）=x（2）-x（1），xy1（2）=y（2）-y（1），xy2（1）=x（3）-x（2），xy2（2）=y（3）-y（2），
然后计算xy1（1）*xy2（2）-xy1（2）*xy2（1），结果为正，节点顺序不变，为负责交换点2 和点3顺序即可满足要求。对三维情况，稍复杂

[此贴子已经被作者于2006-4-8 8:38:25编辑过]