[GiNaC-list] Bug with evalm()?

Martin Sandve Alnæs martinal at simula.no
Mon Feb 12 11:09:23 CET 2007


On 2/12/07, Chris Dams <Chris.Dams at mi.infn.it> wrote:
>
> Dear Martin,
>
> On Sun, 11 Feb 2007, [ISO-8859-1] Martin Sandve Alnæs wrote:
>
> > ex m = series_to_poly(s).evalm();
>
> Can't you simply call your own C++-written tree transversal algorithm at
> this point to check whether there actually are matrix products somewhere
> and then only call evalm to evaluate them?

Hi Chris,
that's a good idea. I'll try that.

Otherwise, I think I can live with a try-catch, and use the original
expression if evalm() fails.

Thanks for the ideas everyone, I'll let you know when I've worked it out.

martin


More information about the GiNaC-list mailing list