home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1993 #1 / NN_1993_1.iso / spool / comp / lang / fortran / 4973 < prev    next >
Encoding:
Internet Message Format  |  1993-01-12  |  2.0 KB
  1. Path: sparky!uunet!pipex!bnr.co.uk!uknet!comlab.ox.ac.uk!oxuniv!nagrwb
  2. From: nagrwb@vax.oxford.ac.uk
  3. Newsgroups: comp.lang.fortran
  4. Subject: Re: Reaching NAG via e-mail?
  5. Message-ID: <1993Jan11.090817.11199@vax.oxford.ac.uk>
  6. Date: 11 Jan 93 09:08:17 GMT
  7. References: <C0K1Hy.BMJ@watserv2.uwaterloo.ca>
  8. Organization: Oxford University VAX 6620
  9. Lines: 43
  10.  
  11. In article <C0K1Hy.BMJ@watserv2.uwaterloo.ca>, kvetzal@watarts.uwaterloo.ca (Ken Vetzal) writes:
  12. > Does anyone know if NAG has an e-mail address which accepts reports of 
  13. > suspected bugs?  
  14. Certainly:    infodesk@nag.co.uk
  15.               infodesk@nag.com   for North American customers
  16. > I find that calling one of their random number 
  17. > generator routines works fine if I request an even number of random
  18. > numbers but blows up (produces Infinity - NaN) if an odd number of
  19. > random numbers is requested.  A sample program which does this is
  20. > included below:
  21. > 
  22. > 
  23. >     implicit double precision(a-h,o-z) 
  24. >     parameter(maxobs=365)
  25. >     double precision e(maxobs)
  26. >     external g05cbf, g05fdf
  27. >     call g05cbf(170565)
  28. >     do 100 n=1,2
  29. >     call g05fdf(0.d0,1.d0,maxobs,e)
  30. >     do 101 i=1,maxobs
  31. >     write(8,*)i,e(i)
  32. > 101    continue
  33. > 100    continue
  34. >     stop
  35. >     end
  36. > 
  37. > This will blow up with maxobs set to 365 or 367 or any other odd
  38. > number I've tried > 100, but seems to work fine at 364 or 366 or
  39. > any other even number > 100.  Thanks in advance for any advice.
  40. > 
  41. > Ken Vetzal (kvetzal@watarts.uwaterloo.ca)
  42. This error in G05FDF is known to us. More precisely, the error only occurs
  43. if two or more calls are made to G05FDF without an intervening call of G05CBF
  44. or G05CCF, and if at least one of the values of N in those calls is odd.
  45. The corrected version of the routine will be in the next release of
  46. the NAG Fortran Library and in a future error bulletin.
  47.  
  48. Our apologies for the inconvenience that this has caused.
  49.  
  50. Richard Brankin
  51. -- 
  52. ** R.W. Brankin -- Numerical Algorithms Group Ltd. -- nagrwb@vax.ox.ac.uk **
  53. **   suggestions for sig. welcome   --   will need   include<std.discl>   **
  54.