~peterpall/maxima/+git/maxima.code:rtoy-3968-zeta-neg

Last commit made on 2022-04-27
Get this branch:
git clone -b rtoy-3968-zeta-neg https://git.launchpad.net/~peterpall/maxima/+git/maxima.code

Branch merges

Branch information

Name:
rtoy-3968-zeta-neg
Repository:
lp:~peterpall/maxima/+git/maxima.code

Recent commits

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

Add more tests for zeta(-x)

Verify that float(zeta(-n/2)) is close to zeta(float(-n/2)) for
various values of odd n.

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

Fix #3968: zeta(s) for s < 0 can be expressed in terms of zeta(1-s)

zeta(s) can be written in terms of zeta(1-s), so if s < 0, we can
write it in terms of 1-s, a positive number. The relationship is

  zeta(s) = 2^s*%pi^(s-1)*sin(%pi/2*s)*gamma(1-s)*zeta(1-s))

This makes, for example, zeta(-7/2) return a result in terms of
zeta(9/2), and when you numerically evaluate it, you get the same
answer as you would get for zeta(-3.5).

This is nicer than just returning zeta(-7/2).

Add a couple of tests to rtest_zeta to verify the result.

19f30f8... by Raymond Toy <email address hidden>

Fix #3967: elliptic_e(5*%pi/4,1) incorrect

elliptic_e(5*%pi/4,1) was returning 2-1/sqrt(2), but
elliptic_e(5*%pi/4,1),numer returned 2.707... or basically
2+1/sqrt(2).

This happens because we were simplifying elliptic_e(5*%pi/4,1) to
sin(5*%pi/4) = -1/sqrt(2). But we need to reduce the arg 5*%pi/4 by
%pi so we get sin(%pi/4) = 1/sqrt(2).

Add a test for this and that the symbolic and numerical values are
consistent.

Also needed to update the test for elliptic_e(x,1) to include the
reduction of the arg for the sin part.

1e7cbee... by Raymond Toy <email address hidden>

Add example for zheev

This is based on the example given by Michel Talon on the maxima list,
https://sourceforge.net/p/maxima/mailman/message/37643060/.

I tweaked it slightly in the case of computing the eigenvectors so
that we print out the eigenvalues and the eigenvectors separately so
that we can see the column vectors more clearly.

9992d6a... by Raymond Toy <email address hidden>

Add some docstrings to functions in simp.lisp

Just add some docstrings to make it easier to understand what the
function does without having to read the code. Just to lower to
burden of understanding code.

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

Fix typo in docs for zheev

The text says "Like @code{zheev}...", but it should say dgeev, as it
does later on in this section.

f1c1976... by kjak

chmod takes options before modes

-R was treated as a file instead of an option

f1ccc67... by kjak

Fix bug #608: taylor(x^a,[x],0,1) unsimplified

Also bug #545: multivar taylor gives 1^2

This is an ancient bug that was present in Macsyma.

SETUP-MULTIVAR-DISREP gives genvars that are associated with multivars
the DISREP property so they disrep to 1. These genvars raised to powers
could be visible as unsimplified 1s raised to powers in output:

(%i1) taylor(x^a,[x],0,1);
(%o1) +1^2*x^a

(%i2) taylor((x*y)^a,[x,y],0,1);
(%o2) +1^4*(x*y)^a

(%i3) taylor(x^a+y^a,[x,y],0,1);
(%o3) +(1^2*y^a+1^2*x^a)

Prior to commit 38d4b2fd the internal genvars themselves could be made
visible by using taytorat:

(%i1) taylor(x^a,[x],0,1);
(%o1) +1^2*x^a

(%i2) taytorat(%);
(%o2) (g137^a)^2*x^a

Now we add to PDISREP! the knowledge that 1^x is equal to 1 for any x,
which prevents the construction of factors like 1^2 above. (PDISREP!
already knows that x^0 is equal to 1 and x^1 is equal to x for any x.)

No problems with the test suite or share test suite. New tests have
been added to rtest_taylor.

4522f7d... by PeterPall

Corrected the permissions for the documentation images in the tmp_html folder

If the example images for the documentation are copied using "cp" if the source tree is write-protected
the images will be, too, and cannot be deleted afterwards due to being write-protected.

the "make distcheck" target tests if all temporary files can be deleted, though, and works on a
write-protected source tree.

On my computer this would repair the "make dist" target if there were no test failures.

2892618... by PeterPall

Building the docu: Hide my debug message when not in debug mode.

I added a message to the ./build_html script that tells
the user which path which file is taken from.

If the script isn't run in debug mode we should hide
that message.