Sources Esope | Cast3M

trrilf
C TRRILF    SOURCE    CB215821  17/11/30    21:17:17     9639           
      FUNCTION TRRILF(P1,P2,P3)
C     *****************************************************************
C     QUALITE DU TRIANGLE : RAYON DU CERCLE INSCRIT SUR ARETE LA PLUS
C     LONGUE.
C     RIL = SURFACE / (DEMI PERIMETRE * ARETE LA PLUS LONGUE)
      IMPLICIT INTEGER(I-N)
      REAL*8 P1(*),P2(*),P3(*)
C
      REAL*8 TRRILF
      REAL*8       XV(3),YV(3),S,D,DMAX
      INTEGER    I
C
      TRRILF = 0.0D0
      XV(1) = P2(1) - P1(1)
      YV(1) = P2(2) - P1(2)
      XV(2) = P3(1) - P2(1)
      YV(2) = P3(2) - P2(2)
      XV(3) = P1(1) - P3(1)
      YV(3) = P1(2) - P3(2)
      S = (XV(1) *  YV(2)) - ( XV(2) * YV(1) )
      IF( S.LT.0.0D0 )GOTO 999
      DMAX = 0.0D0
      DO 10 I=1,3
        D = XV(I)**2 + YV(I)**2
        IF( D .GT. DMAX )DMAX = D
   10 CONTINUE
      TRRILF = ( S / DMAX )
  999 END