Sources Esope | Cast3M

kafint
C KAFINT    SOURCE    CB215821  17/11/30    21:16:32     9639           
      SUBROUTINE KAFINT(R1,Z1,R2,Z2,R,Z,A)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C----------------------------------------------------------------------
C SP appele par KAINTE
C     appartenance du point (R,Z) au triangle defini par les
C     points (R1,Z1) et (R2,Z2)
C     resultat A : 2 si exterieur
C                : represente un cosinus sinon
C----------------------------------------------------------------------
 
      DMIN=1.D-3
      ZMIN=MIN(Z1,Z2)
      ZMAX= MAX(Z1,Z2)
      RMAX= MAX(R1,R2)
 
      IF(R.GE.RMAX.OR.Z.LE.ZMIN.OR.Z.GE.ZMAX) THEN
       A=2.
 
      ELSEIF(ABS(Z-Z1).LE.DMIN.OR.ABS(Z2-Z).LE.DMIN) THEN
       A=2.
      ELSE
       A=R*R*(Z2-Z1)*(Z2-Z1)-R2*R2*(Z-Z1)*(Z-Z1)-R1*R1*(Z2-Z)*(Z2-Z)
       A=A/2./R1/R2/(Z-Z1)/(Z2-Z)
      ENDIF
      RETURN
      END