Télécharger d2def.eso

Retour à la liste

Numérotation des lignes :

  1. C D2DEF SOURCE CHAT 06/03/29 21:18:09 5360
  2. SUBROUTINE D2DEF(XP1,XP2,XP3,XPC,VDIA,IDIMC,
  3. > COEF,TS,iarr)
  4. C *****************************************************************
  5. C OBJET : TAILLE SOUHAITE PAR DEFAUT
  6. C (LA TAILLE SOUHAITE EST EVALUEE AU CENTRE D'UNE BOULE)
  7. C EN ENTREE :
  8. C XP1,XP2,XP3 : LES TROIS POINT DU TRIANGLE
  9. C XPC : UN POINT SUR LE CERCLE (INUTILISE)
  10. C VDIA : LE VECTEUR DIAMETRE PARTANT DE XPC
  11. C IDIMC : DIMENSION DE L'ESPACE
  12. C EN SORTIE :
  13. C TS : TAILLE SOUHAITE POUR LE TRIANGLE
  14. C = LONGUEUR DE SA PLUS PETITE ARETE
  15. C COEF : TS / RC
  16. C RC EST LE RAYON DU CERCLE CIRCONSCRIT
  17. C PLUS COEF EST PETIT PLUS ON RAFFINE
  18. C iarr : CODE D'ERREUR 0 SI OK,
  19. C -1 SI LE RAYON CIRCONSCRIT EST NUL
  20. C NIVEAU : FICHIER
  21. C *****************************************************************
  22. IMPLICIT INTEGER(I-N)
  23. INTEGER IDIMC
  24. REAL*8 XP1(*),XP2(*),XP3(*)
  25. REAL*8 XPC(*),VDIA(*),COEF,TS
  26. INTEGER iarr
  27. C
  28. REAL*8 V(3),D,DMIN,DIAM2,RC,XYZEPS
  29. EXTERNAL XNORVE
  30. REAL*8 XNORVE
  31. XYZEPS=1.D-10
  32. C
  33. C ============================
  34. C ------ ARETE LA PLUS COURTE ------
  35. C ============================
  36. CALL DIFFVE(XP2,XP1,IDIMC,V)
  37. DMIN = XNORVE(V,IDIMC)
  38. CALL DIFFVE(XP3,XP2,IDIMC,V)
  39. D = XNORVE(V,IDIMC)
  40. DMIN = MIN( D, DMIN )
  41. CALL DIFFVE(XP1,XP3,IDIMC,V)
  42. D = XNORVE(V,IDIMC)
  43. DMIN = MIN( D, DMIN )
  44. C ============================
  45. C ------ SUR LE RAYON CIRCONSCRIT ------
  46. C ============================
  47. DIAM2 = VDIA(1)**2 + VDIA(2)**2
  48. RC = SQRT( DIAM2 ) / 2.0D0
  49. IF( RC.GT. -XYZEPS .AND.RC .LT. XYZEPS )GOTO 999
  50. COEF = DMIN / RC
  51. TS = DMIN
  52. iarr = 0
  53. 999 END
  54.  
  55.  
  56.  
  57.  
  58.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales