Bug fix in modular_form_kernel::is_numeric and modular_form_kernel::Laurent_series...
authorStefan Weinzierl <weinzierl@uni-mainz.de>
Tue, 27 Jul 2021 18:55:54 +0000 (20:55 +0200)
committerStefan Weinzierl <weinzierl@uni-mainz.de>
Tue, 27 Jul 2021 18:55:54 +0000 (20:55 +0200)
This bug has been reported by Christoph Dlapa.

ginac/integration_kernel.cpp

index 6571ecdf8e39648d67dc6c64da971c22106e65ff..5e4405190cf20cb4bc4c25a8b6b170875b453374 100644 (file)
@@ -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;