Re: Bug in expm

[Jaroslav Hajek]

OK, I see what was the intention. I'll make the changeset tomorrow.

On Thu, Apr 24, 2008 at 4:36 PM, Marco Caliari <address@hidden> wrote:
>
> > I've tried to convert these two diff patches to a hg changeset for
> > John's convenience, but
> > the dMatrix.cc patch failed. Did you make the patches against the most
> > recent sources?
> > If not, can you please refresh?
> >
>
>  I made them against octave-3.0.1. Now I'm away for ten days...
>  The dMatrix.cc.patch should exchange the lines
>
>  // leading permutations in forward order
>  for (octave_idx_type i = nc - 1; i>= ihi; i--)
>
>  with
>
>  // trailing permutations must be done in reverse order
>  for (octave_idx_type i = 0; i < (ilo-1); i++)
>
>  Marco
>



-- 
RNDr. Jaroslav Hajek
computing expert
Aeronautical Research and Test Institute (VZLU)
Prague, Czech Republic
url: www.highegg.matfyz.cz