pspline 5.0.2 → 5.0.3

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: aa8139dbd0860ead8217998e30e2e6a47668139d
-  data.tar.gz: 0da507577b79ca733d4200a9e0a433877d40cf00
+  metadata.gz: cb96468c596ae7908fdd735008f4a5a22a68a56b
+  data.tar.gz: 7239450009f59fcda3844416ec1a7b575e8aa2ad
 SHA512:
-  metadata.gz: 72f9669ca12d4c7b71dbe765a13813403c981668198a2106a40606b138abb6a35fe4b7359be858597abd18b7e99da8ddf22392d0c894407951e38ce651d1d24d
-  data.tar.gz: d183a49bb333a423ebeac8a98c5dce9436192645aa7de31f5f60a8c77f96958e90310629d92d1fcc290f196b8e55b9de71c44a06257e0e3c67d38872821b244a
+  metadata.gz: a4e07bd3ba5932852d9bdcfa5ec22b7b1826a276b5830b3dab7f639f2a6fd51a9dbdf48c6faacd33cc3a7e4c9298777df5159834095b412210988a26510faf14
+  data.tar.gz: 78858c5c8f9ca31ee18551cfee1877a165ba909ec60f5dbedb3c338398f912224cb372fdc6bd998333edc05b2ce851780d0711ea5f657922087385c86ff6d32e

data/README.md

@@ -22,6 +22,13 @@ Or install it yourself as:
 
 See pspline.wiki at https://github.com/vsp2old/pspline .
 
+## Revision history
+
+5.0.0   Upload for test.   
+5.0.1	Upload for test.   
+5.0.2	Upload for release.   
+5.0.3	Add class PSPLINE::Rfft, class PSPLINE::Cfft.   
+
 ## Development
 
 After checking out the repo, run `bin/setup` to install dependencies. You can also run `bin/console` for an interactive prompt that will allow you to experiment.

data/ext/pspline/example/exfspline.rb

@@ -2,53 +2,31 @@
 require 'pspline'
 include PSPLINE
 
-=begin
-
-1 varable M dimension parametric interpolation
-
-【Module】		PSPLINE
-【Module Function】
-
-（1）obj = PSPLINE#fft_complex_transform([[x1,y1],...,[xn,yn]], f)
-	obj => FFT:[[u1,...,un],[v1,...,vn]]
-	:1 list of complex data
-	:2 type 1 => forward, 0 => backword, -1 => inverse
-
-（2）obj = PSPLINE#fft_complex_get([[u1,...,un],[v1,...,vn]], t)
-	obj => [x,y]
-	:1 FFT_complex data
-	:2 data point ( -N/2 <= t <= N/2 + N%2 )
-
-（3）obj = PSPLINE#fft_complex_bspline([[u1,...,un],[v1,...,vn]], j)
-	obj => PSPLINE#Bspline
-	:1 FFT_complex data
-	:2 dimension
-
-=end
-
 puts "# complex FFT data points"
 
 XY = [ [0, 1],[0, 1],[0, 1],[1, 0,],[1, 0,],[1, 0],[1, 0],[0, 1],[0, 1],[0, 1] ]
 
-fs = fft_complex_transform(XY, 1);
-fs[0].each {|x| printf("% .2f ", x) }
+fs = Cfft.new(XY)
+fs.real.each {|x| printf("% .2f ", x) }
 puts
-fs[1].each {|y| printf("% .2f ", y) }
+fs.imag.each {|y| printf("% .2f ", y) }
 puts
 
+printf "Real :"
 -5.upto(5) {|t|
-	w = fft_complex_get(fs, t)
+	w = fs[t]
 	printf "% .2f ", w[0]
 }
 puts
+printf "Imag :"
 -5.upto(5) {|t|
-	w = fft_complex_get(fs, t)
+	w = fs[t]
 	printf "% .2f ", w[1]
 }
 puts
 
-vv = []
-bs = fft_complex_bspline(fs, 1) {|t, z| vv.push t}
+bs = fs.spline(1)
+vv = fs.axis
 
 puts "# Interpolation points"
 
@@ -57,7 +35,7 @@ s = bs.plot(vv, 4) do |a, b|
 end
 #STDERR.puts s
 
-st = fft_complex_transform(fs, -1)
+st = fs.inverse
 printf "["
 st.each {|x| printf "% .1f ", x[0]}
 printf "]\n["

data/ext/pspline/example/exfspline1.rb

@@ -6,25 +6,25 @@ printf "# real FFT data points\n"
 
 X = [ 0, 0, 0, 1, 1, 1, 1, 0, 0, 0 ]
 
-fs = fft_real_transform(X, 1);
-fs.each {|x| printf("% .2f ", x) }
+fs = Rfft.new(X);
+fs.real.each {|x| printf("% .2f ", x) }
 puts
 
 printf "Real :"
 -5.upto(5) {|t|
-	w = fft_real_get(fs, t)
+	w = fs[t]
 	printf "% .2f ", w[0]
 }
 puts
 printf "Imag :"
 -5.upto(5) {|t|
-	w = fft_real_get(fs, t)
+	w = fs[t]
 	printf "% .2f ", w[1]
 }
 puts
 
-vv = []
-bs = fft_real_bspline(fs, 3) {|t, z| vv.push t}
+bs = fs.spline(3)
+vv = fs.axis
 
 puts "# Interpolation points"
 
@@ -33,7 +33,7 @@ s = bs.plot(vv, 4) do |a, b|
 end
 #STDERR.puts s
 
-st = fft_real_transform(fs, -1)
+st = fs.inverse
 printf "["
 st.each {|x| printf "% .1f ", x}
 puts "]"

data/lib/pspline.rb

@@ -2,5 +2,62 @@ require "pspline/version"
 require "pspline.so"
 
 module PSPLINE
-  # Your code goes here...
+
+class Rfft
+
+	def initialize(data)
+		@R = fft_real_transform(data, 1)
+	end
+	def [] (t)
+		fft_real_get(@R, t)
+	end
+	def spline(j)
+		@V = []
+		fft_real_bspline(@R, j) {|t, z| @V.push t}
+	end
+	def axis
+		@V.dup
+	end
+	def inverse
+		fft_real_transform(@R, -1)
+	end
+	def backword
+		fft_real_transform(@R, 0)
+	end
+	def real
+		@R.dup
+	end
+
+end
+
+class Cfft
+
+	def initialize(data)
+		@F = fft_complex_transform(data, 1)
+	end
+	def [] (t)
+		fft_complex_get(@F, t)
+	end
+	def spline(j)
+		@V = []
+		fft_complex_bspline(@F, j) {|t, z| @V.push t}
+	end
+	def inverse
+		fft_complex_transform(@F, -1)
+	end
+	def backword
+		fft_real_transform(@F, 0)
+	end
+	def axis
+		@V.dup
+	end
+	def real
+		@F[0].dup
+	end
+	def imag
+		@F[1].dup
+	end
+
+end
+
 end

data/lib/pspline/version.rb

@@ -1,3 +1,3 @@
 module PSPLINE
-  VERSION = "5.0.2"
+  VERSION = "5.0.3"
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: pspline
 version: !ruby/object:Gem::Version
-  version: 5.0.2
+  version: 5.0.3
 platform: ruby
 authors:
 - vsp2old
 autorequire: 
 bindir: exe
 cert_chain: []
-date: 2016-02-27 00:00:00.000000000 Z
+date: 2016-02-28 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler