dmacon.f 1.0 KB

1234567891011121314151617181920212223242526272829303132333435
  1. *DECK DMACON
  2. SUBROUTINE DMACON
  3. C***BEGIN PROLOGUE DMACON
  4. C***SUBSIDIARY
  5. C***PURPOSE Subsidiary to DBVSUP
  6. C***LIBRARY SLATEC
  7. C***TYPE DOUBLE PRECISION (MACON-S, DMACON-D)
  8. C***AUTHOR (UNKNOWN)
  9. C***SEE ALSO DBVSUP
  10. C***ROUTINES CALLED D1MACH
  11. C***COMMON BLOCKS DML5MC
  12. C***REVISION HISTORY (YYMMDD)
  13. C 750601 DATE WRITTEN
  14. C 890531 Changed all specific intrinsics to generic. (WRB)
  15. C 890921 Realigned order of variables in certain COMMON blocks.
  16. C (WRB)
  17. C 891214 Prologue converted to Version 4.0 format. (BAB)
  18. C 900328 Added TYPE section. (WRB)
  19. C***END PROLOGUE DMACON
  20. DOUBLE PRECISION D1MACH
  21. INTEGER KE, LPAR
  22. DOUBLE PRECISION DD, EPS, FOURU, SQOVFL, SRU, TWOU, URO
  23. COMMON /DML5MC/ URO,SRU,EPS,SQOVFL,TWOU,FOURU,LPAR
  24. C***FIRST EXECUTABLE STATEMENT DMACON
  25. URO = D1MACH(4)
  26. SRU = SQRT(URO)
  27. DD = -LOG10(URO)
  28. LPAR = 0.5D0*DD
  29. KE = 0.5D0 + 0.75D0*DD
  30. EPS = 10.0D0**(-2*KE)
  31. SQOVFL = SQRT(D1MACH(2))
  32. TWOU = 2.0D0*URO
  33. FOURU = 4.0D0*URO
  34. RETURN
  35. END