home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1993 #1 / NN_1993_1.iso / spool / comp / lang / fortran / 4949 < prev    next >
Encoding:
Text File  |  1993-01-08  |  1.3 KB  |  39 lines
  1. Newsgroups: comp.lang.fortran
  2. Path: sparky!uunet!wupost!zaphod.mps.ohio-state.edu!cs.utexas.edu!torn!watserv2.uwaterloo.ca!watarts.uwaterloo.ca!kvetzal
  3. From: kvetzal@watarts.uwaterloo.ca (Ken Vetzal)
  4. Subject: Reaching NAG via e-mail?
  5. Message-ID: <C0K1Hy.BMJ@watserv2.uwaterloo.ca>
  6. Originator: kvetzal@watarts.uwaterloo.ca
  7. Sender: news@watserv2.uwaterloo.ca
  8. Organization: University of Waterloo
  9. Date: Fri, 8 Jan 1993 21:39:34 GMT
  10. Lines: 27
  11.  
  12. Does anyone know if NAG has an e-mail address which accepts reports of 
  13. suspected bugs?  I find that calling one of their random number 
  14. generator routines works fine if I request an even number of random
  15. numbers but blows up (produces Infinity - NaN) if an odd number of
  16. random numbers is requested.  A sample program which does this is
  17. included below:
  18.  
  19.  
  20.     implicit double precision(a-h,o-z) 
  21.     parameter(maxobs=365)
  22.     double precision e(maxobs)
  23.     external g05cbf, g05fdf
  24.     call g05cbf(170565)
  25.     do 100 n=1,2
  26.     call g05fdf(0.d0,1.d0,maxobs,e)
  27.     do 101 i=1,maxobs
  28.     write(8,*)i,e(i)
  29. 101    continue
  30. 100    continue
  31.     stop
  32.     end
  33.  
  34. This will blow up with maxobs set to 365 or 367 or any other odd
  35. number I've tried > 100, but seems to work fine at 364 or 366 or
  36. any other even number > 100.  Thanks in advance for any advice.
  37.  
  38. Ken Vetzal (kvetzal@watarts.uwaterloo.ca)
  39.