Esope source files

   1 : C NORMV2    SOURCE    MAGN      09/07/30    21:19:42     6458           
   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 CCOPTIO
  34 : -INC SMMATRIK
  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 :