[GiNaC-list] pole_error when expanding power into series

Alexei Sheplyakov alexei.sheplyakov at gmail.com
Tue Nov 2 15:35:00 CET 2010


On Tue, Nov 2, 2010 at 1:34 PM, Camille Gillot <k1000.jlo at wanadoo.fr> wrote:

> However, ginac 1.5.8 is unable to expand expressions like :
> cos(x)^(sin(x)/x), even if it expands well sin(x)/x.

> +       if(!is_a<numeric>(exponent)) {
> +               ex new_exponent = (exponent*log(basis)).series(r,order,options);
> +               return exp(new_exponent).series(r,order,options);
> +       }

I think such transformation is way too aggressive (the exponent might not
depend on x, there might be no singularity at all, etc). I'll try to
fix this in a
more careful way.

Best regards,

More information about the GiNaC-list mailing list