From: Stefan Weinzierl Date: Tue, 27 Jul 2021 18:55:54 +0000 (+0200) Subject: Bug fix in modular_form_kernel::is_numeric and modular_form_kernel::Laurent_series... X-Git-Tag: release_1-8-1~3 X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?a=commitdiff_plain;h=d5a419b8517e3103087cceb3f6a9585b60e46543;p=ginac.git Bug fix in modular_form_kernel::is_numeric and modular_form_kernel::Laurent_series: Call the correct method to compute the series expansion. This bug has been reported by Christoph Dlapa. --- diff --git a/ginac/integration_kernel.cpp b/ginac/integration_kernel.cpp index 6571ecdf..5e440519 100644 --- a/ginac/integration_kernel.cpp +++ b/ginac/integration_kernel.cpp @@ -2012,12 +2012,12 @@ bool modular_form_kernel::is_numeric(void) const symbol qbar("qbar"); // test with a random number and random expansion - return series_to_poly(P.series(qbar,18)).subs(qbar==numeric(1,937)).evalf().info(info_flags::numeric); + return series_to_poly(q_expansion_modular_form(qbar,18)).subs(qbar==numeric(1,937)).evalf().info(info_flags::numeric); } ex modular_form_kernel::Laurent_series(const ex & qbar, int order) const { - ex res = series_to_poly(P.series(qbar,order+1)); + ex res = series_to_poly(q_expansion_modular_form(qbar,order+1)); res = C_norm * res/qbar; res = res.series(qbar,order); return res;