home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1993 #3 / NN_1993_3.iso / spool / sci / math / symbolic / 3512 < prev    next >
Encoding:
Internet Message Format  |  1993-01-25  |  2.7 KB
  1. Path: sparky!uunet!olivea!charnel!rat!usc!howland.reston.ans.net!zaphod.mps.ohio-state.edu!darwin.sura.net!newsserver.jvnc.net!netnews.upenn.edu!netnews.cc.lehigh.edu!ns1.cc.lehigh.edu!fc03
  2. From: fc03@ns1.cc.lehigh.edu (Frederick W. Chapman)
  3. Newsgroups: sci.math.symbolic
  4. Subject: Re: Symbolic Fortran : Help Needed
  5. Message-ID: <1993Jan25.154419.75401@ns1.cc.lehigh.edu>
  6. Date: 25 Jan 93 15:44:19 GMT
  7. Organization: Lehigh University
  8. Lines: 52
  9.  
  10. In article <1993Jan22.081540.1700@atri.curtin.edu.au>,
  11. thng@atri.curtin.edu.au (Ian Li-Jin Thng) writes:
  12.  
  13. >G'day to all symbolic netters,
  14. >
  15. >I'm seeking assistance on anyone who has experience with or have
  16. >come across references regarding symbolic programming using Fortran 
  17. >77 . 
  18. >
  19. >I know that a common method of symbolic programming is to use 
  20. >pointer structures which Fortran 77 do not support. 
  21. >
  22. >I have at present a Fortran program which is numerically intensive 
  23. >at the moment. I need to add a section to that program which is 
  24. >unfortunately symbolically intensive. The requirements of the 
  25. >symbolically intensive portion are as follows :
  26. >
  27. >Essentially to find out the symbolic expressions for the partial 
  28. >derivatives [...]
  29.  
  30. There is a package called JAKEF in the NETLIB collection which will
  31. perform "automatic differentiation" on scalar- and vector-valued
  32. functions of several variables, thereby yielding the gradient and
  33. Jacobian, respectively.  The results are exact up to round-off error.
  34. (I will post an "Introduction to NETLIB" document separately to this
  35. news group.)
  36.  
  37. For symbolic differentiation in FORTRAN (often desirable for
  38. optimization routines which need Jacobians and sometimes Hessians),
  39. you might want to search the literature on the topic "automatic
  40. differentiation", which refers to the process of taking a subroutine
  41. that computes a given differentiable function and automatically
  42. producing a second subroutine that computes the derivative of that
  43. function by an exact rather than approximate method.  (There is a
  44. description of this in the Spring '92 issue of the Maple Technical
  45. Newsletter.)  You might also want to query sci.math.num-analysis on
  46. the topic of "automatic differentiation".
  47.  
  48. Good luck,
  49.  
  50. Frederick W. Chapman
  51. Senior User Consultant
  52. Lehigh University Computing Center
  53.  
  54. -- 
  55.  
  56. o ------------------------------------------------------------------------- o
  57. |  Frederick W. Chapman, User Services, Computing Center, Lehigh University |
  58. |    Campus Phone:  8-3218     Preferred E-mail Address:  fc03@Lehigh.Edu   | 
  59. o ------------------------------------------------------------------------- o
  60. |             Ecstasy is transitory, but a theorem is forever!              |
  61. o ------------------------------------------------------------------------- o
  62.