January 2022 - GiNaC-devel

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-727-g9c577d1f
by Stefan Weinzierl 15 Jan '22

15 Jan '22

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The branch, master has been updated via 9c577d1ff8a40e821d214ce875587edfd4325f23 (commit) from f9450a587900dbd7265b7a1b9b3250a6d574bc13 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 9c577d1ff8a40e821d214ce875587edfd4325f23 Author: Stefan Weinzierl <weinzierl(a)uni-mainz.de> Date: Sat Jan 15 10:08:07 2022 +0100 exam_inifcns_nstdsums.cpp: prec set to 10*pow(10,-Digits) in inifcns_test_HLi to avoid sporadic make check failure. ----------------------------------------------------------------------- Summary of changes: check/exam_inifcns_nstdsums.cpp | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

sporadic make check failure in exam_inifcns_nstdsums
by Richard B. Kreckel 15 Jan '22

15 Jan '22

Hi Stefan, A recent FTBFS for GiNaC-1.8.2 on Debian's mipsel port [1] turned up a little issue with the precision goal not being 100% reproducible. The recursive calls in GiNaC::map_trafo_H_1overx::operator() assemble symbolic functions before evalf()'ing and these are not deterministic. This makes numeric results slightly differ between calls. As a result, regression test suite exam_inifcns_nstdsums.cpp fails in about 1.6% of all cases. It only ever seems to affect H(lst{-2,1,3},numeric(245)/100) - -Li(lst{2,1,3},lst{-numeric(245)/100,-1,1}) in inifcns_test_HLi(), exam_inifcns_nstdsums.cpp:148. I tortured my computer with 1365 calls and here is a histogram of numerical differences between the H and the Li results: [0.00000...1.0e-18[ 4 okay [1.0e-18...1.0e-17[ 530 okay [2.0e-17...3.0e-17[ 503 okay [3.0e-17...4.0e-17[ 232 okay [4.0e-17...5.0e-17[ 74 okay [5.0e-17...6.0e-17[ 18 FAILURE! [6.0e-17...7.0e-17[ 2 FAILURE! [7.0e-17...8.0e-17[ 2 FAILURE! >= 8e-17 - N/A Would it make sense to raise prec to 10*pow(10,-Digits)? If so, for which tests? Or do you have a better idea? All my best, -richard. [1] <https://buildd.debian.org/status/fetch.php?pkg=ginac&arch=mipsel&ver=1.8.2-…> -- Richard B. Kreckel <https://in.terlu.de/~kreckel/>

2 2

Re: [GiNaC-devel] [GiNaC-list] to_polynomial function does not return a polynomial
by Feng Feng 14 Jan '22

14 Jan '22

Dear Vladimir, Thanks so much for checking this, and thanks for the quick update on the git! Best regards! Feng 在 2022年1月13日 22:35，Vladimir V. Kisil<V.Kisil(a)leeds.ac.uk> 写道： Dear Feng, Thanks for spotting this typo. Indeed currently the following lines symbol x("x"); ex a = pow(x+1/x,2); cout << a << endl; exmap mp1, mp2; ex p = a.to_polynomial(mp1); cout << " = " << p << "\n with " << mp1 << endl; p = a.expand().to_polynomial(mp2); cout << " = " << p << "\n with " << mp2 << endl; produce the output with not a polynomial in the first instance: (x^(-1)+x)^2 = (x^(-1)+x)^2 with {} = 2+symbol3^2+x^2 with {symbol3==x^(-1)} If to_rational() is replaced by to_polynomial() as you suggested, then we do have (although slightly different) polynomial outputs in both cases: (x+x^(-1))^2 = (x+symbol3)^2 with {symbol3==x^(-1)} = 2+x^2+symbol4^2 with {symbol4==x^(-1)} I am attaching a one-line patch to fix this and copy the message to the developer list. Best wishes, Vladimir -- Vladimir V. Kisil http://www.maths.leeds.ac.uk/~kisilv/ Book: Geometry of Mobius Maps https://doi.org/10.1142/p835 Soft: Geometry of cycles http://moebinv.sourceforge.net/ Jupyter notebooks: https://github.com/vvkisil?tab=repositories >>>>> On Tue, 11 Jan 2022 07:10:40 +0000, Feng Feng <f.feng(a)outlook.com> said: FF> Dear Vladimir V. Kisil, FF> Thanks very much for suggestion to use expand() method on e3 FF> first. FF> And I have given a look at the source code at FF> power::to_polynomial FF> https://www.ginac.de/reference/normal_8cpp_source.html#l02688 FF> ... FF> and I wonfer if there is a typo on line 2691, the FF> basis.to_rational can be replaced with basis.to_polynomial? FF> Best regards! Feng

1 0

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-726-gf9450a58
by Richard B. Kreckel 13 Jan '22

13 Jan '22

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The branch, master has been updated via f9450a587900dbd7265b7a1b9b3250a6d574bc13 (commit) from 052d2da5b41ac82855d3578050125aa21caf1003 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit f9450a587900dbd7265b7a1b9b3250a6d574bc13 Author: Vladimir V. Kisil <kisilv(a)maths.leeds.ac.uk> Date: Thu Jan 13 20:59:21 2022 +0100 Fix power::to_polynomial() for posint exponents. This makes things like pow(x+1/x,2).to_polynomial(repl) actually produce (x+sym)^2 with sym==1/x in repl, as it should. Reported by Feng Feng <f.feng(a)outlook.com>. ----------------------------------------------------------------------- Summary of changes: ginac/normal.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

Re: [GiNaC-devel] [GiNaC-list] to_polynomial function does not return a polynomial
by Vladimir V. Kisil 13 Jan '22

13 Jan '22

Dear Feng, Thanks for spotting this typo. Indeed currently the following lines symbol x("x"); ex a = pow(x+1/x,2); cout << a << endl; exmap mp1, mp2; ex p = a.to_polynomial(mp1); cout << " = " << p << "\n with " << mp1 << endl; p = a.expand().to_polynomial(mp2); cout << " = " << p << "\n with " << mp2 << endl; produce the output with not a polynomial in the first instance: (x^(-1)+x)^2 = (x^(-1)+x)^2 with {} = 2+symbol3^2+x^2 with {symbol3==x^(-1)} If to_rational() is replaced by to_polynomial() as you suggested, then we do have (although slightly different) polynomial outputs in both cases: (x+x^(-1))^2 = (x+symbol3)^2 with {symbol3==x^(-1)} = 2+x^2+symbol4^2 with {symbol4==x^(-1)} I am attaching a one-line patch to fix this and copy the message to the developer list. Best wishes, Vladimir -- Vladimir V. Kisil http://www.maths.leeds.ac.uk/~kisilv/ Book: Geometry of Mobius Maps https://doi.org/10.1142/p835 Soft: Geometry of cycles http://moebinv.sourceforge.net/ Jupyter notebooks: https://github.com/vvkisil?tab=repositories >>>>> On Tue, 11 Jan 2022 07:10:40 +0000, Feng Feng <f.feng(a)outlook.com> said: FF> Dear Vladimir V. Kisil, FF> Thanks very much for suggestion to use expand() method on e3 FF> first. FF> And I have given a look at the source code at FF> power::to_polynomial FF> https://www.ginac.de/reference/normal_8cpp_source.html#l02688 FF> ... FF> and I wonfer if there is a typo on line 2691, the FF> basis.to_rational can be replaced with basis.to_polynomial? FF> Best regards! Feng

1 0

GiNaC 1.8.2 released
by Richard B. Kreckel 01 Jan '22

01 Jan '22

Happy New Year! GiNaC 1.8.2 has just been released. Here's what's new: * Fix elusive bug in comparing relational objects. * Ensure modular_form_kernel::series() includes an Order term. Happy hacking! -richy. -- Richard B. Kreckel <https://in.terlu.de/~kreckel/>

1 0

[SCM] GiNaC -- a C++ library for symbolic computations tag, release_1-8-2, created. release_1-4-0-725-g052d2da5
by Richard B. Kreckel 01 Jan '22

01 Jan '22

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The tag, release_1-8-2 has been created at 052d2da5b41ac82855d3578050125aa21caf1003 (commit) - Log ----------------------------------------------------------------- commit 052d2da5b41ac82855d3578050125aa21caf1003 Author: Richard Kreckel <kreckel(a)ginac.de> Date: Sat Jan 1 17:08:13 2022 +0100 Finalize 1.8.2 release. Oh, and Happy New Year! ----------------------------------------------------------------------- hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-725-g052d2da5
by Richard B. Kreckel 01 Jan '22

01 Jan '22

1 0