add svdvals example

Author: perazz

example/linalg/CMakeLists.txt

@@ -21,5 +21,6 @@ ADD_EXAMPLE(lapack_getrf)
 ADD_EXAMPLE(lstsq1)
 ADD_EXAMPLE(lstsq2)
 ADD_EXAMPLE(svd1)
+ADD_EXAMPLE(svd2)
 ADD_EXAMPLE(determinant)
 ADD_EXAMPLE(determinant2)

example/linalg/example_svd2.f90

@@ -0,0 +1,27 @@
+! Singular Values
+program example_svd2
+  use stdlib_linalg_constants, only: dp
+  use stdlib_linalg, only: svd
+  implicit none
+
+  real(dp), allocatable :: A(:,:),s(:)
+
+  ! We want to find the singular values of matrix: 
+  ! 
+  !    A = [ 3  2  2]
+  !        [ 2  3 -2]  
+  !   
+  A = transpose(reshape([ 3, 2, 2, &
+                          2, 3,-2], [3,2]))
+
+  ! Get singular values
+  s = svdvals(A)
+
+  ! Singular values: [5, 3]
+  print 1, '    '
+  print 1, 'S = ',s
+  print 1, '    '
+  
+  1 format(a,*(1x,f12.8))
+
+end program example_svd2