Esope source files

   1 : C KCHREP    SOURCE    CHAT      05/01/13    00:52:50     5004           
   2 :       SUBROUTINE  KCHREP(NES,NSOM,O,A,U,AH,UH,H)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)
   5 :       DIMENSION O(NES),A(NES,NSOM),U(NES+1)
   6 :       DIMENSION AH(NES,NSOM),UH(NES+1)
   7 :       DIMENSION H(NES,NES),B(3,3)
   8 : C
   9 : C Calcul des facteurs de forme en 3D
  10 : C     calcul des coordonnées d'un point A dans le repere local du HC
  11 : C
  12 : C     NES  : dim espace
  13 : C     NSOM : nombre de sommets
  14 : C     O    : point de vue
  15 : C     A    : sommets du triangle
  16 : C     U    : vecteur normal unitaire
  17 : C     AH,UH : idem repere local du HC
  18 : C     H    : matrice de passage durepere HC au repere global
  19 : C
  20 : 
  21 :       DO 2 L=1,NSOM
  22 :       DO 20 I=1,NES
  23 :       B(I,L)=A(I,L)-O(I)
  24 :  20   CONTINUE
  25 :   2   CONTINUE
  26 : 
  27 :       DO 1 L=1,NSOM
  28 :       AH(1,L) = B(1,L)*H(1,1) + B(2,L)*H(1,2) + B(3,L)*H(1,3)
  29 :       AH(2,L) = B(1,L)*H(2,1) + B(2,L)*H(2,2) + B(3,L)*H(2,3)
  30 :       AH(3,L) = B(1,L)*H(3,1) + B(2,L)*H(3,2) + B(3,L)*H(3,3)
  31 :   1   CONTINUE
  32 : 
  33 :       UH(1) = U(1)*H(1,1) + U(2)*H(1,2) + U(3)*H(1,3)
  34 :       UH(2) = U(1)*H(2,1) + U(2)*H(2,2) + U(3)*H(2,3)
  35 :       UH(3) = U(1)*H(3,1) + U(2)*H(3,2) + U(3)*H(3,3)
  36 : 
  37 :       UH(4) = -(UH(1)*AH(1,1) + UH(2)*AH(2,1) + UH(3)*AH(3,1))
  38 : 
  39 :       RETURN
  40 :       END
  41 :  
  42 :