paddlec 0.0.1

data/ext/paddlec/pulseaudio.h

@@ -0,0 +1,30 @@
+/* Copyright (C) 2019 Théotime Bollengier <theotime.bollengier@gmail.com>
+ *
+ * This file is part of PaddleC
+ * 
+ * PaddleC is free software: you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License as published by
+ * the Free Software Foundation, either version 3 of the License, or
+ * (at your option) any later version.
+ * 
+ * PaddleC is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ * GNU General Public License for more details.
+ * 
+ * You should have received a copy of the GNU General Public License
+ * along with PaddleC. If not, see <https://www.gnu.org/licenses/>.
+ */
+
+#ifndef PADDLEC_PULSEAUDIO_H
+#define PADDLEC_PULSEAUDIO_H
+
+#include <ruby.h>
+#include "libpaddlec.h"
+
+extern VALUE m_PulseAudio;
+extern VALUE c_PAConnection;
+
+void Init_paddlec_pulseaudio();
+
+#endif /* PADDLEC_PULSEAUDIO_H */

data/lib/paddlec.rb

@@ -0,0 +1,26 @@
+## Copyright (C) 2019 Théotime Bollengier <theotime.bollengier@gmail.com>
+##
+## This file is part of PaddleC
+## 
+## PaddleC is free software: you can redistribute it and/or modify
+## it under the terms of the GNU General Public License as published by
+## the Free Software Foundation, either version 3 of the License, or
+## (at your option) any later version.
+## 
+## PaddleC is distributed in the hope that it will be useful,
+## but WITHOUT ANY WARRANTY; without even the implied warranty of
+## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+## GNU General Public License for more details.
+## 
+## You should have received a copy of the GNU General Public License
+## along with PaddleC. If not, see <https://www.gnu.org/licenses/>.
+
+
+module PaddleC
+end
+
+require 'roctave'
+require_relative 'paddlec/version.rb'
+#require_relative 'paddlec/paddlec'
+require_relative '../ext/paddlec/paddlec'
+

data/lib/paddlec/version.rb

@@ -0,0 +1,3 @@
+module PaddleC
+	VERSION = '0.0.1' # @return [String]
+end

data/paddlec.gemspec

@@ -0,0 +1,55 @@
+require File.expand_path('../lib/paddlec/version', __FILE__)
+
+Gem::Specification.new do |s|
+	s.name = "paddlec"
+	s.version = PaddleC::VERSION
+	s.date = Time.now.strftime '%Y-%m-%d'
+	s.summary = "A Ruby C extension for real-time signal processing."
+	s.description = "PaddleC is a Ruby C extension attempting to provide objects and methods to rapidly set up real-time signal processing test benches in Ruby."
+	s.license = 'GPL-3.0+'
+	s.authors = ["Théotime Bollengier"]
+	s.email = 'theotime.bollengier@gmail.com'
+	s.homepage = 'https://gitlab.com/theotime_bollengier/paddlec'
+	s.add_runtime_dependency 'roctave', '~> 0.0', '>= 0.0.1'
+	s.files = [
+		'ext/libpaddlec/arithmetic.c',
+		'ext/libpaddlec/comparison.c',
+		'ext/libpaddlec/complex.c',
+		'ext/libpaddlec/delay.c',
+		'ext/libpaddlec/fir_filter_avx.c',
+		'ext/libpaddlec/fir_filter.c',
+		'ext/libpaddlec/fir_filter_neon.c',
+		'ext/libpaddlec/fir_filter_sse.c',
+		'ext/libpaddlec/libpaddlec.c',
+		'ext/libpaddlec/libpaddlec.h',
+		'ext/libpaddlec/math.c',
+		'ext/libpaddlec/no_fast_math.c',
+		'ext/libpaddlec/rounding.c',
+		'ext/paddlec/extconf.rb',
+		'ext/paddlec/complex_buffer.c',
+		'ext/paddlec/complex_buffer.h',
+		'ext/paddlec/delay.c',
+		'ext/paddlec/delay.h',
+		'ext/paddlec/fir_filter.c',
+		'ext/paddlec/fir_filter.h',
+		'ext/paddlec/float_buffer.c',
+		'ext/paddlec/float_buffer.h',
+		'ext/paddlec/paddlec.c',
+		'ext/paddlec/paddlec.h',
+		'ext/paddlec/pulseaudio.c',
+		'ext/paddlec/pulseaudio.h',
+		'lib/paddlec/version.rb',
+		'lib/paddlec.rb',
+		'LICENSE',
+		'paddlec.gemspec',
+		'README.md',
+		'samples/fmdemod_chunk_buffer.rb',
+		'samples/fmdemod_chunk.rb',
+		'samples/fmdemod.rb',
+		'samples/stereo_chunk.rb'
+	]
+	s.extensions = [
+		'ext/paddlec/extconf.rb'
+	]
+end
+

data/samples/fmdemod.rb

@@ -0,0 +1,121 @@
+#!/usr/bin/env ruby
+
+require 'roctave'
+require '../ext/paddlec/paddlec'
+
+
+show_filters = false
+
+fs = 2.048e6
+
+# b = File.binread(File.expand_path '~/sdrrec/pouet.raw').unpack('C*').each_slice(2).collect{|a| Complex((a.first - 128) / 128.0, (a.last - 128) / 128.0)}.to_complex_buffer
+# puts b.length
+# File.binwrite('pouet_cf32le.raw', b.pack)
+# exit
+
+
+unless File.file? 'coefficients.f32le' then
+	STDOUT.write 'Generating filter coeficients... '; STDOUT.flush
+	now = Time.now
+
+	b_bbdecimator = Roctave.fir1(254, 95e3*2/fs)
+	Roctave.freqz(b_bbdecimator, :magnitude, nb_points: 2048, fs: fs) if show_filters
+
+	b_differentiator = Roctave.fir_differentiator(32)
+	Roctave.freqz(b_differentiator, :magnitude, nb_points: 2048, fs: fs/8) if show_filters
+
+	b_adecimator = Roctave.fir_low_pass(512, 15e3*2/256e3)
+	Roctave.freqz(b_adecimator, :magnitude, nb_points: 2048, fs: fs/8) if show_filters
+
+	tau = 50e-6
+	ts = 8*8/fs
+	b = [1 - Math.exp(-ts/tau)]
+	a = [1, -Math.exp(-ts/tau)]
+	h, w = Roctave.freqz(b, a, nb_points: 1024)
+	h.collect!{|v| v.abs}
+	w.collect!{|v| v/w.last}
+	n = 2048
+	b_deemphasis = Roctave.fir2(n, w, h, :odd_symmetry)
+	Roctave.freqz(b_deemphasis, :magnitude, nb_points: 2048, fs: fs/8/8) if show_filters
+
+	puts (Time.now - now).round(3)
+
+	File.open('coefficients.f32le', 'wb') do |ofile|
+		ofile.write b_bbdecimator.pack('f*')
+		ofile.write b_differentiator.pack('f*')
+		ofile.write b_adecimator.pack('f*')
+		ofile.write b_deemphasis.pack('f*')
+	end
+else
+	ifile = File.open('coefficients.f32le', 'rb')
+	b = ifile.read.unpack('f*')
+	ifile.close
+	b_bbdecimator = b.shift(255)
+	b_differentiator = b.shift(33)
+	b_adecimator = b.shift(513)
+	b_deemphasis = b.shift(2049)
+end
+
+
+STDOUT.write 'Instantiating filters... '; STDOUT.flush
+now = Time.now
+
+bbdecimator = PaddleC::FirDecimator.new b_bbdecimator, 8
+differentiator = PaddleC::FirFilter.new b_differentiator
+adecimator = PaddleC::FirDecimator.new b_adecimator, 8
+deemphasis = PaddleC::FirFilter.new b_deemphasis
+puts (Time.now - now).round(3)
+
+
+start = Time.now
+STDOUT.write 'Reading input... '; STDOUT.flush
+now = Time.now
+#if2MHz_complex = PaddleC::ComplexBuffer.unpack(File.binread(File.expand_path('~/sdrrec/pouet.raw')), :u8)
+if2MHz_complex = PaddleC::ComplexBuffer.unpack(File.binread(File.expand_path('~/sdrrec/95.0MHz_2.048Msps_cu8.raw')), :u8)
+puts (Time.now - now).round(3)
+puts if2MHz_complex.length
+
+STDOUT.write 'Input decimation... '; STDOUT.flush
+now = Time.now
+baseband256kHz_complex = bbdecimator.decimate(if2MHz_complex)
+puts (Time.now - now).round(3)
+
+STDOUT.write 'Differentiation... '; STDOUT.flush
+now = Time.now
+baseband_diff_complex, baseband_delayed_complex = differentiator.filter(baseband256kHz_complex, delayed: true)
+puts (Time.now - now).round(3)
+
+STDOUT.write 'Demodulation... '; STDOUT.flush
+now = Time.now
+demodulated = (baseband_diff_complex.imag*baseband_delayed_complex.real - baseband_diff_complex.real*baseband_delayed_complex.imag) / baseband_delayed_complex.abs2
+puts (Time.now - now).round(3)
+
+STDOUT.write 'Decimation to audio... '; STDOUT.flush
+now = Time.now
+audioemph = adecimator.decimate(demodulated)
+puts (Time.now - now).round(3)
+
+STDOUT.write 'Deemphasis... '; STDOUT.flush
+now = Time.now
+audio = deemphasis.filter(audioemph) * 0.25
+puts (Time.now - now).round(3)
+
+puts "Processing took #{(Time.now - start).round 3} s"
+
+PaddleC::PulseAudio::Simple::Sink.open(sample_spec: 'float32le 1ch 32000Hz') do |sink|
+	sink << audio
+end
+
+=begin
+STDOUT.write 'Writing output to output.f32le... '; STDOUT.flush
+now = Time.now
+File.binwrite('output.f32le', audio.pack)
+puts (Time.now - now).round(3)
+
+puts "Processing took #{(Time.now - start).round 3} s"
+
+system "ffmpeg -v error -ar 32k -ac 1 -f f32le -i output.f32le -y output.wav"
+
+system "mpv output.wav"
+=end
+

data/samples/fmdemod_chunk.rb

@@ -0,0 +1,120 @@
+#!/usr/bin/env ruby
+
+require 'roctave'
+require '../ext/paddlec/paddlec'
+
+def show_bar(f, t = 80.0)
+	p = (f.pos * t / f.size).round
+	STDOUT.write "\r[#{'='*p}>#{'-'*(t - p)}]  "
+	STDOUT.flush
+end
+
+
+show_filters = false
+
+fs = 2.048e6
+
+# b = File.binread(File.expand_path '~/sdrrec/pouet.raw').unpack('C*').each_slice(2).collect{|a| Complex((a.first - 128) / 128.0, (a.last - 128) / 128.0)}.to_complex_buffer
+# puts b.length
+# File.binwrite('pouet_cf32le.raw', b.pack)
+# exit
+
+
+unless File.file? 'coefficients.f32le' then
+	STDOUT.write 'Generating filter coeficients... '; STDOUT.flush
+	now = Time.now
+
+	b_bbdecimator = Roctave.fir1(254, 95e3*2/fs)
+	Roctave.freqz(b_bbdecimator, :magnitude, nb_points: 2048, fs: fs) if show_filters
+
+	b_differentiator = Roctave.fir_differentiator(32)
+	Roctave.freqz(b_differentiator, :magnitude, nb_points: 2048, fs: fs/8) if show_filters
+
+	b_adecimator = Roctave.fir_low_pass(512, 15e3*2/256e3)
+	Roctave.freqz(b_adecimator, :magnitude, nb_points: 2048, fs: fs/8) if show_filters
+
+	tau = 50e-6
+	ts = 8*8/fs
+	b = [1 - Math.exp(-ts/tau)]
+	a = [1, -Math.exp(-ts/tau)]
+	h, w = Roctave.freqz(b, a, nb_points: 1024)
+	h.collect!{|v| v.abs}
+	w.collect!{|v| v/w.last}
+	n = 2048
+	b_deemphasis = Roctave.fir2(n, w, h, :odd_symmetry)
+	Roctave.freqz(b_deemphasis, :magnitude, nb_points: 2048, fs: fs/8/8) if show_filters
+
+	puts (Time.now - now).round(3)
+
+	File.open('coefficients.f32le', 'wb') do |ofile|
+		ofile.write b_bbdecimator.pack('f*')
+		ofile.write b_differentiator.pack('f*')
+		ofile.write b_adecimator.pack('f*')
+		ofile.write b_deemphasis.pack('f*')
+	end
+else
+	ifile = File.open('coefficients.f32le', 'rb')
+	b = ifile.read.unpack('f*')
+	ifile.close
+	b_bbdecimator = b.shift(255)
+	b_differentiator = b.shift(33)
+	b_adecimator = b.shift(513)
+	b_deemphasis = b.shift(2049)
+end
+
+
+
+bbdecimator = PaddleC::FirDecimator.new b_bbdecimator, 8
+differentiator = PaddleC::FirFilter.new b_differentiator
+adecimator = PaddleC::FirDecimator.new b_adecimator, 8
+deemphasis = PaddleC::FirFilter.new b_deemphasis
+
+time_chunk = 100e-3
+len = (32e3 * time_chunk).round
+
+decoding = 0.0
+waiting = 0.0
+fname = 'pouet_cf32le.raw'
+fname = '~/sdrrec/pouet.raw'
+fname = '~/sdrrec/95.0MHz_2.048Msps_cu8.raw'
+PaddleC::PulseAudio::Simple::Sink.open(stream_name: fname, sample_spec: 'float32 1ch 32000Hz') do |sink|
+	File.open(File.expand_path(fname), 'rb') do |f|
+		loop do
+			show_bar(f)
+
+			now = Time.now
+
+			if2MHz_complex = PaddleC::ComplexBuffer.unpack(f.read(8*8*len*2), :u8)
+
+			baseband256kHz_complex = bbdecimator.decimate(if2MHz_complex)
+
+			baseband_diff_complex, baseband_delayed_complex = differentiator.filter(baseband256kHz_complex, delayed: true)
+
+			demodulated = (baseband_diff_complex.imag*baseband_delayed_complex.real - baseband_diff_complex.real*baseband_delayed_complex.imag) / baseband_delayed_complex.abs2
+
+			audioemph = adecimator.decimate(demodulated)
+
+			audio = deemphasis.filter(audioemph) * 0.25
+
+			decoding += Time.now - now
+
+			now = Time.now
+
+			sink << audio
+
+			waiting += Time.now - now
+
+			break if f.eof?
+		end
+		show_bar(f)
+		puts
+	end
+end
+
+
+puts "decoding took #{decoding.round(3)}s"
+puts "waiting took #{waiting.round(3)}s"
+
+
+
+

data/samples/fmdemod_chunk_buffer.rb

@@ -0,0 +1,144 @@
+#!/usr/bin/env ruby
+
+require 'roctave'
+require '../ext/paddlec/paddlec'
+
+def show_bar(f, t = 80.0)
+	p = (f.pos * t / f.size).round
+	STDOUT.write "\r[#{'='*p}>#{'-'*(t - p)}]  "
+	STDOUT.flush
+end
+
+
+show_filters = false
+
+fs = 2.048e6
+
+# b = File.binread(File.expand_path '~/sdrrec/pouet.raw').unpack('C*').each_slice(2).collect{|a| Complex((a.first - 128) / 128.0, (a.last - 128) / 128.0)}.to_complex_buffer
+# puts b.length
+# File.binwrite('pouet_cf32le.raw', b.pack)
+# exit
+
+
+unless File.file? 'coefficients.f32le' then
+	STDOUT.write 'Generating filter coeficients... '; STDOUT.flush
+	now = Time.now
+
+	b_bbdecimator = Roctave.fir1(254, 95e3*2/fs)
+	Roctave.freqz(b_bbdecimator, :magnitude, nb_points: 2048, fs: fs) if show_filters
+
+	b_differentiator = Roctave.fir_differentiator(32)
+	Roctave.freqz(b_differentiator, :magnitude, nb_points: 2048, fs: fs/8) if show_filters
+
+	b_adecimator = Roctave.fir_low_pass(512, 15e3*2/256e3)
+	Roctave.freqz(b_adecimator, :magnitude, nb_points: 2048, fs: fs/8) if show_filters
+
+	tau = 50e-6
+	ts = 8*8/fs
+	b = [1 - Math.exp(-ts/tau)]
+	a = [1, -Math.exp(-ts/tau)]
+	h, w = Roctave.freqz(b, a, nb_points: 1024)
+	h.collect!{|v| v.abs}
+	w.collect!{|v| v/w.last}
+	n = 2048
+	b_deemphasis = Roctave.fir2(n, w, h, :odd_symmetry)
+	Roctave.freqz(b_deemphasis, :magnitude, nb_points: 2048, fs: fs/8/8) if show_filters
+
+	puts (Time.now - now).round(3)
+
+	File.open('coefficients.f32le', 'wb') do |ofile|
+		ofile.write b_bbdecimator.pack('f*')
+		ofile.write b_differentiator.pack('f*')
+		ofile.write b_adecimator.pack('f*')
+		ofile.write b_deemphasis.pack('f*')
+	end
+else
+	ifile = File.open('coefficients.f32le', 'rb')
+	b = ifile.read.unpack('f*')
+	ifile.close
+	b_bbdecimator = b.shift(255)
+	b_differentiator = b.shift(33)
+	b_adecimator = b.shift(513)
+	b_deemphasis = b.shift(2049)
+end
+
+
+
+bbdecimator = PaddleC::FirDecimator.new b_bbdecimator, 8
+differentiator = PaddleC::FirFilter.new b_differentiator
+adecimator = PaddleC::FirDecimator.new b_adecimator, 8
+deemphasis = PaddleC::FirFilter.new b_deemphasis
+
+time_chunk = 100e-3
+len = (32e3 * time_chunk).round
+audio = PaddleC::FloatBuffer.new
+audioemph = PaddleC::FloatBuffer.new
+demodulated = PaddleC::FloatBuffer.new
+baseband_diff_i = PaddleC::FloatBuffer.new
+baseband_diff_r = PaddleC::FloatBuffer.new
+baseband_del_i = PaddleC::FloatBuffer.new
+baseband_del_r = PaddleC::FloatBuffer.new
+baseband_del_abs2 = PaddleC::FloatBuffer.new
+baseband_diff_complex = PaddleC::ComplexBuffer.new
+baseband_delayed_complex = PaddleC::ComplexBuffer.new
+baseband256kHz_complex = PaddleC::ComplexBuffer.new
+if2MHz_complex = PaddleC::ComplexBuffer.new
+strbuf = ''
+
+decoding = 0.0
+waiting = 0.0
+fname = 'pouet_cf32le.raw'
+fname = '~/sdrrec/95.0MHz_2.048Msps_cu8.raw'
+PaddleC::PulseAudio::Simple::Sink.open(stream_name: fname, sample_spec: 'float32le 1ch 32000Hz') do |sink|
+	File.open(File.expand_path(fname), 'rb') do |f|
+		loop do
+			show_bar(f)
+
+			now = Time.now
+
+			strbuf = f.read(8*8*len*2, strbuf)
+			if2MHz_complex.unpack(strbuf, :u8)
+
+			bbdecimator.decimate(if2MHz_complex, buffer: baseband256kHz_complex)
+			if2MHz_complex.resize 0
+
+			differentiator.filter(baseband256kHz_complex, buffer: baseband_diff_complex, delayed: baseband_delayed_complex)
+
+			baseband_diff_complex.imag(baseband_diff_i)
+			baseband_diff_complex.real(baseband_diff_r)
+			baseband_delayed_complex.imag(baseband_del_i)
+			baseband_delayed_complex.real(baseband_del_r)
+			baseband_delayed_complex.abs2(baseband_del_abs2)
+
+			baseband_diff_i.mult!(baseband_del_r)
+			baseband_diff_r.mult!(baseband_del_i)
+			baseband_diff_i.sub!(baseband_diff_r)
+
+			baseband_diff_i.div(baseband_del_abs2, demodulated)
+
+			adecimator.decimate(demodulated, buffer: audioemph)
+
+			deemphasis.filter(audioemph, buffer: audio).mult!(0.25)
+
+			decoding += Time.now - now
+
+			now = Time.now
+
+			sink << audio
+
+			waiting += Time.now - now
+
+			break if f.eof?
+		end
+		show_bar(f)
+		puts
+	end
+end
+
+
+puts "decoding took #{decoding.round(3)}s"
+puts "waiting took #{waiting.round(3)}s"
+
+
+
+