Typo in the specification of FFPACK's rowEchelonForm() function

[Charles Bouillaguet]

Dear all,

I think that there's a mistake in the comments on lines [761:765] of
ffpack.h that describe the row echelon functions.


    /*
     * If LuTag == FfpackTileRecursive, then after the computation A =
[ L \ M ]
     * such that X A = R is a row echelon decomposition of A,
     * with X =  [ L  0   ] P  and R = M + [Ir 0] Q^T
     *           [    In-r]
     * If LuTag == FfpackTileRecursive then A = [ L \ N ] such that the
same holds with M =  N Q^T
     */

The  "If LuTag == FfpackTileRecursive else if LuTag ==
FfpackTileRecursive" is suspicious. One of these should be
FfpackSlabRecursive.

In addition, the description of the result seems weird. It suggests that
M is not permuted in the first case, but permuted in the second case. Is
that right?

Cheers,

--

Charles Bouillaguet