Esope source files

   1 : C TRANSB    SOURCE    PV        06/04/11    21:16:17     5396           
   2 : C    SOUS PROGRAMME DE VERIFICATION QU'UN SEGMENT POTENTIEL (LL1 LL2)           
   3 : C    EST DU BON COTE DU CONTOUR (LL3 LL1 LL4)                                   
   4 : C                                                                               
   5 :       SUBROUTINE TRANSB(X,IRECL,LL1,LL2,LL3,LL4)                                
   6 :       IMPLICIT INTEGER(I-N)
   7 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
   8 :       DIMENSION X(3,1)                                                          
   9 :       IRECL=1                                                                   
  10 :       X41=X(1,LL1)-X(1,LL4)                                                     
  11 :       Y41=X(2,LL1)-X(2,LL4)                                                     
  12 :       X31=X(1,LL1)-X(1,LL3)                                                     
  13 :       Y31=X(2,LL1)-X(2,LL3)                                                     
  14 :       UV=X41*Y31-X31*Y41                                                        
  15 :       X21=X(1,LL1)-X(1,LL2)                                                     
  16 :       Y21=X(2,LL1)-X(2,LL2)                                                     
  17 :       VX=X41*Y21-Y41*X21                                                        
  18 :       UX=X21*Y31-Y21*X31                                                        
  19 :       IF (UV.GE.0.D0.AND.VX.GT.0.D0.AND.UX.GT.0.D0) IRECL=0                           
  20 :       IF (UV.LT.0.D0.AND.(VX.GT.0.D0.OR.UX.GT.0.D0)) IRECL=0                          
  21 :       RETURN                                                                    
  22 :       END                                                                       
  23 :  
  24 :