home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Professional / OS2PRO194.ISO / os2 / prgramer / adaptor / examples / hpf_gen / pi / main.fcm < prev    next >
Encoding:
Text File  |  1993-03-23  |  1.2 KB  |  44 lines
  1. c-----This program calculates the value of pi, using numerical integration
  2. c-----with parallel processing, and clocks the solution time.
  3. c-----The user selects the number of points of integration and the limits a & b.
  4.  
  5. c-----PI program evaluates the integral of 4/(1+x*x) from a to b
  6. c-----by dividing the area under the curve [(b - a) / #points] slices.
  7. c-----Each processor will evaluate the area of the assigned slices
  8. c-----independently and simultaneously with other processors.
  9. c-----The sum of all the slices is calculated by calling the sum
  10. c-----reduction function
  11.  
  12. c   Variables:
  13. c
  14. c       points          Number of points used in the integration. Input by the user
  15. c       a               Lower limit of integration
  16. c       b               Upper limit of integration
  17. c
  18. c   Calls:
  19. c
  20. c       userinpt        to get number of points of integration from user
  21. c       pi        calculate the value of pi
  22.  
  23.  
  24.     PROGRAM MAIN
  25.  
  26.         integer userinpt
  27.         external userinpt
  28.       
  29.         integer j, points
  30.         double precision a, b
  31.  
  32. c       Receive points of integration from the user
  33. c
  34.  100    j = userinpt(a, b, points)
  35.         if (j.le.0) goto 200
  36.  
  37.     call pi(a, b, points)
  38.     
  39.         goto 100
  40.  
  41.  200    continue
  42.  
  43.         END
  44.