mpmaxr.f 1.1 KB

123456789101112131415161718192021222324252627282930313233343536373839
  1. *DECK MPMAXR
  2. SUBROUTINE MPMAXR (X)
  3. C***BEGIN PROLOGUE MPMAXR
  4. C***SUBSIDIARY
  5. C***PURPOSE Subsidiary to DQDOTA and DQDOTI
  6. C***LIBRARY SLATEC
  7. C***TYPE ALL (MPMAXR-A)
  8. C***AUTHOR (UNKNOWN)
  9. C***DESCRIPTION
  10. C
  11. C Sets X to the largest possible positive 'mp' number.
  12. C
  13. C The argument X(*) is an INTEGER arrays of size 30. See the comments
  14. C in the routine MPBLAS for the reason for this choice.
  15. C
  16. C***SEE ALSO DQDOTA, DQDOTI, MPBLAS
  17. C***ROUTINES CALLED MPCHK
  18. C***COMMON BLOCKS MPCOM
  19. C***REVISION HISTORY (YYMMDD)
  20. C 791001 DATE WRITTEN
  21. C ?????? Modified for use with BLAS. Blank COMMON changed to named
  22. C COMMON. R given dimension 12.
  23. C 891214 Prologue converted to Version 4.0 format. (BAB)
  24. C 900402 Added TYPE section. (WRB)
  25. C 930124 Increased Array size in MPCON for SUN -r8. (RWC)
  26. C***END PROLOGUE MPMAXR
  27. COMMON /MPCOM/ B, T, M, LUN, MXR, R(30)
  28. INTEGER B, T, R, X(*)
  29. C***FIRST EXECUTABLE STATEMENT MPMAXR
  30. CALL MPCHK (1, 4)
  31. IT = B - 1
  32. C SET FRACTION DIGITS TO B-1
  33. DO 10 I = 1, T
  34. 10 X(I+2) = IT
  35. C SET SIGN AND EXPONENT
  36. X(1) = 1
  37. X(2) = M
  38. RETURN
  39. END