wavelet1.rb in gsl-2.1.0

- old
+ new

@@ -3,11 +3,13 @@
 
 n = 256
 nc = 20
 
 data = GSL::Vector.alloc(n)
-data.fscanf("ecg.dat")
+File.new(File.dirname(__FILE__) + '/ecg.dat').each_with_index do |l, i|
+  data[i] = l.to_f
+end
 
 w = GSL::Wavelet.alloc("daubechies", 4)
 work = GSL::Wavelet::Workspace.alloc(n)
 
 # Choose as you like...
@@ -44,7 +46,7 @@
 #data3 = GSL::Wavelet.transform(w, data2, GSL::Wavelet::BACKWARD, work)
 #data3 = GSL::Wavelet.transform(w, data2, GSL::Wavelet::BACKWARD)
 #data3 = GSL::Wavelet.transform_inverse(w, data2, work)
 #data3 = GSL::Wavelet.transform_inverse(w, data2)
 
-GSL::graph(nil, data, data3, "-T X -C -g 3 -x 0 #{data.size} -L 'Red: data, Green: DWT'")
+# GSL::graph(nil, data, data3, "-T X -C -g 3 -x 0 #{data.size} -L 'Red: data, Green: DWT'")
 