Sources Esope | Cast3M

ksubdv
C KSUBDV    SOURCE    CHAT      05/01/13    01:09:11     5004
      SUBROUTINE KSUBDV(A2,AA,NSUB)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C
C Calcul des facteurs de forme en 3D
C Sp appele par KPROJF
C     Subdivision d'un triangle en 4 triangles
C
      DIMENSION A1(3,3)
      DIMENSION A2(3,3),AA(3,3,4)
 
      NSUB = 4
C
C     SOMMETS
C
      DO 10 K=1,3
      A1(K,1) = (A2(K,1)+A2(K,2))/2.
      A1(K,2) = (A2(K,2)+A2(K,3))/2.
      A1(K,3) = (A2(K,3)+A2(K,1))/2.
 10   CONTINUE
C
C     4 ELEMENTS
C
      DO 11 K =1,3
      AA(K,1,1) = A2(K,1)
      AA(K,2,1) = A1(K,1)
      AA(K,3,1) = A1(K,3)
 11   CONTINUE
C
      DO 12 K =1,3
      AA(K,1,2) = A2(K,2)
      AA(K,2,2) = A1(K,2)
      AA(K,3,2) = A1(K,1)
 12   CONTINUE
 
      DO 13 K =1,3
      AA(K,1,3) = A2(K,3)
      AA(K,2,3) = A1(K,3)
      AA(K,3,3) = A1(K,2)
 13   CONTINUE
 
      DO 14 K =1,3
      AA(K,1,4) = A1(K,1)
      AA(K,2,4) = A1(K,2)
      AA(K,3,4) = A1(K,3)
 14   CONTINUE
 
      RETURN
      END