SVD of a matrix based on Lapack interface

Version 1.1 (28.1 KB) by Sung-Eun Jo
SVD_LAPACK computes the singular value decomposition of a matrix by calling LAPACK subroutines.
1.1K Downloads
Updated 28 Oct 2014

View License

SVD_LAPACK provides a SVD routine based on ZGESDD,
which is not supported as Matlab's built-in function, at least, explicitly.
Lapack, a fortran computational library, has two different subroutines
for the Singular Value Decompostion (SVD): xGESVD and xGESDD.
xGESVD is based on an implicit QR iteration and xGESDD uses a
divide-and-conquer approach.
See <http://www.netlib.org/lapack/lug/node32.html> and
<http://www.netlib.org/lapack/lug/node53.html> for Lapack subroutines.
Matlab's built-in function svd seems to use the lapack subroutine xGESVD.
Meanwhile, Octave's built-in function svd support both algorithms
by using svd_driver().

If people want to use xGESDD routines in Matlab,
it is required to write a MEX-file to call the routine.
See <http://www.mathworks.co.kr/matlabcentral/newsreader/view_thread/250196>
But, it is not easy to call the lapack properly, because it needs a good
knowledge of fortran work space.

In this submission, SVD_LAPACK function provides both routines and the lapack interfaces.
It could be an example for calling SVD rountine from Lapack interface.
Although it supports double precision routines, they can be replaced
by single precision routines without any problem.

Cite As

Sung-Eun Jo (2024). SVD of a matrix based on Lapack interface (https://www.mathworks.com/matlabcentral/fileexchange/48288-svd-of-a-matrix-based-on-lapack-interface), MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R2012b
Compatible with any release
Platform Compatibility
Windows macOS Linux
Categories
Find more on Eigenvalues in Help Center and MATLAB Answers
Acknowledgements

Inspired by: lapack

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

html/

Version Published Release Notes
1.1

update examples.

1.0.0.0