src/tool/Dependance_Machine/dmze10.F

   1       subroutine dmze10 ( vinfpo, zero, pd10vi, pd10ze )
   2 c ______________________________________________________________________
   3 c
   4 c                             H O M A R D
   5 c
   6 c Outil de Maillage Adaptatif par Raffinement et Deraffinement d'EDF R&D
   7 c
   8 c Version originale enregistree le 18 juin 1996 sous le numero 96036
   9 c aupres des huissiers de justice Simart et Lavoir a Clamart
  10 c Version 11.2 enregistree le 13 fevrier 2015 sous le numero 2015/014
  11 c aupres des huissiers de justice
  12 c Lavoir, Silinski & Cherqui-Abrahmi a Clamart
  13 c
  14 c    HOMARD est une marque deposee d'Electricite de France
  15 c
  16 c Copyright EDF 1996
  17 c Copyright EDF 1998
  18 c Copyright EDF 2002
  19 c Copyright EDF 2020
  20 c ______________________________________________________________________
  21 c   Dependance Machine : ZEros en puissance de 10
  22 c   -          -         --                    --
  23 c ______________________________________________________________________
  24 c
  25 c
  26 c  retourne la plus petite et la plus grande valeur possible
  27 c  et leurs puissances de 10 associees : x = 0.abc*10**n
  28 c  exemple : avec x = 1.e9,  on a x = 0.1*10**10, donc n = 10
  29 c            avec x = 1.e-40, on a x = 0.1*10**-39, donc n = -39
  30 c ______________________________________________________________________
  31 c .        .     .        .                                            .
  32 c .  nom   . e/s . taille .           description                      .
  33 c .____________________________________________________________________.
  34 c . vinfpo .  s  .  dp    . plus grande valeur positive possible       .
  35 c . zero   .  s  .  dp    . le zero de la machine                      .
  36 c . pd10vi .  s  .   i    . la puissance de 10 associee a vinfpo       .
  37 c . pd10ze .  s  .   i    . la puissance de 10 associee a zero         .
  38 c ______________________________________________________________________
  39 c
  40 c====
  41 c 0. declarations et dimensionnement
  42 c====
  43 c
  44 c 0.1. ==> generalites
  45 c
  46       implicit none
  47       save
  48 c
  49 c 0.2. ==> communs
  50 c
  51 c 0.3. ==> arguments
  52 c
  53       double precision vinfpo, zero
  54       integer pd10vi, pd10ze
  55 c
  56 c 0.4. ==> variables locales
  57 c
  58 c 0.5. ==> initialisations
  59 c
  60 c====
  61 c 1. les valeurs extremes
  62 c===
  63 c
  64 c 1.1. ==> le maximum
  65 c
  66       vinfpo = 1.d29
  67       pd10vi = 30
  68 c
  69 c 1.2. ==> le minimum : l'inverse du maximum
  70 c
  71       zero = 1.d0 / vinfpo
  72 c
  73 c 1.3. ==> les puissances de 10 associees
  74 c
  75       pd10ze = 2 - pd10vi
  76 c
  77       end