Description
For large matrices, this algorithm is often much faster.
i1 : M = random(RR^200, RR^200);
200 200
o1 : Matrix RR <-- RR
53 53
|
i2 : time SVD(M);
-- used 0.0446999s (cpu); 0.0446969s (thread); 0s (gc)
|
i3 : time SVD(M, DivideConquer=>true);
-- used 0.0436347s (cpu); 0.0436391s (thread); 0s (gc)
|