~peterpall/maxima/+git/maxima.code:rtoy-bessel-bfloat

Last commit made on 2022-06-29
Get this branch:
git clone -b rtoy-bessel-bfloat https://git.launchpad.net/~peterpall/maxima/+git/maxima.code

Branch merges

Branch information

Name:
rtoy-bessel-bfloat
Repository:
lp:~peterpall/maxima/+git/maxima.code

Recent commits

dd6fc4c... by Raymond Toy <email address hidden>

Bfloat eval for bessel_y and bessel_k when order is not an integer.

Use the hypergeometric representation to evaluate bessel_y and
bessel_k when the order is not an integer. The representation is not
valid in that case.

b5186e0... by Raymond Toy <email address hidden>

Handle bfloat eval of bessel_j(-n,z) and bessel_i(-n,z)

When the order is a negative integer, these bessel functions can be
represented in terms of the bessel functions of positive order.

ac999b6... by Raymond Toy <email address hidden>

Bigfloat evaluation for bessel_j and bessel_i

The hypergeometric representation for bessel_j and bessel_i is valud
for any complex z and for any complex order (except for negative
integer order). This can be used to evaluate these functions.

a4c09ce... by Robert Dodier <email address hidden>

Clarify description of askequal and fix copy/paste errors.

344ad88... by PeterPall

Restored a comment and a catbox entry I've lost

9c0205f... by PeterPall

Documentation: Added askequal() to the manual

3ee9383... by PeterPall

builtins list: airy() is no more a maxima command

80ff6ec... by Raymond Toy <email address hidden>

Add anchor for guess_exact_value

Makeinfo was complaining that the @mref for guess_exact_value didn't
exist. Therefore add an anchor for it so that the @mref works
correctly.

73105cf... by Marius Gerbershagen

Add some more gcd tests

9e0ed65... by Marius Gerbershagen

fix gcd algorithm (special case for one polynomial linear in one variable)

Example for a wrong result which was returned previously:

(%i1) f : expand((r^2-2*m+a^2)*(a^2*s^2+r^2)^3*(a^2*s^2+r^2-2*m));
(%i2) g : expand((r^2-2*m+a^2)*(a^2*s^2+r^2)^2);
(%i3) gcd(f,g);
                             4 4 2 2 2 4
(%o3) a s + 2 a r s + r

The fix is very simple and just involves checking if rpol1 divides
rpol2, even if rpol1 and rpol2 are polynomial in different variables.