torchaudio 0.1.0 → 0.2.1

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: b527976494325cc12e81342c25d318204d2d7c75bfba7036be4296769cdb30a0
-  data.tar.gz: 2cfde7bd1b0e7a1628818d5bd74657cfbfba6dfa83ef42897f3ad0f98e77f739
+  metadata.gz: 9ed4c14921f1eee18f5e08ddabfae51e09a9b5a7ef408f1dd67fdf7bfe9622fe
+  data.tar.gz: 1e37d5b9abed9cab7bf56a8c30a769bc8ff8f8a3e15e78bbb772847c444571b2
 SHA512:
-  metadata.gz: 8e6f34b014340b5ace3193ab589dae75ed0869ab7606402bd4b09de6042299e6f3a118d439dd381491f489ce9552bca4376a7d5b4693dddc3d1c5f5b26540900
-  data.tar.gz: d651c46f5185ceb70ae3d9c90154c77afe29a5c35854d1a9d98913096b7ab9ba39a745242dd268548ca87f9e109b56c96dee9dc5539cf066f9ad0f773eddbdcd
+  metadata.gz: 9ca5436d7e4309dd9659fdce7ee893b122e9da96e9f7b15bf00de5dea32c635e2828a99939f04ef5bf0d9494ab89957829a65002dc3e855fa8a66f54abbbd181
+  data.tar.gz: d62b2a137c19d3b24facb11eda5c1b81be5841120b505877b8617bee2b9f183dbe4b4d42a95af27447346a3d48476d7faec48b57cd89b88c0ddc9709f1b5d51b

data/CHANGELOG.md

@@ -1,3 +1,24 @@
+## 0.2.1 (2021-07-16)
+
+- Added `create_dct` method
+- Added `ComputeDeltas`, `Fade`, `MFCC`, and `Vol` transforms
+
+## 0.2.0 (2021-05-23)
+
+- Updated to Rice 4
+- Dropped support for Ruby < 2.6
+
+## 0.1.2 (2021-02-06)
+
+- Added `amplitude_to_DB` and `DB_to_amplitude` methods
+- Added `AmplitudeToDB` transform
+- Fixed `save` options
+
+## 0.1.1 (2020-08-26)
+
+- Added `save` method
+- Added transforms
+
 ## 0.1.0 (2020-08-24)
 
 - First release

data/LICENSE.txt

@@ -1,7 +1,7 @@
 BSD 2-Clause License
 
 Copyright (c) 2017 Facebook Inc. (Soumith Chintala),
-Copyright (c) 2020 Andrew Kane,
+Copyright (c) 2020-2021 Andrew Kane,
 All rights reserved.
 
 Redistribution and use in source and binary forms, with or without

data/README.md

@@ -2,6 +2,8 @@
 
 :fire: An audio library for Torch.rb
 
+[![Build Status](https://github.com/ankane/torchaudio/workflows/build/badge.svg?branch=master)](https://github.com/ankane/torchaudio/actions)
+
 ## Installation
 
 First, [install SoX](#sox-installation). For Homebrew, use:
@@ -20,6 +22,67 @@ gem 'torchaudio'
 
 This library follows the [Python API](https://pytorch.org/audio/). Many methods and options are missing at the moment. PRs welcome!
 
+## Tutorial
+
+- [PyTorch tutorial](https://pytorch.org/tutorials/beginner/audio_preprocessing_tutorial.html)
+- [Ruby code](examples/tutorial.rb)
+
+Download the [audio file](https://github.com/pytorch/tutorials/raw/master/_static/img/steam-train-whistle-daniel_simon-converted-from-mp3.wav) and install the [matplotlib](https://github.com/mrkn/matplotlib.rb) gem first.
+
+## Basics
+
+Load a file
+
+```ruby
+waveform, sample_rate = TorchAudio.load("file.wav")
+```
+
+Save a file
+
+```ruby
+TorchAudio.save("new.wave", waveform, sample_rate)
+```
+
+## Transforms
+
+```ruby
+TorchAudio::Transforms::Spectrogram.new.call(waveform)
+```
+
+Supported transforms are:
+
+- AmplitudeToDB
+- ComputeDeltas
+- Fade
+- MelScale
+- MelSpectrogram
+- MFCC
+- MuLawDecoding
+- MuLawEncoding
+- Spectrogram
+- Vol
+
+## Functional
+
+```ruby
+TorchAudio::Functional.lowpass_biquad(waveform, sample_rate, cutoff_freq)
+```
+
+Supported functions are:
+
+- amplitude_to_DB
+- compute_deltas
+- create_dct
+- create_fb_matrix
+- DB_to_amplitude
+- dither
+- gain
+- highpass_biquad
+- lowpass_biquad
+- mu_law_decoding
+- mu_law_encoding
+- spectrogram
+
 ## Datasets
 
 Load a dataset

data/ext/torchaudio/ext.cpp

@@ -1,33 +1,33 @@
 #include <torchaudio/csrc/sox.h>
 
-#include <rice/Module.hpp>
-
-using namespace Rice;
-
-template<>
-inline
-sox_signalinfo_t* from_ruby<sox_signalinfo_t*>(Object x)
-{
-  if (x.is_nil()) {
-    return nullptr;
-  }
-  throw std::runtime_error("Unsupported signalinfo");
-}
-
-template<>
-inline
-sox_encodinginfo_t* from_ruby<sox_encodinginfo_t*>(Object x)
-{
-  if (x.is_nil()) {
-    return nullptr;
-  }
-  throw std::runtime_error("Unsupported encodinginfo");
-}
+#include <rice/rice.hpp>
+#include <rice/stl.hpp>
 
 extern "C"
 void Init_ext()
 {
-  Module rb_mTorchAudio = define_module("TorchAudio");
-  Module rb_mNN = define_module_under(rb_mTorchAudio, "Ext")
-    .define_singleton_method("read_audio_file", &torch::audio::read_audio_file);
+  auto rb_mTorchAudio = Rice::define_module("TorchAudio");
+
+  auto rb_mExt = Rice::define_module_under(rb_mTorchAudio, "Ext")
+    .define_singleton_function(
+      "read_audio_file",
+      [](const std::string& file_name, at::Tensor output, bool ch_first, int64_t nframes, int64_t offset, sox_signalinfo_t* si, sox_encodinginfo_t* ei, const char* ft) {
+        return torch::audio::read_audio_file(file_name, output, ch_first, nframes, offset, si, ei, ft);
+      })
+    .define_singleton_function(
+      "write_audio_file",
+      [](const std::string& file_name, const at::Tensor& tensor, sox_signalinfo_t* si, sox_encodinginfo_t* ei, const char* file_type) {
+        return torch::audio::write_audio_file(file_name, tensor, si, ei, file_type);
+      });
+
+  auto rb_cSignalInfo = Rice::define_class_under<sox_signalinfo_t>(rb_mExt, "SignalInfo")
+    .define_constructor(Rice::Constructor<sox_signalinfo_t>())
+    .define_method("rate", [](sox_signalinfo_t& self) { return self.rate; })
+    .define_method("channels", [](sox_signalinfo_t& self) { return self.channels; })
+    .define_method("precision", [](sox_signalinfo_t& self) { return self.precision; })
+    .define_method("length", [](sox_signalinfo_t& self) { return self.length; })
+    .define_method("rate=", [](sox_signalinfo_t& self, sox_rate_t rate) { self.rate = rate; })
+    .define_method("channels=", [](sox_signalinfo_t& self, unsigned channels) { self.channels = channels; })
+    .define_method("precision=", [](sox_signalinfo_t& self, unsigned precision) { self.precision = precision; })
+    .define_method("length=", [](sox_signalinfo_t& self, sox_uint64_t length) { self.length = length; });
 }

data/ext/torchaudio/extconf.rb

@@ -1,8 +1,6 @@
 require "mkmf-rice"
 
-abort "Missing stdc++" unless have_library("stdc++")
-
-$CXXFLAGS += " -std=c++14"
+$CXXFLAGS += " -std=c++17 $(optflags)"
 
 abort "SoX not found" unless have_library("sox")
 
@@ -24,7 +22,6 @@ apple_clang = RbConfig::CONFIG["CC_VERSION_MESSAGE"] =~ /apple clang/i
 
 # check omp first
 if have_library("omp") || have_library("gomp")
-  $CXXFLAGS += " -DAT_PARALLEL_OPENMP=1"
   $CXXFLAGS += " -Xclang" if apple_clang
   $CXXFLAGS += " -fopenmp"
 end

data/lib/torchaudio.rb

@@ -14,6 +14,17 @@ require "set"
 # modules
 require "torchaudio/datasets/utils"
 require "torchaudio/datasets/yesno"
+require "torchaudio/functional"
+require "torchaudio/transforms/compute_deltas"
+require "torchaudio/transforms/fade"
+require "torchaudio/transforms/mel_scale"
+require "torchaudio/transforms/mel_spectrogram"
+require "torchaudio/transforms/mu_law_encoding"
+require "torchaudio/transforms/mu_law_decoding"
+require "torchaudio/transforms/spectrogram"
+require "torchaudio/transforms/amplitude_to_db"
+require "torchaudio/transforms/mfcc"
+require "torchaudio/transforms/vol"
 require "torchaudio/version"
 
 module TorchAudio
@@ -73,6 +84,70 @@ module TorchAudio
       load(filepath, **kwargs)
     end
 
+    def save(filepath, src, sample_rate, precision: 16, channels_first: true)
+      si = Ext::SignalInfo.new
+      ch_idx = channels_first ? 0 : 1
+      si.rate = sample_rate
+      si.channels = src.dim == 1 ? 1 : src.size(ch_idx)
+      si.length = src.numel
+      si.precision = precision
+      save_encinfo(filepath, src, channels_first: channels_first, signalinfo: si)
+    end
+
+    def save_encinfo(filepath, src, channels_first: true, signalinfo: nil, encodinginfo: nil, filetype: nil)
+      ch_idx, len_idx = channels_first ? [0, 1] : [1, 0]
+
+      # check if save directory exists
+      abs_dirpath = File.dirname(File.expand_path(filepath))
+      unless Dir.exist?(abs_dirpath)
+        raise "Directory does not exist: #{abs_dirpath}"
+      end
+      # check that src is a CPU tensor
+      check_input(src)
+      # Check/Fix shape of source data
+      if src.dim == 1
+        # 1d tensors as assumed to be mono signals
+        src.unsqueeze!(ch_idx)
+      elsif src.dim > 2 || src.size(ch_idx) > 16
+        # assumes num_channels < 16
+        raise ArgumentError, "Expected format where C < 16, but found #{src.size}"
+      end
+      # sox stores the sample rate as a float, though practically sample rates are almost always integers
+      # convert integers to floats
+      if signalinfo
+        if signalinfo.rate && !signalinfo.rate.is_a?(Float)
+          if signalinfo.rate.to_f == signalinfo.rate
+            signalinfo.rate = signalinfo.rate.to_f
+          else
+            raise ArgumentError, "Sample rate should be a float or int"
+          end
+        end
+        # check if the bit precision (i.e. bits per sample) is an integer
+        if signalinfo.precision && ! signalinfo.precision.is_a?(Integer)
+          if signalinfo.precision.to_i == signalinfo.precision
+            signalinfo.precision = signalinfo.precision.to_i
+          else
+            raise ArgumentError, "Bit precision should be an integer"
+          end
+        end
+      end
+      # programs such as librosa normalize the signal, unnormalize if detected
+      if src.min >= -1.0 && src.max <= 1.0
+        src = src * (1 << 31)
+        src = src.long
+      end
+      # set filetype and allow for files with no extensions
+      extension = File.extname(filepath)
+      filetype = extension.length > 0 ? extension[1..-1] : filetype
+      # transpose from C x L -> L x C
+      if channels_first
+        src = src.transpose(1, 0)
+      end
+      # save data to file
+      src = src.contiguous
+      Ext.write_audio_file(filepath, src, signalinfo, encodinginfo, filetype)
+    end
+
     private
 
     def check_input(src)

data/lib/torchaudio/functional.rb

@@ -0,0 +1,324 @@
+module TorchAudio
+  module Functional
+    class << self
+      def spectrogram(waveform, pad, window, n_fft, hop_length, win_length, power, normalized)
+        if pad > 0
+          # TODO add "with torch.no_grad():" back when JIT supports it
+          waveform = Torch::NN::Functional.pad(waveform, [pad, pad], "constant")
+        end
+
+        # pack batch
+        shape = waveform.size
+        waveform = waveform.reshape(-1, shape[-1])
+
+        # default values are consistent with librosa.core.spectrum._spectrogram
+        spec_f =
+          Torch.stft(
+            waveform,
+            n_fft,
+            hop_length: hop_length,
+            win_length: win_length,
+            window: window,
+            center: true,
+            pad_mode: "reflect",
+            normalized: false,
+            onesided: true
+          )
+
+        # unpack batch
+        spec_f = spec_f.reshape(shape[0..-2] + spec_f.shape[-3..-1])
+
+        if normalized
+          spec_f.div!(window.pow(2.0).sum.sqrt)
+        end
+        if power
+          spec_f = complex_norm(spec_f, power: power)
+        end
+
+        spec_f
+      end
+
+      def mu_law_encoding(x, quantization_channels)
+        mu = quantization_channels - 1.0
+        if !x.floating_point?
+          x = x.to(dtype: :float)
+        end
+        mu = Torch.tensor(mu, dtype: x.dtype)
+        x_mu = Torch.sign(x) * Torch.log1p(mu * Torch.abs(x)) / Torch.log1p(mu)
+        x_mu = ((x_mu + 1) / 2 * mu + 0.5).to(dtype: :int64)
+        x_mu
+      end
+
+      def mu_law_decoding(x_mu, quantization_channels)
+        mu = quantization_channels - 1.0
+        if !x_mu.floating_point?
+          x_mu = x_mu.to(dtype: :float)
+        end
+        mu = Torch.tensor(mu, dtype: x_mu.dtype)
+        x = ((x_mu) / mu) * 2 - 1.0
+        x = Torch.sign(x) * (Torch.exp(Torch.abs(x) * Torch.log1p(mu)) - 1.0) / mu
+        x
+      end
+
+      def complex_norm(complex_tensor, power: 1.0)
+        complex_tensor.pow(2.0).sum(-1).pow(0.5 * power)
+      end
+
+      def create_fb_matrix(n_freqs, f_min, f_max, n_mels, sample_rate, norm: nil)
+        if norm && norm != "slaney"
+          raise ArgumentError, "norm must be one of None or 'slaney'"
+        end
+
+        # freq bins
+        # Equivalent filterbank construction by Librosa
+        all_freqs = Torch.linspace(0, sample_rate.div(2), n_freqs)
+
+        # calculate mel freq bins
+        # hertz to mel(f) is 2595. * math.log10(1. + (f / 700.))
+        m_min = 2595.0 * Math.log10(1.0 + (f_min / 700.0))
+        m_max = 2595.0 * Math.log10(1.0 + (f_max / 700.0))
+        m_pts = Torch.linspace(m_min, m_max, n_mels + 2)
+        # mel to hertz(mel) is 700. * (10**(mel / 2595.) - 1.)
+        f_pts = (Torch.pow(10, m_pts / 2595.0) - 1.0) * 700.0
+        # calculate the difference between each mel point and each stft freq point in hertz
+        f_diff = f_pts[1..-1] - f_pts[0...-1]  # (n_mels + 1)
+        slopes = f_pts.unsqueeze(0) - all_freqs.unsqueeze(1)  # (n_freqs, n_mels + 2)
+        # create overlapping triangles
+        zero = Torch.zeros(1)
+        down_slopes = (slopes[0..-1, 0...-2] * -1.0) / f_diff[0...-1]  # (n_freqs, n_mels)
+        up_slopes = slopes[0..-1, 2..-1] / f_diff[1..-1]  # (n_freqs, n_mels)
+        fb = Torch.max(zero, Torch.min(down_slopes, up_slopes))
+
+        if norm && norm == "slaney"
+          # Slaney-style mel is scaled to be approx constant energy per channel
+          enorm = 2.0 / (f_pts[2...(n_mels + 2)] - f_pts[:n_mels])
+          fb *= enorm.unsqueeze(0)
+        end
+
+        fb
+      end
+
+      def compute_deltas(specgram, win_length: 5, mode: "replicate")
+        device = specgram.device
+        dtype = specgram.dtype
+
+        # pack batch
+        shape = specgram.size
+        specgram = specgram.reshape(1, -1, shape[-1])
+
+        raise ArgumentError, "win_length must be >= 3" unless win_length >= 3
+
+        n = (win_length - 1).div(2)
+
+        # twice sum of integer squared
+        denom = n * (n + 1) * (2 * n + 1) / 3
+
+        specgram = Torch::NN::Functional.pad(specgram, [n, n], mode: mode)
+
+        kernel = Torch.arange(-n, n + 1, 1, device: device, dtype: dtype).repeat([specgram.shape[1], 1, 1])
+
+        output = Torch::NN::Functional.conv1d(specgram, kernel, groups: specgram.shape[1]) / denom
+
+        # unpack batch
+        output = output.reshape(shape)
+      end
+
+      def gain(waveform, gain_db: 1.0)
+        return waveform if gain_db == 0
+
+        ratio = 10 ** (gain_db / 20)
+
+        waveform * ratio
+      end
+
+      def dither(waveform, density_function: "TPDF", noise_shaping: false)
+        dithered = _apply_probability_distribution(waveform, density_function: density_function)
+
+        if noise_shaping
+          raise "Not implemented yet"
+          # _add_noise_shaping(dithered, waveform)
+        else
+          dithered
+        end
+      end
+
+      def biquad(waveform, b0, b1, b2, a0, a1, a2)
+        device = waveform.device
+        dtype = waveform.dtype
+
+        output_waveform = lfilter(
+            waveform,
+            Torch.tensor([a0, a1, a2], dtype: dtype, device: device),
+            Torch.tensor([b0, b1, b2], dtype: dtype, device: device)
+        )
+        output_waveform
+      end
+
+      def highpass_biquad(waveform, sample_rate, cutoff_freq, q: 0.707)
+        w0 = 2 * Math::PI * cutoff_freq / sample_rate
+        alpha = Math.sin(w0) / 2.0 / q
+
+        b0 = (1 + Math.cos(w0)) / 2
+        b1 = -1 - Math.cos(w0)
+        b2 = b0
+        a0 = 1 + alpha
+        a1 = -2 * Math.cos(w0)
+        a2 = 1 - alpha
+        biquad(waveform, b0, b1, b2, a0, a1, a2)
+      end
+
+      def lowpass_biquad(waveform, sample_rate, cutoff_freq, q: 0.707)
+        w0 = 2 * Math::PI * cutoff_freq / sample_rate
+        alpha = Math.sin(w0) / 2 / q
+
+        b0 = (1 - Math.cos(w0)) / 2
+        b1 = 1 - Math.cos(w0)
+        b2 = b0
+        a0 = 1 + alpha
+        a1 = -2 * Math.cos(w0)
+        a2 = 1 - alpha
+        biquad(waveform, b0, b1, b2, a0, a1, a2)
+      end
+
+      def lfilter(waveform, a_coeffs, b_coeffs, clamp: true)
+        # pack batch
+        shape = waveform.size
+        waveform = waveform.reshape(-1, shape[-1])
+
+        raise ArgumentError unless (a_coeffs.size(0) == b_coeffs.size(0))
+        raise ArgumentError unless (waveform.size.length == 2)
+        raise ArgumentError unless (waveform.device == a_coeffs.device)
+        raise ArgumentError unless (b_coeffs.device == a_coeffs.device)
+
+        device = waveform.device
+        dtype = waveform.dtype
+        n_channel, n_sample = waveform.size
+        n_order = a_coeffs.size(0)
+        n_sample_padded = n_sample + n_order - 1
+        raise ArgumentError unless (n_order > 0)
+
+        # Pad the input and create output
+        padded_waveform = Torch.zeros(n_channel, n_sample_padded, dtype: dtype, device: device)
+        padded_waveform[0..-1, (n_order - 1)..-1] = waveform
+        padded_output_waveform = Torch.zeros(n_channel, n_sample_padded, dtype: dtype, device: device)
+
+        # Set up the coefficients matrix
+        # Flip coefficients' order
+        a_coeffs_flipped = a_coeffs.flip([0])
+        b_coeffs_flipped = b_coeffs.flip([0])
+
+        # calculate windowed_input_signal in parallel
+        # create indices of original with shape (n_channel, n_order, n_sample)
+        window_idxs = Torch.arange(n_sample, device: device).unsqueeze(0) + Torch.arange(n_order, device: device).unsqueeze(1)
+        window_idxs = window_idxs.repeat([n_channel, 1, 1])
+        window_idxs += (Torch.arange(n_channel, device: device).unsqueeze(-1).unsqueeze(-1) * n_sample_padded)
+        window_idxs = window_idxs.long
+        # (n_order, ) matmul (n_channel, n_order, n_sample) -> (n_channel, n_sample)
+        input_signal_windows = Torch.matmul(b_coeffs_flipped, Torch.take(padded_waveform, window_idxs))
+
+        input_signal_windows.div!(a_coeffs[0])
+        a_coeffs_flipped.div!(a_coeffs[0])
+        input_signal_windows.t.each_with_index do |o0, i_sample|
+          windowed_output_signal = padded_output_waveform[0..-1, i_sample...(i_sample + n_order)]
+          o0.addmv!(windowed_output_signal, a_coeffs_flipped, alpha: -1)
+          padded_output_waveform[0..-1, i_sample + n_order - 1] = o0
+        end
+
+        output = padded_output_waveform[0..-1, (n_order - 1)..-1]
+
+        if clamp
+          output = Torch.clamp(output, -1.0, 1.0)
+        end
+
+        # unpack batch
+        output = output.reshape(shape[0...-1] + output.shape[-1..-1])
+
+        output
+      end
+
+      def amplitude_to_DB(amp, multiplier, amin, db_multiplier, top_db: nil)
+        db = Torch.log10(Torch.clamp(amp, min: amin)) * multiplier
+        db -= multiplier * db_multiplier
+
+        db = db.clamp(min: db.max.item - top_db) if top_db
+
+        db
+      end
+
+      def DB_to_amplitude(db, ref, power)
+        Torch.pow(Torch.pow(10.0, db * 0.1), power) * ref
+      end
+
+      def create_dct(n_mfcc, n_mels, norm: nil)
+        n = Torch.arange(n_mels.to_f)
+        k = Torch.arange(n_mfcc.to_f).unsqueeze!(1)
+        dct = Torch.cos((n + 0.5) * k * Math::PI / n_mels.to_f)
+
+        if norm.nil?
+          dct *= 2.0
+        else
+          raise ArgumentError, "Invalid DCT norm value" unless norm == :ortho
+
+          dct[0] *= 1.0 / Math.sqrt(2.0)
+          dct *= Math.sqrt(2.0 / n_mels)
+        end
+
+        dct.t
+      end
+
+      private
+
+      def _apply_probability_distribution(waveform, density_function: "TPDF")
+        # pack batch
+        shape = waveform.size
+        waveform = waveform.reshape(-1, shape[-1])
+
+        channel_size = waveform.size[0] - 1
+        time_size = waveform.size[-1] - 1
+
+        random_channel = channel_size > 0 ? Torch.randint(channel_size, [1]).item.to_i : 0
+        random_time = time_size > 0 ? Torch.randint(time_size, [1]).item.to_i : 0
+
+        number_of_bits = 16
+        up_scaling = 2 ** (number_of_bits - 1) - 2
+        signal_scaled = waveform * up_scaling
+        down_scaling = 2 ** (number_of_bits - 1)
+
+        signal_scaled_dis = waveform
+        if density_function == "RPDF"
+          rpdf = waveform[random_channel][random_time] - 0.5
+
+          signal_scaled_dis = signal_scaled + rpdf
+        elsif density_function == "GPDF"
+          # TODO Replace by distribution code once
+          # https://github.com/pytorch/pytorch/issues/29843 is resolved
+          # gaussian = torch.distributions.normal.Normal(torch.mean(waveform, -1), 1).sample()
+
+          num_rand_variables = 6
+
+          gaussian = waveform[random_channel][random_time]
+          (num_rand_variables * [time_size]).each do |ws|
+            rand_chan = Torch.randint(channel_size, [1]).item.to_i
+            gaussian += waveform[rand_chan][Torch.randint(ws, [1]).item.to_i]
+          end
+
+          signal_scaled_dis = signal_scaled + gaussian
+        else
+          # dtype needed for https://github.com/pytorch/pytorch/issues/32358
+          # TODO add support for dtype and device to bartlett_window
+          tpdf = Torch.bartlett_window(time_size + 1).to(signal_scaled.device, dtype: signal_scaled.dtype)
+          tpdf = tpdf.repeat([channel_size + 1, 1])
+          signal_scaled_dis = signal_scaled + tpdf
+        end
+
+        quantised_signal_scaled = Torch.round(signal_scaled_dis)
+        quantised_signal = quantised_signal_scaled / down_scaling
+
+        # unpack batch
+        quantised_signal.reshape(shape[0...-1] + quantised_signal.shape[-1..-1])
+      end
+    end
+  end
+
+  F = Functional
+end

data/lib/torchaudio/transforms/amplitude_to_db.rb

@@ -0,0 +1,27 @@
+module TorchAudio
+  module Transforms
+    class AmplitudeToDB < Torch::NN::Module
+      def initialize(stype: :power, top_db: nil)
+        super()
+
+        @stype = stype
+
+        raise ArgumentError, 'top_db must be a positive numerical' if top_db && top_db.negative?
+
+        @top_db = top_db
+        @multiplier = stype == :power ? 10.0 : 20.0
+        @amin = 1e-10
+        @ref_value = 1.0
+        @db_multiplier = Math.log10([@amin, @ref_value].max)
+      end
+
+      def forward(amplitude_spectrogram)
+        F.amplitude_to_DB(
+          amplitude_spectrogram,
+          @multiplier, @amin, @db_multiplier,
+          top_db: @top_db
+        )
+      end
+    end
+  end
+end

data/lib/torchaudio/transforms/compute_deltas.rb

@@ -0,0 +1,15 @@
+module TorchAudio
+  module Transforms
+    class ComputeDeltas < Torch::NN::Module
+      def initialize(win_length: 5, mode: "replicate")
+        super()
+        @win_length = win_length
+        @mode = mode
+      end
+
+      def forward(specgram)
+        F.compute_deltas(specgram, win_length: @win_length, mode: @mode)
+      end
+    end
+  end
+end

data/lib/torchaudio/transforms/fade.rb

@@ -0,0 +1,74 @@
+module TorchAudio
+  module Transforms
+    class Fade < Torch::NN::Module
+      def initialize(fade_in_len: 0, fade_out_len: 0, fade_shape: "linear")
+        super()
+        @fade_in_len = fade_in_len
+        @fade_out_len = fade_out_len
+        @fade_shape = fade_shape
+      end
+
+      def forward(waveform)
+        waveform_length = waveform.size[-1]
+        device = waveform.device
+        fade_in(waveform_length).to(device) * fade_out(waveform_length).to(device) * waveform
+      end
+
+      private
+
+      def fade_in(waveform_length)
+        fade = Torch.linspace(0, 1, @fade_in_len)
+        ones = Torch.ones(waveform_length - @fade_in_len)
+
+        if @fade_shape == "linear"
+          fade = fade
+        end
+
+        if @fade_shape == "exponential"
+          fade = Torch.pow(2, (fade - 1)) * fade
+        end
+
+        if @fade_shape == "logarithmic"
+          fade = Torch.log10(0.1 + fade) + 1
+        end
+
+        if @fade_shape == "quarter_sine"
+          fade = Torch.sin(fade * Math::PI / 2)
+        end
+
+        if @fade_shape == "half_sine"
+          fade = Torch.sin(fade * Math::PI - Math::PI / 2) / 2 + 0.5
+        end
+
+        Torch.cat([fade, ones]).clamp!(0, 1)
+      end
+
+      def fade_out(waveform_length)
+        fade = Torch.linspace(0, 1, @fade_out_len)
+        ones = Torch.ones(waveform_length - @fade_out_len)
+
+        if @fade_shape == "linear"
+          fade = - fade + 1
+        end
+
+        if @fade_shape == "exponential"
+          fade = Torch.pow(2, - fade) * (1 - fade)
+        end
+
+        if @fade_shape == "logarithmic"
+          fade = Torch.log10(1.1 - fade) + 1
+        end
+
+        if @fade_shape == "quarter_sine"
+          fade = Torch.sin(fade * Math::PI / 2 + Math::PI / 2)
+        end
+
+        if @fade_shape == "half_sine"
+          fade = Torch.sin(fade * Math::PI + Math::PI / 2) / 2 + 0.5
+        end
+
+        Torch.cat([ones, fade]).clamp!(0, 1)
+      end
+    end
+  end
+end

data/lib/torchaudio/transforms/mel_scale.rb

@@ -0,0 +1,39 @@
+module TorchAudio
+  module Transforms
+    class MelScale < Torch::NN::Module
+      def initialize(n_mels: 128, sample_rate: 16000, f_min: 0.0, f_max: nil, n_stft: nil)
+        super()
+        @n_mels = n_mels
+        @sample_rate = sample_rate
+        @f_max = f_max || sample_rate.div(2).to_f
+        @f_min = f_min
+
+        raise ArgumentError, "Require f_min: %f < f_max: %f" % [f_min, @f_max] unless f_min <= @f_max
+
+        fb = n_stft.nil? ? Torch.empty(0) : F.create_fb_matrix(n_stft, @f_min, @f_max, @n_mels, @sample_rate)
+        register_buffer("fb", fb)
+      end
+
+      def forward(specgram)
+        shape = specgram.size
+        specgram = specgram.reshape(-1, shape[-2], shape[-1])
+
+        if @fb.numel == 0
+          tmp_fb = F.create_fb_matrix(specgram.size(1), @f_min, @f_max, @n_mels, @sample_rate)
+          # Attributes cannot be reassigned outside __init__ so workaround
+          @fb.resize!(tmp_fb.size)
+          @fb.copy!(tmp_fb)
+        end
+
+        # (channel, frequency, time).transpose(...) dot (frequency, n_mels)
+        # -> (channel, time, n_mels).transpose(...)
+        mel_specgram = Torch.matmul(specgram.transpose(1, 2), @fb).transpose(1, 2)
+
+        # unpack batch
+        mel_specgram = mel_specgram.reshape(shape[0...-2] + mel_specgram.shape[-2..-1])
+
+        mel_specgram
+      end
+    end
+  end
+end

data/lib/torchaudio/transforms/mel_spectrogram.rb

@@ -0,0 +1,37 @@
+module TorchAudio
+  module Transforms
+    class MelSpectrogram < Torch::NN::Module
+      attr_reader :n_mels
+
+      def initialize(
+        sample_rate: 16000, n_fft: 400, win_length: nil, hop_length: nil, f_min: 0.0,
+        f_max: nil, pad: 0, n_mels: 128, window_fn: Torch.method(:hann_window),
+        power: 2.0, normalized: false, wkwargs: nil
+      )
+
+        super()
+        @sample_rate = sample_rate
+        @n_fft = n_fft
+        @win_length = win_length || n_fft
+        @hop_length = hop_length || @win_length.div(2)
+        @pad = pad
+        @power = power
+        @normalized = normalized
+        @n_mels = n_mels  # number of mel frequency bins
+        @f_max = f_max
+        @f_min = f_min
+        @spectrogram =
+          Spectrogram.new(
+            n_fft: @n_fft, win_length: @win_length, hop_length: @hop_length, pad: @pad,
+            window_fn: window_fn, power: @power, normalized: @normalized, wkwargs: wkwargs
+          )
+        @mel_scale = MelScale.new(n_mels: @n_mels, sample_rate: @sample_rate, f_min: @f_min, f_max: @f_max, n_stft: @n_fft.div(2) + 1)
+      end
+
+      def forward(waveform)
+        specgram = @spectrogram.call(waveform)
+        @mel_scale.call(specgram)
+      end
+    end
+  end
+end

data/lib/torchaudio/transforms/mfcc.rb

@@ -0,0 +1,43 @@
+module TorchAudio
+  module Transforms
+    class MFCC < Torch::NN::Module
+
+      SUPPORTED_DCT_TYPES = [2]
+
+      def initialize(sample_rate: 16000, n_mfcc: 40, dct_type: 2, norm: :ortho, log_mels: false, melkwargs: {})
+        super()
+
+        raise ArgumentError, "DCT type not supported: #{dct_type}" unless SUPPORTED_DCT_TYPES.include?(dct_type)
+
+        @sample_rate = sample_rate
+        @n_mfcc = n_mfcc
+        @dct_type = dct_type
+        @norm = norm
+        @top_db = 80.0
+        @amplitude_to_db = TorchAudio::Transforms::AmplitudeToDB.new(stype: :power, top_db: @top_db)
+
+        @melspectrogram = TorchAudio::Transforms::MelSpectrogram.new(sample_rate: @sample_rate, **melkwargs)
+
+        raise ArgumentError, "Cannot select more MFCC coefficients than # mel bins" if @n_mfcc > @melspectrogram.n_mels
+
+        dct_mat = F.create_dct(@n_mfcc, @melspectrogram.n_mels, norm: @norm)
+        register_buffer('dct_mat', dct_mat)
+
+        @log_mels = log_mels
+      end
+
+      def forward(waveform)
+        mel_specgram = @melspectrogram.(waveform)
+        if @log_mels
+          mel_specgram = Torch.log(mel_specgram + 1e-6)
+        else
+          mel_specgram = @amplitude_to_db.(mel_specgram)
+        end
+
+        Torch
+          .matmul(mel_specgram.transpose(-2, -1), @dct_mat)
+          .transpose(-2, -1)
+      end
+    end
+  end
+end

data/lib/torchaudio/transforms/mu_law_decoding.rb

@@ -0,0 +1,14 @@
+module TorchAudio
+  module Transforms
+    class MuLawDecoding < Torch::NN::Module
+      def initialize(quantization_channels: 256)
+        super()
+        @quantization_channels = quantization_channels
+      end
+
+      def forward(x_mu)
+        F.mu_law_decoding(x_mu, @quantization_channels)
+      end
+    end
+  end
+end

data/lib/torchaudio/transforms/mu_law_encoding.rb

@@ -0,0 +1,14 @@
+module TorchAudio
+  module Transforms
+    class MuLawEncoding < Torch::NN::Module
+      def initialize(quantization_channels: 256)
+        super()
+        @quantization_channels = quantization_channels
+      end
+
+      def forward(x)
+        F.mu_law_encoding(x, @quantization_channels)
+      end
+    end
+  end
+end

data/lib/torchaudio/transforms/spectrogram.rb

@@ -0,0 +1,27 @@
+module TorchAudio
+  module Transforms
+    class Spectrogram < Torch::NN::Module
+      def initialize(
+        n_fft: 400, win_length: nil, hop_length: nil, pad: 0,
+        window_fn: Torch.method(:hann_window), power: 2.0, normalized: false, wkwargs: nil
+      )
+
+        super()
+        @n_fft = n_fft
+        # number of FFT bins. the returned STFT result will have n_fft // 2 + 1
+        # number of frequecies due to onesided=True in torch.stft
+        @win_length = win_length || n_fft
+        @hop_length = hop_length || @win_length.div(2) # floor division
+        window = wkwargs.nil? ? window_fn.call(@win_length) : window_fn.call(@win_length, **wkwargs)
+        register_buffer("window", window)
+        @pad = pad
+        @power = power
+        @normalized = normalized
+      end
+
+      def forward(waveform)
+        F.spectrogram(waveform, @pad, @window, @n_fft, @hop_length, @win_length, @power, @normalized)
+      end
+    end
+  end
+end

data/lib/torchaudio/transforms/vol.rb

@@ -0,0 +1,31 @@
+module TorchAudio
+  module Transforms
+    class Vol < Torch::NN::Module
+      def initialize(gain, gain_type: "amplitude")
+        super()
+        @gain = gain
+        @gain_type = gain_type
+
+        if ["amplitude", "power"].include?(gain_type) && gain < 0
+          raise ArgumentError, "If gain_type = amplitude or power, gain must be positive."
+        end
+      end
+
+      def forward(waveform)
+        if @gain_type == "amplitude"
+          waveform = waveform * @gain
+        end
+
+        if @gain_type == "db"
+          waveform = F.gain(waveform, @gain)
+        end
+
+        if @gain_type == "power"
+          waveform = F.gain(waveform, 10 * Math.log10(@gain))
+        end
+
+        Torch.clamp(waveform, -1, 1)
+      end
+    end
+  end
+end

data/lib/torchaudio/version.rb

@@ -1,3 +1,3 @@
 module TorchAudio
-  VERSION = "0.1.0"
+  VERSION = "0.2.1"
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: torchaudio
 version: !ruby/object:Gem::Version
-  version: 0.1.0
+  version: 0.2.1
 platform: ruby
 authors:
 - Andrew Kane
-autorequire: 
+autorequire:
 bindir: bin
 cert_chain: []
-date: 2020-08-25 00:00:00.000000000 Z
+date: 2021-07-16 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: torch-rb
@@ -16,86 +16,30 @@ dependencies:
     requirements:
     - - ">="
       - !ruby/object:Gem::Version
-        version: 0.3.2
+        version: 0.7.0
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
     - - ">="
       - !ruby/object:Gem::Version
-        version: 0.3.2
+        version: 0.7.0
 - !ruby/object:Gem::Dependency
   name: rice
   requirement: !ruby/object:Gem::Requirement
     requirements:
     - - ">="
       - !ruby/object:Gem::Version
-        version: '2.2'
+        version: 4.0.2
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
     - - ">="
       - !ruby/object:Gem::Version
-        version: '2.2'
-- !ruby/object:Gem::Dependency
-  name: bundler
-  requirement: !ruby/object:Gem::Requirement
-    requirements:
-    - - ">="
-      - !ruby/object:Gem::Version
-        version: '0'
-  type: :development
-  prerelease: false
-  version_requirements: !ruby/object:Gem::Requirement
-    requirements:
-    - - ">="
-      - !ruby/object:Gem::Version
-        version: '0'
-- !ruby/object:Gem::Dependency
-  name: rake
-  requirement: !ruby/object:Gem::Requirement
-    requirements:
-    - - ">="
-      - !ruby/object:Gem::Version
-        version: '0'
-  type: :development
-  prerelease: false
-  version_requirements: !ruby/object:Gem::Requirement
-    requirements:
-    - - ">="
-      - !ruby/object:Gem::Version
-        version: '0'
-- !ruby/object:Gem::Dependency
-  name: rake-compiler
-  requirement: !ruby/object:Gem::Requirement
-    requirements:
-    - - ">="
-      - !ruby/object:Gem::Version
-        version: '0'
-  type: :development
-  prerelease: false
-  version_requirements: !ruby/object:Gem::Requirement
-    requirements:
-    - - ">="
-      - !ruby/object:Gem::Version
-        version: '0'
-- !ruby/object:Gem::Dependency
-  name: minitest
-  requirement: !ruby/object:Gem::Requirement
-    requirements:
-    - - ">="
-      - !ruby/object:Gem::Version
-        version: '5'
-  type: :development
-  prerelease: false
-  version_requirements: !ruby/object:Gem::Requirement
-    requirements:
-    - - ">="
-      - !ruby/object:Gem::Version
-        version: '5'
-description: 
-email: andrew@chartkick.com
+        version: 4.0.2
+description:
+email: andrew@ankane.org
 executables: []
 extensions:
 - ext/torchaudio/extconf.rb
@@ -118,12 +62,23 @@ files:
 - lib/torchaudio.rb
 - lib/torchaudio/datasets/utils.rb
 - lib/torchaudio/datasets/yesno.rb
+- lib/torchaudio/functional.rb
+- lib/torchaudio/transforms/amplitude_to_db.rb
+- lib/torchaudio/transforms/compute_deltas.rb
+- lib/torchaudio/transforms/fade.rb
+- lib/torchaudio/transforms/mel_scale.rb
+- lib/torchaudio/transforms/mel_spectrogram.rb
+- lib/torchaudio/transforms/mfcc.rb
+- lib/torchaudio/transforms/mu_law_decoding.rb
+- lib/torchaudio/transforms/mu_law_encoding.rb
+- lib/torchaudio/transforms/spectrogram.rb
+- lib/torchaudio/transforms/vol.rb
 - lib/torchaudio/version.rb
 homepage: https://github.com/ankane/torchaudio
 licenses:
 - BSD-2-Clause
 metadata: {}
-post_install_message: 
+post_install_message:
 rdoc_options: []
 require_paths:
 - lib
@@ -131,15 +86,15 @@ required_ruby_version: !ruby/object:Gem::Requirement
   requirements:
   - - ">="
     - !ruby/object:Gem::Version
-      version: '2.5'
+      version: '2.6'
 required_rubygems_version: !ruby/object:Gem::Requirement
   requirements:
   - - ">="
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubygems_version: 3.1.2
-signing_key: 
+rubygems_version: 3.2.22
+signing_key:
 specification_version: 4
 summary: Data manipulation and transformation for audio signal processing
 test_files: []