From: Dieter K. <cra...@us...> - 2009-11-01 00:19:07
|
Update of /cvsroot/maxima/maxima/src In directory 23jxhf1.ch3.sourceforge.com:/tmp/cvs-serv2170/src Modified Files: hypgeo.lisp Log Message: Cutting out the Laplace transform of the symbol %d. %d has been replaced by $parabolic_cylinder_d, that is the new symbol for the Parabolic Cylinder D function. No problems with the testsuite. Index: hypgeo.lisp =================================================================== RCS file: /cvsroot/maxima/maxima/src/hypgeo.lisp,v retrieving revision 1.71 retrieving revision 1.72 diff -u -d -r1.71 -r1.72 --- hypgeo.lisp 29 Oct 2009 20:48:33 -0000 1.71 +++ hypgeo.lisp 1 Nov 2009 00:15:56 -0000 1.72 @@ -139,10 +139,6 @@ (defun pjac (x n a b) (list '(mqapply) (list '($%p array) n a b) x)) -;; Parabolic cylinder function D -(defun parcyl (x n) - (list '(mqapply) (list '($%d array) n) x)) - ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;; Two general pattern for the routine lt-sf-log. @@ -771,16 +767,6 @@ ((coeffpp)(a zerp))) nil)) -;; Recognize %d[v](w), Parabolic Cylinder function -(defun oned (exp) - (m2 exp - '((mplus) - ((coeffpt) - (u nonzerp) - ((mqapply)(($%d array) (v true)) (w true))) - ((coeffpp) (a zerp))) - nil)) - ;; Recognize parabolic_cylinder_d function (defun m2-parabolic_cylinder_d (expr) (m2 expr @@ -1709,7 +1695,11 @@ (let ((inv4 (inv 4))) (cond ((or $prefer_d (whittindtest (add (div v 2) inv4) inv4)) - (parcyl z v)) + ;; At this time the Parabolic Cylinder D function is not implemented + ;; as a simplifying function. We call nevertheless the simplifer + ;; to simplify the arguments. When we implement the function + ;; The symbol has to be changed to the noun form. + (simplify (list '($parabolic_cylinder_d) v z))) (t (simpdtf z v))))) (defun whittindtest (i1 i2) @@ -2163,13 +2153,6 @@ (return (fractest2 rest arg1 index1 nil 'kti)))))) ;; Laplace transform of Parabolic Cylinder function - (cond ((setq l (oned u)) - (setq index1 (cdras 'v l) - arg1 (cdras 'w l) - rest (cdras 'u l)) - (return (fractest2 rest arg1 index1 nil 'd)))) - - ;; Laplace transform of Parabolic Cylinder function (cond ((setq l (m2-parabolic_cylinder_d u)) (setq index1 (cdras 'v l) arg1 (cdras 'w l) @@ -3020,7 +3003,11 @@ (defun hetd (x n) (mul* (power '$%e (mul* x x (inv 4))) - (parcyl x n))) + ;; At this time the Parabolic Cylinder D function is not implemented + ;; as a simplifying function. We call nevertheless the simplifer + ;; to simplify the arguments. When we implement the function + ;; The symbol has to be changed to the noun form. + (simplify (list '($parabolic_cylinder_d) n x)))) ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; @@ -3112,7 +3099,12 @@ (mul* (power 2 inv2) ; Should this be 2^(1/4)? (power '$%pi (mul* -1 inv2)) (power '$%e (mul* -1 inv2 x x)) - (parcyl (mul* (power 2 inv2) x) -1)))) + ;; At this time the Parabolic Cylinder D function is not implemented + ;; as a simplifying function. We call nevertheless the simplifer + ;; to simplify the arguments. When we implement the function + ;; The symbol has to be changed to the noun form. + (simplify + (list '($parabolic_cylinder_d) -1 (mul (power 2 inv2) x)))))) ;; Lommel S function in terms of Bessel J and Bessel Y. ;; Luke gives |