mak

8 Reputation

2 Badges

15 years, 241 days

MaplePrimes Activity


These are questions asked by mak

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
Page 1 of 1