Télécharger la plaquette Cast3M


fr en

Télécharger perpen.eso

Retour à la liste

Numérotation des lignes : 

perpen

  1. C PERPEN    SOURCE    JB251061  22/04/01    21:16:01     11327          
  2.       SUBROUTINE PERPEN(entree, sortie)
  3. C ====================================================================
  4. C Sous-programme PERPEN : calcul d'un vecteur unitaire perpendiculaire
  5. C                         au vecteur fourni en entree en dimension 3
  6. C ====================================================================
  7.  
  8.       IMPLICIT INTEGER(I-N)
  9.       IMPLICIT REAL*8(A-H,O-Z)
  10.  
  11. -INC PPARAM
  12. -INC CCOPTIO
  13. -INC CCREEL
  14.  
  15.       real*8 entree(3), sortie(3), axe(3), essai(3),
  16.      $        norme_max, norme_essai
  17.  
  18. C     On construit une direction perpendiculaire par produit
  19. C     vectoriel avec la meilleure des 3 directions (1 0 0),
  20. C     (0 1 0) et (0 0 1)
  21.       norme_max = 0.D0
  22.       axe = 0.D0
  23.       do i=1,3
  24.         axe(i) = 1.D0
  25.         call provec(entree, axe, essai)
  26.         call norme(essai, norme_essai)
  27.         norme_essai = norme_essai + xpetit
  28.         if(norme_essai.gt.norme_max) then
  29.           sortie = essai
  30.           norme_max = norme_essai
  31.         endif
  32.         axe(i) = 0.D0
  33.       enddo
  34.       call normer(sortie)
  35.  
  36.       END
  37.  
  38.  

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