1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
|
#!/usr/bin/env ruby
require 'rubygems'
require 'narray'
require 'gsl'
require '../gsl_test.rb'
include GSL::Test
dbleps = 1e-6
expected = Math.sqrt((0..4).inject {|m,x| m+=x*x})
v = GSL::Vector.indgen(5)
v_dnrm2 = GSL::Blas.dnrm2(v)
GSL::Test.test_rel(v_dnrm2, expected, dbleps, "GSL::Blas.dnrm2(GSL::Vector)")
na = NArray.float(5).indgen!
na_dnrm2 = GSL::Blas.dnrm2(na)
GSL::Test.test_rel(na_dnrm2, expected, dbleps, "GSL::Blas.dnrm2(NArray)")
GSL::Test.test_rel(na_dnrm2, v_dnrm2, 0, "GSL::Blas.dnrm2(NArray) == GSL::Blas.dnrm2(GSL::Vector)")
|