mical form of SVD:
[U,S,V] = svd_lapack(A,0)


%%
% Select the lapack subroutine ZGESVD:
[U,S,V,err] = svd_lapack(A,0,'gesvd')

%%
% Select the lapack subroutine ZGESDD:
[U,S,V,err] = svd_lapack(A,0,'gesdd')

%% Diagnostics
% 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/lu