> Maxima 5.17post http://maxima.sourceforge.net
> Using Lisp Clozure Common Lisp Version 1.3-dev-r11537M-trunk
> (WindowsX8632)
> (%i1) to_lisp();
> MAXIMA> (load "bessel.lisp")
> MAXIMA> (setq jvals (make-array 4 :element-type
> 'double-float :initial-element 0.0d0))
>
> #(0.0 0.0 0.0 0.0)
>
> MAXIMA> (slatec:dbesj 2.0 0.0 4 jvals 0)
>
> NIL
> NIL
> NIL
> NIL
> 0
> MAXIMA> jvals
> #(0.0 0.0 0.0 0.0)
>
> Is this a Clozure bug? I still do not understand how calling slatec:dbesj
> has the side-effect of
> modifying the array jvals. But from what I understand, the complex case
> (zbesj) also works by
> side-effects on arrays--that code works at least sometimes with Clozure CL.
This seems to work on a Mac (I downloaded a binary):
Maxima 5.17post http://maxima.sourceforge.net
Using Lisp Clozure Common Lisp Version 1.2-r11300M (DarwinX8664)
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) to_lisp();
Type (to-maxima) to restart, ($quit) to quit Maxima.
MAXIMA> (setq jvals (make-array 4 :element-type
'double-float :initial-element 0.0d0))
#(0.0 0.0 0.0 0.0)
MAXIMA> (slatec:dbesj 2.0 0.0 4 jvals 0)
NIL
NIL
NIL
NIL
0
MAXIMA> jvals
#(0.22389077914123579 0.5767248077568736 0.35283402861563784
0.12894324947440208)
I only get one failed test in the testsuite (rtest_gamma):
********************** Problem 658 ***************
Input:
5 - 3
closeto(float(rectform(beta_incomplete(2, -, %i + ---)))
2 2
- beta_incomplete(2.0, 2.5, %i - 1.5), 0.0, 2.8E-15)
Result:
2.908381734874007E-15
Andrej