Changes commited to PARI GIT during the last two weeks

log(gtodouble()) -> dbllog2

1	1	src/basemath/mftrace.c

gtodouble(glog()) -> dbllog2

4	4	src/basemath/bibli1.c

improve comment

2	2	src/basemath/aprcl.c

gtodouble(glog()) -> dbllog2

2	2	src/basemath/QX_factor.c

gisdouble: return 0 instead of raising an error when input not of real type

conforms to the documentation

1	1	src/headers/pariinl.h

fix gisdouble documentation

4	3	doc/usersch5.tex

cleanup Fincke_Pohst_bound

5	4	src/basemath/buch2.c

heuristic initial precision increase in bnfinit to help small_snorm [#2632]

Commit 784f1ae504 removed this precision estimate
1) it had become possibly wrong when Fincke_Pohst_bound was introduced
2) it was used to reduce the precision of M to a "sufficient" value,
   potentially causing problems for a minor efficiency gain

Before that commit, the bound was valid to compute norms reliably from
floating point embeddings and nothing else. But as a side effect, it
could increase the initial nfinit precision if a t_POL was input. Which
would help in other places, such as small_norm's gaussred_from_QR which
is breaking here.

1) is not a serious concern, the value is of the right order of magnitude
and this is a one-sided error: at worse we would fail to compute largish
vectors which are unlikely to be smooth anyway.
2) can be circumvented: never reduce the precision, but possibly let the
   bound increase it.

So reinstate the bound as a better heuristic guess for the initial
precision. It won't hurt, it may help.

17	0	src/basemath/buch2.c
2	1	src/test/32/bnf
3	0	src/test/in/bnf

cxEk: use E4,E6 instead of G4,G6 in wp series + add 64 bits if k> 14

avoid useless multiplications by powers of 2pi

8	6	src/basemath/elltrans.c
2	2	src/test/32/ell

rewrite and test cxEk

use Weierstrass \wp series expansion for k> 14

27	34	src/basemath/elltrans.c
13	1	src/test/32/ell
1	0	src/test/in/ell

reorder elltrans.c + minor comments updates

1501	1498	src/basemath/elltrans.c

inline ellsigma_theta in ellsigma

10	16	src/basemath/elltrans.c

make elljacobi bench less cryptic

we test v/w against A

3	3	src/test/in/elljacobi

fix elljacobi bench

1	1	src/test/in/elljacobi

elljacobi cleanups

11	3	src/basemath/elltrans.c

restore elljacobi test that fails (not sure why)

1	1	src/test/32/elljacobi

delete useless comment

0	2	src/basemath/elltrans.c

rewrite elljacobi_pol + fix BIBs

48	43	src/basemath/elltrans.c
8	1	src/test/32/elljacobi
5	0	src/test/in/elljacobi

cleanup sigma/theta code, integrate old ellsigma_cx in ellsigma

86	156	src/basemath/elltrans.c
0	1	src/headers/paridecl.h
13	13	src/test/32/ell
3	3	src/test/32/ellnf
3	3	src/test/32/theta

thetaall_ii: don't negate thetanullwerte when z is reduced

This change fixes the following bug in ellsigma_cx:

? install(ellsigma_cx,"GGD0,L,p")
? ellsіgma_cx([2*I, 1], 4/3*I)
%1 = -0.55506326344653125121966647544195623187*I

? ellsigma([2*I,1], 4/3*I)
%2 = 0.55506326344653125121966647544195623186*I

For unknown reasons, it introduces a problem in elljacobi bench.
To be fixed later.

This commit is a prerequisite for the theta / elliptic functions cleanup
patch which follows (otherwise ellsigma becomes completely broken).

1	5	src/basemath/elltrans.c

cleanup elljacobi bench

8	10	src/test/in/elljacobi

ellmanintable typo: u may be a t_FRAC

1	1	src/modules/ellmanin.c

ellmanintable: avoid useless ellchangecurve

6	3	src/modules/ellmanin.c

ellmaninconstant: add flag: use Cremona table

1	0	CHANGES
17	1	src/basemath/ellisog.c
8	4	src/functions/elliptic_curves/ellmaninconstant
2	1	src/headers/paridecl.h
4	0	src/headers/paripriv.h
146	0	src/modules/ellmanin.c
3	39	src/test/32/ellmanin
23	2	src/test/in/ellmanin

ZV_lcm: use gen_product directly

7	1	src/basemath/polarit2.c

test-help: update

1	0	src/test/32/help

doc: ?\: add \qf

1	0	src/language/gplib.c

issquare: use Zn_issquare

1	1	src/basemath/ispower.c

ellwp_cx, ellzeta_cx: simplify by assuming input is reduced

+ remove useless GC

18	19	src/basemath/elltrans.c

doc typo

1	1	doc/usersch1.tex

30- bnfsunit: could fail to detect non-S-units

flawed logic in sunit_val when S contained some of the primes above p
but not all

1	0	CHANGES
24	22	src/basemath/bnfunits.c
2	1	src/test/32/bnfsunit
3	1	src/test/in/bnfsunit

checkom: remove unused parameter z

5	5	src/basemath/elltrans.c

update bench

5	5	src/test/32/compat

matker flag is obsolete

Don't use matker0: either ker or ZM_ker directly

1	3	src/basemath/alglin1.c
4	6	src/functions/linear_algebra/matker