Ray,
I have been working on Airy functions for complex numerical values
using the SLATEC routines in maxima. Some of the translated code
looks wrong, as the subroutines return (values nil nil nil ... nil).
(This is with the translated code in maxima CVS, but the latest
f2cl gives the same result).
For example zseri.lisp and zbinu.lisp need to following changes to
work. There are some other cases as well, but it is bed time here.
David
Index: zseri.lisp
===================================================================
RCS file: /cvsroot/maxima/maxima/src/numerical/slatec/zseri.lisp,v
retrieving revision 1.2
diff -u -r1.2 zseri.lisp
--- zseri.lisp 8 May 2002 13:34:34 -0000 1.2
+++ zseri.lisp 11 May 2005 14:18:54 -0000
@@ -237,5 +237,5 @@
(setf nz (f2cl-lib:int-sub nz))
(go end_label)
end_label
- (return (values nil nil nil nil nil nil nil nz nil nil nil)))))
+ (return (values zr zi fnu kode n yr yi nz tol elim alim)))))
DABilling@calttws92 /usr/local/src/maxima/src/numerical/slatec
$ cvs diff -u zbinu.lisp
Index: zbinu.lisp
===================================================================
RCS file: /cvsroot/maxima/maxima/src/numerical/slatec/zbinu.lisp,v
retrieving revision 1.2
diff -u -r1.2 zbinu.lisp
--- zbinu.lisp 8 May 2002 13:34:34 -0000 1.2
+++ zbinu.lisp 11 May 2005 14:19:10 -0000
@@ -136,5 +136,4 @@
(if (= nw -2) (setf nz -2))
(go end_label)
end_label
- (return (values nil nil nil nil nil nil nil nz nil nil nil nil nil)))))
-
+ (return (values zr zi fnu kode n cyr cyi nz rl fnul tol elim alim)))))
NOTICE
This e-mail and any attachments are private and confidential and may contain privileged information. If you are not an authorised recipient, the copying or distribution of this e-mail and any attachments is prohibited and you must not read, print or act in reliance on this e-mail or attachments.
This notice should not be removed.