Sha256: f0d3d83685017606a7a660581f691fbb42ed247a9f52a04c79b9df78dc9ce52e

Contents?: true

Size: 930 Bytes

Versions: 16

Compression:

Stored size: 930 Bytes

Contents

#!/usr/bin/env ruby
# Turn on warnings
$-w = true

require 'irb/xmp'
require 'gsl'

# Apparently, IRB::Frame has a bug that prevents the defaults from working, so
# an XMP instance must be created explicitly this way instead of using the
# otherwise convenient xmp method.
XMP.new(IRB::Frame.top(-1)).puts <<END
# Setup constants
N = 1000
DECIMATE1 = 10
DECIMATE2 = 100

# Setup random number generator
r = GSL::Rng.alloc

# Create Vector of x values
x0 = GSL::Vector.linspace(0, 20, N)

# Data: Bessel function + noise
y0 = GSL::Sf::bessel_J0(x0) + GSL::Ran::gaussian(r, 0.1, N)

# Decimate y0 by DECIMATE1
y1 = y0.decimate(DECIMATE1)

# Decimate y0 by DECIMATE2
y2 = y0.decimate(DECIMATE2)

# Create Vectors of decimated x values
x1 = GSL::Vector.linspace(0, 20, N/DECIMATE1)
x2 = GSL::Vector.linspace(0, 20, N/DECIMATE2)

# y1 and y2 are shifted vertically for visual purpose
GSL::graph([x0, y0], [x1, y1-1], [x2, y2-2])  
END

Version data entries

16 entries across 16 versions & 4 rubygems

Version Path
rb-gsl-1.16.0.4 examples/vector/decimate.rb
rb-gsl-1.16.0.3 examples/vector/decimate.rb
rb-gsl-1.16.0.3.rc1 examples/vector/decimate.rb
gsl-nmatrix-1.15.3.2 examples/vector/decimate.rb
gsl-nmatrix-1.15.3.1 examples/vector/decimate.rb
rb-gsl-1.16.0.2 examples/vector/decimate.rb
rb-gsl-1.16.0.1 examples/vector/decimate.rb
rb-gsl-1.16.0 examples/vector/decimate.rb
rb-gsl-1.15.3.2 examples/vector/decimate.rb
rb-gsl-1.15.3.1 examples/vector/decimate.rb
blackwinter-gsl-1.15.3.2 examples/vector/decimate.rb
gsl-nmatrix-1.15.3.0 examples/vector/decimate.rb
gsl-1.15.3 examples/vector/decimate.rb
gsl-1.14.7 examples/vector/decimate.rb
gsl-1.14.6 examples/vector/decimate.rb
gsl-1.14.5 examples/vector/decimate.rb