Esope source files

   1 : C TRJ403    SOURCE    CHAT      05/01/13    03:48:11     5004           
   2 :       SUBROUTINE TRJ403 (XNOEU,X,G)
   3 : C                                                                               
   4 : C                                                                               
   5 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC          
   6 : C                                                                               
   7 : C     CALCUL DES COORDONNEES DE REFERENCE DANS UN TET4
   8 : C     XNOEU COORDONNEES DES SOMMETS                                             
   9 : C     X                 REELLES
  10 : C     G                 DE REFERENCE
  11 : C                                                                               
  12 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC         
  13 : C   si X1 X2 X3 X4 sont les noeuds du tetraedre et L1 L2 L3 L4 les
  14 : C  coordonnées barycentriques d' un point M . L2 L3 L4 sont egaux aux
  15 : C  coordonnées de référence de M. On va donc resoudre le systeme:
  16 : C        ->      -->      -->      -->      -->
  17 : C      | OM = L1 OX1 + L2 OX2 + L3 OX3 + L4 OX4
  18 : C      |  1 = L1 + L2 + L3 + L4
  19 : C                                                                               
  20 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  21 : C                                                                               
  22 :       IMPLICIT INTEGER(I-N)
  23 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  24 : C                                                                               
  25 :       DIMENSION XNOEU(3,4),X(3),G(3)
  26 :       A1= XNOEU(1,2)-XNOEU(1,1)
  27 :       A2= XNOEU(1,3)-XNOEU(1,1)
  28 :       A3= XNOEU(1,4)-XNOEU(1,1)
  29 :       A4= XNOEU(2,2)-XNOEU(2,1)
  30 :       A5= XNOEU(2,3)-XNOEU(2,1)
  31 :       A6= XNOEU(2,4)-XNOEU(2,1)
  32 :       A7= XNOEU(3,2)-XNOEU(3,1)
  33 :       A8= XNOEU(3,3)-XNOEU(3,1)
  34 :       A9= XNOEU(3,4)-XNOEU(3,1)
  35 :       B1= X(1)-XNOEU(1,1)
  36 :       B2= X(2)-XNOEU(2,1)
  37 :       B3= X(3)-XNOEU(3,1)
  38 :       DETA=A1*A5*A9+A2*A6*A7+A4*A8*A3-A7*A5*A3-A8*A6*A1-A9*A4*A2
  39 :       G(1)=B1*A5*A9+A2*A6*B3+B2*A8*A3-B3*A5*A3-A8*A6*B1-A9*B2*A2
  40 :       G(2)=A1*B2*A9+B1*A6*A7+A4*B3*A3-A7*B2*A3-B3*A6*A1-A9*A4*B1
  41 :       G(3)=A1*A5*B3+A2*B2*A7+A4*A8*B1-A7*A5*B1-A8*B2*A1-B3*A4*A2
  42 :       G(1)=G(1)/DETA
  43 :       G(2)=G(2)/DETA
  44 :       G(3)=G(3)/DETA
  45 :       RETURN                                                                    
  46 :       END                                                                       
  47 :  
  48 :