Esope source files

   1 : C VPAST     SOURCE    CHAT      05/01/13    04:09:46     5004           
   2 :       SUBROUTINE VPAST(VCORE,B)                                                 
   3 : C=======================================================================        
   4 : C                                                                               
   5 : C     CONSTRUCTION DE LA MATRICE DE PASSAGE                                     
   6 : C     FAISANT PASSER DES COORDONNEES GLOBALES AUX COORDONNEES                   
   7 : C     LOCALES POUR UN TRIANGLE OU UN QUADRANGLE( DANS CE CAS ON PREND           
   8 : C       LE PLAN S APPUYANT SUR LES 3 PREMIERS NOEUDS                            
   9 : C ENTREES                                                                       
  10 : C    VCORE(3,3)=COORD DU TRIANGLE OU DU QUADRANGLE (1 COL = 1 PT)               
  11 : C SORTIES                                                                       
  12 : C    B(3,3)    =MATRICE DE PASSAGE                                              
  13 : C=======================================================================        
  14 :       IMPLICIT INTEGER(I-N)
  15 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  16 :       DIMENSION VCORE(3,*),B(3,*)                                               
  17 : C                                                                               
  18 : C                                                                               
  19 :       X21=VCORE(1,2)-VCORE(1,1)                                                 
  20 :       Y21=VCORE(2,2)-VCORE(2,1)                                                 
  21 :       Z21=VCORE(3,2)-VCORE(3,1)                                                 
  22 :       X31=VCORE(1,3)-VCORE(1,1)                                                 
  23 :       Y31=VCORE(2,3)-VCORE(2,1)                                                 
  24 :       Z31=VCORE(3,3)-VCORE(3,1)                                                 
  25 :       VL12=SQRT(X21*X21+Y21*Y21+Z21*Z21)                                        
  26 :       A1=Y21*Z31-Y31*Z21                                                        
  27 :       A2=X31*Z21-X21*Z31                                                        
  28 :       A3=X21*Y31-X31*Y21                                                        
  29 :       TWOA=SQRT(A1*A1+A2*A2+A3*A3)                                              
  30 : C                                                                               
  31 :       VL12=1.D0/VL12                                                            
  32 :       TWOA=1.D0/TWOA                                                            
  33 :       B(1,1)=X21*VL12                                                           
  34 :       B(1,2)=Y21*VL12                                                           
  35 :       B(1,3)=Z21*VL12                                                           
  36 :       B(3,1)=A1*TWOA                                                            
  37 :       B(3,2)=A2*TWOA                                                            
  38 :       B(3,3)=A3*TWOA                                                            
  39 :       B(2,1)=B(3,2)*B(1,3)-B(3,3)*B(1,2)                                        
  40 :       B(2,2)=B(3,3)*B(1,1)-B(1,3)*B(3,1)                                        
  41 :       B(2,3)=B(3,1)*B(1,2)-B(3,2)*B(1,1)                                        
  42 :       RETURN                                                                    
  43 :       END                                                                       
  44 :