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

w = GSL::Eigen::Gensymm.alloc(4)


a = GSL::Matrix.alloc([0.18, 0.60, 0.57, 0.96], [0.41, 0.24, 0.99, 0.58],
                  [0.14, 0.30, 0.97, 0.66], [0.51, 0.13, 0.19, 0.85])
b = a.clone
v = GSL::Vector.alloc(4)

w.gensymm(a, b)
w.gensymm(a, b, v)
GSL::eigen_gensymm(a, b)
GSL::Eigen::Gensymm::gensymm(a, b, w)