Maxima -- GPL CAS based on DOE-MACSYMA

Update of /cvsroot/maxima/maxima/share/lapack
In directory sfp-cvsdas-1.v30.ch3.sourceforge.com:/tmp/cvs-serv11370/share/lapack

Modified Files:
	lapack-lisp.system 
Log Message:
Fix some missing dependencies.  (They were correct in lapack.system,
but not updated for lapack-lisp.system.)


Index: lapack-lisp.system
===================================================================
RCS file: /cvsroot/maxima/maxima/share/lapack/lapack-lisp.system,v
retrieving revision 1.1
retrieving revision 1.2
diff -u -d -r1.1 -r1.2
--- lapack-lisp.system	8 Jan 2009 18:25:23 -0000	1.1
+++ lapack-lisp.system	14 Dec 2009 22:48:47 -0000	1.2
@@ -52,7 +52,15 @@
 					  "dlaset"))
  	     (:file "dlabad")
  	     (:file "dlacpy")
- 	     (:file "dlamch")
+ 	     (:file "dlamch"
+		    :depends-on ("dlamc2" "dlamc3" "dlamc4" "dlamc5"))
+	     (:file "dlamc1"
+		    :depends-on ("dlamc3"))
+	     (:file "dlamc2"
+		    :depends-on ("dlamc1" "dlamc3" "dlamc4" "dlamc5"))
+	     (:file "dlamc3")
+	     (:file "dlamc4")
+	     (:file "dlamc5")
  	     (:file "dlange"
 		    :depends-on ("dlassq"))
  	     (:file "dlapy2")

Update of /cvsroot/maxima/maxima/share/lapack/lapack
In directory sfp-cvsdas-1.v30.ch3.sourceforge.com:/tmp/cvs-serv11016/share/lapack/lapack

Modified Files:
	dlascl.lisp 
Log Message:
Bug 2913614: LAPACK: dgesvd is broken

Issue caused by incorrect translation by f2cl.  f2cl is fixed, and the
translation updated.


Index: dlascl.lisp
===================================================================
RCS file: /cvsroot/maxima/maxima/share/lapack/lapack/dlascl.lisp,v
retrieving revision 1.5
retrieving revision 1.6
diff -u -d -r1.5 -r1.6
--- dlascl.lisp	8 Apr 2009 00:04:21 -0000	1.5
+++ dlascl.lisp	14 Dec 2009 22:47:00 -0000	1.6
@@ -3,11 +3,11 @@
 ;;;  "f2cl2.l,v 1.37 2008/02/22 22:19:33 rtoy Exp $"
 ;;;  "f2cl3.l,v 1.6 2008/02/22 22:19:33 rtoy Exp $"
 ;;;  "f2cl4.l,v 1.7 2008/02/22 22:19:34 rtoy Exp $"
-;;;  "f2cl5.l,v 1.200 2009/01/19 02:38:17 rtoy Exp $"
+;;;  "f2cl5.l,v 1.201 2009/12/14 22:36:46 rtoy Exp $"
 ;;;  "f2cl6.l,v 1.48 2008/08/24 00:56:27 rtoy Exp $"
 ;;;  "macros.l,v 1.112 2009/01/08 12:57:19 rtoy Exp $")
 
-;;; Using Lisp CMU Common Lisp 19f (19F)
+;;; Using Lisp CMU Common Lisp Snapshot 2009-12 (20A Unicode)
 ;;; 
 ;;; Options: ((:prune-labels nil) (:auto-save t) (:relaxed-array-decls t)
 ;;;           (:coerce-assigns :as-needed) (:array-type ':array)
@@ -27,7 +27,7 @@
              (type (f2cl-lib:integer4) info lda n m ku kl)
              (type (simple-array character (*)) type))
     (f2cl-lib:with-multi-array-data
-        ((type double-float type-%data% type-%offset%)
+        ((type character type-%data% type-%offset%)
          (a double-float a-%data% a-%offset%))
       (prog ((bignum 0.0) (cfrom1 0.0) (cfromc 0.0) (cto1 0.0) (ctoc 0.0)
              (mul 0.0) (smlnum 0.0) (i 0) (itype 0) (j 0) (k1 0) (k2 0) (k3 0)

Update of /cvsroot/maxima/maxima/src
In directory sfp-cvsdas-1.v30.ch3.sourceforge.com:/tmp/cvs-serv8932/src

Modified Files:
	bessel.lisp 
Log Message:
Reformatting the file to prepare further work on the Bessel functions. The functions are put in more logical order and the code is more readable. 
No more changes.

No problems with the testsuite.


Index: bessel.lisp
===================================================================
RCS file: /cvsroot/maxima/maxima/src/bessel.lisp,v
retrieving revision 1.79
retrieving revision 1.80
diff -u -d -r1.79 -r1.80
--- bessel.lisp	18 Nov 2009 22:57:35 -0000	1.79
+++ bessel.lisp	14 Dec 2009 22:38:07 -0000	1.80
@@ -13,485 +13,37 @@
 
 (defmvar $besselexpand nil)
 
-;; Compute value of Modified Bessel function of the first kind of order n
-(defun bessel-i (order arg)
-  (cond ((zerop (imagpart arg))
-	 ;; We have numeric args and the first arg is purely
-	 ;; real. Call the real-valued Bessel function.  Use special
-	 ;; routines for order 0 and 1, when possible
-	 (let ((arg (realpart arg)))
-	   (cond ((zerop order)
[...2324 lines suppressed...]
@@ -1725,8 +1664,6 @@
 (defun $struve_l (v z)
   (simplify (list '(%struve_l) v z)))
 
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-
 (defprop $struve_l %struve_l alias)
 (defprop $struve_l %struve_l verb)
 (defprop %struve_l $struve_l reversealias)
@@ -1965,3 +1902,10 @@
        (eqtest (list '(%struve_l) order arg) expr)))))
 
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+
+(defmspec $gauss (form)
+  (format t
+"NOTE: The gauss function is superseded by random_normal in the `distrib' package.
+Perhaps you meant to enter `~a'.~%"
+    (print-invert-case (implode (mstring `(($random_normal) ,@ (cdr form))))))
+  '$done)

Update of /cvsroot/maxima/maxima/src
In directory sfp-cvsdas-1.v30.ch3.sourceforge.com:/tmp/cvs-serv24198/src

Modified Files:
	nparse.lisp 
Log Message:
Fix a few minor issues in test-make-number.  Print out each test too,
so we can see what's happening since the tests are rather slow.

No functional changes anywhere else.



Index: nparse.lisp
===================================================================
RCS file: /cvsroot/maxima/maxima/src/nparse.lisp,v
retrieving revision 1.54
retrieving revision 1.55
diff -u -d -r1.54 -r1.55
--- nparse.lisp	10 Dec 2009 18:36:46 -0000	1.54
+++ nparse.lisp	14 Dec 2009 14:41:04 -0000	1.55
@@ -277,18 +277,18 @@
 (defvar *fast-bfloat-extra-bits* 0)
 
 ;; Here is a test routine to test the fast bfloat conversion
-#||
+#+nil
 (defun test-make-number (&optional (n 1000))
   (let ((failures 0))
     (dotimes (k n)
       (flet ((digit-list (n)
-	       (coerce (format nil "~D" (random n)) 'list)))
+	       (coerce (format nil "~D" n) 'list)))
 	(let ((numlist nil))
 	  ;; Generate a random number with 30 fraction digits and an
 	  ;; large exponent.
-	  (push (digit-list 10) numlist)
+	  (push (digit-list (random 10)) numlist)
 	  (push '(#\.) numlist)
-	  (push (digit-list (expt 10 30)) numlist)
+	  (push (digit-list (random (expt 10 30))) numlist)
 	  (push '(#\B) numlist)
 	  (push (if (zerop (random 2)) '(#\+) '(#\-)) numlist)
 	  (push (digit-list (+ $fast_bfloat_threshold
@@ -300,13 +300,17 @@
 			(make-number (copy-list numlist))))
 		(fast (let (($fast_bfloat_conversion t))
 			(make-number (copy-list numlist)))))
+	    (format t "Test ~3A: " k)
+	    (map nil #'(lambda (x)
+			 (map nil #'princ x))
+		 (reverse numlist))
+	    (terpri)
 	    (unless (equalp true fast)
 	      (incf failures)
 	      (format t "NUM:  ~A~%  TRUE: ~S~%  FAST: ~S~%"
 		      (reverse numlist) true fast))))))
     (format t "~D failures in ~D tests (~F%)~%"
 	    failures n (* 100 failures (/ (float n))))))
-||#
 
 
 ;; WARNING: MAKE-NUMBER destructively modifies it argument!  Should we

Update of /cvsroot/maxima/maxima/tests
In directory sfp-cvsdas-1.v30.ch3.sourceforge.com:/tmp/cvs-serv22203/tests

Modified Files:
	rtest16.mac 
Log Message:
Bug : 2914176 Conversion of rational to bfloat is inaccurate

Implement method to convert rationals to bfloats using only one
rounding.   (This might slow down conversion of large rationals to
bfloats.) 

src/float.lisp:
o Modify $bfloat and add special case to handle Maxima rationals,
  instead of using the default floatprog.
o Modify RATBIGFLOAT to call the new routine
o Add FLOAT-RATIO routine to convert rationals to bfloats with just
  one rounding operation.

tests/rtest16.mac:
o Add test.


Index: rtest16.mac
===================================================================
RCS file: /cvsroot/maxima/maxima/tests/rtest16.mac,v
retrieving revision 1.106
retrieving revision 1.107
diff -u -d -r1.106 -r1.107
--- rtest16.mac	13 Dec 2009 22:50:45 -0000	1.106
+++ rtest16.mac	14 Dec 2009 14:21:43 -0000	1.107
@@ -1156,3 +1156,16 @@
 0;
 limit(cos(1/x)*sin(x)-sin(x)+a,x,inf);
 a;
+
+/* Bug ID: 2914176  	 Conversion of rational to bfloat is inaccurate
+ *
+ * The difference should be 1/262144, but we don't check for that.
+ */
+(oldfpprec:fpprec, fpprec:5, done);
+done;
+is(bfloat((2^20+1)/(2^20-1)) - 1b0 > 0);
+true;
+
+(fpprec:oldfpprec, done);
+done;
+

Update of /cvsroot/maxima/maxima/src
In directory sfp-cvsdas-1.v30.ch3.sourceforge.com:/tmp/cvs-serv22203/src

Modified Files:
	float.lisp 
Log Message:
Bug : 2914176 Conversion of rational to bfloat is inaccurate

Implement method to convert rationals to bfloats using only one
rounding.   (This might slow down conversion of large rationals to
bfloats.) 

src/float.lisp:
o Modify $bfloat and add special case to handle Maxima rationals,
  instead of using the default floatprog.
o Modify RATBIGFLOAT to call the new routine
o Add FLOAT-RATIO routine to convert rationals to bfloats with just
  one rounding operation.

tests/rtest16.mac:
o Add test.


Index: float.lisp
===================================================================
RCS file: /cvsroot/maxima/maxima/src/float.lisp,v
retrieving revision 1.58
retrieving revision 1.59
diff -u -d -r1.58 -r1.59
--- float.lisp	14 Dec 2009 04:55:49 -0000	1.58
+++ float.lisp	14 Dec 2009 14:21:42 -0000	1.59
@@ -390,6 +390,8 @@
 	       (exptbigfloat ($bfloat (cadr x)) (caddr x))))
 	  ((eq (caar x) 'mncexpt)
 	   (list '(mncexpt) ($bfloat (cadr x)) (caddr x)))
+	  ((eq (caar x) 'rat)
+	   (ratbigfloat (cdr x)))
 	  ((setq y (safe-get (caar x) 'floatprog))
 	   (funcall y (mapcar #'$bfloat (cdr x))))
 	  ((or (trigp (caar x)) (arcp (caar x)) (eq (caar x) '$entier))
@@ -429,8 +431,71 @@
 		   ((equal fans tst) nfans)
 		   (t (simplify (list '(mplus) fans nfans)))))))
 
-(defmfun ratbigfloat (l)
-  (bcons (fpquotient (cdar l) (cdadr l))))
+(defmfun ratbigfloat (r)
+  (bcons (float-ratio r)))
+
+;; This is borrowed from CMUCL (float-ratio-float), and modified for
+;; converting ratios to Maxima's bfloat numbers.
+(defun float-ratio (x)
+  (let* ((signed-num (first x))
+	 (plusp (plusp signed-num))
+	 (num (if plusp signed-num (- signed-num)))
+	 (den (second x))
+	 (digits fpprec)
+	 (scale 0))
+    (declare (fixnum digits scale))
+    ;;
+    ;; Strip any trailing zeros from the denominator and move it into the scale
+    ;; factor (to minimize the size of the operands.)
+    (let ((den-twos (1- (integer-length (logxor den (1- den))))))
+      (declare (fixnum den-twos))
+      (decf scale den-twos)
+      (setq den (ash den (- den-twos))))
+    ;;
+    ;; Guess how much we need to scale by from the magnitudes of the numerator
+    ;; and denominator.  We want one extra bit for a guard bit.
+    (let* ((num-len (integer-length num))
+	   (den-len (integer-length den))
+	   (delta (- den-len num-len))
+	   (shift (1+ (the fixnum (+ delta digits))))
+	   (shifted-num (ash num shift)))
+      (declare (fixnum delta shift))
+      (decf scale delta)
+      (labels ((float-and-scale (bits)
+		 (let* ((bits (ash bits -1))
+			(len (integer-length bits)))
+		   (cond ((> len digits)
+			  (assert (= len (the fixnum (1+ digits))))
+			  (multiple-value-bind (f0)
+			      (floatit (ash bits -1))
+			    (list (first f0) (+ (second f0)
+						(1+ scale)))))
+			 (t
+			  (multiple-value-bind (f0)
+			      (floatit bits)
+			    (list (first f0) (+ (second f0) scale)))))))
+	       (floatit (bits)
+		 (let ((sign (if plusp 1 -1)))
+		   (list (* sign bits) 0))))
+	(loop
+	  (multiple-value-bind (fraction-and-guard rem)
+	      (truncate shifted-num den)
+	    (let ((extra (- (integer-length fraction-and-guard) digits)))
+	      (declare (fixnum extra))
+	      (cond ((/= extra 1)
+		     (assert (> extra 1)))
+		    ((oddp fraction-and-guard)
+		     (return
+		       (if (zerop rem)
+			   (float-and-scale
+			    (if (zerop (logand fraction-and-guard 2))
+				fraction-and-guard
+				(1+ fraction-and-guard)))
+			   (float-and-scale (1+ fraction-and-guard)))))
+		    (t
+		     (return (float-and-scale fraction-and-guard)))))
+	    (setq shifted-num (ash shifted-num -1))
+	    (incf scale)))))))
 
 (defun decimalsin (x)
   (do ((i (quotient (* 59. x) 196.) (1+ i))) ;log[10](2)=.301029

Update of /cvsroot/maxima/maxima/tests
In directory sfp-cvsdas-1.v30.ch3.sourceforge.com:/tmp/cvs-serv30905/tests

Modified Files:
	rtest_expintegral.mac rtest_gamma.mac 
Log Message:
Bug 2909980: optimization for log of bfloats

Optimize computation of log of bfloats by using the fact that each
bfloat is frac*2^exp, where 1/2 < frac <= 1.

src/float.lisp:
o Add new variable bigfloat_log2 to hold the current bigfloat value of
  log(2b0).
o Add function fplog2 to return log(2b0) of the correct precision.
o Modify big-float-log to compute logs using the formula
  log(frac*2^exp) = log(frac) + exp*log(2).

tests/rtest_expintegral.mac:
o Test 73 (expintegral_e(1,.25b0)) no longer is accurate to 64 digits;
  reduce accuracy to 63.

tests/rtest_gamma.mac:
o Test 636 (beta_incomplete(1b0,1,z)) is no longer exactly
  1b0*z. Change test to allow a small error.



Index: rtest_expintegral.mac
===================================================================
RCS file: /cvsroot/maxima/maxima/tests/rtest_expintegral.mac,v
retrieving revision 1.7
retrieving revision 1.8
diff -u -d -r1.7 -r1.8
--- rtest_expintegral.mac	1 Aug 2009 13:07:02 -0000	1.7
+++ rtest_expintegral.mac	14 Dec 2009 04:55:50 -0000	1.8
@@ -261,7 +261,7 @@
 test_value(
   expintegral_e(1,0.25b0),
   1.044282634443738194536438161232282251891528374744802718635140468b0,
-  64);
+  63);
 true;
 
 test_value(

Index: rtest_gamma.mac
===================================================================
RCS file: /cvsroot/maxima/maxima/tests/rtest_gamma.mac,v
retrieving revision 1.39
retrieving revision 1.40
diff -u -d -r1.39 -r1.40
--- rtest_gamma.mac	29 Oct 2009 17:49:06 -0000	1.39
+++ rtest_gamma.mac	14 Dec 2009 04:55:50 -0000	1.40
@@ -3165,8 +3165,9 @@
 
 beta_incomplete(1.0,1,z);
 1.0*z;
-beta_incomplete(1.0b0,1,z);
-1.0b0*z;
+/* Unfortunate, but we don't get exactly 1.0b0 */
+(closeto(beta_incomplete(1.0b0,1,z)/z,1b0,1b-15));
+true;
 beta_incomplete(1.0,1,1/2);
 0.5;
 beta_incomplete(1.0b0,1,1/2);

Update of /cvsroot/maxima/maxima/src
In directory sfp-cvsdas-1.v30.ch3.sourceforge.com:/tmp/cvs-serv30905/src

Modified Files:
	float.lisp 
Log Message:
Bug 2909980: optimization for log of bfloats

Optimize computation of log of bfloats by using the fact that each
bfloat is frac*2^exp, where 1/2 < frac <= 1.

src/float.lisp:
o Add new variable bigfloat_log2 to hold the current bigfloat value of
  log(2b0).
o Add function fplog2 to return log(2b0) of the correct precision.
o Modify big-float-log to compute logs using the formula
  log(frac*2^exp) = log(frac) + exp*log(2).

tests/rtest_expintegral.mac:
o Test 73 (expintegral_e(1,.25b0)) no longer is accurate to 64 digits;
  reduce accuracy to 63.

tests/rtest_gamma.mac:
o Test 636 (beta_incomplete(1b0,1,z)) is no longer exactly
  1b0*z. Change test to allow a small error.



Index: float.lisp
===================================================================
RCS file: /cvsroot/maxima/maxima/src/float.lisp,v
retrieving revision 1.57
retrieving revision 1.58
diff -u -d -r1.57 -r1.58
--- float.lisp	18 Nov 2009 19:46:42 -0000	1.57
+++ float.lisp	14 Dec 2009 04:55:49 -0000	1.58
@@ -81,6 +81,9 @@
 (defmvar bigfloat%gamma '((bigfloat simp 56.) 41592772053807304. 0)
   "Bigfloat representation of %gamma")
 
+(defmvar bigfloat_log2 '((bigfloat simp 56.) 49946518145322874. 0)
+  "Bigfloat representation of log(2)")
+
 ;; Internal specials
 
 ;; Number of bits of precision in the mantissa of newly created bigfloats.
@@ -104,6 +107,8 @@
 (defvar max-bfloat-%pi bigfloat%pi)
 (defvar max-bfloat-%e  bigfloat%e)
 (defvar max-bfloat-%gamma bigfloat%gamma)
+(defvar max-bfloat-log2 bigfloat_log2)
+
 
 (declare-top (special *cancelled $float $bfloat $ratprint $ratepsilon $domain $m1pbranch adjust))
 
@@ -714,6 +719,16 @@
 	(t
 	 (cdr (setq max-bfloat-%gamma (setq bigfloat%gamma (fpgamma1)))))))
 
+(defun fplog2 ()
+  (cond ((= fpprec (caddar bigfloat_log2))
+	 (cdr bigfloat_log2))
+	((< fpprec (caddar bigfloat_log2))
+	 (cdr (setq bigfloat_log2 (bigfloatp bigfloat_log2))))
+	((< fpprec (caddar max-bfloat-log2))
+	 (cdr (setq bigfloat_log2 (bigfloatp max-bfloat-log2))))
+	(t
+	 (cdr (setq max-bfloat-log2 (setq bigfloat_log2 (bcons (fplog (intofp 2)))))))))
+  
 (defun fpone ()
   (cond (*decfp (intofp 1))
 	((= fpprec (caddar bigfloatone)) (cdr bigfloatone))
@@ -1755,13 +1770,30 @@
   (if y
       (multiple-value-bind (u v) (complex-log x y)
 	(add (bcons u) (mul '$%i (bcons v))))
-      (let ((fp-x (cdr (bigfloatp x))))
-	(if (fplessp fp-x (intofp 0))
-	    ;; ??? Do we want to return an exact %i*%pi or a float
-	    ;; approximation?
-	    (add (bcons (fplog (fpminus fp-x)))
-		 (mul '$%i (bcons (fppi))))
-	    (bcons (fplog fp-x))))))
+      (flet ((%log (x)
+	       ;; x is (mantissa exp), where mantissa = frac*2^fpprec,
+	       ;; with 1/2 < frac <= 1 and x is frac*2^exp.  To
+	       ;; compute log(x), use log(x) = log(frac)+ exp*log(2).
+	       (cdr
+		(let* ((extra 8)
+		       (fpprec (+ fpprec extra))
+		       (log-frac
+			(fplog #+nil
+			       (cdr ($bfloat
+				     (cl-rat-to-maxima (/ (car x)
+							  (ash 1 (- fpprec 8))))))
+			       (list (ash (car x) extra) 0)))
+		       (log-exp (fptimes* (intofp (second x)) (fplog2)))
+		       (result (bcons (fpplus log-frac log-exp))))
+		  (let ((fpprec (- fpprec extra)))
+		    (bigfloatp result))))))
+	(let ((fp-x (cdr (bigfloatp x))))
+	  (if (fplessp fp-x (intofp 0))
+	      ;; ??? Do we want to return an exact %i*%pi or a float
+	      ;; approximation?
+	      (add (bcons (%log (fpminus fp-x)))
+		   (mul '$%i (bcons (fppi))))
+	      (bcons (%log fp-x)))))))
 
 (defun big-float-sqrt (x &optional y)
   (if y

Update of /cvsroot/maxima/maxima
In directory sfp-cvsdas-1.v30.ch3.sourceforge.com:/tmp/cvs-serv23591

Modified Files:
      Tag: RELEASE-5_20-BRANCH
	configure.in 
Log Message:
Version number 5.20.1.

Index: configure.in
===================================================================
RCS file: /cvsroot/maxima/maxima/configure.in,v
retrieving revision 1.117.2.1
retrieving revision 1.117.2.2
diff -u -d -r1.117.2.1 -r1.117.2.2
--- configure.in	1 Dec 2009 20:02:58 -0000	1.117.2.1
+++ configure.in	14 Dec 2009 03:16:12 -0000	1.117.2.2
@@ -2,7 +2,7 @@
 
 AC_INIT(src/macsys.lisp)
 
-AM_INIT_AUTOMAKE(maxima,5.20.0)
+AM_INIT_AUTOMAKE(maxima,5.20.1)
 
 AM_PATH_LISPDIR

Update of /cvsroot/maxima/maxima
In directory sfp-cvsdas-1.v30.ch3.sourceforge.com:/tmp/cvs-serv19431

Modified Files:
      Tag: RELEASE-5_20-BRANCH
	Makefile.am 
Log Message:
Put ChangeLog-5.20 on list of files.

Index: Makefile.am
===================================================================
RCS file: /cvsroot/maxima/maxima/Makefile.am,v
retrieving revision 1.32
retrieving revision 1.32.4.1
diff -u -d -r1.32 -r1.32.4.1
--- Makefile.am	4 Aug 2009 03:14:29 -0000	1.32
+++ Makefile.am	14 Dec 2009 01:57:27 -0000	1.32.4.1
@@ -12,7 +12,7 @@
   ChangeLog-5.9.2 ChangeLog-5.9.3 ChangeLog-5.10.0 ChangeLog-5.11.0 \
   ChangeLog-5.12.0 ChangeLog-5.13.0 ChangeLog-5.14.0 ChangeLog-5.15.0 \
   ChangeLog-5.16 ChangeLog-5.17 ChangeLog-5.17-special-functions \
-  ChangeLog-5.18 ChangeLog-5.19
+  ChangeLog-5.18 ChangeLog-5.19 ChangeLog-5.20
 
 include common.mk

Update of /cvsroot/maxima/maxima
In directory sfp-cvsdas-1.v30.ch3.sourceforge.com:/tmp/cvs-serv19002

Added Files:
      Tag: RELEASE-5_20-BRANCH
	ChangeLog-5.20 
Log Message:
Change log for Maxima 5.20 branch.

--- NEW FILE: ChangeLog-5.20 ---
   Maxima 5.20 change log
   Compiled 2009-12-13
   Thanks to Dieter Kaiser for most of this text


   Backwards-incompatible changes:

 * removed package jtroot3 (superseded by bfallroots)

 * internal symbols of specint replaced by more descriptive names:

    %f         --> hypergeometric
    %e         --> elliptic_ec
    %kelliptic --> elliptic_kc
    %p[a,b,c]  --> jacobi_p
    %p[a,b]    --> assoc_legendre_p
    %p[a]      --> legendre_p
    %q[a,b]    --> assoc_legendre_q
    %q[a]      --> legendre_q
    hstruve    --> struve_h
    lstruve    --> struve_l
    %d         --> parabolic_cylinder_d
    %m         --> whittaker_m
    %w         --> whittaker_w
    %ei        --> expintegral_ei
    %h[n,x]    --> hankel_1 and hankel_2
    %he        --> hermite


   New items in core:

 * general mechanism for functions to distribute over operators
   (distribute_over declaration)


   Other revisions:

 * improved Laplace transforms of special functions

 * assume database:

   The inferences between the declarations of variables as integer,
   real, complex, .. have been cut out. The constant %i is declared
   to be imaginary. The constants %pi, %e, %gamma, %phi are declared
   to be real.
  
   The database handles expressions with constants more completely.
   Facts like assume(x < %pi/2) are now handled correctly.
  
   The sign of sums of constant expressions is determined more
   correctly.
  
   Askinteger take into account more facts from the database.

 * improvements of special functions:

   More functionality for the Hankel 1 and 2 functions:
   Expansion for a half integral order has been implemented.

   More consistent numerical evaluation of Bessel and Hankel functions.

   More general integrals for Bessel functions and implementation of
   integrals for the Airy functions.

   Gamma function accepts a Taylor expansion as an argument.

   Implementation of Struve H and L functions:
   Numerical evaluation for real and complex numbers in float and
   bigfloat precision. Handling of specific values. Derivatives of the
   functions and expansion for a half integral order is implemented.

   Adding mirror symmetry for the Airy functions.

   Improving the limits of the log and factorial functions.

 * package fft: revised for greater speed
 * package to_poly_solver: improvements
 * regularize error and warning messages, and translate via gettext (ongoing)
 * bug fixes and other improvements in plotting code
 * bug fixes in translator


   Bug fixes:

   2906049: integration failure with option integrate_use_rootsof :true
   2901855: limit(sqrt(x),x,minf) not fully evaluated
   2886564: conjugate(atan2(y,x)) not simplified
   2882408: solve_rat_ineq returns false answers
   2876387: ode2 asks if sqrt(3) is an integer
   2876284: nicedummies of subscripted
   2876277: hgfred([3/2,-2],[5/2],-x) not fully simplified
   2873057: incorrect version is reported in manual front page
   2872605: abs_integrate bug
   2872505: Assume database inconsistent after reset()
   2847387: hgfred([3/2,-b],[5/2],-1) bogus
   2843705: limit of psi[i]
   2842198: next_prime speed improvement
   2842060: unsimplified result from integrate
   2841504: Limit of the factorial function - 4 problems
   2840566: defint fails to determine if one of its limit is real
   2836339: unsimplified negation
   2835634: logcontract broken
   2835098: SIGN-PREP strangeness
   2834336: ratsimp vs facsum
   2824360: missing bug reporting node in manual's detailed menu
   2808568: Reserved words that aren't reserved
   2802006: integrate(1/(sqrt(x)+1), x, 0, 1);
   2770575: rtestsum test 226
   2636628: solve_rat_ineq doesnt find all soloutions
   2609426: integrate(cos(a)/sqrt((tan(a))^2+1), a,-%pi/2,%pi/2);
   2148461: docs for factor refer to "the field of integers"
   2018842: unsimplified result from jacobi_p
   1977146: radexpand does not work as explained in documentation
   1725951: unsimplified boolean
   1646397: listofvars includes true
   1551310: mod and floor should distribute over matrix and list
   1038624: askinteger ignores asksign database
    840360: qunit(4) internal errors
    831163: part(x) should give warning
    816797: limit(%i*log(a),a,0) nounform (%i*und problem)
    788892: zeroa handled inconsistently
    752067: Can't untrace ?meval
    706455: Should powerseries do Laurent expansions?
    652470: Pickapart error (due to MEMSIMILAR)
    640332: Need to specdisrep more systematically
    627759: Ratdisrep of aggregates

Update of /cvsroot/maxima/maxima
In directory sfp-cvsdas-1.v30.ch3.sourceforge.com:/tmp/cvs-serv18780

Added Files:
	ChangeLog-5.20 
Log Message:
Change log for Maxima 5.20 branch.

--- NEW FILE: ChangeLog-5.20 ---
   Maxima 5.20 change log
   Compiled 2009-12-13
   Thanks to Dieter Kaiser for most of this text


   Backwards-incompatible changes:

 * removed package jtroot3 (superseded by bfallroots)

 * internal symbols of specint replaced by more descriptive names:

    %f         --> hypergeometric
    %e         --> elliptic_ec
    %kelliptic --> elliptic_kc
    %p[a,b,c]  --> jacobi_p
    %p[a,b]    --> assoc_legendre_p
    %p[a]      --> legendre_p
    %q[a,b]    --> assoc_legendre_q
    %q[a]      --> legendre_q
    hstruve    --> struve_h
    lstruve    --> struve_l
    %d         --> parabolic_cylinder_d
    %m         --> whittaker_m
    %w         --> whittaker_w
    %ei        --> expintegral_ei
    %h[n,x]    --> hankel_1 and hankel_2
    %he        --> hermite


   New items in core:

 * general mechanism for functions to distribute over operators
   (distribute_over declaration)


   Other revisions:

 * improved Laplace transforms of special functions

 * assume database:

   The inferences between the declarations of variables as integer,
   real, complex, .. have been cut out. The constant %i is declared
   to be imaginary. The constants %pi, %e, %gamma, %phi are declared
   to be real.
  
   The database handles expressions with constants more completely.
   Facts like assume(x < %pi/2) are now handled correctly.
  
   The sign of sums of constant expressions is determined more
   correctly.
  
   Askinteger take into account more facts from the database.

 * improvements of special functions:

   More functionality for the Hankel 1 and 2 functions:
   Expansion for a half integral order has been implemented.

   More consistent numerical evaluation of Bessel and Hankel functions.

   More general integrals for Bessel functions and implementation of
   integrals for the Airy functions.

   Gamma function accepts a Taylor expansion as an argument.

   Implementation of Struve H and L functions:
   Numerical evaluation for real and complex numbers in float and
   bigfloat precision. Handling of specific values. Derivatives of the
   functions and expansion for a half integral order is implemented.

   Adding mirror symmetry for the Airy functions.

   Improving the limits of the log and factorial functions.

 * package fft: revised for greater speed
 * package to_poly_solver: improvements
 * regularize error and warning messages, and translate via gettext (ongoing)
 * bug fixes and other improvements in plotting code
 * bug fixes in translator


   Bug fixes:

   2906049: integration failure with option integrate_use_rootsof :true
   2901855: limit(sqrt(x),x,minf) not fully evaluated
   2886564: conjugate(atan2(y,x)) not simplified
   2882408: solve_rat_ineq returns false answers
   2876387: ode2 asks if sqrt(3) is an integer
   2876284: nicedummies of subscripted
   2876277: hgfred([3/2,-2],[5/2],-x) not fully simplified
   2873057: incorrect version is reported in manual front page
   2872605: abs_integrate bug
   2872505: Assume database inconsistent after reset()
   2847387: hgfred([3/2,-b],[5/2],-1) bogus
   2843705: limit of psi[i]
   2842198: next_prime speed improvement
   2842060: unsimplified result from integrate
   2841504: Limit of the factorial function - 4 problems
   2840566: defint fails to determine if one of its limit is real
   2836339: unsimplified negation
   2835634: logcontract broken
   2835098: SIGN-PREP strangeness
   2834336: ratsimp vs facsum
   2824360: missing bug reporting node in manual's detailed menu
   2808568: Reserved words that aren't reserved
   2802006: integrate(1/(sqrt(x)+1), x, 0, 1);
   2770575: rtestsum test 226
   2636628: solve_rat_ineq doesnt find all soloutions
   2609426: integrate(cos(a)/sqrt((tan(a))^2+1), a,-%pi/2,%pi/2);
   2148461: docs for factor refer to "the field of integers"
   2018842: unsimplified result from jacobi_p
   1977146: radexpand does not work as explained in documentation
   1725951: unsimplified boolean
   1646397: listofvars includes true
   1551310: mod and floor should distribute over matrix and list
   1038624: askinteger ignores asksign database
    840360: qunit(4) internal errors
    831163: part(x) should give warning
    816797: limit(%i*log(a),a,0) nounform (%i*und problem)
    788892: zeroa handled inconsistently
    752067: Can't untrace ?meval
    706455: Should powerseries do Laurent expansions?
    652470: Pickapart error (due to MEMSIMILAR)
    640332: Need to specdisrep more systematically
    627759: Ratdisrep of aggregates

Update of /cvsroot/maxima/maxima/src
In directory sfp-cvsdas-1.v30.ch3.sourceforge.com:/tmp/cvs-serv14920/src

Modified Files:
	defint.lisp 
Log Message:
src/defint.lisp:
o mtosc: handle terms w/o exp(m*%i*x) produced by call to $expand
   fixes fourier integral incorrect - ID: 2875784

tests/rtestint.mac:
o add tests



Index: defint.lisp
===================================================================
RCS file: /cvsroot/maxima/maxima/src/defint.lisp,v
retrieving revision 1.68
retrieving revision 1.69
diff -u -d -r1.68 -r1.69
--- defint.lisp	14 Nov 2009 20:51:01 -0000	1.68
+++ defint.lisp	14 Dec 2009 00:57:03 -0000	1.69
@@ -779,7 +779,7 @@
   (setq *global-defint-assumptions*
 	(cons (assume '((mlessp) epsilon 1.0e-8))
 	      *global-defint-assumptions*))
-  ;; EPSILON is used in principal vaule code to denote the familiar
+  ;; EPSILON is used in principal value code to denote the familiar
   ;; mathematical entity.
   (setq *global-defint-assumptions*
 	(cons (assume '((mgreaterp) prin-inf 1.0e+8))
@@ -999,6 +999,9 @@
 	 (format t "Principal Value~%")
 	 (setq pcprntd t))))
 
+;; e is of form poly(x)*exp(m*%i*x)
+;; s is degree of denominator
+;; adds e to bptu or bptd according to sign of m
 (defun rib (e s)
   (let (*updn c)
     (cond ((or (mnump e) (constant e))
@@ -1013,9 +1016,10 @@
 		      (cond (*updn (setq bptu (cons e bptu)))
 			    (t (setq bptd (cons e bptd))))))))))
 
+;; check term is of form poly(x)*exp(m*%i*x)
+;; n is degree of denominator
 (defun ptimes%e (term n)
-  (cond ((polyinx term var nil) term)
-	((and (mexptp term)
+  (cond ((and (mexptp term)
 	      (eq (cadr term) '$%e)
 	      (polyinx (caddr term) var nil)
 	      (eq ($sign (m+ (deg ($realpart (caddr term))) -1))
@@ -1065,7 +1069,7 @@
 					   (cadr n))
 					  (t 0.))))))))
 
-
+;; exponentialize sin and cos
 (defun sconvert (e)
   (cond ((atom e) e)
 	((polyinx e var nil) e)
@@ -1107,10 +1111,19 @@
 	   (t (return (simplify (cons (list (caar e))
 				      (mapcar #'esap (cdr e)))))))))
 
+;; computes integral from minf to inf for expressions of the form
+;; exp(%i*m*x)*r(x) by residues on either the upper half
+;;		  plane or the lower half plane, depending on whether
+;;		  m is positive or negative.  [wang p. 77]
+;;
+;; exponentializes sin and cos before applying residue method.
+;; can handle some expressions with poles on real line, such as
+;; sin(x)*cos(x)/x.
 (defun mtosc (grand)
   (numden grand)
   (let ((n nn*)
 	(d dn*)
+	ratterms ratans
 	plf bptu bptd s upans downans)
     (cond ((not (or (polyinx d var nil)
 		    (and (setq grand (%einvolve d))
@@ -1127,25 +1140,32 @@
 				       (t (sconvert n)))))
 		(cond ((mplusp n)  (setq n (cdr n)))
 		      (t (setq n (list n))))
-		(do ((do-var n (cdr do-var)))
-		    ((null do-var) t)
-		  (cond ((rib (car do-var) s))
+		(dolist (term n t)
+		  (cond ((polyinx term var nil)
+			 ;; call to $expand can create rational terms
+			 ;; with no exp(m*%i*x)
+			 (setq ratterms (cons term ratterms)))
+			((rib term s))
 			(t (return nil))))
 ;;;Function RIB sets up the values of BPTU and BPTD
 		(cond ((car plf)
 		       (setq bptu (subst (car plf) 'x* bptu))
 		       (setq bptd (subst (car plf) 'x* bptd))
+		       (setq ratterms (subst (car plf) 'x* ratterms))
 		       t)	 ;CROCK, CROCK. This is TERRIBLE code.
 		      (t t))
 ;;;If there is BPTU then CSEMIUP must succeed.
 ;;;Likewise for BPTD.
-		(cond (bptu (cond ((setq upans (csemiup (m+l bptu) d var)))
-				  (t nil)))
+		(cond (ratterms (setq ratans 
+				      (defint (m// (m+l ratterms) d) var '$minf '$inf)))
+		      (t (setq ratans 0)))
+		(cond (bptu (setq upans (csemiup (m+l bptu) d var)))
 		      (t (setq upans 0)))
-		(cond (bptd (cond ((setq downans (csemidown (m+l bptd) d var)))
-				  (t nil)))
+		(cond (bptd (setq downans (csemidown (m+l bptd) d var)))
 		      (t (setq downans 0))))
-	   (sratsimp (m* '$%pi (m+ upans (m- downans))))))))
+	   
+	   (sratsimp (m+ ratans
+			 (m* '$%pi (m+ upans (m- downans)))))))))
 
 
 (defun evenfn (e var)

Update of /cvsroot/maxima/maxima/tests
In directory sfp-cvsdas-1.v30.ch3.sourceforge.com:/tmp/cvs-serv14920/tests

Modified Files:
	rtestint.mac 
Log Message:
src/defint.lisp:
o mtosc: handle terms w/o exp(m*%i*x) produced by call to $expand
   fixes fourier integral incorrect - ID: 2875784

tests/rtestint.mac:
o add tests



Index: rtestint.mac
===================================================================
RCS file: /cvsroot/maxima/maxima/tests/rtestint.mac,v
retrieving revision 1.47
retrieving revision 1.48
diff -u -d -r1.47 -r1.48
--- rtestint.mac	1 Dec 2009 17:09:50 -0000	1.47
+++ rtestint.mac	14 Dec 2009 00:57:03 -0000	1.48
@@ -1248,3 +1248,14 @@
 integrate((a*x^2+b*x+c)/(d*x^3+a*x^2+b),x),integrate_use_rootsof:true;
 'lsum((c+%r1*b+%r1^2*a)*d*log(x-%r1)/(3*%r1^2*d+2*%r1*a),%r1,
             rootsof(d*x^3+a*x^2+b));
+
+/* fourier integral incorrect - ID: 2875784 */
+integrate(cos(x)*sin(x)/x, x, minf, inf);
+%pi/2;
+
+integrate(exp(%i*x)*sin(x)/x, x, minf, inf);
+%pi/2;
+
+integrate(exp(-%i*x)*sin(x)/x, x, minf, inf);
+%pi/2;
+

Update of /cvsroot/maxima/maxima/src
In directory sfp-cvsdas-1.v30.ch3.sourceforge.com:/tmp/cvs-serv13878/src

Modified Files:
	comm.lisp 
Log Message:
Reworking the implementation of the function i-$dependicies

- Check that all arguments are symbols
- Do not add duplicates to the list of dependencies
- Overwrite existing dependencies with new dependencies

Related bug report:
 Bug ID: 893633 - depends(a,[b,b,b])
 
 No problems with the testsuite and share_testsuite

Index: comm.lisp
===================================================================
RCS file: /cvsroot/maxima/maxima/src/comm.lisp,v
retrieving revision 1.38
retrieving revision 1.39
diff -u -d -r1.38 -r1.39
--- comm.lisp	3 Nov 2009 22:34:46 -0000	1.38
+++ comm.lisp	14 Dec 2009 00:41:15 -0000	1.39
@@ -251,7 +251,8 @@
        (l))
       ((null args) (i-$dependencies (nreverse l)))
     (if ($listp (first args))
-	(mapc #'(lambda (e) (push (depends1 e (second args)) l)) (cdr (first args)))
+	(mapc #'(lambda (e) (push (depends1 e (second args)) l))
+	      (cdr (first args)))
 	(push (depends1 (first args) (second args)) l))))
 
 (defun depends1 (x y)
@@ -261,19 +262,35 @@
 (defmspec $dependencies (form)
   (i-$dependencies (cdr form)))
 
-(defmfun i-$dependencies (l)
+(defun i-$dependencies (l &aux res)
   (dolist (z l)
-    (cond ((atom z) (merror (intl:gettext "depends: argument must be a non-atomic expression; found ~M") z))
-	  ((or (eq (caar z) 'mqapply) (member 'array (cdar z) :test #'eq))
-	   (merror (intl:gettext "depends: argument cannot be a subscripted expression; found ~M") z))
-	  (t (let ((y (mget (caar z) 'depends)))
-	       (mputprop (caar z)
-			 (setq y (union* (reverse (cdr z)) y))
-			 'depends)
-	       (unless (cdr $dependencies)
-		 (setq $dependencies (copy-list '((mlist simp)))))
-	       (add2lnc (cons (cons (caar z) nil) y) $dependencies)))))
-  (cons '(mlist simp) l))
+    (cond
+      ((atom z)
+       (merror
+         (intl:gettext
+           "depends: argument must be a non-atomic expression; found ~M") z))
+      ((or (eq (caar z) 'mqapply)
+           (member 'array (cdar z) :test #'eq))
+       (merror
+         (intl:gettext
+           "depends: argument cannot be a subscripted expression; found ~M") z))
+      (t
+       (do ((zz z (cdr zz))
+            (y nil))
+           ((null zz)
+            (mputprop (caar z) (setq y (reverse y)) 'depends)
+            (setq res (push (cons (ncons (caar z)) y) res))
+            (unless (cdr $dependencies)
+              (setq $dependencies '((mlist simp))))
+            (add2lnc (cons (cons (caar z) nil) y) $dependencies))
+         (cond ((not (symbolp (cadr zz)))
+                (merror
+                  (intl:gettext "depends: argument must be a symbol; found ~M")
+                  (cadr zz)))
+               ((and (cadr zz)
+                     (not (member (cadr zz) y)))
+                (setq y (push (cadr zz) y))))))))
+  (cons '(mlist simp) (reverse res)))
 
 (defmspec $gradef (l)
   (setq l (cdr l))

Update of /cvsroot/maxima/maxima/src
In directory sfp-cvsdas-1.v30.ch3.sourceforge.com:/tmp/cvs-serv12818/src

Modified Files:
	simp.lisp 
Log Message:
src/simp.lisp:
o ordfn: don't mix comparisons of mtimes lists and mplus lists
o ordexpt: remove clause that compares expt of x directly to y
   fixes orderlessp not transitive - ID: 776441

src/rtest_lambert_w.mac:
integrate(lambert_w(1/x),x) works now



Index: simp.lisp
===================================================================
RCS file: /cvsroot/maxima/maxima/src/simp.lisp,v
retrieving revision 1.93
retrieving revision 1.94
diff -u -d -r1.93 -r1.94
--- simp.lisp	29 Nov 2009 20:35:07 -0000	1.93
+++ simp.lisp	14 Dec 2009 00:28:02 -0000	1.94
@@ -2447,6 +2447,7 @@
 	((eq e a))
 	(t (great e a))))
 
+;; compare lists a and b elementwise from back to front
 (defun ordlist (a b cx cy)
   (prog (l1 l2 c d)
      (setq l1 (length a) l2 (length b))
@@ -2460,19 +2461,28 @@
      (setq l1 (1- l1) l2 (1- l2))
      (go loop)))
 
+(defun term-list (x)
+  (if (mplusp x)
+      (cdr x)
+    (list x)))
+
+(defun factor-list (x)
+  (if (mtimesp x)
+      (cdr x)
+    (list x)))
+
+;; one of the exprs x or y should be one of:
+;; %del, mexpt, mplus, mtimes
 (defun ordfn (x y)
   (let ((cx (caar x)) (cy (caar y)) u)
     (cond ((eq cx '%del) (if (eq cy '%del) (great (cadr x) (cadr y)) t))
 	  ((eq cy '%del) nil)
-	  ((member cx '(mplus mtimes) :test #'eq)
-	   (cond ((member cy '(mplus mtimes) :test #'eq) (ordlist (cdr x) (cdr y) cx cy))
-		 ((alike1 (setq u (car (last x))) y) (not (ordhack x)))
-		 ((and (eq cx 'mplus) (eq cy 'mexpt) (mplusp (cadr y)))
-		  (not (ordmexpt y x)))
-		 (t (great u y))))
-	  ((member cy '(mplus mtimes) :test #'eq) (not (ordfn y x)))
+	  ((or (eq cx 'mtimes) (eq cy 'mtimes))
+	   (ordlist (factor-list x) (factor-list y) 'mtimes 'mtimes))
+	  ((or (eq cx 'mplus) (eq cy 'mplus))
+	   (ordlist (term-list x) (term-list y) 'mplus 'mplus))
 	  ((eq cx 'mexpt) (ordmexpt x y))
-	  (t (not (ordmexpt y x))))))	; (EQ CY 'MEXPT)
+	  ((eq cy 'mexpt) (not (ordmexpt y x))))))
 
 (defun ordhack (x)
   (if (and (cddr x) (null (cdddr x)))
@@ -2495,8 +2505,6 @@
 	       (t (let ((x1 (simpln1 x)) (y1 (simpln1 y)))
 		    (if (alike1 x1 y1) (great (cadr x) (cadr y))
 			(great x1 y1))))))
-	((maxima-constantp (cadr x))
-	 (if (alike1 (caddr x) y) t (great (caddr x) y)))
 	((alike1 (cadr x) y) (great (caddr x) 1))
 	((mnump (caddr x)) (great (cadr x) y))
 	(t (great (simpln1 x) (simpln (list '(%log) y) 1 t)))))

Update of /cvsroot/maxima/maxima/tests
In directory sfp-cvsdas-1.v30.ch3.sourceforge.com:/tmp/cvs-serv12818/tests

Modified Files:
	rtest_lambert_w.mac 
Log Message:
src/simp.lisp:
o ordfn: don't mix comparisons of mtimes lists and mplus lists
o ordexpt: remove clause that compares expt of x directly to y
   fixes orderlessp not transitive - ID: 776441

src/rtest_lambert_w.mac:
integrate(lambert_w(1/x),x) works now



Index: rtest_lambert_w.mac
===================================================================
RCS file: /cvsroot/maxima/maxima/tests/rtest_lambert_w.mac,v
retrieving revision 1.7
retrieving revision 1.8
diff -u -d -r1.7 -r1.8
--- rtest_lambert_w.mac	26 Dec 2008 22:15:20 -0000	1.7
+++ rtest_lambert_w.mac	14 Dec 2009 00:28:03 -0000	1.8
@@ -124,8 +124,8 @@
 
 /* SF bug report 2468610 the integrator loops endlessly 
    errcatch() doesn't catch the error for clisp-2.46 on cygwin */
-/* errcatch(integrate(lambert_w(1/x),x));
-[integrate(lambert_w(1/x),x)]; */
+integrate(lambert_w(1/x),x);
+'integrate(lambert_w(1/x),x);
 
 taylor(lambert_w(x),x,0,6);
 x-x^2+3*x^3/2-(8*x^4/3)+125*x^5/24-(54*x^6/5);