diff options
| -rw-r--r-- | src/interp/compat.boot (renamed from src/interp/compat.boot.pamphlet) | 43 | ||||
| -rw-r--r-- | src/interp/sfsfun.boot (renamed from src/interp/sfsfun.boot.pamphlet) | 159 | ||||
| -rw-r--r-- | src/interp/word.boot (renamed from src/interp/word.boot.pamphlet) | 22 |
3 files changed, 77 insertions, 147 deletions
diff --git a/src/interp/compat.boot.pamphlet b/src/interp/compat.boot index 2db27b0a..87e1227c 100644 --- a/src/interp/compat.boot.pamphlet +++ b/src/interp/compat.boot @@ -1,30 +1,7 @@ -\documentclass{article} -\usepackage{axiom} -\begin{document} -\title{\$SPAD/src/interp compat.boot} -\author{The Axiom Team} -\maketitle -\begin{abstract} -\end{abstract} -\eject -\tableofcontents -\eject -\section{readline} -Reads a line, defaulting stream to *standard-input* - -This was commented out in the NAG version, probably due to the -use of the saturn GUI interface but the real reason is uncertain. -<<readline>>= -READLINE(:s) == - s => read_-line(first s) - read_-line(_*STANDARD_-INPUT_*) - -@ -\section{compat.boot} -\section{License} -<<license>>= -- Copyright (c) 1991-2002, The Numerical ALgorithms Group Ltd. -- All rights reserved. +-- Copyright (C) 2007, Gabriel Dos Reis. +-- All rights reserved. -- -- Redistribution and use in source and binary forms, with or without -- modification, are permitted provided that the following conditions are @@ -54,9 +31,6 @@ READLINE(:s) == -- NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS -- SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. -@ -<<*>>= -<<license>> import '"pathname" )package "BOOT" @@ -82,7 +56,7 @@ isSharpVarWithNum x == -- RREAD which takes erroval to return if key is missing rread(key,rstream,errorval) == if IDENTP key then key := PNAME key - RREAD(key,rstream,errorval) + RREAD(key,rstream,errorval) rwrite(key,val,stream) == if IDENTP key then key := PNAME key @@ -105,10 +79,7 @@ editFile file == makeBigFloat(mantissa,expon) == [$BFtag,mantissa,:expon] -<<readline>> -@ -\eject -\begin{thebibliography}{99} -\bibitem{1} nothing -\end{thebibliography} -\end{document} +READLINE(:s) == + s => read_-line(first s) + read_-line(_*STANDARD_-INPUT_*) + diff --git a/src/interp/sfsfun.boot.pamphlet b/src/interp/sfsfun.boot index 50bd7b5b..c6c28253 100644 --- a/src/interp/sfsfun.boot.pamphlet +++ b/src/interp/sfsfun.boot @@ -1,52 +1,7 @@ -\documentclass{article} -\usepackage{axiom} -\begin{document} -\title{\$SPAD/src/interp sfsfun.boot} -\author{The Axiom Team} -\maketitle -\begin{abstract} -\end{abstract} -\eject -\tableofcontents -\eject -\begin{verbatim} -NOTEfrom TTT: at least BesselJAsymptOrder needs work - -1. This file contains the contents of BWC's original files: - floaterrors.boot - floatutils.boot - rgamma.boot - cgamma.boot - rpsi.boot - cpsi.boot - f01.boot - chebf01cmake.boot - chebevalsf.boot - besselIJ.boot - -2. All declarations have been commented out with "--@@" - since the boot translator is generating bad lisp code from them. - -3. The functions PsiAsymptotic, PsiEps and PsiAsymptoticOrder - had inconpatible definitions in rpsi.boot and cpsi.boot -- - the local variables were declared float in one file and COMPLEX in - the other. The type declarations have been commented out and the - duplicate definitions have been deleted. - -4. BesselIJ was not compiling. I have modified the code from that - file to make it compile. It should be checked for correctness. - -SMW June 25, 1991 - -"Fixes" to BesselJ, B. Char June 14, 1992. Needs extensive testing and - further fixes to BesselI and BesselJ. -More fixes to BesselJ, T. Tsikas 24 Feb, 1995. - -\end{verbatim} -\section{License} -<<license>>= -- Copyright (c) 1991-2002, The Numerical ALgorithms Group Ltd. -- All rights reserved. +-- Copyright (C) 2007, Gabriel Dos Reis. +-- All rights reserved. -- -- Redistribution and use in source and binary forms, with or without -- modification, are permitted provided that the following conditions are @@ -76,9 +31,41 @@ More fixes to BesselJ, T. Tsikas 24 Feb, 1995. -- NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS -- SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. -@ -<<*>>= -<<license>> + +-- NOTEfrom TTT: at least BesselJAsymptOrder needs work + +-- 1. This file contains the contents of BWC's original files: +-- floaterrors.boot +-- floatutils.boot +-- rgamma.boot +-- cgamma.boot +-- rpsi.boot +-- cpsi.boot +-- f01.boot +-- chebf01cmake.boot +-- chebevalsf.boot +-- besselIJ.boot + +-- 2. All declarations have been commented out with "--@@" +-- since the boot translator is generating bad lisp code from them. + +-- 3. The functions PsiAsymptotic, PsiEps and PsiAsymptoticOrder +-- had inconpatible definitions in rpsi.boot and cpsi.boot -- +-- the local variables were declared float in one file and COMPLEX in +-- the other. The type declarations have been commented out and the +-- duplicate definitions have been deleted. + +-- 4. BesselIJ was not compiling. I have modified the code from that +-- file to make it compile. It should be checked for correctness. + +-- SMW June 25, 1991 + +-- "Fixes" to BesselJ, B. Char June 14, 1992. Needs extensive testing and +-- further fixes to BesselI and BesselJ. +-- More fixes to BesselJ, T. Tsikas 24 Feb, 1995. + + + -- Used to be SPECFNSF )package "BOOT" @@ -119,7 +106,7 @@ horner(l,x) == rgamma (x) == if COMPLEXP(x) then FloatError('"Gamma not implemented for complex value ~D",x) - ZEROP (x-1.0) => 1.0 + ZEROP (x-1.0) => 1.0 if x>20 then gammaStirling(x) else gammaRatapprox(x) lnrgamma (x) == @@ -742,12 +729,12 @@ BesselJ(v,z) == n := 50 --- number of terms in Chebychev series. --- tests for negative integer order (FLOATP(v) and ZEROP fracpart(v) and (v<0)) or (COMPLEXP(v) and ZEROP IMAGPART(v) and ZEROP fracpart(REALPART(v)) and REALPART(v)<0.0) => - --- odd or even according to v (9.1.5 A&S) - --- $J_{-n}(z)=(-1)^{n} J_{n}(z)$ + --- odd or even according to v (9.1.5 A&S) + --- $J_{-n}(z)=(-1)^{n} J_{n}(z)$ BesselJ(-v,z)*EXPT(-1.0,v) (FLOATP(z) and (z<0)) or (COMPLEXP(z) and REALPART(z)<0.0) => --- negative argument (9.1.35 A&S) - --- $J_{\nu}(z e^{m \pi i}) = e^{m \nu \pi i} J_{\nu}(z)$ + --- $J_{\nu}(z e^{m \pi i}) = e^{m \nu \pi i} J_{\nu}(z)$ BesselJ(v,-z)*EXPT(-1.0,v) ZEROP z and ((FLOATP(v) and (v>=0.0)) or (COMPLEXP(v) and ZEROP IMAGPART(v) and REALPART(v)>=0.0)) => --- zero arg, pos. real order @@ -764,29 +751,29 @@ BesselJ(v,z) == w := 2.0*arg vp1 := v+1.0 [sum,arr] := chebf01coefmake(vp1,w,n) - ---if we get NaNs then half n - while not _=(sum,sum) repeat - n:=FLOOR(n/2) + ---if we get NaNs then half n + while not _=(sum,sum) repeat + n:=FLOOR(n/2) [sum,arr] := chebf01coefmake(vp1,w,n) - ---now n is safe, can we increase it (we know that 2*n is bad)? + ---now n is safe, can we increase it (we know that 2*n is bad)? chebstarevalarr(arr,arg/w,n)/cgamma(vp1)*EXPT(z/2.0,v) true => BesselJRecur(v,z) FloatError('"BesselJ not implemented for ~S", [v,z]) BesselJRecur(v,z) == - -- boost order + -- boost order --Numerical.Recipes. suggest so:=v+sqrt(n.s.f.^2*v) - so:=15.0*z - -- reduce order until non-zero + so:=15.0*z + -- reduce order until non-zero while ZEROP ABS(BesselJAsymptOrder(so,z)) repeat so:=so/2.0 - if ABS(so)<ABS(z) then so:=v+18.*SQRT(v) - m:= FLOOR(ABS(so-v))+1 - w:=MAKE_-ARRAY(m) - SETF(AREF(w,m-1),BesselJAsymptOrder(v+m-1,z)) + if ABS(so)<ABS(z) then so:=v+18.*SQRT(v) + m:= FLOOR(ABS(so-v))+1 + w:=MAKE_-ARRAY(m) + SETF(AREF(w,m-1),BesselJAsymptOrder(v+m-1,z)) SETF(AREF(w,m-2),BesselJAsymptOrder(v+m-2,z)) for i in m-3 .. 0 by -1 repeat SETF(AREF(w,i), 2.0 * (v+i+1.0) * AREF(w,i+1) /z -AREF(w,i+2)) - AREF(w,0) + AREF(w,0) BesselI(v,z) == B1 := 15.0 @@ -952,24 +939,24 @@ BesselJAsymptOrder(v,z) == ca := 1.0/tanhalpha Pi := PI - ca2:=ca*ca - ca4:=ca2*ca2 - ca8:=ca4*ca4 + ca2:=ca*ca + ca4:=ca2*ca2 + ca8:=ca4*ca4 EXP(-v*(alpha-tanhalpha))/SQRT(2.0*Pi*v*tanhalpha)*_ - (1.0+_ - horner([ -5.0, 3.0],_ - ca2)*ca/(v*24.0)+_ - horner([ 385.0, -462.0, 81.0],_ - ca2)*ca2/(1152.0*v*v)+_ - horner([ -425425.0, 765765.0, -369603.0, 30375.0],_ - ca2)*ca2*ca/(414720.0*v*v*v)+_ - horner([ 185910725.0, -446185740.0, 349922430.0, -94121676.0, 4465125.0],_ - ca2)*ca4/(39813120.0*v*v*v*v)+_ - horner([ -188699385875.0, 566098157625.0, -614135872350.0, 284499769554.0, -49286948607.0, 1519035525.0],_ - ca2)*ca4*ca/(6688604160.0*v*v*v*v*v)+_ - horner([1023694168371875.0,-3685299006138750.0,5104696716244125.0,-3369032068261860.0,1050760774457901.0,-127577298354750.0,2757049477875.0],_ - ca2)*ca4*ca2/(4815794995200.0*v*v*v*v*v*v)) - + (1.0+_ + horner([ -5.0, 3.0],_ + ca2)*ca/(v*24.0)+_ + horner([ 385.0, -462.0, 81.0],_ + ca2)*ca2/(1152.0*v*v)+_ + horner([ -425425.0, 765765.0, -369603.0, 30375.0],_ + ca2)*ca2*ca/(414720.0*v*v*v)+_ + horner([ 185910725.0, -446185740.0, 349922430.0, -94121676.0, 4465125.0],_ + ca2)*ca4/(39813120.0*v*v*v*v)+_ + horner([ -188699385875.0, 566098157625.0, -614135872350.0, 284499769554.0, -49286948607.0, 1519035525.0],_ + ca2)*ca4*ca/(6688604160.0*v*v*v*v*v)+_ + horner([1023694168371875.0,-3685299006138750.0,5104696716244125.0,-3369032068261860.0,1050760774457901.0,-127577298354750.0,2757049477875.0],_ + ca2)*ca4*ca2/(4815794995200.0*v*v*v*v*v*v)) + --- See Olver, p. 376-382. BesselIAsymptOrder(v,vz) == @@ -1023,9 +1010,3 @@ BesselKAsymptOrder (v,vz) == -@ -\eject -\begin{thebibliography}{99} -\bibitem{1} nothing -\end{thebibliography} -\end{document} diff --git a/src/interp/word.boot.pamphlet b/src/interp/word.boot index ac76dca3..95dfc7a1 100644 --- a/src/interp/word.boot.pamphlet +++ b/src/interp/word.boot @@ -1,16 +1,3 @@ -\documentclass{article} -\usepackage{axiom} -\begin{document} -\title{\$SPAD/src/interp word.boot} -\author{The Axiom Team} -\maketitle -\begin{abstract} -\end{abstract} -\eject -\tableofcontents -\eject -\section{License} -<<license>>= -- Copyright (c) 1991-2002, The Numerical ALgorithms Group Ltd. -- All rights reserved. -- @@ -42,9 +29,6 @@ -- NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS -- SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. -@ -<<*>>= -<<license>> --======================================================================= -- Build Directories @@ -414,9 +398,3 @@ obSearch x == [y for i in 0..MAXINDEX OBARRAY() | (IDENTP (y := vec.i) or CVEC y) and match?(pattern,COPY y)] -@ -\eject -\begin{thebibliography}{99} -\bibitem{1} nothing -\end{thebibliography} -\end{document} |