Article View: sci.math.symbolic
Article #954Symbolic Fortran : Help Needed
From: Ian Li-Jin Thng
Date: Fri, 22 Jan 1993 08:15
Date: Fri, 22 Jan 1993 08:15
61 lines
2707 bytes
2707 bytes
G'day to all symbolic netters,
I'm seeking assistance on anyone who has experience with or have
come across references regarding symbolic programming using Fortran
77 .
I know that a common method of symbolic programming is to use
pointer structures which Fortran 77 do not support.
I have at present a Fortran program which is numerically intensive
at the moment. I need to add a section to that program which is
unfortunately symbolically intensive. The requirements of the
symbolically intensive portion are as follows :
Essentially to find out the symbolic expressions for the partial
derivatives of
w(T) . Q . w
where
(T) denotes transpose of a vector or matrix
w is a vector whose elements consists of different variables
eg : w = [2 , 2a+b, a^2, b^4, 3a , 4c-a, 4, 5.6 ](T)
Q is a matrix whose elements are similar to the elements in w,
ie , they can be symbolic elements or purely numerical entries
Help is needed on how to program in Fortran 77 to solve the following
problems :
a) How to represent these vector structures and matrix structures
elegantly.
b) Perform matrix operation (ie multiplications, additions, subtraction)
properly on these structures
c) Sort the variables out ,say, in ascending orders of the variable 'a',
followed by variable 'b' and finally in variable 'c' .
d) Obtaining the partial derivatives of w(T) . Q . w wrt variables 'a'
'b' and 'c' .
I would really appreciate anyone who can provide info on the above
problems . You can contact me at ...
------------------------------------------------------------------------------
Ian Thng - Forever Trying to be Young
Email:thng@atri.curtin.edu.au
Adaptive Signal Processing Laboratory
Australian Telecommunications Research Institute
Curtin University of Technology
Ph : +61-9-3513270
Fax : +61-9-3513244
Bothered with seat belts ? Ride a Harley Mate !
------------------------------------------------------------------------------
A thousand thanx in advance to all who can help me out !
Technical Details
Message-ID:
<1993Jan22.081540.1700@atri.curtin.edu.au>Path:
arc1-ca.pugleaf.net!sparky!uunet!munnari.oz.au!uniwa!cujo!atri.curtin.edu.au!Giuseppe!thngShow Full Headers
Complete Article Headers
Newsgroups: sci.math.symbolic Path: arc1-ca.pugleaf.net!sparky!uunet!munnari.oz.au!uniwa!cujo!atri.curtin.edu.au!Giuseppe!thng From: thng@atri.curtin.edu.au (Ian Li-Jin Thng) Subject: Symbolic Fortran : Help Needed Message-ID: <1993Jan22.081540.1700@atri.curtin.edu.au> Sender: news@atri.curtin.edu.au (news admin directory) Reply-To: thng@atri.curtin.edu.au Organization: Australian Telecommunication Research Institute Date: Fri, 22 Jan 93 08:15:40 GMT Lines: 61