Télécharger normv2.eso

Retour à la liste

Numérotation des lignes :

normv2
  1. C NORMV2 SOURCE PV 20/09/26 21:19:03 10724
  2. SUBROUTINE NORMV2(VECT,NORMV,
  3. $ IMPR,IRET)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. IMPLICIT INTEGER (I-N)
  6. C***********************************************************************
  7. C NOM : NORMV2
  8. C DESCRIPTION : Equilibrage d'un vecteur dual
  9. C <=> Déséquilibrage d'un vecteur primal
  10. C
  11. C
  12. C
  13. C LANGAGE : ESOPE
  14. C AUTEUR : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  15. C mél : gounand@semt2.smts.cea.fr
  16. C********************************************************
  17. C***********************************************************************
  18. C SYNTAXE GIBIANE :
  19. C ENTREES :
  20. C ENTREES/SORTIES :
  21. C SORTIES :
  22. C CODE RETOUR (IRET) : = 0 si tout s'est bien passé
  23. C***********************************************************************
  24. C VERSION : v1, 29/04/2003, version initiale
  25. C HISTORIQUE : v1, 29/04/2003, création
  26. C HISTORIQUE :
  27. C HISTORIQUE :
  28. C***********************************************************************
  29. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  30. C en cas de modification de ce sous-programme afin de faciliter
  31. C la maintenance !
  32. C***********************************************************************
  33. -INC PPARAM
  34. -INC CCOPTIO
  35. POINTEUR VECT.IZA
  36. POINTEUR NORMV.IZA
  37. *
  38. INTEGER IMPR,IRET
  39. *
  40. * Executable statements
  41. *
  42. IF (IMPR.GT.5) WRITE(IOIMP,*) 'Entrée dans normv2.eso'
  43. SEGACT NORMV
  44. NTTDDL=NORMV.A(/1)
  45. SEGACT VECT*MOD
  46. NTTDD2=VECT.A(/1)
  47. IF (NTTDD2.NE.NTTDDL) THEN
  48. WRITE(IOIMP,*) 'Erreur grave'
  49. GOTO 9999
  50. ENDIF
  51. DO ITTDDL=1,NTTDDL
  52. VAL=VECT.A(ITTDDL)
  53. VAL=VAL/SQRT(NORMV.A(ITTDDL))
  54. VECT.A(ITTDDL)=VAL
  55. ENDDO
  56. SEGDES VECT
  57. SEGDES NORMV
  58. *
  59. * Normal termination
  60. *
  61. IRET=0
  62. RETURN
  63. *
  64. * Format handling
  65. *
  66. *
  67. * Error handling
  68. *
  69. 9999 CONTINUE
  70. IRET=1
  71. WRITE(IOIMP,*) 'An error was detected in subroutine normv2'
  72. RETURN
  73. *
  74. * End of subroutine NORMV2
  75. *
  76. END
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  

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