next | previous | forward | backward | up | top | index | toc | Macaulay2 web site
Macaulay2Doc :: SVD(..., DivideConquer => ...)

SVD(..., DivideConquer => ...) -- Use the lapack divide and conquer SVD algorithm

Synopsis

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.220478 seconds
i3 : time SVD(M, DivideConquer=>true);
GC Warning: Repeated allocation of very large block (appr. size 1298432):
        May lead to memory leak and poor performance.
     -- used 0.18531 seconds

Further information