home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD1.mdf / fortran / library / ssp / matops / recp.for < prev    next >
Encoding:
Text File  |  1986-01-01  |  896 b   |  47 lines
  1. C
  2. C     ..................................................................
  3. C
  4. C     FUNCTION RECP
  5. C
  6. C     PURPOSE
  7. C        CALCULATE RECIPROCAL OF AN ELEMENT. THIS IS A FORTRAN
  8. C        FUNCTION SUBPROGRAM WHICH MAY BE USED AS AN ARGUMENT BY
  9. C        SUBROUTINE MFUN.
  10. C
  11. C     USAGE
  12. C        RECP(E)
  13. C
  14. C     DESCRIPTION OF PARAMETERS
  15. C        E - MATRIX ELEMENT
  16. C
  17. C     REMARKS
  18. C        RECIPROCAL OF ZERO IS TAKEN TO BE 1.0E75
  19. C
  20. C     SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
  21. C        NONE
  22. C
  23. C     METHOD
  24. C        RECIPROCAL OF ELEMENT E IS PLACED IN RECP
  25. C
  26. C     ..................................................................
  27. C
  28.       FUNCTION RECP(E)
  29. C
  30.       BIG=1.0E38
  31. C
  32. C     TEST ELEMENT FOR ZERO
  33. C
  34.       IF(E) 1,2,1
  35. C
  36. C        IF NON-ZERO, CALCULATE RECIPROCAL
  37. C
  38.     1 RECP=1.0/E
  39.       RETURN
  40. C
  41. C        IF ZERO, SET EQUAL TO INFINITY
  42. C
  43.     2 RECP=SIGN(BIG,E)
  44.       RETURN
  45.       END
  46. 
  47.