Fast Reduced Row Echelon Form

A much faster version of MATLAB's RREF

Description

The algorithm is exactly the same as MATLAB's RREF, only difference being that vectorization is used instead of for-loop in the original code to subtract multiples of the pivot row
from all other rows. A typical speed-up range is about 2-4 times of the original code. However, the actual speed-up depends on the size of A. The speed-up is quite considerable (in some cases more than 1000 times) if the number of columns in A far surpass its number of rows.

• MATLAB release: MATLAB 7.4 (R2007a)

[Yi Cao]

It is indead very fast compared to native MATLAB rref. However, I just wonder when the function is most useful.

[Nikola Toljic]

Very good.

[Babak]

Very fast and robust. It may not be a bad idea to preallocate 'jb' and then delete the empty elements. Otherwise, awesome!