it supports double precision routines, they can be replaced
% by single precision routines without any problem.

%% xGESVD vs. xGESDD
%
% Speaking of performance, xGESDD is faster than xGESVD.
% If only the singular values are needed to calculate,
% both show similar speed. See the following examples:

%%
% Generate a big matrix:
A=randn(1000,1000)+i*randn(1000,1000);

%%
% Measure times when only sinular values are calculated:
tic, S=svd(A,0); toc
tic, S=svd_lapack(A,0,'gesvd'); toc
tic, S