home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #20 / NN_1992_20.iso / spool / comp / lang / fortran / 3467 < prev    next >
Encoding:
Internet Message Format  |  1992-09-08  |  2.6 KB
  1. Path: sparky!uunet!usc!wupost!waikato.ac.nz!comp.vuw.ac.nz!windy!lhn!mdlcpgs
  2. Newsgroups: comp.lang.fortran
  3. Subject: Re: Permutation through recursion.
  4. Message-ID: <1992Sep8.134515.8196@lhn.gns.cri.nz>
  5. From: mdlcpgs@lhn.gns.cri.nz (Phil Scadden)
  6. Date: 8 Sep 92 13:45:15 +1200
  7. References: <998@engcon.marshall.ltv.com>
  8. Organization: Institute of Geological and Nuclear Sciences, Lower Hutt, New Zealand
  9. Lines: 66
  10.  
  11. In article <998@engcon.marshall.ltv.com>, ropella@engcon.marshall.ltv.com (GEROPELLA) writes:
  12. > 
  13. > The following is a VAX fortran polynomial evaluator.  But,
  14. > it only works for up to 3 variables.  What I would like from
  15. > ya'll is an algorithm that will evaluate any polynomial in
  16. > any number of dimensions for any order.  I feel that recursion
  17. > is the only answer and, of course, VAX fortran doesn't have
  18. > recursion built in; but I'm really going to implement it on
  19. > a DEC station, anyway, and DEC fortran does have recursion.
  20. > 
  21. > If somebody could tell me how to modify this, or if you have
  22. > a better way of doing it (in fortran, though); please post
  23. > the algorithm.  (Since our company was just bought, I'm 
  24. > afraid our email is getting screwed up.)  Or try to  email
  25. > me.
  26. > 
  27. > 
  28. > CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  29. > C...
  30. >       REAL*4 FUNCTION Poly(X,Eqn,var,order)
  31. >       IMPLICIT NONE
  32. >       REAL*4  X(4),TFLYOUT,mult,Eqn(*)
  33. >       INTEGER var,order,i,j,k,l,numcoef,p,ndx(4)
  34. > 
  35. >       TFLYOUT = 0.0
  36. >       numcoef = 0
  37. >       mult    =1.0
  38. > 
  39. > C  The following loop creates the alpha-index in this order:
  40. > C  000,001,002,010,011,020,100,101,110,200
  41. > C  where the first digit represents the power of the x variable,
  42. > C  the second represents the power of the y variable, and the third
  43. > C  represents the power of the z variable in the terms of the polynomial.
  44. > 
  45. >       do i=0,order
  46. >          ndx(1) = i
  47. >          do j=0,order-i
  48. >             ndx(2) = j
  49. >             do k=0,order-i-j
  50. >                ndx(3) = k
  51. >                mult = 1.0
  52. >                numcoef=numcoef+1
  53. >                do l=1,var
  54. >                   if ((x(l) .eq. 0.0) .and. (ndx(l) .eq. 0)) then
  55. >                     continue
  56. >                   else 
  57. >                     mult=mult*(X(l)**ndx(l))
  58. >                   end if
  59. >                end do
  60. > c                 write(*,*) (ndx(p),p=1,3)
  61. > c                 write(*,*) tflyout,' + ',Eqn(numcoef),' * ',mult
  62. >                TFLYOUT=TFLYOUT+Eqn(numcoef)*mult
  63. >             end do
  64. >          end do
  65. >       end do
  66. > 
  67. >       Poly = TFLYOUT
  68. > 
  69. >       RETURN
  70. >       END
  71. > 
  72. > 
  73. > glen e. ropella
  74. > ropella@engcon.marshall.ltv.com
  75. > ropella@engcon.uucp
  76. > 
  77.