You can subscribe to this list here.
2001 
_{Jan}

_{Feb}

_{Mar}

_{Apr}

_{May}

_{Jun}

_{Jul}

_{Aug}

_{Sep}

_{Oct}

_{Nov}
(20) 
_{Dec}
(17) 

2002 
_{Jan}
(39) 
_{Feb}
(21) 
_{Mar}
(33) 
_{Apr}
(135) 
_{May}
(53) 
_{Jun}
(88) 
_{Jul}
(47) 
_{Aug}
(59) 
_{Sep}
(207) 
_{Oct}
(40) 
_{Nov}
(7) 
_{Dec}
(26) 
2003 
_{Jan}
(49) 
_{Feb}
(39) 
_{Mar}
(117) 
_{Apr}
(50) 
_{May}
(62) 
_{Jun}
(6) 
_{Jul}
(19) 
_{Aug}
(24) 
_{Sep}
(11) 
_{Oct}
(11) 
_{Nov}
(49) 
_{Dec}
(9) 
2004 
_{Jan}
(29) 
_{Feb}
(123) 
_{Mar}
(32) 
_{Apr}
(53) 
_{May}
(52) 
_{Jun}
(19) 
_{Jul}
(33) 
_{Aug}
(10) 
_{Sep}
(76) 
_{Oct}
(86) 
_{Nov}
(171) 
_{Dec}
(163) 
2005 
_{Jan}
(147) 
_{Feb}
(121) 
_{Mar}
(120) 
_{Apr}
(126) 
_{May}
(120) 
_{Jun}
(213) 
_{Jul}
(76) 
_{Aug}
(79) 
_{Sep}
(140) 
_{Oct}
(83) 
_{Nov}
(156) 
_{Dec}
(202) 
2006 
_{Jan}
(181) 
_{Feb}
(171) 
_{Mar}
(157) 
_{Apr}
(98) 
_{May}
(96) 
_{Jun}
(97) 
_{Jul}
(193) 
_{Aug}
(76) 
_{Sep}
(130) 
_{Oct}
(63) 
_{Nov}
(196) 
_{Dec}
(253) 
2007 
_{Jan}
(256) 
_{Feb}
(293) 
_{Mar}
(276) 
_{Apr}
(258) 
_{May}
(181) 
_{Jun}
(91) 
_{Jul}
(108) 
_{Aug}
(69) 
_{Sep}
(107) 
_{Oct}
(179) 
_{Nov}
(137) 
_{Dec}
(121) 
2008 
_{Jan}
(124) 
_{Feb}
(129) 
_{Mar}
(192) 
_{Apr}
(201) 
_{May}
(90) 
_{Jun}
(86) 
_{Jul}
(115) 
_{Aug}
(142) 
_{Sep}
(49) 
_{Oct}
(91) 
_{Nov}
(95) 
_{Dec}
(218) 
2009 
_{Jan}
(230) 
_{Feb}
(149) 
_{Mar}
(118) 
_{Apr}
(72) 
_{May}
(77) 
_{Jun}
(68) 
_{Jul}
(102) 
_{Aug}
(72) 
_{Sep}
(89) 
_{Oct}
(76) 
_{Nov}
(125) 
_{Dec}
(86) 
2010 
_{Jan}
(75) 
_{Feb}
(90) 
_{Mar}
(89) 
_{Apr}
(121) 
_{May}
(111) 
_{Jun}
(66) 
_{Jul}
(75) 
_{Aug}
(66) 
_{Sep}
(66) 
_{Oct}
(166) 
_{Nov}
(121) 
_{Dec}
(73) 
2011 
_{Jan}
(74) 
_{Feb}

_{Mar}

_{Apr}
(14) 
_{May}
(22) 
_{Jun}
(31) 
_{Jul}
(53) 
_{Aug}
(37) 
_{Sep}
(23) 
_{Oct}
(25) 
_{Nov}
(31) 
_{Dec}
(28) 
2012 
_{Jan}
(18) 
_{Feb}
(11) 
_{Mar}
(32) 
_{Apr}
(17) 
_{May}
(48) 
_{Jun}
(37) 
_{Jul}
(23) 
_{Aug}
(54) 
_{Sep}
(15) 
_{Oct}
(11) 
_{Nov}
(19) 
_{Dec}
(22) 
2013 
_{Jan}
(11) 
_{Feb}
(32) 
_{Mar}
(24) 
_{Apr}
(37) 
_{May}
(31) 
_{Jun}
(14) 
_{Jul}
(26) 
_{Aug}
(33) 
_{Sep}
(40) 
_{Oct}
(21) 
_{Nov}
(36) 
_{Dec}
(84) 
2014 
_{Jan}
(23) 
_{Feb}
(20) 
_{Mar}
(27) 
_{Apr}
(24) 
_{May}
(31) 
_{Jun}
(27) 
_{Jul}
(3) 
_{Aug}

_{Sep}

_{Oct}

_{Nov}

_{Dec}

S  M  T  W  T  F  S 



1

2

3

4

5

6
(5) 
7
(5) 
8
(1) 
9

10
(1) 
11

12

13

14
(1) 
15
(2) 
16
(1) 
17
(3) 
18

19

20

21

22

23
(1) 
24
(1) 
25
(1) 
26

27
(1) 
28
(12) 
29
(5) 
30

31



From: Raymond Toy <rtoy@us...>  20021008 17:06:29

Update of /cvsroot/maxima/maxima/src In directory uswprcvs1:/tmp/cvsserv12070/src Modified Files: float.lisp Log Message: o Fix bug in fpshift wherein right shifts of negative numbers were not truncated properly o Add some comments on the algorithm used by fplog. o Move the multiplication by 2 that was inside the loop to outside the loop. Index: float.lisp =================================================================== RCS file: /cvsroot/maxima/maxima/src/float.lisp,v retrieving revision 1.3 retrieving revision 1.4 diff u d r1.3 r1.4  float.lisp 27 Jun 2002 17:37:20 0000 1.3 +++ float.lisp 8 Oct 2002 17:06:25 0000 1.4 @@ 575,13 +575,22 @@ FPPREC)) (RETURN (FPSHIFT L (MINUS *M))))))) +;; Compute (* L (expt d n)) where D is 2 or 10 depending on +;; *decfp. Throw away an fractional part by truncating to zero. (DEFUN FPSHIFT (L N)  ;; It seems that BIGLSH is broken so we replace it by ASH  ;(COND ((NULL *DECFP) (BIGLSH L N))  (COND ((NULL *DECFP) (ASH L N))  ((GREATERP N 0.) (TIMES L (EXPT 10. N)))  ((LESSP N 0.) (QUOTIENT L (EXPT 10. (MINUS N))))  (T L))) + (COND ((NULL *DECFP) + (cond ((and (minusp n) (minusp l)) + ;; Left shift of negative number requires some + ;; care. (That is, (truncate l (expt 2 n)), but use + ;; shifts instead.) + ( (ash ( l) n))) + (t + (ash l n)))) + ((GREATERP N 0.) + (TIMES L (EXPT 10. N))) + ((LESSP N 0.) + (QUOTIENT L (EXPT 10. (MINUS N)))) + (T L))) ;; Bignum LSH  N is assumed (and declared above) to be a fixnum. ;; This isn't really LSH, since the sign bit isn't propagated when @@ 1005,11 +1014,38 @@ (COND (MINUS (ADD A (MUL '$%I ($BFLOAT '$%PI)))) (T A))) NIL)) +;;; Computes the log of a bigfloat number. +;;; +;;; Uses the series +;;; +;;; log(1+x) = sum((x/(x+2))^(2*n+1)/(2*n+1),n,0,inf); +;;; +;;; +;;; INF x 2 n + 1 +;;; ==== () +;;; \ x + 2 +;;; = 2 >  +;;; / 2 n + 1 +;;; ==== +;;; n = 0 +;;; +;;; +;;; which converges for x > 0. +;;; +;;; Note that FPLOG is given 1+X, not X. +;;; +;;; However, to aid convergence of the series, we scale 1+x until 1/e +;;; < 1+x <= e. +;;; (DEFUN FPLOG (X) (PROG (OVER TWO ANS OLDANS TERM E) (IF (NOT (GREATERP (CAR X) 0)) (MERROR "Nonpositive argument to FPLOG"))  (SETQ E (FPE) OVER (FPQUOTIENT (FPONE) E) ANS 0) + (SETQ E (FPE) + OVER (FPQUOTIENT (FPONE) E) + ANS 0) + ;; Scale X until 1/e < X <= E. ANS keeps track of how + ;; many factors of E were used. Set X to NIL if X is E. (DO () (NIL) (COND ((EQUAL X E) (SETQ X NIL) (RETURN NIL)) ((AND (FPLESSP X E) (FPLESSP OVER X)) @@ 1020,21 +1056,27 @@ (T (SETQ ANS (ADD1 ANS)) (SETQ X (FPQUOTIENT X e))))) (COND ((NULL X) (RETURN (INTOFP (ADD1 ANS)))))  (SETQ X (FPDIFFERENCE X (FPONE)) ANS (INTOFP ANS)) + ;; Prepare X for the series. The series is for 1 + x, so + ;; get x from our X. TERM is (x/(x+2)). X becomes + ;; (x/(x+2))^2. + (SETQ X (FPDIFFERENCE X (FPONE)) + ANS (INTOFP ANS)) (SETQ X (FPEXPT (SETQ TERM (FPQUOTIENT X (FPPLUS X (SETQ TWO (INTOFP 2))))) 2))  (DO ((N 0 (f1+ N))) ((EQUAL ANS OLDANS)) + ;; Sum the series until the sum (in ANS) doesn't change + ;; anymore. + (DO ((N 1 (+ N 2))) + ((EQUAL ANS OLDANS)) (SETQ OLDANS ANS)  (SETQ  ANS  (FPPLUS  ANS  (FPTIMES* TWO (FPQUOTIENT TERM (INTOFP (f1+ (f* 2 N))))))) + (SETQ ANS + (FPPLUS + ANS + (FPQUOTIENT TERM (INTOFP N)))) (SETQ TERM (FPTIMES* TERM X)))  (RETURN ANS))) + (RETURN (fptimes* two ANS)))) (DEFUN MABSBIGFLOAT (L) (PROG (R) 