Esope source files

   1 : C PSDIAG    SOURCE    GOUNAND   06/03/06    21:19:20     5319           
   2 :       SUBROUTINE PSDIAG(N,INVDIA,Z,R)
   3 :       IMPLICIT REAL*8 (A-H,O-Z)
   4 :       IMPLICIT INTEGER (I-N)
   5 : C***********************************************************************
   6 : C NOM         : PSDIAG
   7 : C DESCRIPTION : Preconditionner solve diagonal
   8 : C
   9 : C
  10 : C
  11 : C LANGAGE     : ESOPE
  12 : C AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  13 : C               mél : gounand@semt2.smts.cea.fr
  14 : C***********************************************************************
  15 : C APPELES          : 
  16 : C APPELES (E/S)    : 
  17 : C APPELES (BLAS)   :
  18 : C APPELES (CALCUL) :
  19 : C APPELE PAR       :
  20 : C***********************************************************************
  21 : C SYNTAXE GIBIANE    :
  22 : C ENTREES            :
  23 : C ENTREES/SORTIES    :
  24 : C SORTIES            :
  25 : C CODE RETOUR (IRET) : = 0 si tout s'est bien passé
  26 : C***********************************************************************
  27 : C VERSION    : v1, 09/02/2006, version initiale
  28 : C HISTORIQUE : v1, 09/02/2006, création
  29 : C HISTORIQUE :
  30 : C HISTORIQUE :
  31 : C***********************************************************************
  32 : C Prière de PRENDRE LE TEMPS de compléter les commentaires
  33 : C en cas de modification de ce sous-programme afin de faciliter
  34 : C la maintenance !
  35 : C***********************************************************************
  36 :       INTEGER N
  37 :       REAL*8 INVDIA(N)
  38 :       REAL*8 Z(N)
  39 :       REAL*8 R(N)
  40 : *     
  41 : * Executable statements
  42 : *
  43 :       DO I=1,N
  44 :          Z(I)=R(I)*INVDIA(I)
  45 :       ENDDO
  46 : *
  47 : * Normal termination
  48 : *
  49 :       RETURN
  50 : *
  51 : * End of subroutine PSDIAG
  52 : *
  53 :       END
  54 : 
  55 :