Esope source files

   1 : C VLH1      SOURCE    CHAT      05/01/13    04:07:41     5004           
   2 :       SUBROUTINE VLH1(H1,H1P,X,G)
   3 : C************************************************************************
   4 : C
   5 : C PROJET            :  CASTEM 2000
   6 : C
   7 : C NOM               :  VLH1 (subroutine du groupe RIECOM.ESO)
   8 : C
   9 : C DESCRIPTION       :  Fonction H1 de SMOLLER (voir RIEMAN)
  10 : C
  11 : C LANGAGE           :  FORTRAN 77 
  12 : C      
  13 : C AUTEUR            :  P. GALON DRN/DMT/SEMT/TTMF  
  14 : C
  15 : C************************************************************************
  16 : C
  17 : C APPELES           : aucun
  18 : C
  19 : C************************************************************************
  20 : C 
  21 : C Entrées : X, G
  22 : C
  23 : C Sorties :    H1   FONCTION H1
  24 : C              H1P  FONCTION DERIVEE DE H1
  25 : C
  26 : C************************************************************************
  27 : C
  28 : C HISTORIQUE (Anomalies et modifications éventuelles)
  29 : C
  30 : C HISTORIQUE :  créé le 1.10.96
  31 : C
  32 : C HISTORIQUE :  modifié par A. BECCANTINI le 7.1.98 pour
  33 : C               ajouter des commentaires et eliminer la
  34 : C               partié relative au cas X = 0.      
  35 : C                          
  36 : C 
  37 : C************************************************************************
  38 : C
  39 : C N.B.: Toutes le variables sont DECLAREES
  40 : C      
  41 :       IMPLICIT INTEGER(I-N)
  42 :       IMPLICIT REAL*8 (A-H,O-Z)
  43 :       REAL*8 H1, H1P, X, G, GM1, UNSGM1, T, ETX, AUX1, AUX2, B, EX
  44 : C
  45 :       GM1 = G-1.0D0
  46 :       UNSGM1 = 1.0D0 / GM1
  47 :       T = 0.5D0*GM1/G
  48 :       IF(X .GE. 0.D0) THEN
  49 : C
  50 : C******** X>0
  51 : C
  52 :          ETX = EXP(-T*X)
  53 :          H1 = 2.D0*UNSGM1*(1.D0-ETX)
  54 :          H1P = 2.D0*UNSGM1*T*ETX
  55 : C
  56 : C***** NB Function continue et derivable in X = 0
  57 : C        
  58 : C      X=0 
  59 : C
  60 : C      H1 = 0.0D0
  61 : C      H1P = 1.D0/G
  62 : C         
  63 :       ELSE
  64 : C
  65 : C******** X<0
  66 : C
  67 :          B = (G+1.0D0)*UNSGM1
  68 :          EX = EXP(-X)
  69 :          AUX1 = 1.0D0 / SQRT(1.D0+B*EX)
  70 :          AUX2 = 2.0D0*SQRT(T)*UNSGM1*AUX1
  71 :          H1 = AUX2*(1.0D0-EX)
  72 :          H1P = AUX2*EX* ( 1.D0 + 0.5D0*(1.D0-EX)*B*AUX1*AUX1 )
  73 :       ENDIF
  74 :       RETURN
  75 :       END
  76 :  
  77 : 
  78 : 
  79 : 
  80 : 
  81 : 
  82 : 
  83 : 
  84 : 
  85 : 
  86 :  
  87 :