Sources Esope | Cast3M

sppoin
C SPPOIN    SOURCE    CHAT      05/01/13    03:24:08     5004
      INTEGER FUNCTION  SPPOIN(IDIMC,U,POINT,SPHERE,ZERO)
C     **********************************************************************
C     OBJET : LE POINT EST-IL DANS LA SPHERE ?
C     EN ENTREE:
C       IDIMC  : DIMENSION DE L'ESPACE
C       POINT : COORDONNEES DU POINT A TESTER
C       SPHERE : VECTEUR DIAMETRE DE LA SPHERE (CF SPCIRC)
C       U     : LE POINT DE LA SPHERE QUI A SERVIT A SON CALCUL(CF SPCIRC)
C     EN SORTIE : 1 SI "POINT" EST DANS "SPHERE", O SINON
C     **********************************************************************
      IMPLICIT INTEGER(I-N)
      INTEGER IDIMC
      REAL*8    U(*),POINT(*),SPHERE(*),ZERO
C
      REAL*8 V,FAC,S
      INTEGER I
C      DATA FAC /.999D0/
      DATA FAC /.999999999D0/
C
      S = 0.0D0
      SPPOIN=0
      DO 10 I=1,IDIMC
        V = POINT(I)-U(I)
        S = S + V * ( (FAC*SPHERE(I)) - V )
   10 CONTINUE
      IF( S.LT.ZERO )RETURN
C     --- DANS SPHERE ---
      SPPOIN=1
      END