numo-pocketfft 0.2.2 → 0.4.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 4bd70edeb6b4247e01c8cc7cf439830a158518550f3c95c4ca5b502ae6ef7c2d
-  data.tar.gz: 2ff73865a15805b7e8e69d55ac87f892d3d0e4aea9b0e2af7cd129b3fa489d63
+  metadata.gz: 9869e7964f7c3ace43a835f43318400b01051a2087f5392445ba22371ca16ef5
+  data.tar.gz: 4b2703db3b805fd135514b651a4d1d09f85dc0cede5333ae0f82faee74a72873
 SHA512:
-  metadata.gz: ef1cfc56b929763495695a28855217292a65126cdbf7134aaa208334722b75456724e0196719932c7fedbcad2dc010e8d34d937c65554429992b623ed6bfdf9f
-  data.tar.gz: '08917a93b37b1c81519a5ad11fde74585077eee3d80761b17a2cbac24b8a67b01ae48f1afdedf8aee8a872af28eaf3df48a6c06006a6b53f2ab0381941bd5573'
+  metadata.gz: dadb8c33ea4eab319eea76546116735be99558dee9c70777c76a77099e45b1689d7293903b875032782942f8e1ec427aa25a655a2a39e2fec2aabcd6ca9ae4f3
+  data.tar.gz: 8d8c2d1b153d84ca3dd03233d612ce28393f76b164ad42b968cbb705083eb0861fcfeb3ab56e09b2fab55ff137dd764745bbb91de0005fc42ce4a9a30b1dde3e

data/CHANGELOG.md

@@ -1,3 +1,19 @@
+# 0.4.0
+- Refactor native extension codes with C99 style.
+- Change not ot use git submodule for pocketfft codes bundle.
+- Introduce conventional commits.
+
+# 0.3.2
+- Update type declaration file.
+- Remove dependent gem's type declaration file from installation files.
+
+# 0.3.1
+- Fix version specifier of runtime dependencies.
+
+# 0.3.0
+- Add type declaration file: sig/numo/pocketfft.rbs
+- Refactor to avoid generating unnecessary arrays.
+
 # 0.2.2
 - Fix bug that caused segmentation fault due to garbage collection ([#4](https://github.com/yoshoku/numo-pocketfft/pull/4)).
 - Fix some configuration files.

data/LICENSE.txt

@@ -1,4 +1,4 @@
-Copyright (c) 2019-2021 Atsushi Tatsuma
+Copyright (c) 2019-2022 Atsushi Tatsuma
 All rights reserved.
 
 Redistribution and use in source and binary forms, with or without

data/README.md

@@ -1,7 +1,7 @@
 # Numo::Pocketfft
 
 [![Build Status](https://github.com/yoshoku/numo-pocketfft/workflows/build/badge.svg)](https://github.com/yoshoku/numo-pocketfft/actions?query=workflow%3Abuild)
-[![Coverage Status](https://coveralls.io/repos/github/yoshoku/numo-pocketfft/badge.svg?branch=master)](https://coveralls.io/github/yoshoku/numo-pocketfft?branch=master)
+[![Coverage Status](https://coveralls.io/repos/github/yoshoku/numo-pocketfft/badge.svg?branch=main)](https://coveralls.io/github/yoshoku/numo-pocketfft?branch=main)
 [![Gem Version](https://badge.fury.io/rb/numo-pocketfft.svg)](https://badge.fury.io/rb/numo-pocketfft)
 [![BSD 3-Clause License](https://img.shields.io/badge/License-BSD%203--Clause-orange.svg)](https://github.com/yoshoku/numo-pocketfft/blob/main/LICENSE.txt)
 [![Documentation](http://img.shields.io/badge/docs-rdoc.info-blue.svg)](https://yoshoku.github.io/numo-pocketfft/doc/)
@@ -59,12 +59,6 @@ Numo::Pocketfft.ifft2(Numo::Pocketfft.fft2(c))
 # [0.794815+0.539948i, 0.201042+0.737815i]]
 ```
 
-## Development
-
-After checking out the repo, run `bin/setup` to install dependencies. Then, run `rake spec` to run the tests. You can also run `bin/console` for an interactive prompt that will allow you to experiment.
-
-To install this gem onto your local machine, run `bundle exec rake install`. To release a new version, update the version number in `version.rb`, and then run `bundle exec rake release`, which will create a git tag for the version, push git commits and tags, and push the `.gem` file to [rubygems.org](https://rubygems.org).
-
 ## Contributing
 
 Bug reports and pull requests are welcome on GitHub at https://github.com/yoshoku/numo-pocketfft. This project is intended to be a safe, welcoming space for collaboration, and contributors are expected to adhere to the [Contributor Covenant](http://contributor-covenant.org) code of conduct.
@@ -72,7 +66,3 @@ Bug reports and pull requests are welcome on GitHub at https://github.com/yoshok
 ## License
 
 The gem is available as open source under the terms of the [BSD-3-Clause License](https://opensource.org/licenses/BSD-3-Clause).
-
-## Code of Conduct
-
-Everyone interacting in the Numo::Pocketfft project’s codebases, issue trackers, chat rooms and mailing lists is expected to follow the [code of conduct](https://github.com/yoshoku/numo-pocketfft/blob/main/CODE_OF_CONDUCT.md).

data/ext/numo/pocketfft/extconf.rb

@@ -28,12 +28,7 @@ end
 
 $CFLAGS = "#{$CFLAGS} -std=c99"
 
-$srcs = Dir.glob("#{$srcdir}/*.c").map { |path| File.basename(path) }
-$srcs << 'pocketfft.c'
-Dir.glob("#{$srcdir}/*/") do |path|
-  dir = File.basename(path)
-  $INCFLAGS << " -I$(srcdir)/#{dir}"
-  $VPATH << "$(srcdir)/#{dir}"
-end
+$srcs = Dir.glob("#{$srcdir}/**/*.c").map { |path| File.basename(path) }
+$VPATH << "$(srcdir)/src"
 
 create_makefile('numo/pocketfft/pocketfftext')

data/ext/numo/pocketfft/pocketfftext.c

@@ -3,22 +3,7 @@
 VALUE mNumo;
 VALUE mPocketfft;
 
-VALUE numo_pocketfft_fft(VALUE x_val, int is_forward)
-{
-  narray_t* x_nary;
-  double* x_pt;
-  size_t length;
-  int n_dims;
-  int n_repeats;
-  int i;
-  int res;
-  int fail;
-  double fct;
-  VALUE z_val;
-  double* z_pt;
-  narray_t* z_nary;
-  cfft_plan plan = NULL;
-
+VALUE numo_pocketfft_fft(VALUE x_val, int is_forward) {
   if (CLASS_OF(x_val) != numo_cDComplex) {
     x_val = rb_funcall(numo_cDComplex, rb_intern("cast"), 1, x_val);
   }
@@ -26,29 +11,30 @@ VALUE numo_pocketfft_fft(VALUE x_val, int is_forward)
     x_val = nary_dup(x_val);
   }
 
+  narray_t* x_nary;
   GetNArray(x_val, x_nary);
-  n_dims = NA_NDIM(x_nary);
-  length = NA_SHAPE(x_nary)[n_dims - 1];
-  x_pt = (double*)na_get_pointer_for_read(x_val);
+  int n_dims = NA_NDIM(x_nary);
+  size_t length = NA_SHAPE(x_nary)[n_dims - 1];
+  double* x_pt = (double*)na_get_pointer_for_read(x_val);
 
-  plan = make_cfft_plan(length);
+  cfft_plan plan = make_cfft_plan(length);
   if (!plan) {
     rb_raise(rb_eNoMemError, "Failed to allocate memory for plan of pocketfft.");
     return Qnil;
   }
 
-  z_val = nary_s_new_like(numo_cDComplex, x_val);
-  z_pt = (double*)na_get_pointer_for_write(z_val);
+  VALUE z_val = nary_s_new_like(numo_cDComplex, x_val);
+  double* z_pt = (double*)na_get_pointer_for_write(z_val);
+  narray_t* z_nary;
   GetNArray(z_val, z_nary);
-  for (i = 0; i < (int)(NA_SIZE(z_nary) * 2); z_pt[i++] = 0.0);
+  memset(z_pt, 0, (NA_SIZE(z_nary) * 2) * sizeof(double));
 
-  fail = 0;
-  fct = is_forward == 1 ? 1.0 : 1.0 / length;
-  n_repeats = (int)(NA_SIZE(x_nary)) / length;
-  for (i = 0; i < n_repeats; i++) {
+  int fail = 0;
+  double fct = is_forward == 1 ? 1.0 : 1.0 / length;
+  int n_repeats = (int)(NA_SIZE(x_nary)) / length;
+  for (int i = 0; i < n_repeats; i++) {
     memcpy(z_pt, x_pt, 2 * length * sizeof(double));
-    res = is_forward == 1 ? cfft_forward(plan, z_pt, fct) : cfft_backward(plan, z_pt, fct);
-    if (res != 0) {
+    if ((is_forward == 1 ? cfft_forward(plan, z_pt, fct) : cfft_backward(plan, z_pt, fct)) != 0) {
       fail = 1;
       break;
     }
@@ -63,8 +49,8 @@ VALUE numo_pocketfft_fft(VALUE x_val, int is_forward)
   RB_GC_GUARD(x_val);
 
   if (fail) {
-    rb_raise(rb_eNoMemError, "Failed to allocate memory in function of pocketfft.");
     rb_funcall(z_val, rb_intern("free"), 0);
+    rb_raise(rb_eNoMemError, "Failed to allocate memory in function of pocketfft.");
     return Qnil;
   }
 
@@ -74,38 +60,17 @@ VALUE numo_pocketfft_fft(VALUE x_val, int is_forward)
 /**
  * @!visibility private
  */
-static VALUE numo_pocketfft_cfft(VALUE self, VALUE x_val)
-{
-  return numo_pocketfft_fft(x_val, 1);
-}
+static VALUE numo_pocketfft_cfft(VALUE self, VALUE x_val) { return numo_pocketfft_fft(x_val, 1); }
 
 /**
  * @!visibility private
  */
-static VALUE numo_pocketfft_icfft(VALUE self, VALUE x_val)
-{
-  return numo_pocketfft_fft(x_val, 0);
-}
+static VALUE numo_pocketfft_icfft(VALUE self, VALUE x_val) { return numo_pocketfft_fft(x_val, 0); }
 
 /**
  * @!visibility private
  */
-static VALUE numo_pocketfft_rfft(VALUE self, VALUE x_val)
-{
-  narray_t* x_nary;
-  double* x_pt;
-  int n_dims;
-  size_t length;
-  int n_repeats;
-  int i;
-  int fail;
-  size_t* z_shape;
-  VALUE z_val;
-  narray_t* z_nary;
-  double* z_pt;
-  int z_step;
-  rfft_plan plan = NULL;
-
+static VALUE numo_pocketfft_rfft(VALUE self, VALUE x_val) {
   if (CLASS_OF(x_val) != numo_cDFloat) {
     x_val = rb_funcall(numo_cDFloat, rb_intern("cast"), 1, x_val);
   }
@@ -113,31 +78,32 @@ static VALUE numo_pocketfft_rfft(VALUE self, VALUE x_val)
     x_val = nary_dup(x_val);
   }
 
+  narray_t* x_nary;
   GetNArray(x_val, x_nary);
-  n_dims = NA_NDIM(x_nary);
-  length = NA_SHAPE(x_nary)[n_dims - 1];
-  x_pt = (double*)na_get_pointer_for_read(x_val);
+  int n_dims = NA_NDIM(x_nary);
+  size_t length = NA_SHAPE(x_nary)[n_dims - 1];
+  double* x_pt = (double*)na_get_pointer_for_read(x_val);
 
-  plan = make_rfft_plan(length);
+  rfft_plan plan = make_rfft_plan(length);
   if (!plan) {
     rb_raise(rb_eNoMemError, "Failed to allocate memory for plan of pocketfft.");
     return Qnil;
   }
 
-  z_shape = ALLOCA_N(size_t, n_dims);
-  for (i = 0; i < n_dims - 1; i++) {
+  size_t* z_shape = ALLOCA_N(size_t, n_dims);
+  for (int i = 0; i < n_dims - 1; i++)
     z_shape[i] = NA_SHAPE(x_nary)[i];
-  }
   z_shape[n_dims - 1] = length / 2 + 1;
-  z_val = rb_narray_new(numo_cDComplex, n_dims, z_shape);
-  z_pt = (double*)na_get_pointer_for_write(z_val);
+  VALUE z_val = rb_narray_new(numo_cDComplex, n_dims, z_shape);
+  double* z_pt = (double*)na_get_pointer_for_write(z_val);
+  narray_t* z_nary;
   GetNArray(z_val, z_nary);
-  for (i = 0; i < (int)(NA_SIZE(z_nary) * 2); z_pt[i++] = 0.0);
+  memset(z_pt, 0, (NA_SIZE(z_nary) * 2) * sizeof(double));
 
-  fail = 0;
-  z_step = (int)(NA_SHAPE(z_nary)[n_dims - 1]) * 2;
-  n_repeats = (int)(NA_SIZE(x_nary)) / length;
-  for (i = 0; i < n_repeats; i++) {
+  int fail = 0;
+  int z_step = (int)(NA_SHAPE(z_nary)[n_dims - 1]) * 2;
+  int n_repeats = (int)(NA_SIZE(x_nary)) / length;
+  for (int i = 0; i < n_repeats; i++) {
     z_pt[z_step - 1] = 0.0;
     memcpy(z_pt + 1, x_pt, length * sizeof(double));
     if (rfft_forward(plan, z_pt + 1, 1.0) != 0) {
@@ -157,8 +123,8 @@ static VALUE numo_pocketfft_rfft(VALUE self, VALUE x_val)
   RB_GC_GUARD(x_val);
 
   if (fail) {
-    rb_raise(rb_eNoMemError, "Failed to allocate memory in function of pocketfft.");
     rb_funcall(z_val, rb_intern("free"), 0);
+    rb_raise(rb_eNoMemError, "Failed to allocate memory in function of pocketfft.");
     return Qnil;
   }
 
@@ -168,22 +134,7 @@ static VALUE numo_pocketfft_rfft(VALUE self, VALUE x_val)
 /**
  * @!visibility private
  */
-static VALUE numo_pocketfft_irfft(VALUE self, VALUE x_val)
-{
-  narray_t* x_nary;
-  double* x_pt;
-  size_t length;
-  int n_dims;
-  int n_repeats;
-  int i;
-  int fail;
-  double fct;
-  size_t* z_shape;
-  VALUE z_val;
-  double* z_pt;
-  narray_t* z_nary;
-  rfft_plan plan = NULL;
-
+static VALUE numo_pocketfft_irfft(VALUE self, VALUE x_val) {
   if (CLASS_OF(x_val) != numo_cDComplex) {
     x_val = rb_funcall(numo_cDComplex, rb_intern("cast"), 1, x_val);
   }
@@ -191,31 +142,33 @@ static VALUE numo_pocketfft_irfft(VALUE self, VALUE x_val)
     x_val = nary_dup(x_val);
   }
 
+  narray_t* x_nary;
   GetNArray(x_val, x_nary);
-  n_dims = NA_NDIM(x_nary);
-  length = NA_SHAPE(x_nary)[n_dims - 1];
-  x_pt = (double*)na_get_pointer_for_read(x_val);
+  int n_dims = NA_NDIM(x_nary);
+  size_t length = NA_SHAPE(x_nary)[n_dims - 1];
+  double* x_pt = (double*)na_get_pointer_for_read(x_val);
 
-  plan = make_rfft_plan(length);
+  rfft_plan plan = make_rfft_plan(length);
   if (!plan) {
     rb_raise(rb_eNoMemError, "Failed to allocate memory for plan of pocketfft.");
     return Qnil;
   }
 
-  z_shape = ALLOCA_N(size_t, n_dims);
-  for (i = 0; i < n_dims - 1; i++) {
+  size_t* z_shape = ALLOCA_N(size_t, n_dims);
+  for (int i = 0; i < n_dims - 1; i++) {
     z_shape[i] = NA_SHAPE(x_nary)[i];
   }
   z_shape[n_dims - 1] = length;
-  z_val = rb_narray_new(numo_cDFloat, n_dims, z_shape);
-  z_pt = (double*)na_get_pointer_for_write(z_val);
+  VALUE z_val = rb_narray_new(numo_cDFloat, n_dims, z_shape);
+  double* z_pt = (double*)na_get_pointer_for_write(z_val);
+  narray_t* z_nary;
   GetNArray(z_val, z_nary);
-  for (i = 0; i < (int)NA_SIZE(z_nary); z_pt[i++] = 0.0);
+  memset(z_pt, 0, NA_SIZE(z_nary) * sizeof(double));
 
-  fail = 0;
-  fct = 1.0 / length;
-  n_repeats = (int)(NA_SIZE(z_nary)) / length;
-  for (i = 0; i < n_repeats; i++) {
+  int fail = 0;
+  double fct = 1.0 / length;
+  int n_repeats = (int)(NA_SIZE(z_nary)) / length;
+  for (int i = 0; i < n_repeats; i++) {
     memcpy(z_pt + 1, x_pt + 2, (length - 1) * sizeof(double));
     z_pt[0] = x_pt[0];
     if (rfft_backward(plan, z_pt, fct) != 0) {
@@ -233,16 +186,15 @@ static VALUE numo_pocketfft_irfft(VALUE self, VALUE x_val)
   RB_GC_GUARD(x_val);
 
   if (fail) {
-    rb_raise(rb_eNoMemError, "Failed to allocate memory in function of pocketfft.");
     rb_funcall(z_val, rb_intern("free"), 0);
+    rb_raise(rb_eNoMemError, "Failed to allocate memory in function of pocketfft.");
     return Qnil;
   }
 
   return z_val;
 }
 
-void Init_pocketfftext()
-{
+void Init_pocketfftext() {
   rb_require("numo/narray");
 
   mNumo = rb_define_module("Numo");

data/ext/numo/pocketfft/pocketfftext.h

@@ -5,9 +5,10 @@
 #include <string.h>
 
 #include <ruby.h>
+
 #include <numo/narray.h>
 #include <numo/template.h>
 
-#include "pocketfft/pocketfft.h"
+#include "src/pocketfft.h"
 
 #endif /* NUMO_POCKETFFT_H */

data/lib/numo/pocketfft/version.rb

@@ -5,6 +5,6 @@ module Numo
   # Numo::Pocketfft is the module that has functions for Fourier transform.
   module Pocketfft
     # The version of Numo::Pocketfft you are using.
-    VERSION = '0.2.2'
+    VERSION = '0.4.0'
   end
 end

data/lib/numo/pocketfft.rb

@@ -64,8 +64,11 @@ module Numo
       raise ArgumentError, 'Expect class of input array to be Numo::NArray.' unless a.is_a?(Numo::NArray)
       raise ArgumentError, 'Expect input array to be non-empty.' if a.empty?
 
-      b = a.dup
-      (0...b.ndim).to_a.reverse.each { |ax_id| b = raw_fft(b, ax_id, inverse: false, real: false) }
+      return raw_fft(a, 0, inverse: false, real: false) if a.ndim == 1
+
+      last_axis_id = a.ndim - 1
+      b = raw_fft(a, last_axis_id, inverse: false, real: false)
+      (last_axis_id - 1).downto(0) { |ax_id| b = raw_fft(b, ax_id, inverse: false, real: false) }
       b
     end
 
@@ -77,8 +80,11 @@ module Numo
       raise ArgumentError, 'Expect class of input array to be Numo::NArray.' unless a.is_a?(Numo::NArray)
       raise ArgumentError, 'Expect input array to be non-empty.' if a.empty?
 
-      b = a.dup
-      (0...b.ndim).to_a.each { |ax_id| b = raw_fft(b, ax_id, inverse: true, real: false) }
+      return raw_fft(a, 0, inverse: true, real: false) if a.ndim == 1
+
+      last_axis_id = a.ndim - 1
+      b = raw_fft(a, 0, inverse: true, real: false)
+      1.upto(last_axis_id) { |ax_id| b = raw_fft(b, ax_id, inverse: true, real: false) }
       b
     end
 
@@ -138,9 +144,11 @@ module Numo
       raise ArgumentError, 'Expect class of input array to be Numo::NArray.' unless a.is_a?(Numo::NArray)
       raise ArgumentError, 'Expect input array to be non-empty.' if a.empty?
 
+      return raw_fft(a, 0, inverse: false, real: true) if a.ndim == 1
+
       last_axis_id = a.ndim - 1
       b = raw_fft(a, last_axis_id, inverse: false, real: true)
-      (0...last_axis_id).to_a.reverse.each { |ax_id| b = raw_fft(b, ax_id, inverse: false, real: false) }
+      (last_axis_id - 1).downto(0) { |ax_id| b = raw_fft(b, ax_id, inverse: false, real: false) }
       b
     end
 
@@ -152,9 +160,11 @@ module Numo
       raise ArgumentError, 'Expect class of input array to be Numo::NArray.' unless a.is_a?(Numo::NArray)
       raise ArgumentError, 'Expect input array to be non-empty.' if a.empty?
 
+      return raw_fft(a, 0, inverse: true, real: true) if a.ndim == 1
+
       last_axis_id = a.ndim - 1
-      b = a.dup
-      (0...last_axis_id).to_a.each { |ax_id| b = raw_fft(b, ax_id, inverse: true, real: false) }
+      b = raw_fft(a, 0, inverse: true, real: false)
+      1.upto(last_axis_id - 1) { |ax_id| b = raw_fft(b, ax_id, inverse: true, real: false) }
       raw_fft(b, last_axis_id, inverse: true, real: true)
     end
 
@@ -209,8 +219,7 @@ module Numo
             b_shape = a.shape
             b_shape[-1] = n
             b = Numo::DComplex.zeros(*b_shape)
-            b_range = [true] * b.ndim
-            b_range[-1] = 0...a.shape[-1]
+            b_range = Array.new(b.ndim) { |idx| idx < b.ndim - 1 ? true : 0...a.shape[-1] }
             b[*b_range] = a
             # inverse of dft for real data
             ext_irfft(b)

data/sig/numo/pocketfft.rbs

@@ -0,0 +1,28 @@
+module Numo
+  module Pocketfft
+    VERSION: String
+
+    def self?.fft: (Numo::DFloat | Numo::DComplex a) -> Numo::DComplex
+    def self?.fft2: (Numo::DFloat | Numo::DComplex a) -> Numo::DComplex
+    def self?.fftn: (Numo::DFloat | Numo::DComplex a) -> Numo::DComplex
+    def self?.ifft: (Numo::DComplex a) -> Numo::DComplex
+    def self?.ifft2: (Numo::DComplex a) -> Numo::DComplex
+    def self?.ifftn: (Numo::DComplex a) -> Numo::DComplex
+    def self?.rfft: (Numo::DFloat a) -> Numo::DComplex
+    def self?.rfft2: (Numo::DFloat a) -> Numo::DComplex
+    def self?.rfftn: (Numo::DFloat a) -> Numo::DComplex
+    def self?.irfft: (Numo::DComplex a) -> Numo::DFloat
+    def self?.irfft2: (Numo::DComplex a) -> Numo::DFloat
+    def self?.irfftn: (Numo::DComplex a) -> Numo::DFloat
+    def self?.fftconvolve: (Numo::DFloat | Numo::DComplex a, Numo::DFloat | Numo::DComplex b) -> Numo::DComplex
+                         | (Numo::DFloat a, Numo::DFloat b) -> Numo::DFloat
+
+    private
+
+    def self?.ext_cfft: (Numo::DComplex a) -> Numo::DComplex
+    def self?.ext_icfft: (Numo::DComplex a) -> Numo::DComplex
+    def self?.ext_rfft: (Numo::DFloat a) -> Numo::DComplex
+    def self?.ext_irfft: (Numo::DComplex a) -> Numo::DFloat
+    def self?.raw_fft: (untyped a, Integer axis_id, inverse: bool inverse, real: bool real) -> untyped
+  end
+end

metadata

@@ -1,27 +1,27 @@
 --- !ruby/object:Gem::Specification
 name: numo-pocketfft
 version: !ruby/object:Gem::Version
-  version: 0.2.2
+  version: 0.4.0
 platform: ruby
 authors:
 - yoshoku
 autorequire:
 bindir: bin
 cert_chain: []
-date: 2021-01-16 00:00:00.000000000 Z
+date: 2022-01-01 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: numo-narray
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - "~>"
+    - - ">="
       - !ruby/object:Gem::Version
         version: 0.9.1
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - "~>"
+    - - ">="
       - !ruby/object:Gem::Version
         version: 0.9.1
 description: 'Numo::Pocketfft provides functions for descrete Fourier Transform based
@@ -35,31 +35,18 @@ extensions:
 - ext/numo/pocketfft/extconf.rb
 extra_rdoc_files: []
 files:
-- ".coveralls.yml"
-- ".github/workflows/build.yml"
-- ".github/workflows/coverage.yml"
-- ".gitignore"
-- ".gitmodules"
-- ".rspec"
 - CHANGELOG.md
-- CODE_OF_CONDUCT.md
-- Gemfile
 - LICENSE.txt
 - README.md
-- Rakefile
 - ext/numo/pocketfft/extconf.rb
-- ext/numo/pocketfft/pocketfft/.gitlab-ci.yml
-- ext/numo/pocketfft/pocketfft/LICENSE.md
-- ext/numo/pocketfft/pocketfft/README.md
-- ext/numo/pocketfft/pocketfft/TESTING
-- ext/numo/pocketfft/pocketfft/ffttest.c
-- ext/numo/pocketfft/pocketfft/pocketfft.c
-- ext/numo/pocketfft/pocketfft/pocketfft.h
 - ext/numo/pocketfft/pocketfftext.c
 - ext/numo/pocketfft/pocketfftext.h
+- ext/numo/pocketfft/src/LICENSE.md
+- ext/numo/pocketfft/src/pocketfft.c
+- ext/numo/pocketfft/src/pocketfft.h
 - lib/numo/pocketfft.rb
 - lib/numo/pocketfft/version.rb
-- numo-pocketfft.gemspec
+- sig/numo/pocketfft.rbs
 homepage: https://github.com/yoshoku/numo-pocketfft
 licenses:
 - BSD-3-Clause
@@ -84,7 +71,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubygems_version: 3.1.4
+rubygems_version: 3.2.32
 signing_key:
 specification_version: 4
 summary: Numo::Pocketfft provides functions for descrete Fourier Transform based on

data/.coveralls.yml

@@ -1 +0,0 @@
-service_name: github-ci

data/.github/workflows/build.yml

@@ -1,27 +0,0 @@
-name: build
-
-on: [push, pull_request]
-
-jobs:
-  build:
-    runs-on: ubuntu-20.04
-    strategy:
-      matrix:
-        ruby: [ '2.5', '2.6', '2.7' ]
-    steps:
-      - uses: actions/checkout@v2
-      - name: Checkout submodule
-        shell: bash
-        run: |
-          auth_header="$(git config --local --get http.https://github.com/.extraheader)"
-          git submodule sync --recursive
-          git submodule update --init --force --recursive
-      - name: Set up Ruby ${{ matrix.ruby }}
-        uses: actions/setup-ruby@v1
-        with:
-          ruby-version: ${{ matrix.ruby }}
-      - name: Build and test with Rake
-        run: |
-          gem install bundler
-          bundle install --jobs 4 --retry 3
-          bundle exec rake

data/.github/workflows/coverage.yml

@@ -1,32 +0,0 @@
-name: coverage
-
-on:
-  push:
-    branches: [ main ]
-  pull_request:
-    branches: [ main ]
-
-jobs:
-  coverage:
-    runs-on: ubuntu-20.04
-    steps:
-      - uses: actions/checkout@v2
-      - name: Checkout submodule
-        shell: bash
-        run: |
-          auth_header="$(git config --local --get http.https://github.com/.extraheader)"
-          git submodule sync --recursive
-          git submodule update --init --force --recursive
-      - name: Set up Ruby 2.7
-        uses: actions/setup-ruby@v1
-        with:
-          ruby-version: '2.7'
-      - name: Build and test with Rake
-        run: |
-          gem install bundler
-          bundle install
-          bundle exec rake
-      - name: Coveralls GitHub Action
-        uses: coverallsapp/github-action@v1.1.2
-        with:
-          github-token: ${{ secrets.GITHUB_TOKEN }}

data/.gitignore

@@ -1,20 +0,0 @@
-/.bundle/
-/.yardoc
-/Gemfile.lock
-/_yardoc/
-/coverage/
-/doc/
-/pkg/
-/spec/reports/
-/tmp/
-/bin/
-
-# rspec failure tracking
-.rspec_status
-
-*.swp
-*.bundle
-mkmf.log
-tags
-.DS_Store
-.ruby-version

data/.gitmodules

@@ -1,4 +0,0 @@
-[submodule "ext/numo/pocketfft/pocketfft"]
-	path = ext/numo/pocketfft/pocketfft
-	url = https://gitlab.mpcdf.mpg.de/mtr/pocketfft.git
-  ignore = dirty

data/.rspec

@@ -1,3 +0,0 @@
---format documentation
---color
---require spec_helper

data/CODE_OF_CONDUCT.md

@@ -1,74 +0,0 @@
-# Contributor Covenant Code of Conduct
-
-## Our Pledge
-
-In the interest of fostering an open and welcoming environment, we as
-contributors and maintainers pledge to making participation in our project and
-our community a harassment-free experience for everyone, regardless of age, body
-size, disability, ethnicity, gender identity and expression, level of experience,
-nationality, personal appearance, race, religion, or sexual identity and
-orientation.
-
-## Our Standards
-
-Examples of behavior that contributes to creating a positive environment
-include:
-
-* Using welcoming and inclusive language
-* Being respectful of differing viewpoints and experiences
-* Gracefully accepting constructive criticism
-* Focusing on what is best for the community
-* Showing empathy towards other community members
-
-Examples of unacceptable behavior by participants include:
-
-* The use of sexualized language or imagery and unwelcome sexual attention or
-advances
-* Trolling, insulting/derogatory comments, and personal or political attacks
-* Public or private harassment
-* Publishing others' private information, such as a physical or electronic
-  address, without explicit permission
-* Other conduct which could reasonably be considered inappropriate in a
-  professional setting
-
-## Our Responsibilities
-
-Project maintainers are responsible for clarifying the standards of acceptable
-behavior and are expected to take appropriate and fair corrective action in
-response to any instances of unacceptable behavior.
-
-Project maintainers have the right and responsibility to remove, edit, or
-reject comments, commits, code, wiki edits, issues, and other contributions
-that are not aligned to this Code of Conduct, or to ban temporarily or
-permanently any contributor for other behaviors that they deem inappropriate,
-threatening, offensive, or harmful.
-
-## Scope
-
-This Code of Conduct applies both within project spaces and in public spaces
-when an individual is representing the project or its community. Examples of
-representing a project or community include using an official project e-mail
-address, posting via an official social media account, or acting as an appointed
-representative at an online or offline event. Representation of a project may be
-further defined and clarified by project maintainers.
-
-## Enforcement
-
-Instances of abusive, harassing, or otherwise unacceptable behavior may be
-reported by contacting the project team at yoshoku@outlook.com. All
-complaints will be reviewed and investigated and will result in a response that
-is deemed necessary and appropriate to the circumstances. The project team is
-obligated to maintain confidentiality with regard to the reporter of an incident.
-Further details of specific enforcement policies may be posted separately.
-
-Project maintainers who do not follow or enforce the Code of Conduct in good
-faith may face temporary or permanent repercussions as determined by other
-members of the project's leadership.
-
-## Attribution
-
-This Code of Conduct is adapted from the [Contributor Covenant][homepage], version 1.4,
-available at [http://contributor-covenant.org/version/1/4][version]
-
-[homepage]: http://contributor-covenant.org
-[version]: http://contributor-covenant.org/version/1/4/

data/Gemfile

@@ -1,10 +0,0 @@
-source "https://rubygems.org"
-
-# Specify your gem's dependencies in numo-pocketfft.gemspec
-gemspec
-
-gem 'rake', '~> 12.0'
-gem 'rake-compiler', '~> 1.0'
-gem 'rspec', '~> 3.0'
-gem 'simplecov', '~> 0.19'
-gem 'simplecov-lcov', '~> 0.8'

data/Rakefile

@@ -1,16 +0,0 @@
-require 'bundler/gem_tasks'
-require 'rspec/core/rake_task'
-
-RSpec::Core::RakeTask.new(:spec)
-
-require 'rake/extensiontask'
-
-task :build => :compile
-
-Rake::ExtensionTask.new('pocketfftext') do |ext|
-  ext.ext_dir = 'ext/numo/pocketfft'
-  ext.lib_dir = 'lib/numo/pocketfft'
-end
-
-task :default => [:clobber, :compile, :spec]
-

data/ext/numo/pocketfft/pocketfft/.gitlab-ci.yml

@@ -1,12 +0,0 @@
-image: debian:testing-slim
-
-stages:
-  - testing
-
-test:
-  stage: testing
-  script:
-    - apt-get update
-    - apt-get install -y gcc
-    - gcc -O3 -W -Wall -std=c99 -march=native *.c -lm
-    - ./a.out

data/ext/numo/pocketfft/pocketfft/README.md

@@ -1,55 +0,0 @@
-PocketFFT
----------
-
-This is a heavily modified implementation of FFTPack [1,2], with the following
-advantages:
-
-- strictly C99 compliant
-- more accurate twiddle factor computation
-- very fast plan generation
-- worst case complexity for transform sizes with large prime factors is
-  `N*log(N)`, because Bluestein's algorithm [3] is used for these cases.
-
-License
--------
-
-3-clause BSD (see LICENSE.md)
-
-
-Some code details
------------------
-
-Twiddle factor computation:
-
-- making use of symmetries to reduce number of sin/cos evaluations
-- all angles are reduced to the range `[0; pi/4]` for higher accuracy
-- an adapted implementation of `sincospi()` is used, which actually computes
-  `sin(x)` and `(cos(x)-1)`.
-- if `n` sin/cos pairs are required, the adjusted `sincospi()` is only called
-  `2*sqrt(n)` times; the remaining values are obtained by evaluating the
-  angle addition theorems in a numerically accurate way.
-
-Parallel invocation:
-
-- Plans only contain read-only data; all temporary arrays are allocated and
-  deallocated during an individual FFT execution. This means that a single plan
-  can be used in several threads at the same time.
-
-Efficient codelets are available for the factors:
-
-- 2, 3, 4, 5, 7, 11 for complex-valued FFTs
-- 2, 3, 4, 5 for real-valued FFTs
-
-Larger prime factors are handled by somewhat less efficient, generic routines.
-
-For lengths with very large prime factors, Bluestein's algorithm is used, and
-instead of an FFT of length `n`, a convolution of length `n2 >= 2*n-1`
-is performed, where `n2` is chosen to be highly composite.
-
-
-[1] Swarztrauber, P. 1982, Vectorizing the Fast Fourier Transforms
-    (New York: Academic Press), 51
-
-[2] https://www.netlib.org/fftpack/
-
-[3] https://en.wikipedia.org/wiki/Chirp_Z-transform

data/ext/numo/pocketfft/pocketfft/TESTING

@@ -1,6 +0,0 @@
-To run a quick test for all FFT sizes up to 8192, do
-
-gcc -std=c99 -O2 -W -Wall *.c -lm
-./a.out
-
-If there is no output, the tests have succeeded.

data/ext/numo/pocketfft/pocketfft/ffttest.c

@@ -1,96 +0,0 @@
-/*
- * This file is part of pocketfft.
- * Licensed under a 3-clause BSD style license - see LICENSE.md
- */
-
-/*
- *  Test codes for pocketfft.
- *
- *  Copyright (C) 2004-2018 Max-Planck-Society
- *  \author Martin Reinecke
- */
-
-#include <stdio.h>
-#include <stdlib.h>
-#include <math.h>
-#include <string.h>
-#include "pocketfft.h"
-
-#define maxlen 8192
-
-static void fill_random (double *data, size_t length)
-  {
-  for (size_t m=0; m<length; ++m)
-    data[m] = rand()/(RAND_MAX+1.0)-0.5;
-  }
-
-static double errcalc (double *data, double *odata, size_t length)
-  {
-  double sum = 0, errsum = 0;
-  for (size_t m=0; m<length; ++m)
-    {
-    errsum += (data[m]-odata[m])*(data[m]-odata[m]);
-    sum += odata[m]*odata[m];
-    }
-  return sqrt(errsum/sum);
-  }
-
-static int test_real(void)
-  {
-  double data[maxlen], odata[maxlen];
-  const double epsilon=2e-15;
-  int ret = 0;
-  fill_random (odata, maxlen);
-  double errsum=0;
-  for (int length=1; length<=maxlen; ++length)
-    {
-    memcpy (data,odata,length*sizeof(double));
-    rfft_plan plan = make_rfft_plan (length);
-    rfft_forward (plan, data, 1.);
-    rfft_backward (plan, data, 1./length);
-    destroy_rfft_plan (plan);
-    double err = errcalc (data, odata, length);
-    if (err>epsilon)
-      {
-      printf("problem at real length %i: %e\n",length,err);
-      ret = 1;
-      }
-    errsum+=err;
-    }
-  printf("errsum: %e\n",errsum);
-  return ret;
-  }
-
-static int test_complex(void)
-  {
-  double data[2*maxlen], odata[2*maxlen];
-  fill_random (odata, 2*maxlen);
-  const double epsilon=2e-15;
-  int ret = 0;
-  double errsum=0;
-  for (int length=1; length<=maxlen; ++length)
-    {
-    memcpy (data,odata,2*length*sizeof(double));
-    cfft_plan plan = make_cfft_plan (length);
-    cfft_forward(plan, data, 1.);
-    cfft_backward(plan, data, 1./length);
-    destroy_cfft_plan (plan);
-    double err = errcalc (data, odata, 2*length);
-    if (err>epsilon)
-      {
-      printf("problem at complex length %i: %e\n",length,err);
-      ret = 1;
-      }
-    errsum+=err;
-    }
-  printf("errsum: %e\n",errsum);
-  return ret;
-  }
-
-int main(void)
-  {
-  int ret = 0;
-  ret = test_real();
-  ret += test_complex();
-  return ret;
-  }

data/numo-pocketfft.gemspec

@@ -1,49 +0,0 @@
-lib = File.expand_path('lib', __dir__)
-$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-require 'numo/pocketfft/version'
-
-Gem::Specification.new do |spec|
-  spec.name          = 'numo-pocketfft'
-  spec.version       = Numo::Pocketfft::VERSION
-  spec.authors       = ['yoshoku']
-  spec.email         = ['yoshoku@outlook.com']
-
-  spec.summary       = <<~MSG
-    Numo::Pocketfft provides functions for descrete Fourier Transform based on pocketfft.
-  MSG
-  spec.description = <<~MSG
-    Numo::Pocketfft provides functions for descrete Fourier Transform based on pocketfft.
-  MSG
-  spec.homepage      = 'https://github.com/yoshoku/numo-pocketfft'
-  spec.license       = 'BSD-3-Clause'
-
-  # Specify which files should be added to the gem when it is released.
-  # The `git ls-files -z` loads the files in the RubyGem that have been added into git.
-  spec.files         = Dir.chdir(File.expand_path(__dir__)) do
-    `git ls-files -z`.split("\x0").reject { |f| f.match(%r{^(test|spec|features)/}) }
-  end
-
-  # Add files in submodule: https://gist.github.com/mattconnolly/5875987
-  gem_dir = __dir__ + '/'
-  `git submodule --quiet foreach pwd`.split($OUTPUT_RECORD_SEPARATOR).each do |submodule_path|
-    Dir.chdir(submodule_path) do
-      submodule_relative_path = submodule_path.sub gem_dir, ''
-      `git ls-files`.split($OUTPUT_RECORD_SEPARATOR).each do |filename|
-        spec.files << "#{submodule_relative_path}/#{filename}"
-      end
-    end
-  end
-
-  spec.require_paths = ['lib']
-  spec.extensions    = ['ext/numo/pocketfft/extconf.rb']
-
-  spec.metadata      = {
-    'homepage_uri' => 'https://github.com/yoshoku/numo-pocketfft',
-    'changelog_uri' => 'https://github.com/yoshoku/numo-pocketfft/blob/main/CHANGELOG.md',
-    'source_code_uri' => 'https://github.com/yoshoku/numo-pocketfft',
-    'documentation_uri' => 'https://yoshoku.github.io/numo-pocketfft/doc/',
-    'bug_tracker_uri' => 'https://github.com/yoshoku/numo-pocketfft/issues'
-  }
-
-  spec.add_runtime_dependency 'numo-narray', '~> 0.9.1'
-end