verify avantage between fortran basic function and loop implementation
do with singular value (=s[min(m,n)-1])