numo-narray-alt 0.9.3

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA256:
+  metadata.gz: 97dffa7f83cd2acb55f74fead971a6d0d76bdf70b9c1aa17f3e1799bd80cd3cf
+  data.tar.gz: 1bc2bb54d66c86942cffefa80e5ba6bc930453a8b1d4bf6f7925cf82d3964e2a
+SHA512:
+  metadata.gz: d5be63d5efb393e5a9830aca9557b0c3433ff9b8ded112447e7f05d1f71da1d8342032316c3697afc640882d00a094fa25d2d8a15cc2e5e6131c4125b7e40a74
+  data.tar.gz: 11f664a3464b58d9ea29e7c34f34b4e03d1280c8f2228546bba3d9420571c6769e3161f6b969d5f04117fa7079308208a1bec0f573c85ae24b91a7864182770f

data/Gemfile

@@ -0,0 +1,14 @@
+# frozen_string_literal: true
+
+source 'https://rubygems.org'
+
+# Specify your gem's dependencies in narray-devel.gemspec
+gemspec
+
+gem 'rake', '~> 13.3'
+gem 'rake-compiler', '~> 1.3'
+gem 'rubocop', '~> 1.80'
+gem 'rubocop-performance', '~> 1.26'
+gem 'rubocop-rake', '~> 0.7.1'
+gem 'simplecov', '~> 0.22.0'
+gem 'test-unit', '~> 3.7'

data/LICENSE

@@ -0,0 +1,30 @@
+BSD 3-Clause License
+
+Copyright (c) 1999-2025, Masahiro TANAKA
+Copyright (c) 2025, Atsushi TATSUMA
+All rights reserved.
+
+Redistribution and use in source and binary forms, with or without
+modification, are permitted provided that the following conditions are met:
+
+* Redistributions of source code must retain the above copyright notice, this
+  list of conditions and the following disclaimer.
+
+* Redistributions in binary form must reproduce the above copyright notice,
+  this list of conditions and the following disclaimer in the documentation
+  and/or other materials provided with the distribution.
+
+* Neither the name of the copyright holder nor the names of its
+  contributors may be used to endorse or promote products derived from
+  this software without specific prior written permission.
+
+THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
+DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE
+FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
+SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
+CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
+OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.

data/README.md

@@ -0,0 +1,71 @@
+# Numo::NArray Alternative
+
+Numo::NArray Alternative is an experimental project forked from [Numo::NArray](https://github.com/ruby-numo/numo-narray).
+Its goal is to adopt more modern Ruby features and development practices.
+The project owner is also the developer of [Rumale](https://github.com/yoshoku/rumale), a machine learning library,
+and may add new features to this project to support the needs of Rumale.
+In addition, since Rumale depends on Numo::NArray, namespaces, class names, and method names will be kept unchanged
+to allow this library to be used as a drop-in replacement.
+
+The project owner has the utmost respect for Numo::NArray and its creator, Prof. Masahiro Tanaka.
+This project is in no way intended to adversely affect the development of the original Numo::NArray.
+
+## Installation
+
+```shell
+$ gem install numo-narray-alt
+```
+
+## Usage
+
+The same as Numo::NArray.
+
+```irb
+irb(main):001> require 'numo/narray'
+=> true
+irb(main):002> x = Numo::DFloat[[1, 2, 3], [4, 5, 6]]
+=>
+Numo::DFloat#shape=[2,3]
+...
+irb(main):003> x.dot(x.transpose)
+=>
+Numo::DFloat#shape=[2,2]
+[[14, 32],
+ [32, 77]]
+irb(main):004>
+```
+
+## Development
+
+preparation:
+
+```shell
+$ git clone https://github.com/yoshoku/numo-narray-alt
+$ cd numo-narray-alt
+$ bundle install
+```
+
+build and test:
+
+```
+$ bundle exec rake compile
+$ bundle exec rake test
+```
+
+linter:
+
+```shell
+$ bundle exec rubocop
+$ clang-format --dry-run --Werror --style=file ext/**/*.h ext/**/*.c
+```
+
+## Contributing
+
+This project is in an early experimental stage.
+Therefore, we are not accepting any Issues or Pull Requests at this time.
+Thank you for your understanding.
+
+## License
+
+This project was originally forked from [Numo::NArray](https://github.com/ruby-numo/numo-narray).
+It continues to be licensed under the [BSD-3-Clause License](https://github.com/yoshoku/numo-narray-alt/blob/main/LICENSE).

data/Rakefile

@@ -0,0 +1,24 @@
+# frozen_string_literal: true
+
+require 'bundler/gem_tasks'
+
+task :doc do
+  dir = 'ext/numo/narray'
+  src = %w[array.c data.c index.c math.c narray.c rand.c struct.c]
+        .map { |s| File.join(dir, s) } +
+        [File.join(dir, 't_*.c'), 'lib/numo/narray/extra.rb']
+  sh 'cd ext/numo/narray; ruby extconf.rb; make src'
+  sh "rm -rf yard .yardoc; yard doc -o yard -m markdown -r README.md #{src.join(' ')}"
+end
+
+require 'rake/testtask'
+Rake::TestTask.new(:test) do |t|
+  t.libs << 'test'
+  t.libs << 'lib'
+  t.verbose = false
+  t.warning = false
+  t.test_files = FileList['test/**/*_test.rb']
+end
+
+require 'rake/extensiontask'
+Rake::ExtensionTask.new('numo/narray')

data/ext/numo/narray/SFMT-params.h

@@ -0,0 +1,97 @@
+#ifndef SFMT_PARAMS_H
+#define SFMT_PARAMS_H
+
+#if !defined(MEXP)
+// #ifdef __GNUC__
+//   #warning "MEXP is not defined. I assume MEXP is 19937."
+// #endif
+#define MEXP 19937
+#endif
+/*-----------------
+  BASIC DEFINITIONS
+  -----------------*/
+/** Mersenne Exponent. The period of the sequence
+ *  is a multiple of 2^MEXP-1.
+ * #define MEXP 19937 */
+/** SFMT generator has an internal state array of 128-bit integers,
+ * and N is its size. */
+#define N (MEXP / 128 + 1)
+/** N32 is the size of internal state array when regarded as an array
+ * of 32-bit integers.*/
+#define N32 (N * 4)
+/** N64 is the size of internal state array when regarded as an array
+ * of 64-bit integers.*/
+#define N64 (N * 2)
+
+/*----------------------
+  the parameters of SFMT
+  following definitions are in paramsXXXX.h file.
+  ----------------------*/
+/** the pick up position of the array.
+#define POS1 122
+*/
+
+/** the parameter of shift left as four 32-bit registers.
+#define SL1 18
+ */
+
+/** the parameter of shift left as one 128-bit register.
+ * The 128-bit integer is shifted by (SL2 * 8) bits.
+#define SL2 1
+*/
+
+/** the parameter of shift right as four 32-bit registers.
+#define SR1 11
+*/
+
+/** the parameter of shift right as one 128-bit register.
+ * The 128-bit integer is shifted by (SL2 * 8) bits.
+#define SR2 1
+*/
+
+/** A bitmask, used in the recursion.  These parameters are introduced
+ * to break symmetry of SIMD.
+#define MSK1 0xdfffffefU
+#define MSK2 0xddfecb7fU
+#define MSK3 0xbffaffffU
+#define MSK4 0xbffffff6U
+*/
+
+/** These definitions are part of a 128-bit period certification vector.
+#define PARITY1	0x00000001U
+#define PARITY2	0x00000000U
+#define PARITY3	0x00000000U
+#define PARITY4	0xc98e126aU
+*/
+
+#if MEXP == 607
+#include "SFMT-params607.h"
+#elif MEXP == 1279
+#include "SFMT-params1279.h"
+#elif MEXP == 2281
+#include "SFMT-params2281.h"
+#elif MEXP == 4253
+#include "SFMT-params4253.h"
+#elif MEXP == 11213
+#include "SFMT-params11213.h"
+#elif MEXP == 19937
+#include "SFMT-params19937.h"
+#elif MEXP == 44497
+#include "SFMT-params44497.h"
+#elif MEXP == 86243
+#include "SFMT-params86243.h"
+#elif MEXP == 132049
+#include "SFMT-params132049.h"
+#elif MEXP == 216091
+#include "SFMT-params216091.h"
+#else
+#ifdef __GNUC__
+#error "MEXP is not valid."
+#undef MEXP
+#else
+#undef MEXP
+#endif
+
+#endif
+
+#endif /* SFMT_PARAMS_H */

data/ext/numo/narray/SFMT-params19937.h

@@ -0,0 +1,48 @@
+#ifndef SFMT_PARAMS19937_H
+#define SFMT_PARAMS19937_H
+
+#define POS1 122
+#define SL1 18
+#define SL2 1
+#define SR1 11
+#define SR2 1
+#define MSK1 0xdfffffefU
+#define MSK2 0xddfecb7fU
+#define MSK3 0xbffaffffU
+#define MSK4 0xbffffff6U
+#define PARITY1 0x00000001U
+#define PARITY2 0x00000000U
+#define PARITY3 0x00000000U
+#define PARITY4 0x13c9e684U
+
+/* PARAMETERS FOR ALTIVEC */
+#if defined(__APPLE__) /* For OSX */
+#define ALTI_SL1 (vector unsigned int)(SL1, SL1, SL1, SL1)
+#define ALTI_SR1 (vector unsigned int)(SR1, SR1, SR1, SR1)
+#define ALTI_MSK (vector unsigned int)(MSK1, MSK2, MSK3, MSK4)
+#define ALTI_MSK64 (vector unsigned int)(MSK2, MSK1, MSK4, MSK3)
+#define ALTI_SL2_PERM (vector unsigned char)(1, 2, 3, 23, 5, 6, 7, 0, 9, 10, 11, 4, 13, 14, 15, 8)
+#define ALTI_SL2_PERM64 (vector unsigned char)(1, 2, 3, 4, 5, 6, 7, 31, 9, 10, 11, 12, 13, 14, 15, 0)
+#define ALTI_SR2_PERM (vector unsigned char)(7, 0, 1, 2, 11, 4, 5, 6, 15, 8, 9, 10, 17, 12, 13, 14)
+#define ALTI_SR2_PERM64 (vector unsigned char)(15, 0, 1, 2, 3, 4, 5, 6, 17, 8, 9, 10, 11, 12, 13, 14)
+#else /* For OTHER OSs(Linux?) */
+#define ALTI_SL1                                                                                                               \
+  { SL1, SL1, SL1, SL1 }
+#define ALTI_SR1                                                                                                               \
+  { SR1, SR1, SR1, SR1 }
+#define ALTI_MSK                                                                                                               \
+  { MSK1, MSK2, MSK3, MSK4 }
+#define ALTI_MSK64                                                                                                             \
+  { MSK2, MSK1, MSK4, MSK3 }
+#define ALTI_SL2_PERM                                                                                                          \
+  { 1, 2, 3, 23, 5, 6, 7, 0, 9, 10, 11, 4, 13, 14, 15, 8 }
+#define ALTI_SL2_PERM64                                                                                                        \
+  { 1, 2, 3, 4, 5, 6, 7, 31, 9, 10, 11, 12, 13, 14, 15, 0 }
+#define ALTI_SR2_PERM                                                                                                          \
+  { 7, 0, 1, 2, 11, 4, 5, 6, 15, 8, 9, 10, 17, 12, 13, 14 }
+#define ALTI_SR2_PERM64                                                                                                        \
+  { 15, 0, 1, 2, 3, 4, 5, 6, 17, 8, 9, 10, 11, 12, 13, 14 }
+#endif /* For OSX */
+#define IDSTR "SFMT-19937:122-18-1-11-1:dfffffef-ddfecb7f-bffaffff-bffffff6"
+
+#endif /* SFMT_PARAMS19937_H */