proposal to move bffac.lisp to archive (keeping bffac.mac)

Hello,

I would like to suggest that share/numeric/bffac.lisp 
be moved to the archive. 

There are two files with nearly identical content,
.../bffac.lisp and .../bffac.mac. It appears that bffac.lisp
was translated from a .mac file. A comment at the top of 
bffac.mac states that it was created by executing fundef on
the functions loaded from bffac.lisp. 

It appears that bffac.mac has been slightly modified since
it was dumped from the bffac.lisp functions. (1) bffac.lisp
contains an apparently-obsolete Maxima construct, dolist.
(2) bffac.lisp  refers to an undefined function bfac in 
some places, where bffac.mac refers to bffac, which is defined.

It appears, then, that bffac.mac is the more up-to-date of
the two files. To avoid confusion and focus attention, 
I propose that bffac.lisp be moved to the archive. I'll
carry out that task unless there are objections.

regards,
Robert Dodier

PS. I compared the bffac.lisp and bffac.mac by loading
the functions from each, simplifying the function bodies,
dumping the functions (as displayed) to a file, and 
diff'ing the files. Aside from (1) and (2) above, there
are some expressions which are equivalent but did not
simplify to the same thing for some reason. Here is the diff.

--- dispfun-bffac.lisp  2004-12-30 11:07:16.000000000 -0700
+++ dispfun-bffac.mac   2004-12-30 11:05:59.000000000 -0700
@@ -43,15 +43,15 @@
  else (if s = 0 then - - else (bzeta(1 - s, fpprec) %pi      2
                        2
 
-                       %pi s
- bfac(- s, fpprec) sin(-----), bfloat(%%)))) $ 
-                         2
+                        %pi s
+ bffac(- s, fpprec) sin(-----), bfloat(%%)))) $ 
+                          2
 vonschtoonk(p) := apply(*, subst([* = lambda([[l]], 1), 
 
                ^ = lambda([[l]], 1)], factor(divrlst(p) + 1))) $ 
 divrlst(p) := (args(expand(ratsimp(subst(^ = 
 
-                          b/2 + 1
+                          1 + b/2
                1 - part(a)                  2
 lambda([a, b], ------------------), factor(p ))))), ev(%%, part))
                   1 - part(a)
@@ -97,13 +97,13 @@
   s          s
  p          k
 
-                       1
- else (if s = 0 then - - else dolist(bzeta(1 - s, fpprec)
+                       1                               s - 1  s
+ else (if s = 0 then - - else (bzeta(1 - s, fpprec) %pi      2
                        2
 
-    s - 1  s                       %pi s
- %pi      2  bfac(- s, fpprec) sin(-----), bfloat(%%))) $ 
-                                     2
+                        %pi s
+ bffac(- s, fpprec) sin(-----), bfloat(%%))) $ 
+                          2
 bzeta(s, fpp) := if s > 0 then block([fpprec : 7, m : 0, 
 
                                             1        s
@@ -158,7 +158,7 @@
             p
 
                                s - 1  s
- else (bzeta(1 - s, fpprec) %pi      2  bfac(- s, fpprec)
+ else (bzeta(1 - s, fpprec) %pi      2  bffac(- s, fpprec)
 
      %pi s
  sin(-----), bfloat(%%))) $ 
@@ -266,9 +266,9 @@
  - m))) $ 
 bfzeta(s, fpp) := if s > 0 then block([fpprec : 7, m : 0, 
 
-                                            1        s
-t : 0.0B0, p, k : 2 max(1, entier(ev(- (s - -) log(-----)
-                                            2      2 %pi
+                                      1            s
+t : 0.0B0, p, k : 2 max(1, entier(ev((- - s) log(-----)
+                                      2          2 %pi
 
                             3
  + log(10) fpp - log(%pi) + -, numer)/2))], 
@@ -318,7 +318,7 @@
             p
 
                                 s - 1  s
- else (bfzeta(1 - s, fpprec) %pi      2  bfac(- s, fpprec)
+ else (bfzeta(1 - s, fpprec) %pi      2  bffac(- s, fpprec)
 
      %pi s
  sin(-----), bfloat(%%))) $





		
__________________________________ 
Do you Yahoo!? 
Meet the all-new My Yahoo! - Try it today! 
http://my.yahoo.com