X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=check%2Fseries_expansion.cpp;h=42170c0b1104323621bb357c0c6deb763d623642;hp=bd7cbbb348e9da37928f6ef194fef6046d927f38;hb=1a21a1a5d97950126a8df787302286b6cc1a902c;hpb=84806997c064f7429cbec8b446b91445430fb1fb diff --git a/check/series_expansion.cpp b/check/series_expansion.cpp index bd7cbbb3..42170c0b 100644 --- a/check/series_expansion.cpp +++ b/check/series_expansion.cpp @@ -31,7 +31,7 @@ static symbol x("x"); static unsigned check_series(const ex &e, const ex &point, const ex &d, int order = 8) { ex es = e.series(x, point, order); - ex ep = static_cast(es.bp)->convert_to_poly(); + ex ep = ex_to_pseries(es).convert_to_poly(); if (!(ep - d).is_zero()) { clog << "series expansion of " << e << " at " << point << " erroneously returned " << ep << " (instead of " << d @@ -131,8 +131,28 @@ static unsigned series3(void) return result; } -// Series of special functions +// Order term handling static unsigned series4(void) +{ + unsigned result = 0; + ex e, d; + + e = 1 + x + pow(x, 2) + pow(x, 3); + d = Order(1); + result += check_series(e, 0, d, 0); + d = 1 + Order(x); + result += check_series(e, 0, d, 1); + d = 1 + x + Order(pow(x, 2)); + result += check_series(e, 0, d, 2); + d = 1 + x + pow(x, 2) + Order(pow(x, 3)); + result += check_series(e, 0, d, 3); + d = 1 + x + pow(x, 2) + pow(x, 3); + result += check_series(e, 0, d, 4); + return result; +} + +// Series of special functions +static unsigned series5(void) { unsigned result = 0; ex e, d; @@ -190,6 +210,7 @@ unsigned series_expansion(void) result += series2(); result += series3(); result += series4(); + result += series5(); if (!result) { cout << " passed ";