SVD of a matrix based on Lapack interface
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.
Citar como
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. Recuperado .
Compatibilidad con la versión de MATLAB
Compatibilidad con las plataformas
Windows macOS LinuxCategorías
- MATLAB > Mathematics > Linear Algebra > Eigenvalues >
Etiquetas
Agradecimientos
Inspirado por: lapack
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!Descubra Live Editor
Cree scripts con código, salida y texto formateado en un documento ejecutable.
html/
Versión | Publicado | Notas de la versión | |
---|---|---|---|
1.1 | update examples. |
|
|
1.0.0.0 |