ead/250196</a> But, it is not easy to call the lapack properly, because it needs a good knowledge of fortran work space.</p><p><b>SVD_LAPACK</b> 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.</p><h2>xGESVD vs. xGESDD<a name="17"></a></h2><p>Speaking of performance, xGESDD is faster than xGESVD. If only the si