home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #18 / NN_1992_18.iso / spool / comp / lang / apl / 950 < prev    next >
Encoding:
Internet Message Format  |  1992-08-15  |  3.3 KB
  1. Path: sparky!uunet!psinntp!kepler1!andrew
  2. From: andrew@rentec.com (Andrew Mullhaupt)
  3. Newsgroups: comp.lang.apl
  4. Subject: Re: Promoting APL2
  5. Message-ID: <1176@kepler1.rentec.com>
  6. Date: 15 Aug 92 20:48:21 GMT
  7. References: <168417F7D.JOHNSON2@vm.cc.purdue.edu> <1992Aug13.171145.24691@csi.jpl.nasa.gov>
  8. Organization: Renaissance Technologies Corp., Setauket, NY.
  9. Lines: 70
  10.  
  11. In article <1992Aug13.171145.24691@csi.jpl.nasa.gov> sam@csi.jpl.nasa.gov (Sam Sirlin) writes:
  12. >
  13. >In article <168417F7D.JOHNSON2@vm.cc.purdue.edu>, JOHNSON@vm.cc.purdue.edu writes:
  14. >I don't know what would end up being most popular, but for me there
  15. >are a few neccessities:
  16. >
  17. >1. linear algebra - as much of eispack/linpack/lapack as you can get.
  18. >Perhaps you could link them in via some sort of quad na? I've
  19. >translated the qr and qz eigensystem paths into APL if you want them.
  20. >I have them in readable ASCII form.
  21.  
  22. NO! NO! A thousand times NO!. These guys have a 3090, and they should use
  23. ESSL, which is provided by IBM. Although ESSL has some trivial flaws, it
  24. is probably the only commercial package which is as good or better than
  25. LINPACK/EISPACK/LAPACK. ESSL has the world's _best_ FFT's bar none. (what
  26. else would you expect from Winograd et. al.?) Last time I used it, IBM
  27. had gone to the trouble to write the required interface for all of ESSL
  28. routines to quad-NA. Use the IBM provided ESSL.
  29.  
  30. Also: there was talk that OSL would be done in the same way - if so use
  31. that for all of your constraint-bound optimizations.
  32.  
  33. >2. ode solvers - the more the better. Some sort of rk45 for example.
  34. >This also needs to include random colored noise generators (for
  35. >example the OU of ACSL). A stiff method would be nice.
  36.  
  37. A painful issue with APL, since you have to call the function defining
  38. the DE from the solver. I've been told that APL2 can _now_ be called from
  39. FORTRAN, but I've never seen it. If you can't, it's too bad. If you can,
  40. your best bet is to use FORTRAN versions.
  41.  
  42. >4. random number generators - for the various processes. At least
  43. >Gaussian and Poisson.
  44.  
  45. ESSL provides a good uniform generator.
  46.  
  47. >5. integration methods - I don't really need to much here, trapezoidal
  48. >rules and Simpson seem to work. 
  49.  
  50. ESSL provides several excellent methods. Linking in quadpack would
  51. be even better, but we run into the "can FORTRAN call APL" problem
  52. again.
  53.  
  54. >6. special functions - Bessel, erf, etc.
  55.  
  56. Although you can profitably construct these using ESSL, you might
  57. be better off linking in PD FORTRAN code.
  58.  
  59. >7. controls analysis - frequency response calulation etc. Well that's
  60. >what I want anyway... 
  61.  
  62. A good chunk of the foundationns you need for this is provided by
  63. ESSL's rich and "world's best" FT/Correlations/etc. functions.
  64.  
  65. >8. data translation - functions for reading/writing data from/for
  66. >other languages (fortran). This also includes text manipulation - one
  67. >of APL's strong points since (at least STSC/J) can read exactly what's
  68. >in a file.
  69.  
  70. But they live on a 3090, which might be running MVS. In which case
  71. the file system is a disaster beyond comprehension.
  72.  
  73. A few more things about ESSL - ESSL is finely tuned to take advantage
  74. of the resources the 3090 offers, even across different models of the
  75. 3090. It can take advantage of vector heads, as well as run on multiple
  76. heads, for some functions. APL2 on the other hand only makes use of
  77. vectorization, at least last time I looked.
  78.  
  79. Later,
  80. Andrew Mullhaupt
  81.