C KAFINT SOURCE CHAT 05/01/13 00:50:28 5004 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=DMIN1(Z1,Z2) ZMAX=DMAX1(Z1,Z2) RMAX=DMAX1(R1,R2) IF(R.GE.RMAX.OR.Z.LE.ZMIN.OR.Z.GE.ZMAX) THEN A=2. ELSEIF(DABS(Z-Z1).LE.DMIN.OR.DABS(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