Advanced Search

Question: Row echelon form of sparse matrices

Question:Row echelon form of sparse matrices

mak 8 Maple
November 05 2006
0
I have a 661x662 entry matrix with <= 3 non-zero entries per row. The entries are small integers, non-zero ones being stored in a list s. I initialize the matrix by

m:=Matrix(661,662,s,storage=sparse,datatype=integer[1]);

and then ask for

a:=ReducedRowEchelonFrom(m);

This operation takes half an hour. Is there any way to speed it up, eg. using the NAG routines?

Thanks,
Mak
Please Wait...