Contents

#!/usr/bin/env ruby
require("gsl")

A = GSL::Matrix[[3, 5, 2], [6, 2, 1], [4, 7, 3]]
p A
u, v, s = A.SV_decomp
p u.class
p v.class
p s.class

# u and v are orthonormal
p u*u.trans
p v*v.trans

# Reconstruct the matrix A from u, v, and s.
p u*GSL::Matrix.diagonal(s)*v.trans

Version data entries

22 entries across 22 versions & 4 rubygems

Version	Path
gsl-2.1.0.3	examples/linalg/SV.rb
gsl-2.1.0.2	examples/linalg/SV.rb
gsl-2.1.0.1	examples/linalg/SV.rb
gsl-2.1.0	examples/linalg/SV.rb
gsl-1.16.0.6	examples/linalg/SV.rb
rb-gsl-1.16.0.5	examples/linalg/SV.rb
rb-gsl-1.16.0.4	examples/linalg/SV.rb
rb-gsl-1.16.0.3	examples/linalg/SV.rb
rb-gsl-1.16.0.3.rc1	examples/linalg/SV.rb
gsl-nmatrix-1.15.3.2	examples/linalg/SV.rb
gsl-nmatrix-1.15.3.1	examples/linalg/SV.rb
rb-gsl-1.16.0.2	examples/linalg/SV.rb
rb-gsl-1.16.0.1	examples/linalg/SV.rb
rb-gsl-1.16.0	examples/linalg/SV.rb
rb-gsl-1.15.3.2	examples/linalg/SV.rb
rb-gsl-1.15.3.1	examples/linalg/SV.rb
blackwinter-gsl-1.15.3.2	examples/linalg/SV.rb
gsl-nmatrix-1.15.3.0	examples/linalg/SV.rb
gsl-1.15.3	examples/linalg/SV.rb
gsl-1.14.7	examples/linalg/SV.rb