Esope source files

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