Interface to SVD computation subroutine

subroutine svdcomp(A,S,V,svord,m,n,isw,error)

A :

Input matrix stored as linear array in row major order. On output contains matrix U in row major order

S :

Output array of singular values

V :

Output matrix V as linear array in row major order

svord :

Indices pointing to elements of S in decreasing order, and to corresponding columns of A and V

isw :

Number of sweeps needed to achieve convergence

error :

0 if convergence is achieved, 1 otherwise