From 7503302eee01f14848843c27b5afd54a19517f67 Mon Sep 17 00:00:00 2001 From: headlessNode Date: Tue, 29 Sep 2026 15:10:18 +0500 Subject: [PATCH 1/3] feat: add ndarray/base/kernels/generic/nullary-strided1d-by/blocked --- .../nullary-strided1d-by/blocked/README.md | 581 ++++++++++++++++++ .../blocked/docs/repl.txt | 148 +++++ .../blocked/examples/index.js | 120 ++++ .../nullary-strided1d-by/blocked/lib/2d.js | 256 ++++++++ .../nullary-strided1d-by/blocked/lib/3d.js | 279 +++++++++ .../blocked/lib/callback_wrapper.js | 72 +++ .../nullary-strided1d-by/blocked/lib/index.js | 144 +++++ .../nullary-strided1d-by/blocked/lib/main.js | 154 +++++ .../nullary-strided1d-by/blocked/package.json | 65 ++ .../nullary-strided1d-by/blocked/test/test.js | 35 ++ 10 files changed, 1854 insertions(+) create mode 100644 lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/README.md create mode 100644 lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/docs/repl.txt create mode 100644 lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/examples/index.js create mode 100644 lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/2d.js create mode 100644 lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/3d.js create mode 100644 lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/callback_wrapper.js create mode 100644 lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/index.js create mode 100644 lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/main.js create mode 100644 lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/package.json create mode 100644 lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/test/test.js diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/README.md b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/README.md new file mode 100644 index 000000000000..b94588359475 --- /dev/null +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/README.md @@ -0,0 +1,581 @@ + + +# kernel + +> Return a kernel for applying a one-dimensional strided array function to an input ndarray according to a callback function using loop blocking. + +
+ +
+ + + +
+ +## Usage + +```javascript +var kernel = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked' ); +``` + +#### kernel( ndims ) + +Returns a kernel for applying a one-dimensional strided array function to an input ndarray according to a callback function using loop blocking. + + + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); +var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +var strategy = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/strategy' ); + +function clbk( value, idx ) { + return ( idx[ 1 ]*100 ) + ( idx[ 2 ]*10 ) + idx[ 3 ]; +} + +// Create a data buffer: +var xbuf = new Float64Array( 12 ); + +// Define the array shape: +var xsh = [ 1, 3, 2, 2 ]; + +// Define the array strides: +var sx = [ 12, 4, 2, 1 ]; + +// Define the index offset: +var ox = 0; + +// Create an input ndarray descriptor: +var x = { + 'dtype': 'float64', + 'data': xbuf, + 'shape': xsh, + 'strides': sx, + 'offset': ox, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the starting index (inclusive): +var start = { + 'dtype': 'float64', + 'data': new Float64Array( [ 0.0 ] ), + 'shape': [ 1, 3 ], + 'strides': [ 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the ending index (exclusive): +var end = { + 'dtype': 'float64', + 'data': new Float64Array( [ 4.0 ] ), + 'shape': [ 1, 3 ], + 'strides': [ 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Initialize ndarray descriptors representing subarray views: +var views = [ + { + 'dtype': x.dtype, + 'data': x.data, + 'shape': [ 2, 2 ], + 'strides': [ 2, 1 ], + 'offset': x.offset, + 'order': x.order + }, + { + 'dtype': start.dtype, + 'data': start.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': start.offset, + 'order': start.order + }, + { + 'dtype': end.dtype, + 'data': end.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': end.offset, + 'order': end.order + } +]; + +// Resolve an input strategy when iterating over subarray views: +var strategyX = strategy( views[ 0 ] ); + +// Create a workspace array for storing iteration indices: +var ibuf = zeros( xsh.length ); + +// Define the loop and core dimensions: +var ldims = [ 0, 1 ]; +var cdims = [ 2, 3 ]; + +// Resolve a kernel: +var f = kernel( 2 ); + +// Apply strided function: +f( gfillBy, [ x, start, end ], x, views, [ 1, 3 ], [ 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); + +var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +// returns [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] +``` + +The function accepts the following arguments: + +- **ndims**: number of loop dimensions. + +If the function is provided `ndims < 2` or a value greater than the maximum number of supported loop dimensions, the function returns `null`. + +```javascript +var f = kernel( 100000 ); +// returns null +``` + +The returned function accepts the following arguments: + +- **fcn**: function which will be applied to a one-dimensional input subarray according to a callback function. +- **arrays**: array containing an input ndarray [descriptor][@stdlib/ndarray/base/descriptor], followed by any additional ndarray arguments. +- **ref**: original input ndarray provided to a callback function. +- **views**: initialized ndarray [descriptors][@stdlib/ndarray/base/descriptor] representing subarray views. +- **shape**: loop dimensions. +- **stridesX**: loop dimension strides for the input ndarray. +- **strategyX**: strategy for marshaling data to and from an input ndarray view. +- **ibuf**: workspace for storing iteration indices. +- **ldims**: list of loop dimensions. +- **cdims**: list of "core" dimensions. +- **options**: function options which are passed through to `fcn`. +- **hasOpts**: boolean indicating whether to pass an options argument to a strided array function. +- **clbk**: callback function. +- **thisArg**: callback execution context. + + + +#### kernel.kernel2d( fcn, arrays, ref, views, shape, stridesX, strategyX, ibuf, ldims, cdims, options, hasOpts, clbk, thisArg ) + + + +Applies a one-dimensional strided array function to an input ndarray according to a callback function via loop blocking. + + + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); +var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +var strategy = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/strategy' ); + +function clbk( value, idx ) { + return ( idx[ 1 ]*100 ) + ( idx[ 2 ]*10 ) + idx[ 3 ]; +} + +// Create a data buffer: +var xbuf = new Float64Array( 12 ); + +// Define the array shape: +var xsh = [ 1, 3, 2, 2 ]; + +// Define the array strides: +var sx = [ 12, 4, 2, 1 ]; + +// Define the index offset: +var ox = 0; + +// Create an input ndarray descriptor: +var x = { + 'dtype': 'float64', + 'data': xbuf, + 'shape': xsh, + 'strides': sx, + 'offset': ox, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the starting index (inclusive): +var start = { + 'dtype': 'float64', + 'data': new Float64Array( [ 0.0 ] ), + 'shape': [ 1, 3 ], + 'strides': [ 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the ending index (exclusive): +var end = { + 'dtype': 'float64', + 'data': new Float64Array( [ 4.0 ] ), + 'shape': [ 1, 3 ], + 'strides': [ 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Initialize ndarray descriptors representing subarray views: +var views = [ + { + 'dtype': x.dtype, + 'data': x.data, + 'shape': [ 2, 2 ], + 'strides': [ 2, 1 ], + 'offset': x.offset, + 'order': x.order + }, + { + 'dtype': start.dtype, + 'data': start.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': start.offset, + 'order': start.order + }, + { + 'dtype': end.dtype, + 'data': end.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': end.offset, + 'order': end.order + } +]; + +// Resolve an input strategy when iterating over subarray views: +var strategyX = strategy( views[ 0 ] ); + +// Create a workspace array for storing iteration indices: +var ibuf = zeros( xsh.length ); + +// Define the loop and core dimensions: +var ldims = [ 0, 1 ]; +var cdims = [ 2, 3 ]; + +// Apply strided function: +kernel.kernel2d( gfillBy, [ x, start, end ], x, views, [ 1, 3 ], [ 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); + +var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +// returns [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] +``` + +The function has the following parameters: + +- **fcn**: function which will be applied to a one-dimensional input subarray according to a callback function. +- **arrays**: array containing an input ndarray [descriptor][@stdlib/ndarray/base/descriptor], followed by any additional ndarray arguments. +- **ref**: original input ndarray provided to a callback function. +- **views**: initialized ndarray [descriptors][@stdlib/ndarray/base/descriptor] representing subarray views. +- **shape**: loop dimensions. Should have two elements. +- **stridesX**: loop dimension strides for the input ndarray. +- **strategyX**: strategy for marshaling data to and from an input ndarray view. +- **ibuf**: workspace for storing iteration indices. +- **ldims**: list of loop dimensions. +- **cdims**: list of "core" dimensions. +- **options**: function options which are passed through to `fcn`. +- **hasOpts**: boolean indicating whether to pass an options argument to a strided array function. +- **clbk**: callback function. +- **thisArg**: callback execution context. + + + +#### kernel.kernel3d( fcn, arrays, ref, views, shape, stridesX, strategyX, ibuf, ldims, cdims, options, hasOpts, clbk, thisArg ) + + + +Applies a one-dimensional strided array function to an input ndarray according to a callback function via loop blocking. + + + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); +var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +var strategy = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/strategy' ); + +function clbk( value, idx ) { + return ( idx[ 2 ]*100 ) + ( idx[ 3 ]*10 ) + idx[ 4 ]; +} + +// Create a data buffer: +var xbuf = new Float64Array( 12 ); + +// Define the array shape: +var xsh = [ 1, 1, 3, 2, 2 ]; + +// Define the array strides: +var sx = [ 12, 12, 4, 2, 1 ]; + +// Define the index offset: +var ox = 0; + +// Create an input ndarray descriptor: +var x = { + 'dtype': 'float64', + 'data': xbuf, + 'shape': xsh, + 'strides': sx, + 'offset': ox, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the starting index (inclusive): +var start = { + 'dtype': 'float64', + 'data': new Float64Array( [ 0.0 ] ), + 'shape': [ 1, 1, 3 ], + 'strides': [ 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the ending index (exclusive): +var end = { + 'dtype': 'float64', + 'data': new Float64Array( [ 4.0 ] ), + 'shape': [ 1, 1, 3 ], + 'strides': [ 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Initialize ndarray descriptors representing subarray views: +var views = [ + { + 'dtype': x.dtype, + 'data': x.data, + 'shape': [ 2, 2 ], + 'strides': [ 2, 1 ], + 'offset': x.offset, + 'order': x.order + }, + { + 'dtype': start.dtype, + 'data': start.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': start.offset, + 'order': start.order + }, + { + 'dtype': end.dtype, + 'data': end.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': end.offset, + 'order': end.order + } +]; + +// Resolve an input strategy when iterating over subarray views: +var strategyX = strategy( views[ 0 ] ); + +// Create a workspace array for storing iteration indices: +var ibuf = zeros( xsh.length ); + +// Define the loop and core dimensions: +var ldims = [ 0, 1, 2 ]; +var cdims = [ 3, 4 ]; + +// Apply strided function: +kernel.kernel3d( gfillBy, [ x, start, end ], x, views, [ 1, 1, 3 ], [ 12, 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); + +var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +// returns [ [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] ] +``` + +The function has the following parameters: + +- **fcn**: function which will be applied to a one-dimensional input subarray according to a callback function. +- **arrays**: array containing an input ndarray [descriptor][@stdlib/ndarray/base/descriptor], followed by any additional ndarray arguments. +- **ref**: original input ndarray provided to a callback function. +- **views**: initialized ndarray [descriptors][@stdlib/ndarray/base/descriptor] representing subarray views. +- **shape**: loop dimensions. Should have three elements. +- **stridesX**: loop dimension strides for the input ndarray. +- **strategyX**: strategy for marshaling data to and from an input ndarray view. +- **ibuf**: workspace for storing iteration indices. +- **ldims**: list of loop dimensions. +- **cdims**: list of "core" dimensions. +- **options**: function options which are passed through to `fcn`. +- **hasOpts**: boolean indicating whether to pass an options argument to a strided array function. +- **clbk**: callback function. +- **thisArg**: callback execution context. + +
+ + + +
+ +## Notes + +- The strided array function is expected to have the following signature: + + ```text + fcn( arrays[, options], clbk ) + ``` + + where + + - **arrays**: array containing a one-dimensional subarray of the input ndarray and any additional ndarray arguments as subarrays. + - **options**: function options (_optional_). + - **clbk**: callback function. + +- The callback function is provided the following arguments: + + - **value**: current array element. + - **idx**: current array element indices. + - **arr**: the original input ndarray. + +
+ + + +
+ +## Examples + + + + + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); +var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +var strategy = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/strategy' ); +var kernel = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked' ); + +function clbk( value, idx ) { + return ( idx[ 2 ]*100 ) + ( idx[ 3 ]*10 ) + idx[ 4 ]; +} + +// Create a data buffer: +var xbuf = new Float64Array( 12 ); + +// Define the array shape: +var xsh = [ 1, 1, 3, 2, 2 ]; + +// Define the array strides: +var sx = [ 12, 12, 4, 2, 1 ]; + +// Define the index offset: +var ox = 0; + +// Create an input ndarray descriptor: +var x = { + 'dtype': 'float64', + 'data': xbuf, + 'shape': xsh, + 'strides': sx, + 'offset': ox, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the starting index (inclusive): +var start = { + 'dtype': 'float64', + 'data': new Float64Array( [ 0.0 ] ), + 'shape': [ 1, 1, 3 ], + 'strides': [ 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the ending index (exclusive): +var end = { + 'dtype': 'float64', + 'data': new Float64Array( [ 4.0 ] ), + 'shape': [ 1, 1, 3 ], + 'strides': [ 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Initialize ndarray descriptors representing subarray views: +var views = [ + { + 'dtype': x.dtype, + 'data': x.data, + 'shape': [ 2, 2 ], + 'strides': [ 2, 1 ], + 'offset': x.offset, + 'order': x.order + }, + { + 'dtype': start.dtype, + 'data': start.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': start.offset, + 'order': start.order + }, + { + 'dtype': end.dtype, + 'data': end.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': end.offset, + 'order': end.order + } +]; + +// Resolve an input strategy when iterating over subarray views: +var strategyX = strategy( views[ 0 ] ); + +// Create a workspace array for storing iteration indices: +var ibuf = zeros( xsh.length ); + +// Define the loop and core dimensions: +var ldims = [ 0, 1, 2 ]; +var cdims = [ 3, 4 ]; + +// Resolve a kernel: +var f = kernel( 3 ); + +// Apply strided function: +f( gfillBy, [ x, start, end ], x, views, [ 1, 1, 3 ], [ 12, 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); + +console.log( ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ) ); +``` + +
+ + + + + + + + + + + + diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/docs/repl.txt b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/docs/repl.txt new file mode 100644 index 000000000000..7cb907bbedb9 --- /dev/null +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/docs/repl.txt @@ -0,0 +1,148 @@ + +{{alias}}( ndims ) + Returns a kernel for applying a one-dimensional strided array function to an + input ndarray according to a callback function using loop blocking. + + The returned function has the following parameters: + + - fcn: function which will be applied to a one-dimensional input subarray + according to a callback function. + - arrays: array containing an input ndarray descriptor, followed by any + additional ndarray arguments. + - ref: original input ndarray provided to a callback function. + - views: initialized ndarray descriptors representing subarray views. + - shape: loop dimensions. + - stridesX: loop dimension strides for the input ndarray. + - strategyX: strategy for marshaling data to and from an input ndarray view. + - ibuf: workspace for storing iteration indices. + - ldims: list of loop dimensions. + - cdims: list of "core" dimensions. + - options: function options which are passed through to `fcn`. + - hasOpts: boolean indicating whether to pass an options argument to a + strided array function. + - clbk: callback function. + - thisArg: callback execution context. + + If `ndims < 2` or greater than the maximum supported number of loop + dimensions, the function returns `null`. + + Parameters + ---------- + ndims: integer + Number of loop dimensions. + + Returns + ------- + fcn: Function|null + Kernel function. + + Examples + -------- + > var f = {{alias}}( 2 ) + + + +{{alias}}.kernel2d( fcn, arr, ref, v, sh, sx, st, ib, ld, cd, o, po, cb, ctx ) + Applies a one-dimensional strided array function to an input ndarray + according to a callback function via loop blocking. + + Parameters + ---------- + fcn: Function + Wrapper for a one-dimensional strided array function. + + arr: Array + Array containing an input ndarray descriptor, followed by any additional + ndarray arguments. + + ref: ndarray + Original input ndarray provided to a callback function. + + v: Array + Initialized ndarray descriptors representing subarray views. + + sh: Array + Loop dimensions. Should have two elements. + + sx: Array + Loop dimension strides for the input ndarray. + + st: Object + Strategy for marshaling data to and from an input ndarray view. + + ib: Array + Workspace for storing iteration indices. + + ld: Array + List of loop dimensions. + + cd: Array + List of "core" dimensions. + + o: Object + Function options which are passed through to `fcn`. + + po: boolean + Boolean indicating whether to pass an options argument to a strided + array function. + + cb: Function + Callback function. + + ctx: any + Callback execution context. + + +{{alias}}.kernel3d( fcn, arr, ref, v, sh, sx, st, ib, ld, cd, o, po, cb, ctx ) + Applies a one-dimensional strided array function to an input ndarray + according to a callback function via loop blocking. + + Parameters + ---------- + fcn: Function + Wrapper for a one-dimensional strided array function. + + arr: Array + Array containing an input ndarray descriptor, followed by any additional + ndarray arguments. + + ref: ndarray + Original input ndarray provided to a callback function. + + v: Array + Initialized ndarray descriptors representing subarray views. + + sh: Array + Loop dimensions. Should have three elements. + + sx: Array + Loop dimension strides for the input ndarray. + + st: Object + Strategy for marshaling data to and from an input ndarray view. + + ib: Array + Workspace for storing iteration indices. + + ld: Array + List of loop dimensions. + + cd: Array + List of "core" dimensions. + + o: Object + Function options which are passed through to `fcn`. + + po: boolean + Boolean indicating whether to pass an options argument to a strided + array function. + + cb: Function + Callback function. + + ctx: any + Callback execution context. + + See Also + -------- + diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/examples/index.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/examples/index.js new file mode 100644 index 000000000000..1c93f3555ca7 --- /dev/null +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/examples/index.js @@ -0,0 +1,120 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +/* eslint-disable max-len */ + +'use strict'; + +var Float64Array = require( '@stdlib/array/float64' ); +var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +var strategy = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/strategy' ); +var kernel = require( './../lib' ); + +function clbk( value, idx ) { + return ( idx[ 2 ]*100 ) + ( idx[ 3 ]*10 ) + idx[ 4 ]; +} + +// Create a data buffer: +var xbuf = new Float64Array( 12 ); + +// Define the array shape: +var xsh = [ 1, 1, 3, 2, 2 ]; + +// Define the array strides: +var sx = [ 12, 12, 4, 2, 1 ]; + +// Define the index offset: +var ox = 0; + +// Create an input ndarray descriptor: +var x = { + 'dtype': 'float64', + 'data': xbuf, + 'shape': xsh, + 'strides': sx, + 'offset': ox, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the starting index (inclusive): +var start = { + 'dtype': 'float64', + 'data': new Float64Array( [ 0.0 ] ), + 'shape': [ 1, 1, 3 ], + 'strides': [ 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the ending index (exclusive): +var end = { + 'dtype': 'float64', + 'data': new Float64Array( [ 4.0 ] ), + 'shape': [ 1, 1, 3 ], + 'strides': [ 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Initialize ndarray descriptors representing subarray views: +var views = [ + { + 'dtype': x.dtype, + 'data': x.data, + 'shape': [ 2, 2 ], + 'strides': [ 2, 1 ], + 'offset': x.offset, + 'order': x.order + }, + { + 'dtype': start.dtype, + 'data': start.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': start.offset, + 'order': start.order + }, + { + 'dtype': end.dtype, + 'data': end.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': end.offset, + 'order': end.order + } +]; + +// Resolve an input strategy when iterating over subarray views: +var strategyX = strategy( views[ 0 ] ); + +// Create a workspace array for storing iteration indices: +var ibuf = zeros( xsh.length ); + +// Define the loop and core dimensions: +var ldims = [ 0, 1, 2 ]; +var cdims = [ 3, 4 ]; + +// Resolve a kernel: +var f = kernel( 3 ); + +// Apply strided function: +f( gfillBy, [ x, start, end ], x, views, [ 1, 1, 3 ], [ 12, 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); + +console.log( ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ) ); diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/2d.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/2d.js new file mode 100644 index 000000000000..e9e05ffb1859 --- /dev/null +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/2d.js @@ -0,0 +1,256 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +/* eslint-disable max-len, max-params */ + +'use strict'; + +// MODULES // + +var loopOrder = require( '@stdlib/ndarray/base/nullary-loop-interchange-order' ); +var blockSize = require( '@stdlib/ndarray/base/nullary-tiling-block-size' ); +var takeIndexed = require( '@stdlib/array/base/take-indexed' ); +var copyIndexed = require( '@stdlib/array/base/copy-indexed' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var incrementOffsets = require( '@stdlib/ndarray/base/kernels/utils/increment-offsets' ); +var setViewOffsets = require( '@stdlib/ndarray/base/set-descriptor-offsets' ); +var offsets = require( '@stdlib/ndarray/base/offsets' ); +var wrap = require( './callback_wrapper.js' ); + + +// MAIN // + +/** +* Applies a one-dimensional strided array function to an input ndarray according to a callback function via loop blocking. +* +* @private +* @param {Function} fcn - wrapper for a one-dimensional strided array function +* @param {Array} arrays - ndarray descriptors +* @param {ndarrayLike} ref - original input ndarray provided to a callback function +* @param {Array} views - initialized ndarray descriptors representing subarray views +* @param {Array} shape - loop dimensions +* @param {Array} stridesX - loop dimension strides for the input ndarray +* @param {Object} strategyX - strategy for marshaling data to and from an input ndarray view +* @param {NonNegativeIntegerArray} ibuf - workspace for storing iteration indices +* @param {NonNegativeIntegerArray} ldims - list of loop dimensions +* @param {NonNegativeIntegerArray} cdims - list of "core" dimensions +* @param {Options} opts - function options +* @param {boolean} hasOpts - boolean indicating whether to pass an options argument to a strided array function +* @param {Function} clbk - callback function +* @param {*} thisArg - callback execution context +* @returns {void} +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +* var zeros = require( '@stdlib/array/base/zeros' ); +* var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +* var strategy = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/strategy' ); +* +* function clbk( value, idx ) { +* return ( idx[ 1 ]*100 ) + ( idx[ 2 ]*10 ) + idx[ 3 ]; +* } +* +* // Create a data buffer: +* var xbuf = new Float64Array( 12 ); +* +* // Define the array shape: +* var xsh = [ 1, 3, 2, 2 ]; +* +* // Define the array strides: +* var sx = [ 12, 4, 2, 1 ]; +* +* // Define the index offset: +* var ox = 0; +* +* // Create an input ndarray descriptor: +* var x = { +* 'dtype': 'float64', +* 'data': xbuf, +* 'shape': xsh, +* 'strides': sx, +* 'offset': ox, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the starting index (inclusive): +* var start = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 0.0 ] ), +* 'shape': [ 1, 3 ], +* 'strides': [ 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the ending index (exclusive): +* var end = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 4.0 ] ), +* 'shape': [ 1, 3 ], +* 'strides': [ 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Initialize ndarray descriptors representing subarray views: +* var views = [ +* { +* 'dtype': x.dtype, +* 'data': x.data, +* 'shape': [ 2, 2 ], +* 'strides': [ 2, 1 ], +* 'offset': x.offset, +* 'order': x.order +* }, +* { +* 'dtype': start.dtype, +* 'data': start.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': start.offset, +* 'order': start.order +* }, +* { +* 'dtype': end.dtype, +* 'data': end.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': end.offset, +* 'order': end.order +* } +* ]; +* +* // Resolve an input strategy when iterating over subarray views: +* var strategyX = strategy( views[ 0 ] ); +* +* // Create a workspace array for storing iteration indices: +* var ibuf = zeros( xsh.length ); +* +* // Define the loop and core dimensions: +* var ldims = [ 0, 1 ]; +* var cdims = [ 2, 3 ]; +* +* // Apply strided function: +* kernel( gfillBy, [ x, start, end ], x, views, [ 1, 3 ], [ 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); +* +* var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +* // returns [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] +*/ +function kernel( fcn, arrays, ref, views, shape, stridesX, strategyX, ibuf, ldims, cdims, opts, hasOpts, clbk, thisArg ) { + var bsize; + var dv0; + var dv1; + var ov1; + var sh; + var s0; + var s1; + var sv; + var ov; + var iv; + var i0; + var i1; + var j0; + var j1; + var N; + var x; + var v; + var o; + var k; + var f; + + // Note on variable naming convention: S#, dv#, i#, j# where # corresponds to the loop number, with `0` being the innermost loop... + + N = arrays.length; + x = arrays[ 0 ]; + + // Resolve the loop interchange order: + o = loopOrder( shape, stridesX ); + sh = o.sh; + sv = [ o.sx ]; + for ( k = 1; k < N; k++ ) { + sv.push( takeIndexed( arrays[k].strides, o.idx ) ); + } + // Determine the block size: + bsize = blockSize( x.dtype ); + + // Resolve a list of pointers to the first indexed elements in the respective ndarrays: + ov = offsets( arrays ); + + // Cache offset increments for the innermost loop... + dv0 = []; + for ( k = 0; k < N; k++ ) { + dv0.push( sv[k][0] ); + } + // Initialize loop variables... + ov1 = zeros( N ); + dv1 = zeros( N ); + iv = zeros( N ); + + // Shallow copy the list of views to an internal array so that we can update with reshaped views without impacting the original list of views: + v = copyIndexed( views ); + + // Iterate over blocks... + for ( j1 = sh[1]; j1 > 0; ) { + if ( j1 < bsize ) { + s1 = j1; + j1 = 0; + } else { + s1 = bsize; + j1 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov1[ k ] = ov[k] + ( j1*sv[k][1] ); + } + for ( j0 = sh[0]; j0 > 0; ) { + if ( j0 < bsize ) { + s0 = j0; + j0 = 0; + } else { + s0 = bsize; + j0 -= bsize; + } + // Compute index offsets and loop offset increments for the first ndarray elements in the current block... + for ( k = 0; k < N; k++ ) { + iv[ k ] = ov1[k] + ( j0*sv[k][0] ); + dv1[ k ] = sv[k][1] - ( s0*sv[k][0] ); + } + // Iterate over the loop dimensions... + for ( i1 = 0; i1 < s1; i1++ ) { + for ( i0 = 0; i0 < s0; i0++ ) { + setViewOffsets( views, iv ); + v[ 0 ] = strategyX.input( views[ 0 ] ); + f = wrap( ref, views[ 0 ], ibuf, ldims, [ j1+i1, j0+i0 ], cdims, clbk, thisArg ); + if ( hasOpts ) { + fcn( v, opts, f ); + } else { + fcn( v, f ); + } + strategyX.output( views[ 0 ] ); + incrementOffsets( iv, dv0 ); + } + incrementOffsets( iv, dv1 ); + } + } + } +} + + +// EXPORTS // + +module.exports = kernel; diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/3d.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/3d.js new file mode 100644 index 000000000000..d6509e37e361 --- /dev/null +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/3d.js @@ -0,0 +1,279 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +/* eslint-disable max-len, max-params, max-depth */ + +'use strict'; + +// MODULES // + +var loopOrder = require( '@stdlib/ndarray/base/nullary-loop-interchange-order' ); +var blockSize = require( '@stdlib/ndarray/base/nullary-tiling-block-size' ); +var takeIndexed = require( '@stdlib/array/base/take-indexed' ); +var copyIndexed = require( '@stdlib/array/base/copy-indexed' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var incrementOffsets = require( '@stdlib/ndarray/base/kernels/utils/increment-offsets' ); +var setViewOffsets = require( '@stdlib/ndarray/base/set-descriptor-offsets' ); +var offsets = require( '@stdlib/ndarray/base/offsets' ); +var wrap = require( './callback_wrapper.js' ); + + +// MAIN // + +/** +* Applies a one-dimensional strided array function to an input ndarray according to a callback function via loop blocking. +* +* @private +* @param {Function} fcn - wrapper for a one-dimensional strided array function +* @param {Array} arrays - ndarray descriptors +* @param {ndarrayLike} ref - original input ndarray provided to a callback function +* @param {Array} views - initialized ndarray descriptors representing subarray views +* @param {Array} shape - loop dimensions +* @param {Array} stridesX - loop dimension strides for the input ndarray +* @param {Object} strategyX - strategy for marshaling data to and from an input ndarray view +* @param {NonNegativeIntegerArray} ibuf - workspace for storing iteration indices +* @param {NonNegativeIntegerArray} ldims - list of loop dimensions +* @param {NonNegativeIntegerArray} cdims - list of "core" dimensions +* @param {Options} opts - function options +* @param {boolean} hasOpts - boolean indicating whether to pass an options argument to a strided array function +* @param {Function} clbk - callback function +* @param {*} thisArg - callback execution context +* @returns {void} +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +* var zeros = require( '@stdlib/array/base/zeros' ); +* var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +* var strategy = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/strategy' ); +* +* function clbk( value, idx ) { +* return ( idx[ 2 ]*100 ) + ( idx[ 3 ]*10 ) + idx[ 4 ]; +* } +* +* // Create a data buffer: +* var xbuf = new Float64Array( 12 ); +* +* // Define the array shape: +* var xsh = [ 1, 1, 3, 2, 2 ]; +* +* // Define the array strides: +* var sx = [ 12, 12, 4, 2, 1 ]; +* +* // Define the index offset: +* var ox = 0; +* +* // Create an input ndarray descriptor: +* var x = { +* 'dtype': 'float64', +* 'data': xbuf, +* 'shape': xsh, +* 'strides': sx, +* 'offset': ox, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the starting index (inclusive): +* var start = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 0.0 ] ), +* 'shape': [ 1, 1, 3 ], +* 'strides': [ 0, 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the ending index (exclusive): +* var end = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 4.0 ] ), +* 'shape': [ 1, 1, 3 ], +* 'strides': [ 0, 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Initialize ndarray descriptors representing subarray views: +* var views = [ +* { +* 'dtype': x.dtype, +* 'data': x.data, +* 'shape': [ 2, 2 ], +* 'strides': [ 2, 1 ], +* 'offset': x.offset, +* 'order': x.order +* }, +* { +* 'dtype': start.dtype, +* 'data': start.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': start.offset, +* 'order': start.order +* }, +* { +* 'dtype': end.dtype, +* 'data': end.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': end.offset, +* 'order': end.order +* } +* ]; +* +* // Resolve an input strategy when iterating over subarray views: +* var strategyX = strategy( views[ 0 ] ); +* +* // Create a workspace array for storing iteration indices: +* var ibuf = zeros( xsh.length ); +* +* // Define the loop and core dimensions: +* var ldims = [ 0, 1, 2 ]; +* var cdims = [ 3, 4 ]; +* +* // Apply strided function: +* kernel( gfillBy, [ x, start, end ], x, views, [ 1, 1, 3 ], [ 12, 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); +* +* var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +* // returns [ [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] ] +*/ +function kernel( fcn, arrays, ref, views, shape, stridesX, strategyX, ibuf, ldims, cdims, opts, hasOpts, clbk, thisArg ) { + var bsize; + var dv0; + var dv1; + var dv2; + var ov1; + var ov2; + var sh; + var s0; + var s1; + var s2; + var sv; + var ov; + var iv; + var i0; + var i1; + var i2; + var j0; + var j1; + var j2; + var N; + var x; + var v; + var o; + var k; + var f; + + // Note on variable naming convention: S#, dv#, i#, j# where # corresponds to the loop number, with `0` being the innermost loop... + + N = arrays.length; + x = arrays[ 0 ]; + + // Resolve the loop interchange order: + o = loopOrder( shape, stridesX ); + sh = o.sh; + sv = [ o.sx ]; + for ( k = 1; k < N; k++ ) { + sv.push( takeIndexed( arrays[k].strides, o.idx ) ); + } + // Determine the block size: + bsize = blockSize( x.dtype ); + + // Resolve a list of pointers to the first indexed elements in the respective ndarrays: + ov = offsets( arrays ); + + // Cache offset increments for the innermost loop... + dv0 = []; + for ( k = 0; k < N; k++ ) { + dv0.push( sv[k][0] ); + } + // Initialize loop variables... + ov1 = zeros( N ); + ov2 = zeros( N ); + dv1 = zeros( N ); + dv2 = zeros( N ); + iv = zeros( N ); + + // Shallow copy the list of views to an internal array so that we can update with reshaped views without impacting the original list of views: + v = copyIndexed( views ); + + // Iterate over blocks... + for ( j2 = sh[2]; j2 > 0; ) { + if ( j2 < bsize ) { + s2 = j2; + j2 = 0; + } else { + s2 = bsize; + j2 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov2[ k ] = ov[k] + ( j2*sv[k][2] ); + } + for ( j1 = sh[1]; j1 > 0; ) { + if ( j1 < bsize ) { + s1 = j1; + j1 = 0; + } else { + s1 = bsize; + j1 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov1[ k ] = ov2[k] + ( j1*sv[k][1] ); + dv2[ k ] = sv[k][2] - ( s1*sv[k][1] ); + } + for ( j0 = sh[0]; j0 > 0; ) { + if ( j0 < bsize ) { + s0 = j0; + j0 = 0; + } else { + s0 = bsize; + j0 -= bsize; + } + // Compute index offsets and loop offset increments for the first ndarray elements in the current block... + for ( k = 0; k < N; k++ ) { + iv[ k ] = ov1[k] + ( j0*sv[k][0] ); + dv1[ k ] = sv[k][1] - ( s0*sv[k][0] ); + } + // Iterate over the loop dimensions... + for ( i2 = 0; i2 < s2; i2++ ) { + for ( i1 = 0; i1 < s1; i1++ ) { + for ( i0 = 0; i0 < s0; i0++ ) { + setViewOffsets( views, iv ); + v[ 0 ] = strategyX.input( views[ 0 ] ); + f = wrap( ref, views[ 0 ], ibuf, ldims, [ j2+i2, j1+i1, j0+i0 ], cdims, clbk, thisArg ); + if ( hasOpts ) { + fcn( v, opts, f ); + } else { + fcn( v, f ); + } + strategyX.output( views[ 0 ] ); + incrementOffsets( iv, dv0 ); + } + incrementOffsets( iv, dv1 ); + } + incrementOffsets( iv, dv2 ); + } + } + } + } +} + + +// EXPORTS // + +module.exports = kernel; diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/callback_wrapper.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/callback_wrapper.js new file mode 100644 index 000000000000..ad25b5d7bbcf --- /dev/null +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/callback_wrapper.js @@ -0,0 +1,72 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var put = require( '@stdlib/array/base/put' ); +var ind2sub = require( '@stdlib/ndarray/base/ind2sub' ).assign; +var zeros = require( '@stdlib/array/base/zeros' ); + + +// VARIABLES // + +var MODE = 'throw'; + + +// MAIN // + +/** +* Wraps a provided callback function. +* +* @private +* @param {ndarray} arr - input ndarray +* @param {ndarray} view - subarray view of the input ndarray +* @param {NonNegativeIntegerArray} idx - workspace for storing iteration indices +* @param {NonNegativeIntegerArray} ldims - loop dimensions +* @param {NonNegativeIntegerArray} lidx - current loop iteration indices +* @param {NonNegativeIntegerArray} cdims - core dimensions +* @param {Function} clbk - callback function +* @param {thisArg} thisArg - callback execution context +* @returns {Function} callback wrapper +*/ +function wrap( arr, view, idx, ldims, lidx, cdims, clbk, thisArg ) { + var cidx = zeros( cdims.length ); // workspace for storing core iteration indices + put( idx, ldims, lidx, MODE ); + return wrapper; + + /** + * Invokes a callback function. + * + * @private + * @param {*} v - value + * @param {NonNegativeInteger} aidx - current array element index + * @returns {*} result + */ + function wrapper( v, aidx ) { + ind2sub( view.shape, view.strides, view.offset, view.order, aidx, MODE, cidx ); // eslint-disable-line max-len + put( idx, cdims, cidx, MODE ); + return clbk.call( thisArg, v, idx.slice(), arr ); + } +} + + +// EXPORTS // + +module.exports = wrap; diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/index.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/index.js new file mode 100644 index 000000000000..dfdf4d55e482 --- /dev/null +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/index.js @@ -0,0 +1,144 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +/** +* Return a kernel for applying a one-dimensional strided array function to an input ndarray according to a callback function using loop blocking. +* +* @module @stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +* var zeros = require( '@stdlib/array/base/zeros' ); +* var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +* var strategy = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/strategy' ); +* var kernel = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked' ); +* +* function clbk( value, idx ) { +* return ( idx[ 1 ]*100 ) + ( idx[ 2 ]*10 ) + idx[ 3 ]; +* } +* +* // Create a data buffer: +* var xbuf = new Float64Array( 12 ); +* +* // Define the array shape: +* var xsh = [ 1, 3, 2, 2 ]; +* +* // Define the array strides: +* var sx = [ 12, 4, 2, 1 ]; +* +* // Define the index offset: +* var ox = 0; +* +* // Create an input ndarray descriptor: +* var x = { +* 'dtype': 'float64', +* 'data': xbuf, +* 'shape': xsh, +* 'strides': sx, +* 'offset': ox, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the starting index (inclusive): +* var start = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 0.0 ] ), +* 'shape': [ 1, 3 ], +* 'strides': [ 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the ending index (exclusive): +* var end = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 4.0 ] ), +* 'shape': [ 1, 3 ], +* 'strides': [ 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Initialize ndarray descriptors representing subarray views: +* var views = [ +* { +* 'dtype': x.dtype, +* 'data': x.data, +* 'shape': [ 2, 2 ], +* 'strides': [ 2, 1 ], +* 'offset': x.offset, +* 'order': x.order +* }, +* { +* 'dtype': start.dtype, +* 'data': start.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': start.offset, +* 'order': start.order +* }, +* { +* 'dtype': end.dtype, +* 'data': end.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': end.offset, +* 'order': end.order +* } +* ]; +* +* // Resolve an input strategy when iterating over subarray views: +* var strategyX = strategy( views[ 0 ] ); +* +* // Create a workspace array for storing iteration indices: +* var ibuf = zeros( xsh.length ); +* +* // Define the loop and core dimensions: +* var ldims = [ 0, 1 ]; +* var cdims = [ 2, 3 ]; +* +* // Resolve a kernel: +* var f = kernel( 2 ); +* +* // Apply strided function: +* f( gfillBy, [ x, start, end ], x, views, [ 1, 3 ], [ 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); +* +* var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +* // returns [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] +*/ + +// MODULES // + +var setReadOnly = require( '@stdlib/utils/define-nonenumerable-read-only-property' ); +var kernel2d = require( './2d.js' ); +var kernel3d = require( './3d.js' ); +var main = require( './main.js' ); + + +// MAIN // + +setReadOnly( main, 'kernel2d', kernel2d ); +setReadOnly( main, 'kernel3d', kernel3d ); + + +// EXPORTS // + +module.exports = main; diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/main.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/main.js new file mode 100644 index 000000000000..074af133012f --- /dev/null +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/main.js @@ -0,0 +1,154 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var kernel2d = require( './2d.js' ); +var kernel3d = require( './3d.js' ); + + +// VARIABLES // + +var KERNELS = [ + kernel2d, + kernel3d +]; +var MAX_DIMS = KERNELS.length + 1; + + +// MAIN // + +/** +* Returns a kernel for applying a one-dimensional strided array function to an input ndarray according to a callback function using loop blocking. +* +* @param {integer} ndims - number of loop dimensions +* @returns {(Function|null)} kernel function or null +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +* var zeros = require( '@stdlib/array/base/zeros' ); +* var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +* var strategy = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/strategy' ); +* +* function clbk( value, idx ) { +* return ( idx[ 1 ]*100 ) + ( idx[ 2 ]*10 ) + idx[ 3 ]; +* } +* +* // Create a data buffer: +* var xbuf = new Float64Array( 12 ); +* +* // Define the array shape: +* var xsh = [ 1, 3, 2, 2 ]; +* +* // Define the array strides: +* var sx = [ 12, 4, 2, 1 ]; +* +* // Define the index offset: +* var ox = 0; +* +* // Create an input ndarray descriptor: +* var x = { +* 'dtype': 'float64', +* 'data': xbuf, +* 'shape': xsh, +* 'strides': sx, +* 'offset': ox, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the starting index (inclusive): +* var start = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 0.0 ] ), +* 'shape': [ 1, 3 ], +* 'strides': [ 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the ending index (exclusive): +* var end = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 4.0 ] ), +* 'shape': [ 1, 3 ], +* 'strides': [ 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Initialize ndarray descriptors representing subarray views: +* var views = [ +* { +* 'dtype': x.dtype, +* 'data': x.data, +* 'shape': [ 2, 2 ], +* 'strides': [ 2, 1 ], +* 'offset': x.offset, +* 'order': x.order +* }, +* { +* 'dtype': start.dtype, +* 'data': start.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': start.offset, +* 'order': start.order +* }, +* { +* 'dtype': end.dtype, +* 'data': end.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': end.offset, +* 'order': end.order +* } +* ]; +* +* // Resolve an input strategy when iterating over subarray views: +* var strategyX = strategy( views[ 0 ] ); +* +* // Create a workspace array for storing iteration indices: +* var ibuf = zeros( xsh.length ); +* +* // Define the loop and core dimensions: +* var ldims = [ 0, 1 ]; +* var cdims = [ 2, 3 ]; +* +* // Resolve a kernel: +* var f = kernel( 2 ); +* +* // Apply strided function: +* f( gfillBy, [ x, start, end ], x, views, [ 1, 3 ], [ 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); +* +* var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +* // returns [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] +*/ +function kernel( ndims ) { + if ( ndims < 2 || ndims > MAX_DIMS ) { + return null; + } + return KERNELS[ ndims-2 ]; +} + + +// EXPORTS // + +module.exports = kernel; diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/package.json b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/package.json new file mode 100644 index 000000000000..5b9f3be3c9bc --- /dev/null +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/package.json @@ -0,0 +1,65 @@ +{ + "name": "@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked", + "version": "0.0.0", + "description": "Return a kernel for applying a one-dimensional strided array function to an input ndarray according to a callback function using loop blocking.", + "license": "Apache-2.0", + "author": { + "name": "The Stdlib Authors", + "url": "https://github.com/stdlib-js/stdlib/graphs/contributors" + }, + "contributors": [ + { + "name": "The Stdlib Authors", + "url": "https://github.com/stdlib-js/stdlib/graphs/contributors" + } + ], + "main": "./lib", + "directories": { + "doc": "./docs", + "example": "./examples", + "lib": "./lib", + "test": "./test" + }, + "types": "./docs/types", + "scripts": {}, + "homepage": "https://github.com/stdlib-js/stdlib", + "repository": { + "type": "git", + "url": "git://github.com/stdlib-js/stdlib.git" + }, + "bugs": { + "url": "https://github.com/stdlib-js/stdlib/issues" + }, + "dependencies": {}, + "devDependencies": {}, + "engines": { + "node": ">=0.10.0", + "npm": ">2.7.0" + }, + "os": [ + "aix", + "darwin", + "freebsd", + "linux", + "macos", + "openbsd", + "sunos", + "win32", + "windows" + ], + "keywords": [ + "stdlib", + "base", + "strided", + "array", + "ndarray", + "nullary", + "apply", + "fill", + "callback", + "vector", + "kernel", + "blocked" + ], + "__stdlib__": {} +} diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/test/test.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/test/test.js new file mode 100644 index 000000000000..053851cfb7f7 --- /dev/null +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/test/test.js @@ -0,0 +1,35 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var tape = require( 'tape' ); +var kernel = require( './../lib' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof kernel, 'function', 'main export is a function' ); + t.end(); +}); + +// FIXME: add tests From 2b6cb63f9cb4057e8ec60af11737b164d6526f8b Mon Sep 17 00:00:00 2001 From: headlessNode Date: Tue, 29 Sep 2026 15:18:54 +0500 Subject: [PATCH 2/3] fix: repalce strategy --- .../generic/nullary-strided1d-by/blocked/README.md | 8 ++++---- .../nullary-strided1d-by/blocked/examples/index.js | 2 +- .../generic/nullary-strided1d-by/blocked/lib/2d.js | 2 +- .../generic/nullary-strided1d-by/blocked/lib/3d.js | 2 +- .../generic/nullary-strided1d-by/blocked/lib/index.js | 2 +- .../generic/nullary-strided1d-by/blocked/lib/main.js | 2 +- 6 files changed, 9 insertions(+), 9 deletions(-) diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/README.md b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/README.md index b94588359475..290fdbd46356 100644 --- a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/README.md +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/README.md @@ -47,7 +47,7 @@ var Float64Array = require( '@stdlib/array/float64' ); var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); var zeros = require( '@stdlib/array/base/zeros' ); var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); -var strategy = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/strategy' ); +var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); function clbk( value, idx ) { return ( idx[ 1 ]*100 ) + ( idx[ 2 ]*10 ) + idx[ 3 ]; @@ -186,7 +186,7 @@ var Float64Array = require( '@stdlib/array/float64' ); var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); var zeros = require( '@stdlib/array/base/zeros' ); var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); -var strategy = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/strategy' ); +var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); function clbk( value, idx ) { return ( idx[ 1 ]*100 ) + ( idx[ 2 ]*10 ) + idx[ 3 ]; @@ -311,7 +311,7 @@ var Float64Array = require( '@stdlib/array/float64' ); var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); var zeros = require( '@stdlib/array/base/zeros' ); var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); -var strategy = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/strategy' ); +var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); function clbk( value, idx ) { return ( idx[ 2 ]*100 ) + ( idx[ 3 ]*10 ) + idx[ 4 ]; @@ -464,7 +464,7 @@ var Float64Array = require( '@stdlib/array/float64' ); var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); var zeros = require( '@stdlib/array/base/zeros' ); var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); -var strategy = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/strategy' ); +var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); var kernel = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked' ); function clbk( value, idx ) { diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/examples/index.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/examples/index.js index 1c93f3555ca7..dd301b39f5bf 100644 --- a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/examples/index.js +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/examples/index.js @@ -24,7 +24,7 @@ var Float64Array = require( '@stdlib/array/float64' ); var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); var zeros = require( '@stdlib/array/base/zeros' ); var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); -var strategy = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/strategy' ); +var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); var kernel = require( './../lib' ); function clbk( value, idx ) { diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/2d.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/2d.js index e9e05ffb1859..b6fa65946513 100644 --- a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/2d.js +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/2d.js @@ -60,7 +60,7 @@ var wrap = require( './callback_wrapper.js' ); * var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); * var zeros = require( '@stdlib/array/base/zeros' ); * var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); -* var strategy = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/strategy' ); +* var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); * * function clbk( value, idx ) { * return ( idx[ 1 ]*100 ) + ( idx[ 2 ]*10 ) + idx[ 3 ]; diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/3d.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/3d.js index d6509e37e361..0abf606b6c46 100644 --- a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/3d.js +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/3d.js @@ -60,7 +60,7 @@ var wrap = require( './callback_wrapper.js' ); * var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); * var zeros = require( '@stdlib/array/base/zeros' ); * var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); -* var strategy = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/strategy' ); +* var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); * * function clbk( value, idx ) { * return ( idx[ 2 ]*100 ) + ( idx[ 3 ]*10 ) + idx[ 4 ]; diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/index.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/index.js index dfdf4d55e482..9a737f43d174 100644 --- a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/index.js +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/index.js @@ -28,7 +28,7 @@ * var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); * var zeros = require( '@stdlib/array/base/zeros' ); * var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); -* var strategy = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/strategy' ); +* var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); * var kernel = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked' ); * * function clbk( value, idx ) { diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/main.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/main.js index 074af133012f..90714b9c2fea 100644 --- a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/main.js +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/main.js @@ -46,7 +46,7 @@ var MAX_DIMS = KERNELS.length + 1; * var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); * var zeros = require( '@stdlib/array/base/zeros' ); * var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); -* var strategy = require( '@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/strategy' ); +* var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); * * function clbk( value, idx ) { * return ( idx[ 1 ]*100 ) + ( idx[ 2 ]*10 ) + idx[ 3 ]; From 3a1df6c1b595036f8f28874cf1e6f70db528e69c Mon Sep 17 00:00:00 2001 From: headlessNode Date: Sat, 3 Oct 2026 21:14:32 +0500 Subject: [PATCH 3/3] feat: add kernels --- .../nullary-strided1d-by/blocked/README.md | 875 ++++++++++++++++++ .../blocked/docs/repl.txt | 357 +++++++ .../nullary-strided1d-by/blocked/lib/10d.js | 440 +++++++++ .../nullary-strided1d-by/blocked/lib/4d.js | 302 ++++++ .../nullary-strided1d-by/blocked/lib/5d.js | 325 +++++++ .../nullary-strided1d-by/blocked/lib/6d.js | 348 +++++++ .../nullary-strided1d-by/blocked/lib/7d.js | 371 ++++++++ .../nullary-strided1d-by/blocked/lib/8d.js | 394 ++++++++ .../nullary-strided1d-by/blocked/lib/9d.js | 417 +++++++++ .../nullary-strided1d-by/blocked/lib/index.js | 14 + .../nullary-strided1d-by/blocked/lib/main.js | 16 +- 11 files changed, 3858 insertions(+), 1 deletion(-) create mode 100644 lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/10d.js create mode 100644 lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/4d.js create mode 100644 lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/5d.js create mode 100644 lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/6d.js create mode 100644 lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/7d.js create mode 100644 lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/8d.js create mode 100644 lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/9d.js diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/README.md b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/README.md index 290fdbd46356..00c63731b14f 100644 --- a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/README.md +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/README.md @@ -421,6 +421,881 @@ The function has the following parameters: - **clbk**: callback function. - **thisArg**: callback execution context. + + +#### kernel.kernel4d( fcn, arrays, ref, views, shape, stridesX, strategyX, ibuf, ldims, cdims, options, hasOpts, clbk, thisArg ) + + + +Applies a one-dimensional strided array function to an input ndarray according to a callback function via loop blocking. + + + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); +var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); + +function clbk( value, idx ) { + return ( idx[ 3 ]*100 ) + ( idx[ 4 ]*10 ) + idx[ 5 ]; +} + +// Create a data buffer: +var xbuf = new Float64Array( 12 ); + +// Define the array shape: +var xsh = [ 1, 1, 1, 3, 2, 2 ]; + +// Define the array strides: +var sx = [ 12, 12, 12, 4, 2, 1 ]; + +// Define the index offset: +var ox = 0; + +// Create an input ndarray descriptor: +var x = { + 'dtype': 'float64', + 'data': xbuf, + 'shape': xsh, + 'strides': sx, + 'offset': ox, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the starting index (inclusive): +var start = { + 'dtype': 'float64', + 'data': new Float64Array( [ 0.0 ] ), + 'shape': [ 1, 1, 1, 3 ], + 'strides': [ 0, 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the ending index (exclusive): +var end = { + 'dtype': 'float64', + 'data': new Float64Array( [ 4.0 ] ), + 'shape': [ 1, 1, 1, 3 ], + 'strides': [ 0, 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Initialize ndarray descriptors representing subarray views: +var views = [ + { + 'dtype': x.dtype, + 'data': x.data, + 'shape': [ 2, 2 ], + 'strides': [ 2, 1 ], + 'offset': x.offset, + 'order': x.order + }, + { + 'dtype': start.dtype, + 'data': start.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': start.offset, + 'order': start.order + }, + { + 'dtype': end.dtype, + 'data': end.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': end.offset, + 'order': end.order + } +]; + +// Resolve an input strategy when iterating over subarray views: +var strategyX = strategy( views[ 0 ] ); + +// Create a workspace array for storing iteration indices: +var ibuf = zeros( xsh.length ); + +// Define the loop and core dimensions: +var ldims = [ 0, 1, 2, 3 ]; +var cdims = [ 4, 5 ]; + +// Apply strided function: +kernel.kernel4d( gfillBy, [ x, start, end ], x, views, [ 1, 1, 1, 3 ], [ 12, 12, 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); + +var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +// returns [ [ [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] ] ] +``` + +The function has the following parameters: + +- **fcn**: function which will be applied to a one-dimensional input subarray according to a callback function. +- **arrays**: array containing an input ndarray [descriptor][@stdlib/ndarray/base/descriptor], followed by any additional ndarray arguments. +- **ref**: original input ndarray provided to a callback function. +- **views**: initialized ndarray [descriptors][@stdlib/ndarray/base/descriptor] representing subarray views. +- **shape**: loop dimensions. Should have four elements. +- **stridesX**: loop dimension strides for the input ndarray. +- **strategyX**: strategy for marshaling data to and from an input ndarray view. +- **ibuf**: workspace for storing iteration indices. +- **ldims**: list of loop dimensions. +- **cdims**: list of "core" dimensions. +- **options**: function options which are passed through to `fcn`. +- **hasOpts**: boolean indicating whether to pass an options argument to a strided array function. +- **clbk**: callback function. +- **thisArg**: callback execution context. + + + +#### kernel.kernel5d( fcn, arrays, ref, views, shape, stridesX, strategyX, ibuf, ldims, cdims, options, hasOpts, clbk, thisArg ) + + + +Applies a one-dimensional strided array function to an input ndarray according to a callback function via loop blocking. + + + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); +var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); + +function clbk( value, idx ) { + return ( idx[ 4 ]*100 ) + ( idx[ 5 ]*10 ) + idx[ 6 ]; +} + +// Create a data buffer: +var xbuf = new Float64Array( 12 ); + +// Define the array shape: +var xsh = [ 1, 1, 1, 1, 3, 2, 2 ]; + +// Define the array strides: +var sx = [ 12, 12, 12, 12, 4, 2, 1 ]; + +// Define the index offset: +var ox = 0; + +// Create an input ndarray descriptor: +var x = { + 'dtype': 'float64', + 'data': xbuf, + 'shape': xsh, + 'strides': sx, + 'offset': ox, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the starting index (inclusive): +var start = { + 'dtype': 'float64', + 'data': new Float64Array( [ 0.0 ] ), + 'shape': [ 1, 1, 1, 1, 3 ], + 'strides': [ 0, 0, 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the ending index (exclusive): +var end = { + 'dtype': 'float64', + 'data': new Float64Array( [ 4.0 ] ), + 'shape': [ 1, 1, 1, 1, 3 ], + 'strides': [ 0, 0, 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Initialize ndarray descriptors representing subarray views: +var views = [ + { + 'dtype': x.dtype, + 'data': x.data, + 'shape': [ 2, 2 ], + 'strides': [ 2, 1 ], + 'offset': x.offset, + 'order': x.order + }, + { + 'dtype': start.dtype, + 'data': start.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': start.offset, + 'order': start.order + }, + { + 'dtype': end.dtype, + 'data': end.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': end.offset, + 'order': end.order + } +]; + +// Resolve an input strategy when iterating over subarray views: +var strategyX = strategy( views[ 0 ] ); + +// Create a workspace array for storing iteration indices: +var ibuf = zeros( xsh.length ); + +// Define the loop and core dimensions: +var ldims = [ 0, 1, 2, 3, 4 ]; +var cdims = [ 5, 6 ]; + +// Apply strided function: +kernel.kernel5d( gfillBy, [ x, start, end ], x, views, [ 1, 1, 1, 1, 3 ], [ 12, 12, 12, 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); + +var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +// returns [ [ [ [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] ] ] ] +``` + +The function has the following parameters: + +- **fcn**: function which will be applied to a one-dimensional input subarray according to a callback function. +- **arrays**: array containing an input ndarray [descriptor][@stdlib/ndarray/base/descriptor], followed by any additional ndarray arguments. +- **ref**: original input ndarray provided to a callback function. +- **views**: initialized ndarray [descriptors][@stdlib/ndarray/base/descriptor] representing subarray views. +- **shape**: loop dimensions. Should have five elements. +- **stridesX**: loop dimension strides for the input ndarray. +- **strategyX**: strategy for marshaling data to and from an input ndarray view. +- **ibuf**: workspace for storing iteration indices. +- **ldims**: list of loop dimensions. +- **cdims**: list of "core" dimensions. +- **options**: function options which are passed through to `fcn`. +- **hasOpts**: boolean indicating whether to pass an options argument to a strided array function. +- **clbk**: callback function. +- **thisArg**: callback execution context. + + + +#### kernel.kernel6d( fcn, arrays, ref, views, shape, stridesX, strategyX, ibuf, ldims, cdims, options, hasOpts, clbk, thisArg ) + + + +Applies a one-dimensional strided array function to an input ndarray according to a callback function via loop blocking. + + + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); +var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); + +function clbk( value, idx ) { + return ( idx[ 5 ]*100 ) + ( idx[ 6 ]*10 ) + idx[ 7 ]; +} + +// Create a data buffer: +var xbuf = new Float64Array( 12 ); + +// Define the array shape: +var xsh = [ 1, 1, 1, 1, 1, 3, 2, 2 ]; + +// Define the array strides: +var sx = [ 12, 12, 12, 12, 12, 4, 2, 1 ]; + +// Define the index offset: +var ox = 0; + +// Create an input ndarray descriptor: +var x = { + 'dtype': 'float64', + 'data': xbuf, + 'shape': xsh, + 'strides': sx, + 'offset': ox, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the starting index (inclusive): +var start = { + 'dtype': 'float64', + 'data': new Float64Array( [ 0.0 ] ), + 'shape': [ 1, 1, 1, 1, 1, 3 ], + 'strides': [ 0, 0, 0, 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the ending index (exclusive): +var end = { + 'dtype': 'float64', + 'data': new Float64Array( [ 4.0 ] ), + 'shape': [ 1, 1, 1, 1, 1, 3 ], + 'strides': [ 0, 0, 0, 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Initialize ndarray descriptors representing subarray views: +var views = [ + { + 'dtype': x.dtype, + 'data': x.data, + 'shape': [ 2, 2 ], + 'strides': [ 2, 1 ], + 'offset': x.offset, + 'order': x.order + }, + { + 'dtype': start.dtype, + 'data': start.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': start.offset, + 'order': start.order + }, + { + 'dtype': end.dtype, + 'data': end.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': end.offset, + 'order': end.order + } +]; + +// Resolve an input strategy when iterating over subarray views: +var strategyX = strategy( views[ 0 ] ); + +// Create a workspace array for storing iteration indices: +var ibuf = zeros( xsh.length ); + +// Define the loop and core dimensions: +var ldims = [ 0, 1, 2, 3, 4, 5 ]; +var cdims = [ 6, 7 ]; + +// Apply strided function: +kernel.kernel6d( gfillBy, [ x, start, end ], x, views, [ 1, 1, 1, 1, 1, 3 ], [ 12, 12, 12, 12, 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); + +var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +// returns [ [ [ [ [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] ] ] ] ] +``` + +The function has the following parameters: + +- **fcn**: function which will be applied to a one-dimensional input subarray according to a callback function. +- **arrays**: array containing an input ndarray [descriptor][@stdlib/ndarray/base/descriptor], followed by any additional ndarray arguments. +- **ref**: original input ndarray provided to a callback function. +- **views**: initialized ndarray [descriptors][@stdlib/ndarray/base/descriptor] representing subarray views. +- **shape**: loop dimensions. Should have six elements. +- **stridesX**: loop dimension strides for the input ndarray. +- **strategyX**: strategy for marshaling data to and from an input ndarray view. +- **ibuf**: workspace for storing iteration indices. +- **ldims**: list of loop dimensions. +- **cdims**: list of "core" dimensions. +- **options**: function options which are passed through to `fcn`. +- **hasOpts**: boolean indicating whether to pass an options argument to a strided array function. +- **clbk**: callback function. +- **thisArg**: callback execution context. + + + +#### kernel.kernel7d( fcn, arrays, ref, views, shape, stridesX, strategyX, ibuf, ldims, cdims, options, hasOpts, clbk, thisArg ) + + + +Applies a one-dimensional strided array function to an input ndarray according to a callback function via loop blocking. + + + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); +var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); + +function clbk( value, idx ) { + return ( idx[ 6 ]*100 ) + ( idx[ 7 ]*10 ) + idx[ 8 ]; +} + +// Create a data buffer: +var xbuf = new Float64Array( 12 ); + +// Define the array shape: +var xsh = [ 1, 1, 1, 1, 1, 1, 3, 2, 2 ]; + +// Define the array strides: +var sx = [ 12, 12, 12, 12, 12, 12, 4, 2, 1 ]; + +// Define the index offset: +var ox = 0; + +// Create an input ndarray descriptor: +var x = { + 'dtype': 'float64', + 'data': xbuf, + 'shape': xsh, + 'strides': sx, + 'offset': ox, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the starting index (inclusive): +var start = { + 'dtype': 'float64', + 'data': new Float64Array( [ 0.0 ] ), + 'shape': [ 1, 1, 1, 1, 1, 1, 3 ], + 'strides': [ 0, 0, 0, 0, 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the ending index (exclusive): +var end = { + 'dtype': 'float64', + 'data': new Float64Array( [ 4.0 ] ), + 'shape': [ 1, 1, 1, 1, 1, 1, 3 ], + 'strides': [ 0, 0, 0, 0, 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Initialize ndarray descriptors representing subarray views: +var views = [ + { + 'dtype': x.dtype, + 'data': x.data, + 'shape': [ 2, 2 ], + 'strides': [ 2, 1 ], + 'offset': x.offset, + 'order': x.order + }, + { + 'dtype': start.dtype, + 'data': start.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': start.offset, + 'order': start.order + }, + { + 'dtype': end.dtype, + 'data': end.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': end.offset, + 'order': end.order + } +]; + +// Resolve an input strategy when iterating over subarray views: +var strategyX = strategy( views[ 0 ] ); + +// Create a workspace array for storing iteration indices: +var ibuf = zeros( xsh.length ); + +// Define the loop and core dimensions: +var ldims = [ 0, 1, 2, 3, 4, 5, 6 ]; +var cdims = [ 7, 8 ]; + +// Apply strided function: +kernel.kernel7d( gfillBy, [ x, start, end ], x, views, [ 1, 1, 1, 1, 1, 1, 3 ], [ 12, 12, 12, 12, 12, 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); + +var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +// returns [ [ [ [ [ [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] ] ] ] ] ] +``` + +The function has the following parameters: + +- **fcn**: function which will be applied to a one-dimensional input subarray according to a callback function. +- **arrays**: array containing an input ndarray [descriptor][@stdlib/ndarray/base/descriptor], followed by any additional ndarray arguments. +- **ref**: original input ndarray provided to a callback function. +- **views**: initialized ndarray [descriptors][@stdlib/ndarray/base/descriptor] representing subarray views. +- **shape**: loop dimensions. Should have seven elements. +- **stridesX**: loop dimension strides for the input ndarray. +- **strategyX**: strategy for marshaling data to and from an input ndarray view. +- **ibuf**: workspace for storing iteration indices. +- **ldims**: list of loop dimensions. +- **cdims**: list of "core" dimensions. +- **options**: function options which are passed through to `fcn`. +- **hasOpts**: boolean indicating whether to pass an options argument to a strided array function. +- **clbk**: callback function. +- **thisArg**: callback execution context. + + + +#### kernel.kernel8d( fcn, arrays, ref, views, shape, stridesX, strategyX, ibuf, ldims, cdims, options, hasOpts, clbk, thisArg ) + + + +Applies a one-dimensional strided array function to an input ndarray according to a callback function via loop blocking. + + + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); +var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); + +function clbk( value, idx ) { + return ( idx[ 7 ]*100 ) + ( idx[ 8 ]*10 ) + idx[ 9 ]; +} + +// Create a data buffer: +var xbuf = new Float64Array( 12 ); + +// Define the array shape: +var xsh = [ 1, 1, 1, 1, 1, 1, 1, 3, 2, 2 ]; + +// Define the array strides: +var sx = [ 12, 12, 12, 12, 12, 12, 12, 4, 2, 1 ]; + +// Define the index offset: +var ox = 0; + +// Create an input ndarray descriptor: +var x = { + 'dtype': 'float64', + 'data': xbuf, + 'shape': xsh, + 'strides': sx, + 'offset': ox, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the starting index (inclusive): +var start = { + 'dtype': 'float64', + 'data': new Float64Array( [ 0.0 ] ), + 'shape': [ 1, 1, 1, 1, 1, 1, 1, 3 ], + 'strides': [ 0, 0, 0, 0, 0, 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the ending index (exclusive): +var end = { + 'dtype': 'float64', + 'data': new Float64Array( [ 4.0 ] ), + 'shape': [ 1, 1, 1, 1, 1, 1, 1, 3 ], + 'strides': [ 0, 0, 0, 0, 0, 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Initialize ndarray descriptors representing subarray views: +var views = [ + { + 'dtype': x.dtype, + 'data': x.data, + 'shape': [ 2, 2 ], + 'strides': [ 2, 1 ], + 'offset': x.offset, + 'order': x.order + }, + { + 'dtype': start.dtype, + 'data': start.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': start.offset, + 'order': start.order + }, + { + 'dtype': end.dtype, + 'data': end.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': end.offset, + 'order': end.order + } +]; + +// Resolve an input strategy when iterating over subarray views: +var strategyX = strategy( views[ 0 ] ); + +// Create a workspace array for storing iteration indices: +var ibuf = zeros( xsh.length ); + +// Define the loop and core dimensions: +var ldims = [ 0, 1, 2, 3, 4, 5, 6, 7 ]; +var cdims = [ 8, 9 ]; + +// Apply strided function: +kernel.kernel8d( gfillBy, [ x, start, end ], x, views, [ 1, 1, 1, 1, 1, 1, 1, 3 ], [ 12, 12, 12, 12, 12, 12, 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); + +var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +// returns [ [ [ [ [ [ [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] ] ] ] ] ] ] +``` + +The function has the following parameters: + +- **fcn**: function which will be applied to a one-dimensional input subarray according to a callback function. +- **arrays**: array containing an input ndarray [descriptor][@stdlib/ndarray/base/descriptor], followed by any additional ndarray arguments. +- **ref**: original input ndarray provided to a callback function. +- **views**: initialized ndarray [descriptors][@stdlib/ndarray/base/descriptor] representing subarray views. +- **shape**: loop dimensions. Should have eight elements. +- **stridesX**: loop dimension strides for the input ndarray. +- **strategyX**: strategy for marshaling data to and from an input ndarray view. +- **ibuf**: workspace for storing iteration indices. +- **ldims**: list of loop dimensions. +- **cdims**: list of "core" dimensions. +- **options**: function options which are passed through to `fcn`. +- **hasOpts**: boolean indicating whether to pass an options argument to a strided array function. +- **clbk**: callback function. +- **thisArg**: callback execution context. + + + +#### kernel.kernel9d( fcn, arrays, ref, views, shape, stridesX, strategyX, ibuf, ldims, cdims, options, hasOpts, clbk, thisArg ) + + + +Applies a one-dimensional strided array function to an input ndarray according to a callback function via loop blocking. + + + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); +var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); + +function clbk( value, idx ) { + return ( idx[ 8 ]*100 ) + ( idx[ 9 ]*10 ) + idx[ 10 ]; +} + +// Create a data buffer: +var xbuf = new Float64Array( 12 ); + +// Define the array shape: +var xsh = [ 1, 1, 1, 1, 1, 1, 1, 1, 3, 2, 2 ]; + +// Define the array strides: +var sx = [ 12, 12, 12, 12, 12, 12, 12, 12, 4, 2, 1 ]; + +// Define the index offset: +var ox = 0; + +// Create an input ndarray descriptor: +var x = { + 'dtype': 'float64', + 'data': xbuf, + 'shape': xsh, + 'strides': sx, + 'offset': ox, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the starting index (inclusive): +var start = { + 'dtype': 'float64', + 'data': new Float64Array( [ 0.0 ] ), + 'shape': [ 1, 1, 1, 1, 1, 1, 1, 1, 3 ], + 'strides': [ 0, 0, 0, 0, 0, 0, 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the ending index (exclusive): +var end = { + 'dtype': 'float64', + 'data': new Float64Array( [ 4.0 ] ), + 'shape': [ 1, 1, 1, 1, 1, 1, 1, 1, 3 ], + 'strides': [ 0, 0, 0, 0, 0, 0, 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Initialize ndarray descriptors representing subarray views: +var views = [ + { + 'dtype': x.dtype, + 'data': x.data, + 'shape': [ 2, 2 ], + 'strides': [ 2, 1 ], + 'offset': x.offset, + 'order': x.order + }, + { + 'dtype': start.dtype, + 'data': start.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': start.offset, + 'order': start.order + }, + { + 'dtype': end.dtype, + 'data': end.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': end.offset, + 'order': end.order + } +]; + +// Resolve an input strategy when iterating over subarray views: +var strategyX = strategy( views[ 0 ] ); + +// Create a workspace array for storing iteration indices: +var ibuf = zeros( xsh.length ); + +// Define the loop and core dimensions: +var ldims = [ 0, 1, 2, 3, 4, 5, 6, 7, 8 ]; +var cdims = [ 9, 10 ]; + +// Apply strided function: +kernel.kernel9d( gfillBy, [ x, start, end ], x, views, [ 1, 1, 1, 1, 1, 1, 1, 1, 3 ], [ 12, 12, 12, 12, 12, 12, 12, 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); + +var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +// returns [ [ [ [ [ [ [ [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] ] ] ] ] ] ] ] +``` + +The function has the following parameters: + +- **fcn**: function which will be applied to a one-dimensional input subarray according to a callback function. +- **arrays**: array containing an input ndarray [descriptor][@stdlib/ndarray/base/descriptor], followed by any additional ndarray arguments. +- **ref**: original input ndarray provided to a callback function. +- **views**: initialized ndarray [descriptors][@stdlib/ndarray/base/descriptor] representing subarray views. +- **shape**: loop dimensions. Should have nine elements. +- **stridesX**: loop dimension strides for the input ndarray. +- **strategyX**: strategy for marshaling data to and from an input ndarray view. +- **ibuf**: workspace for storing iteration indices. +- **ldims**: list of loop dimensions. +- **cdims**: list of "core" dimensions. +- **options**: function options which are passed through to `fcn`. +- **hasOpts**: boolean indicating whether to pass an options argument to a strided array function. +- **clbk**: callback function. +- **thisArg**: callback execution context. + + + +#### kernel.kernel10d( fcn, arrays, ref, views, shape, stridesX, strategyX, ibuf, ldims, cdims, options, hasOpts, clbk, thisArg ) + + + +Applies a one-dimensional strided array function to an input ndarray according to a callback function via loop blocking. + + + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); +var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); + +function clbk( value, idx ) { + return ( idx[ 9 ]*100 ) + ( idx[ 10 ]*10 ) + idx[ 11 ]; +} + +// Create a data buffer: +var xbuf = new Float64Array( 12 ); + +// Define the array shape: +var xsh = [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 2, 2 ]; + +// Define the array strides: +var sx = [ 12, 12, 12, 12, 12, 12, 12, 12, 12, 4, 2, 1 ]; + +// Define the index offset: +var ox = 0; + +// Create an input ndarray descriptor: +var x = { + 'dtype': 'float64', + 'data': xbuf, + 'shape': xsh, + 'strides': sx, + 'offset': ox, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the starting index (inclusive): +var start = { + 'dtype': 'float64', + 'data': new Float64Array( [ 0.0 ] ), + 'shape': [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 3 ], + 'strides': [ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Create an ndarray descriptor for the ending index (exclusive): +var end = { + 'dtype': 'float64', + 'data': new Float64Array( [ 4.0 ] ), + 'shape': [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 3 ], + 'strides': [ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ], + 'offset': 0, + 'order': 'row-major' +}; + +// Initialize ndarray descriptors representing subarray views: +var views = [ + { + 'dtype': x.dtype, + 'data': x.data, + 'shape': [ 2, 2 ], + 'strides': [ 2, 1 ], + 'offset': x.offset, + 'order': x.order + }, + { + 'dtype': start.dtype, + 'data': start.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': start.offset, + 'order': start.order + }, + { + 'dtype': end.dtype, + 'data': end.data, + 'shape': [], + 'strides': [ 0 ], + 'offset': end.offset, + 'order': end.order + } +]; + +// Resolve an input strategy when iterating over subarray views: +var strategyX = strategy( views[ 0 ] ); + +// Create a workspace array for storing iteration indices: +var ibuf = zeros( xsh.length ); + +// Define the loop and core dimensions: +var ldims = [ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 ]; +var cdims = [ 10, 11 ]; + +// Apply strided function: +kernel.kernel10d( gfillBy, [ x, start, end ], x, views, [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 3 ], [ 12, 12, 12, 12, 12, 12, 12, 12, 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); + +var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +// returns [ [ [ [ [ [ [ [ [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] ] ] ] ] ] ] ] ] +``` + +The function has the following parameters: + +- **fcn**: function which will be applied to a one-dimensional input subarray according to a callback function. +- **arrays**: array containing an input ndarray [descriptor][@stdlib/ndarray/base/descriptor], followed by any additional ndarray arguments. +- **ref**: original input ndarray provided to a callback function. +- **views**: initialized ndarray [descriptors][@stdlib/ndarray/base/descriptor] representing subarray views. +- **shape**: loop dimensions. Should have ten elements. +- **stridesX**: loop dimension strides for the input ndarray. +- **strategyX**: strategy for marshaling data to and from an input ndarray view. +- **ibuf**: workspace for storing iteration indices. +- **ldims**: list of loop dimensions. +- **cdims**: list of "core" dimensions. +- **options**: function options which are passed through to `fcn`. +- **hasOpts**: boolean indicating whether to pass an options argument to a strided array function. +- **clbk**: callback function. +- **thisArg**: callback execution context. + diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/docs/repl.txt b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/docs/repl.txt index 7cb907bbedb9..d66f6b222245 100644 --- a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/docs/repl.txt +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/docs/repl.txt @@ -143,6 +143,363 @@ ctx: any Callback execution context. + +{{alias}}.kernel4d( fcn, arr, ref, v, sh, sx, st, ib, ld, cd, o, po, cb, ctx ) + Applies a one-dimensional strided array function to an input ndarray + according to a callback function via loop blocking. + + Parameters + ---------- + fcn: Function + Wrapper for a one-dimensional strided array function. + + arr: Array + Array containing an input ndarray descriptor, followed by any additional + ndarray arguments. + + ref: ndarray + Original input ndarray provided to a callback function. + + v: Array + Initialized ndarray descriptors representing subarray views. + + sh: Array + Loop dimensions. Should have four elements. + + sx: Array + Loop dimension strides for the input ndarray. + + st: Object + Strategy for marshaling data to and from an input ndarray view. + + ib: Array + Workspace for storing iteration indices. + + ld: Array + List of loop dimensions. + + cd: Array + List of "core" dimensions. + + o: Object + Function options which are passed through to `fcn`. + + po: boolean + Boolean indicating whether to pass an options argument to a strided + array function. + + cb: Function + Callback function. + + ctx: any + Callback execution context. + + +{{alias}}.kernel5d( fcn, arr, ref, v, sh, sx, st, ib, ld, cd, o, po, cb, ctx ) + Applies a one-dimensional strided array function to an input ndarray + according to a callback function via loop blocking. + + Parameters + ---------- + fcn: Function + Wrapper for a one-dimensional strided array function. + + arr: Array + Array containing an input ndarray descriptor, followed by any additional + ndarray arguments. + + ref: ndarray + Original input ndarray provided to a callback function. + + v: Array + Initialized ndarray descriptors representing subarray views. + + sh: Array + Loop dimensions. Should have five elements. + + sx: Array + Loop dimension strides for the input ndarray. + + st: Object + Strategy for marshaling data to and from an input ndarray view. + + ib: Array + Workspace for storing iteration indices. + + ld: Array + List of loop dimensions. + + cd: Array + List of "core" dimensions. + + o: Object + Function options which are passed through to `fcn`. + + po: boolean + Boolean indicating whether to pass an options argument to a strided + array function. + + cb: Function + Callback function. + + ctx: any + Callback execution context. + + +{{alias}}.kernel6d( fcn, arr, ref, v, sh, sx, st, ib, ld, cd, o, po, cb, ctx ) + Applies a one-dimensional strided array function to an input ndarray + according to a callback function via loop blocking. + + Parameters + ---------- + fcn: Function + Wrapper for a one-dimensional strided array function. + + arr: Array + Array containing an input ndarray descriptor, followed by any additional + ndarray arguments. + + ref: ndarray + Original input ndarray provided to a callback function. + + v: Array + Initialized ndarray descriptors representing subarray views. + + sh: Array + Loop dimensions. Should have six elements. + + sx: Array + Loop dimension strides for the input ndarray. + + st: Object + Strategy for marshaling data to and from an input ndarray view. + + ib: Array + Workspace for storing iteration indices. + + ld: Array + List of loop dimensions. + + cd: Array + List of "core" dimensions. + + o: Object + Function options which are passed through to `fcn`. + + po: boolean + Boolean indicating whether to pass an options argument to a strided + array function. + + cb: Function + Callback function. + + ctx: any + Callback execution context. + + +{{alias}}.kernel7d( fcn, arr, ref, v, sh, sx, st, ib, ld, cd, o, po, cb, ctx ) + Applies a one-dimensional strided array function to an input ndarray + according to a callback function via loop blocking. + + Parameters + ---------- + fcn: Function + Wrapper for a one-dimensional strided array function. + + arr: Array + Array containing an input ndarray descriptor, followed by any additional + ndarray arguments. + + ref: ndarray + Original input ndarray provided to a callback function. + + v: Array + Initialized ndarray descriptors representing subarray views. + + sh: Array + Loop dimensions. Should have seven elements. + + sx: Array + Loop dimension strides for the input ndarray. + + st: Object + Strategy for marshaling data to and from an input ndarray view. + + ib: Array + Workspace for storing iteration indices. + + ld: Array + List of loop dimensions. + + cd: Array + List of "core" dimensions. + + o: Object + Function options which are passed through to `fcn`. + + po: boolean + Boolean indicating whether to pass an options argument to a strided + array function. + + cb: Function + Callback function. + + ctx: any + Callback execution context. + + +{{alias}}.kernel8d( fcn, arr, ref, v, sh, sx, st, ib, ld, cd, o, po, cb, ctx ) + Applies a one-dimensional strided array function to an input ndarray + according to a callback function via loop blocking. + + Parameters + ---------- + fcn: Function + Wrapper for a one-dimensional strided array function. + + arr: Array + Array containing an input ndarray descriptor, followed by any additional + ndarray arguments. + + ref: ndarray + Original input ndarray provided to a callback function. + + v: Array + Initialized ndarray descriptors representing subarray views. + + sh: Array + Loop dimensions. Should have eight elements. + + sx: Array + Loop dimension strides for the input ndarray. + + st: Object + Strategy for marshaling data to and from an input ndarray view. + + ib: Array + Workspace for storing iteration indices. + + ld: Array + List of loop dimensions. + + cd: Array + List of "core" dimensions. + + o: Object + Function options which are passed through to `fcn`. + + po: boolean + Boolean indicating whether to pass an options argument to a strided + array function. + + cb: Function + Callback function. + + ctx: any + Callback execution context. + + +{{alias}}.kernel9d( fcn, arr, ref, v, sh, sx, st, ib, ld, cd, o, po, cb, ctx ) + Applies a one-dimensional strided array function to an input ndarray + according to a callback function via loop blocking. + + Parameters + ---------- + fcn: Function + Wrapper for a one-dimensional strided array function. + + arr: Array + Array containing an input ndarray descriptor, followed by any additional + ndarray arguments. + + ref: ndarray + Original input ndarray provided to a callback function. + + v: Array + Initialized ndarray descriptors representing subarray views. + + sh: Array + Loop dimensions. Should have nine elements. + + sx: Array + Loop dimension strides for the input ndarray. + + st: Object + Strategy for marshaling data to and from an input ndarray view. + + ib: Array + Workspace for storing iteration indices. + + ld: Array + List of loop dimensions. + + cd: Array + List of "core" dimensions. + + o: Object + Function options which are passed through to `fcn`. + + po: boolean + Boolean indicating whether to pass an options argument to a strided + array function. + + cb: Function + Callback function. + + ctx: any + Callback execution context. + + +{{alias}}.kernel10d( fcn, arr, ref, v, sh, sx, st, ib, ld, cd, o, po, cb, ctx ) + Applies a one-dimensional strided array function to an input ndarray + according to a callback function via loop blocking. + + Parameters + ---------- + fcn: Function + Wrapper for a one-dimensional strided array function. + + arr: Array + Array containing an input ndarray descriptor, followed by any additional + ndarray arguments. + + ref: ndarray + Original input ndarray provided to a callback function. + + v: Array + Initialized ndarray descriptors representing subarray views. + + sh: Array + Loop dimensions. Should have ten elements. + + sx: Array + Loop dimension strides for the input ndarray. + + st: Object + Strategy for marshaling data to and from an input ndarray view. + + ib: Array + Workspace for storing iteration indices. + + ld: Array + List of loop dimensions. + + cd: Array + List of "core" dimensions. + + o: Object + Function options which are passed through to `fcn`. + + po: boolean + Boolean indicating whether to pass an options argument to a strided + array function. + + cb: Function + Callback function. + + ctx: any + Callback execution context. + See Also -------- diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/10d.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/10d.js new file mode 100644 index 000000000000..cddb238b97d4 --- /dev/null +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/10d.js @@ -0,0 +1,440 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +/* eslint-disable max-len, max-params, max-depth */ + +'use strict'; + +// MODULES // + +var loopOrder = require( '@stdlib/ndarray/base/nullary-loop-interchange-order' ); +var blockSize = require( '@stdlib/ndarray/base/nullary-tiling-block-size' ); +var takeIndexed = require( '@stdlib/array/base/take-indexed' ); +var copyIndexed = require( '@stdlib/array/base/copy-indexed' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var incrementOffsets = require( '@stdlib/ndarray/base/kernels/utils/increment-offsets' ); +var setViewOffsets = require( '@stdlib/ndarray/base/set-descriptor-offsets' ); +var offsets = require( '@stdlib/ndarray/base/offsets' ); +var wrap = require( './callback_wrapper.js' ); + + +// MAIN // + +/** +* Applies a one-dimensional strided array function to an input ndarray according to a callback function via loop blocking. +* +* @private +* @param {Function} fcn - wrapper for a one-dimensional strided array function +* @param {Array} arrays - ndarray descriptors +* @param {ndarrayLike} ref - original input ndarray provided to a callback function +* @param {Array} views - initialized ndarray descriptors representing subarray views +* @param {Array} shape - loop dimensions +* @param {Array} stridesX - loop dimension strides for the input ndarray +* @param {Object} strategyX - strategy for marshaling data to and from an input ndarray view +* @param {NonNegativeIntegerArray} ibuf - workspace for storing iteration indices +* @param {NonNegativeIntegerArray} ldims - list of loop dimensions +* @param {NonNegativeIntegerArray} cdims - list of "core" dimensions +* @param {Options} opts - function options +* @param {boolean} hasOpts - boolean indicating whether to pass an options argument to a strided array function +* @param {Function} clbk - callback function +* @param {*} thisArg - callback execution context +* @returns {void} +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +* var zeros = require( '@stdlib/array/base/zeros' ); +* var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +* var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); +* +* function clbk( value, idx ) { +* return ( idx[ 9 ]*100 ) + ( idx[ 10 ]*10 ) + idx[ 11 ]; +* } +* +* // Create a data buffer: +* var xbuf = new Float64Array( 12 ); +* +* // Define the array shape: +* var xsh = [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 2, 2 ]; +* +* // Define the array strides: +* var sx = [ 12, 12, 12, 12, 12, 12, 12, 12, 12, 4, 2, 1 ]; +* +* // Define the index offset: +* var ox = 0; +* +* // Create an input ndarray descriptor: +* var x = { +* 'dtype': 'float64', +* 'data': xbuf, +* 'shape': xsh, +* 'strides': sx, +* 'offset': ox, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the starting index (inclusive): +* var start = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 0.0 ] ), +* 'shape': [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 3 ], +* 'strides': [ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the ending index (exclusive): +* var end = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 4.0 ] ), +* 'shape': [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 3 ], +* 'strides': [ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Initialize ndarray descriptors representing subarray views: +* var views = [ +* { +* 'dtype': x.dtype, +* 'data': x.data, +* 'shape': [ 2, 2 ], +* 'strides': [ 2, 1 ], +* 'offset': x.offset, +* 'order': x.order +* }, +* { +* 'dtype': start.dtype, +* 'data': start.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': start.offset, +* 'order': start.order +* }, +* { +* 'dtype': end.dtype, +* 'data': end.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': end.offset, +* 'order': end.order +* } +* ]; +* +* // Resolve an input strategy when iterating over subarray views: +* var strategyX = strategy( views[ 0 ] ); +* +* // Create a workspace array for storing iteration indices: +* var ibuf = zeros( xsh.length ); +* +* // Define the loop and core dimensions: +* var ldims = [ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 ]; +* var cdims = [ 10, 11 ]; +* +* // Apply strided function: +* kernel( gfillBy, [ x, start, end ], x, views, [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 3 ], [ 12, 12, 12, 12, 12, 12, 12, 12, 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); +* +* var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +* // returns [ [ [ [ [ [ [ [ [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] ] ] ] ] ] ] ] ] +*/ +function kernel( fcn, arrays, ref, views, shape, stridesX, strategyX, ibuf, ldims, cdims, opts, hasOpts, clbk, thisArg ) { // eslint-disable-line max-statements, max-lines-per-function + var bsize; + var dv0; + var dv1; + var dv2; + var dv3; + var dv4; + var dv5; + var dv6; + var dv7; + var dv8; + var dv9; + var ov1; + var ov2; + var ov3; + var ov4; + var ov5; + var ov6; + var ov7; + var ov8; + var ov9; + var sh; + var s0; + var s1; + var s2; + var s3; + var s4; + var s5; + var s6; + var s7; + var s8; + var s9; + var sv; + var ov; + var iv; + var i0; + var i1; + var i2; + var i3; + var i4; + var i5; + var i6; + var i7; + var i8; + var i9; + var j0; + var j1; + var j2; + var j3; + var j4; + var j5; + var j6; + var j7; + var j8; + var j9; + var N; + var x; + var v; + var o; + var k; + var f; + + // Note on variable naming convention: S#, dv#, i#, j# where # corresponds to the loop number, with `0` being the innermost loop... + + N = arrays.length; + x = arrays[ 0 ]; + + // Resolve the loop interchange order: + o = loopOrder( shape, stridesX ); + sh = o.sh; + sv = [ o.sx ]; + for ( k = 1; k < N; k++ ) { + sv.push( takeIndexed( arrays[k].strides, o.idx ) ); + } + // Determine the block size: + bsize = blockSize( x.dtype ); + + // Resolve a list of pointers to the first indexed elements in the respective ndarrays: + ov = offsets( arrays ); + + // Cache offset increments for the innermost loop... + dv0 = []; + for ( k = 0; k < N; k++ ) { + dv0.push( sv[k][0] ); + } + // Initialize loop variables... + ov1 = zeros( N ); + ov2 = zeros( N ); + ov3 = zeros( N ); + ov4 = zeros( N ); + ov5 = zeros( N ); + ov6 = zeros( N ); + ov7 = zeros( N ); + ov8 = zeros( N ); + ov9 = zeros( N ); + dv1 = zeros( N ); + dv2 = zeros( N ); + dv3 = zeros( N ); + dv4 = zeros( N ); + dv5 = zeros( N ); + dv6 = zeros( N ); + dv7 = zeros( N ); + dv8 = zeros( N ); + dv9 = zeros( N ); + iv = zeros( N ); + + // Shallow copy the list of views to an internal array so that we can update with reshaped views without impacting the original list of views: + v = copyIndexed( views ); + + // Iterate over blocks... + for ( j9 = sh[9]; j9 > 0; ) { + if ( j9 < bsize ) { + s9 = j9; + j9 = 0; + } else { + s9 = bsize; + j9 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov9[ k ] = ov[k] + ( j9*sv[k][9] ); + } + for ( j8 = sh[8]; j8 > 0; ) { + if ( j8 < bsize ) { + s8 = j8; + j8 = 0; + } else { + s8 = bsize; + j8 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov8[ k ] = ov9[k] + ( j8*sv[k][8] ); + dv9[ k ] = sv[k][9] - ( s8*sv[k][8] ); + } + for ( j7 = sh[7]; j7 > 0; ) { + if ( j7 < bsize ) { + s7 = j7; + j7 = 0; + } else { + s7 = bsize; + j7 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov7[ k ] = ov8[k] + ( j7*sv[k][7] ); + dv8[ k ] = sv[k][8] - ( s7*sv[k][7] ); + } + for ( j6 = sh[6]; j6 > 0; ) { + if ( j6 < bsize ) { + s6 = j6; + j6 = 0; + } else { + s6 = bsize; + j6 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov6[ k ] = ov7[k] + ( j6*sv[k][6] ); + dv7[ k ] = sv[k][7] - ( s6*sv[k][6] ); + } + for ( j5 = sh[5]; j5 > 0; ) { + if ( j5 < bsize ) { + s5 = j5; + j5 = 0; + } else { + s5 = bsize; + j5 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov5[ k ] = ov6[k] + ( j5*sv[k][5] ); + dv6[ k ] = sv[k][6] - ( s5*sv[k][5] ); + } + for ( j4 = sh[4]; j4 > 0; ) { + if ( j4 < bsize ) { + s4 = j4; + j4 = 0; + } else { + s4 = bsize; + j4 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov4[ k ] = ov5[k] + ( j4*sv[k][4] ); + dv5[ k ] = sv[k][5] - ( s4*sv[k][4] ); + } + for ( j3 = sh[3]; j3 > 0; ) { + if ( j3 < bsize ) { + s3 = j3; + j3 = 0; + } else { + s3 = bsize; + j3 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov3[ k ] = ov4[k] + ( j3*sv[k][3] ); + dv4[ k ] = sv[k][4] - ( s3*sv[k][3] ); + } + for ( j2 = sh[2]; j2 > 0; ) { + if ( j2 < bsize ) { + s2 = j2; + j2 = 0; + } else { + s2 = bsize; + j2 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov2[ k ] = ov3[k] + ( j2*sv[k][2] ); + dv3[ k ] = sv[k][3] - ( s2*sv[k][2] ); + } + for ( j1 = sh[1]; j1 > 0; ) { + if ( j1 < bsize ) { + s1 = j1; + j1 = 0; + } else { + s1 = bsize; + j1 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov1[ k ] = ov2[k] + ( j1*sv[k][1] ); + dv2[ k ] = sv[k][2] - ( s1*sv[k][1] ); + } + for ( j0 = sh[0]; j0 > 0; ) { + if ( j0 < bsize ) { + s0 = j0; + j0 = 0; + } else { + s0 = bsize; + j0 -= bsize; + } + // Compute index offsets and loop offset increments for the first ndarray elements in the current block... + for ( k = 0; k < N; k++ ) { + iv[ k ] = ov1[k] + ( j0*sv[k][0] ); + dv1[ k ] = sv[k][1] - ( s0*sv[k][0] ); + } + // Iterate over the loop dimensions... + for ( i9 = 0; i9 < s9; i9++ ) { + for ( i8 = 0; i8 < s8; i8++ ) { + for ( i7 = 0; i7 < s7; i7++ ) { + for ( i6 = 0; i6 < s6; i6++ ) { + for ( i5 = 0; i5 < s5; i5++ ) { + for ( i4 = 0; i4 < s4; i4++ ) { + for ( i3 = 0; i3 < s3; i3++ ) { + for ( i2 = 0; i2 < s2; i2++ ) { + for ( i1 = 0; i1 < s1; i1++ ) { + for ( i0 = 0; i0 < s0; i0++ ) { + setViewOffsets( views, iv ); + v[ 0 ] = strategyX.input( views[ 0 ] ); + f = wrap( ref, views[ 0 ], ibuf, ldims, [ j9+i9, j8+i8, j7+i7, j6+i6, j5+i5, j4+i4, j3+i3, j2+i2, j1+i1, j0+i0 ], cdims, clbk, thisArg ); + if ( hasOpts ) { + fcn( v, opts, f ); + } else { + fcn( v, f ); + } + strategyX.output( views[ 0 ] ); + incrementOffsets( iv, dv0 ); + } + incrementOffsets( iv, dv1 ); + } + incrementOffsets( iv, dv2 ); + } + incrementOffsets( iv, dv3 ); + } + incrementOffsets( iv, dv4 ); + } + incrementOffsets( iv, dv5 ); + } + incrementOffsets( iv, dv6 ); + } + incrementOffsets( iv, dv7 ); + } + incrementOffsets( iv, dv8 ); + } + incrementOffsets( iv, dv9 ); + } + } + } + } + } + } + } + } + } + } + } +} + + +// EXPORTS // + +module.exports = kernel; diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/4d.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/4d.js new file mode 100644 index 000000000000..46c82cc6f85d --- /dev/null +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/4d.js @@ -0,0 +1,302 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +/* eslint-disable max-len, max-params, max-depth */ + +'use strict'; + +// MODULES // + +var loopOrder = require( '@stdlib/ndarray/base/nullary-loop-interchange-order' ); +var blockSize = require( '@stdlib/ndarray/base/nullary-tiling-block-size' ); +var takeIndexed = require( '@stdlib/array/base/take-indexed' ); +var copyIndexed = require( '@stdlib/array/base/copy-indexed' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var incrementOffsets = require( '@stdlib/ndarray/base/kernels/utils/increment-offsets' ); +var setViewOffsets = require( '@stdlib/ndarray/base/set-descriptor-offsets' ); +var offsets = require( '@stdlib/ndarray/base/offsets' ); +var wrap = require( './callback_wrapper.js' ); + + +// MAIN // + +/** +* Applies a one-dimensional strided array function to an input ndarray according to a callback function via loop blocking. +* +* @private +* @param {Function} fcn - wrapper for a one-dimensional strided array function +* @param {Array} arrays - ndarray descriptors +* @param {ndarrayLike} ref - original input ndarray provided to a callback function +* @param {Array} views - initialized ndarray descriptors representing subarray views +* @param {Array} shape - loop dimensions +* @param {Array} stridesX - loop dimension strides for the input ndarray +* @param {Object} strategyX - strategy for marshaling data to and from an input ndarray view +* @param {NonNegativeIntegerArray} ibuf - workspace for storing iteration indices +* @param {NonNegativeIntegerArray} ldims - list of loop dimensions +* @param {NonNegativeIntegerArray} cdims - list of "core" dimensions +* @param {Options} opts - function options +* @param {boolean} hasOpts - boolean indicating whether to pass an options argument to a strided array function +* @param {Function} clbk - callback function +* @param {*} thisArg - callback execution context +* @returns {void} +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +* var zeros = require( '@stdlib/array/base/zeros' ); +* var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +* var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); +* +* function clbk( value, idx ) { +* return ( idx[ 3 ]*100 ) + ( idx[ 4 ]*10 ) + idx[ 5 ]; +* } +* +* // Create a data buffer: +* var xbuf = new Float64Array( 12 ); +* +* // Define the array shape: +* var xsh = [ 1, 1, 1, 3, 2, 2 ]; +* +* // Define the array strides: +* var sx = [ 12, 12, 12, 4, 2, 1 ]; +* +* // Define the index offset: +* var ox = 0; +* +* // Create an input ndarray descriptor: +* var x = { +* 'dtype': 'float64', +* 'data': xbuf, +* 'shape': xsh, +* 'strides': sx, +* 'offset': ox, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the starting index (inclusive): +* var start = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 0.0 ] ), +* 'shape': [ 1, 1, 1, 3 ], +* 'strides': [ 0, 0, 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the ending index (exclusive): +* var end = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 4.0 ] ), +* 'shape': [ 1, 1, 1, 3 ], +* 'strides': [ 0, 0, 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Initialize ndarray descriptors representing subarray views: +* var views = [ +* { +* 'dtype': x.dtype, +* 'data': x.data, +* 'shape': [ 2, 2 ], +* 'strides': [ 2, 1 ], +* 'offset': x.offset, +* 'order': x.order +* }, +* { +* 'dtype': start.dtype, +* 'data': start.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': start.offset, +* 'order': start.order +* }, +* { +* 'dtype': end.dtype, +* 'data': end.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': end.offset, +* 'order': end.order +* } +* ]; +* +* // Resolve an input strategy when iterating over subarray views: +* var strategyX = strategy( views[ 0 ] ); +* +* // Create a workspace array for storing iteration indices: +* var ibuf = zeros( xsh.length ); +* +* // Define the loop and core dimensions: +* var ldims = [ 0, 1, 2, 3 ]; +* var cdims = [ 4, 5 ]; +* +* // Apply strided function: +* kernel( gfillBy, [ x, start, end ], x, views, [ 1, 1, 1, 3 ], [ 12, 12, 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); +* +* var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +* // returns [ [ [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] ] ] +*/ +function kernel( fcn, arrays, ref, views, shape, stridesX, strategyX, ibuf, ldims, cdims, opts, hasOpts, clbk, thisArg ) { + var bsize; + var dv0; + var dv1; + var dv2; + var dv3; + var ov1; + var ov2; + var ov3; + var sh; + var s0; + var s1; + var s2; + var s3; + var sv; + var ov; + var iv; + var i0; + var i1; + var i2; + var i3; + var j0; + var j1; + var j2; + var j3; + var N; + var x; + var v; + var o; + var k; + var f; + + // Note on variable naming convention: S#, dv#, i#, j# where # corresponds to the loop number, with `0` being the innermost loop... + + N = arrays.length; + x = arrays[ 0 ]; + + // Resolve the loop interchange order: + o = loopOrder( shape, stridesX ); + sh = o.sh; + sv = [ o.sx ]; + for ( k = 1; k < N; k++ ) { + sv.push( takeIndexed( arrays[k].strides, o.idx ) ); + } + // Determine the block size: + bsize = blockSize( x.dtype ); + + // Resolve a list of pointers to the first indexed elements in the respective ndarrays: + ov = offsets( arrays ); + + // Cache offset increments for the innermost loop... + dv0 = []; + for ( k = 0; k < N; k++ ) { + dv0.push( sv[k][0] ); + } + // Initialize loop variables... + ov1 = zeros( N ); + ov2 = zeros( N ); + ov3 = zeros( N ); + dv1 = zeros( N ); + dv2 = zeros( N ); + dv3 = zeros( N ); + iv = zeros( N ); + + // Shallow copy the list of views to an internal array so that we can update with reshaped views without impacting the original list of views: + v = copyIndexed( views ); + + // Iterate over blocks... + for ( j3 = sh[3]; j3 > 0; ) { + if ( j3 < bsize ) { + s3 = j3; + j3 = 0; + } else { + s3 = bsize; + j3 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov3[ k ] = ov[k] + ( j3*sv[k][3] ); + } + for ( j2 = sh[2]; j2 > 0; ) { + if ( j2 < bsize ) { + s2 = j2; + j2 = 0; + } else { + s2 = bsize; + j2 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov2[ k ] = ov3[k] + ( j2*sv[k][2] ); + dv3[ k ] = sv[k][3] - ( s2*sv[k][2] ); + } + for ( j1 = sh[1]; j1 > 0; ) { + if ( j1 < bsize ) { + s1 = j1; + j1 = 0; + } else { + s1 = bsize; + j1 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov1[ k ] = ov2[k] + ( j1*sv[k][1] ); + dv2[ k ] = sv[k][2] - ( s1*sv[k][1] ); + } + for ( j0 = sh[0]; j0 > 0; ) { + if ( j0 < bsize ) { + s0 = j0; + j0 = 0; + } else { + s0 = bsize; + j0 -= bsize; + } + // Compute index offsets and loop offset increments for the first ndarray elements in the current block... + for ( k = 0; k < N; k++ ) { + iv[ k ] = ov1[k] + ( j0*sv[k][0] ); + dv1[ k ] = sv[k][1] - ( s0*sv[k][0] ); + } + // Iterate over the loop dimensions... + for ( i3 = 0; i3 < s3; i3++ ) { + for ( i2 = 0; i2 < s2; i2++ ) { + for ( i1 = 0; i1 < s1; i1++ ) { + for ( i0 = 0; i0 < s0; i0++ ) { + setViewOffsets( views, iv ); + v[ 0 ] = strategyX.input( views[ 0 ] ); + f = wrap( ref, views[ 0 ], ibuf, ldims, [ j3+i3, j2+i2, j1+i1, j0+i0 ], cdims, clbk, thisArg ); + if ( hasOpts ) { + fcn( v, opts, f ); + } else { + fcn( v, f ); + } + strategyX.output( views[ 0 ] ); + incrementOffsets( iv, dv0 ); + } + incrementOffsets( iv, dv1 ); + } + incrementOffsets( iv, dv2 ); + } + incrementOffsets( iv, dv3 ); + } + } + } + } + } +} + + +// EXPORTS // + +module.exports = kernel; diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/5d.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/5d.js new file mode 100644 index 000000000000..c50fe0c6801a --- /dev/null +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/5d.js @@ -0,0 +1,325 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +/* eslint-disable max-len, max-params, max-depth */ + +'use strict'; + +// MODULES // + +var loopOrder = require( '@stdlib/ndarray/base/nullary-loop-interchange-order' ); +var blockSize = require( '@stdlib/ndarray/base/nullary-tiling-block-size' ); +var takeIndexed = require( '@stdlib/array/base/take-indexed' ); +var copyIndexed = require( '@stdlib/array/base/copy-indexed' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var incrementOffsets = require( '@stdlib/ndarray/base/kernels/utils/increment-offsets' ); +var setViewOffsets = require( '@stdlib/ndarray/base/set-descriptor-offsets' ); +var offsets = require( '@stdlib/ndarray/base/offsets' ); +var wrap = require( './callback_wrapper.js' ); + + +// MAIN // + +/** +* Applies a one-dimensional strided array function to an input ndarray according to a callback function via loop blocking. +* +* @private +* @param {Function} fcn - wrapper for a one-dimensional strided array function +* @param {Array} arrays - ndarray descriptors +* @param {ndarrayLike} ref - original input ndarray provided to a callback function +* @param {Array} views - initialized ndarray descriptors representing subarray views +* @param {Array} shape - loop dimensions +* @param {Array} stridesX - loop dimension strides for the input ndarray +* @param {Object} strategyX - strategy for marshaling data to and from an input ndarray view +* @param {NonNegativeIntegerArray} ibuf - workspace for storing iteration indices +* @param {NonNegativeIntegerArray} ldims - list of loop dimensions +* @param {NonNegativeIntegerArray} cdims - list of "core" dimensions +* @param {Options} opts - function options +* @param {boolean} hasOpts - boolean indicating whether to pass an options argument to a strided array function +* @param {Function} clbk - callback function +* @param {*} thisArg - callback execution context +* @returns {void} +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +* var zeros = require( '@stdlib/array/base/zeros' ); +* var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +* var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); +* +* function clbk( value, idx ) { +* return ( idx[ 4 ]*100 ) + ( idx[ 5 ]*10 ) + idx[ 6 ]; +* } +* +* // Create a data buffer: +* var xbuf = new Float64Array( 12 ); +* +* // Define the array shape: +* var xsh = [ 1, 1, 1, 1, 3, 2, 2 ]; +* +* // Define the array strides: +* var sx = [ 12, 12, 12, 12, 4, 2, 1 ]; +* +* // Define the index offset: +* var ox = 0; +* +* // Create an input ndarray descriptor: +* var x = { +* 'dtype': 'float64', +* 'data': xbuf, +* 'shape': xsh, +* 'strides': sx, +* 'offset': ox, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the starting index (inclusive): +* var start = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 0.0 ] ), +* 'shape': [ 1, 1, 1, 1, 3 ], +* 'strides': [ 0, 0, 0, 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the ending index (exclusive): +* var end = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 4.0 ] ), +* 'shape': [ 1, 1, 1, 1, 3 ], +* 'strides': [ 0, 0, 0, 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Initialize ndarray descriptors representing subarray views: +* var views = [ +* { +* 'dtype': x.dtype, +* 'data': x.data, +* 'shape': [ 2, 2 ], +* 'strides': [ 2, 1 ], +* 'offset': x.offset, +* 'order': x.order +* }, +* { +* 'dtype': start.dtype, +* 'data': start.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': start.offset, +* 'order': start.order +* }, +* { +* 'dtype': end.dtype, +* 'data': end.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': end.offset, +* 'order': end.order +* } +* ]; +* +* // Resolve an input strategy when iterating over subarray views: +* var strategyX = strategy( views[ 0 ] ); +* +* // Create a workspace array for storing iteration indices: +* var ibuf = zeros( xsh.length ); +* +* // Define the loop and core dimensions: +* var ldims = [ 0, 1, 2, 3, 4 ]; +* var cdims = [ 5, 6 ]; +* +* // Apply strided function: +* kernel( gfillBy, [ x, start, end ], x, views, [ 1, 1, 1, 1, 3 ], [ 12, 12, 12, 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); +* +* var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +* // returns [ [ [ [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] ] ] ] +*/ +function kernel( fcn, arrays, ref, views, shape, stridesX, strategyX, ibuf, ldims, cdims, opts, hasOpts, clbk, thisArg ) { + var bsize; + var dv0; + var dv1; + var dv2; + var dv3; + var dv4; + var ov1; + var ov2; + var ov3; + var ov4; + var sh; + var s0; + var s1; + var s2; + var s3; + var s4; + var sv; + var ov; + var iv; + var i0; + var i1; + var i2; + var i3; + var i4; + var j0; + var j1; + var j2; + var j3; + var j4; + var N; + var x; + var v; + var o; + var k; + var f; + + // Note on variable naming convention: S#, dv#, i#, j# where # corresponds to the loop number, with `0` being the innermost loop... + + N = arrays.length; + x = arrays[ 0 ]; + + // Resolve the loop interchange order: + o = loopOrder( shape, stridesX ); + sh = o.sh; + sv = [ o.sx ]; + for ( k = 1; k < N; k++ ) { + sv.push( takeIndexed( arrays[k].strides, o.idx ) ); + } + // Determine the block size: + bsize = blockSize( x.dtype ); + + // Resolve a list of pointers to the first indexed elements in the respective ndarrays: + ov = offsets( arrays ); + + // Cache offset increments for the innermost loop... + dv0 = []; + for ( k = 0; k < N; k++ ) { + dv0.push( sv[k][0] ); + } + // Initialize loop variables... + ov1 = zeros( N ); + ov2 = zeros( N ); + ov3 = zeros( N ); + ov4 = zeros( N ); + dv1 = zeros( N ); + dv2 = zeros( N ); + dv3 = zeros( N ); + dv4 = zeros( N ); + iv = zeros( N ); + + // Shallow copy the list of views to an internal array so that we can update with reshaped views without impacting the original list of views: + v = copyIndexed( views ); + + // Iterate over blocks... + for ( j4 = sh[4]; j4 > 0; ) { + if ( j4 < bsize ) { + s4 = j4; + j4 = 0; + } else { + s4 = bsize; + j4 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov4[ k ] = ov[k] + ( j4*sv[k][4] ); + } + for ( j3 = sh[3]; j3 > 0; ) { + if ( j3 < bsize ) { + s3 = j3; + j3 = 0; + } else { + s3 = bsize; + j3 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov3[ k ] = ov4[k] + ( j3*sv[k][3] ); + dv4[ k ] = sv[k][4] - ( s3*sv[k][3] ); + } + for ( j2 = sh[2]; j2 > 0; ) { + if ( j2 < bsize ) { + s2 = j2; + j2 = 0; + } else { + s2 = bsize; + j2 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov2[ k ] = ov3[k] + ( j2*sv[k][2] ); + dv3[ k ] = sv[k][3] - ( s2*sv[k][2] ); + } + for ( j1 = sh[1]; j1 > 0; ) { + if ( j1 < bsize ) { + s1 = j1; + j1 = 0; + } else { + s1 = bsize; + j1 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov1[ k ] = ov2[k] + ( j1*sv[k][1] ); + dv2[ k ] = sv[k][2] - ( s1*sv[k][1] ); + } + for ( j0 = sh[0]; j0 > 0; ) { + if ( j0 < bsize ) { + s0 = j0; + j0 = 0; + } else { + s0 = bsize; + j0 -= bsize; + } + // Compute index offsets and loop offset increments for the first ndarray elements in the current block... + for ( k = 0; k < N; k++ ) { + iv[ k ] = ov1[k] + ( j0*sv[k][0] ); + dv1[ k ] = sv[k][1] - ( s0*sv[k][0] ); + } + // Iterate over the loop dimensions... + for ( i4 = 0; i4 < s4; i4++ ) { + for ( i3 = 0; i3 < s3; i3++ ) { + for ( i2 = 0; i2 < s2; i2++ ) { + for ( i1 = 0; i1 < s1; i1++ ) { + for ( i0 = 0; i0 < s0; i0++ ) { + setViewOffsets( views, iv ); + v[ 0 ] = strategyX.input( views[ 0 ] ); + f = wrap( ref, views[ 0 ], ibuf, ldims, [ j4+i4, j3+i3, j2+i2, j1+i1, j0+i0 ], cdims, clbk, thisArg ); + if ( hasOpts ) { + fcn( v, opts, f ); + } else { + fcn( v, f ); + } + strategyX.output( views[ 0 ] ); + incrementOffsets( iv, dv0 ); + } + incrementOffsets( iv, dv1 ); + } + incrementOffsets( iv, dv2 ); + } + incrementOffsets( iv, dv3 ); + } + incrementOffsets( iv, dv4 ); + } + } + } + } + } + } +} + + +// EXPORTS // + +module.exports = kernel; diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/6d.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/6d.js new file mode 100644 index 000000000000..9bb8e129202c --- /dev/null +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/6d.js @@ -0,0 +1,348 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +/* eslint-disable max-len, max-params, max-depth */ + +'use strict'; + +// MODULES // + +var loopOrder = require( '@stdlib/ndarray/base/nullary-loop-interchange-order' ); +var blockSize = require( '@stdlib/ndarray/base/nullary-tiling-block-size' ); +var takeIndexed = require( '@stdlib/array/base/take-indexed' ); +var copyIndexed = require( '@stdlib/array/base/copy-indexed' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var incrementOffsets = require( '@stdlib/ndarray/base/kernels/utils/increment-offsets' ); +var setViewOffsets = require( '@stdlib/ndarray/base/set-descriptor-offsets' ); +var offsets = require( '@stdlib/ndarray/base/offsets' ); +var wrap = require( './callback_wrapper.js' ); + + +// MAIN // + +/** +* Applies a one-dimensional strided array function to an input ndarray according to a callback function via loop blocking. +* +* @private +* @param {Function} fcn - wrapper for a one-dimensional strided array function +* @param {Array} arrays - ndarray descriptors +* @param {ndarrayLike} ref - original input ndarray provided to a callback function +* @param {Array} views - initialized ndarray descriptors representing subarray views +* @param {Array} shape - loop dimensions +* @param {Array} stridesX - loop dimension strides for the input ndarray +* @param {Object} strategyX - strategy for marshaling data to and from an input ndarray view +* @param {NonNegativeIntegerArray} ibuf - workspace for storing iteration indices +* @param {NonNegativeIntegerArray} ldims - list of loop dimensions +* @param {NonNegativeIntegerArray} cdims - list of "core" dimensions +* @param {Options} opts - function options +* @param {boolean} hasOpts - boolean indicating whether to pass an options argument to a strided array function +* @param {Function} clbk - callback function +* @param {*} thisArg - callback execution context +* @returns {void} +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +* var zeros = require( '@stdlib/array/base/zeros' ); +* var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +* var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); +* +* function clbk( value, idx ) { +* return ( idx[ 5 ]*100 ) + ( idx[ 6 ]*10 ) + idx[ 7 ]; +* } +* +* // Create a data buffer: +* var xbuf = new Float64Array( 12 ); +* +* // Define the array shape: +* var xsh = [ 1, 1, 1, 1, 1, 3, 2, 2 ]; +* +* // Define the array strides: +* var sx = [ 12, 12, 12, 12, 12, 4, 2, 1 ]; +* +* // Define the index offset: +* var ox = 0; +* +* // Create an input ndarray descriptor: +* var x = { +* 'dtype': 'float64', +* 'data': xbuf, +* 'shape': xsh, +* 'strides': sx, +* 'offset': ox, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the starting index (inclusive): +* var start = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 0.0 ] ), +* 'shape': [ 1, 1, 1, 1, 1, 3 ], +* 'strides': [ 0, 0, 0, 0, 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the ending index (exclusive): +* var end = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 4.0 ] ), +* 'shape': [ 1, 1, 1, 1, 1, 3 ], +* 'strides': [ 0, 0, 0, 0, 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Initialize ndarray descriptors representing subarray views: +* var views = [ +* { +* 'dtype': x.dtype, +* 'data': x.data, +* 'shape': [ 2, 2 ], +* 'strides': [ 2, 1 ], +* 'offset': x.offset, +* 'order': x.order +* }, +* { +* 'dtype': start.dtype, +* 'data': start.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': start.offset, +* 'order': start.order +* }, +* { +* 'dtype': end.dtype, +* 'data': end.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': end.offset, +* 'order': end.order +* } +* ]; +* +* // Resolve an input strategy when iterating over subarray views: +* var strategyX = strategy( views[ 0 ] ); +* +* // Create a workspace array for storing iteration indices: +* var ibuf = zeros( xsh.length ); +* +* // Define the loop and core dimensions: +* var ldims = [ 0, 1, 2, 3, 4, 5 ]; +* var cdims = [ 6, 7 ]; +* +* // Apply strided function: +* kernel( gfillBy, [ x, start, end ], x, views, [ 1, 1, 1, 1, 1, 3 ], [ 12, 12, 12, 12, 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); +* +* var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +* // returns [ [ [ [ [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] ] ] ] ] +*/ +function kernel( fcn, arrays, ref, views, shape, stridesX, strategyX, ibuf, ldims, cdims, opts, hasOpts, clbk, thisArg ) { + var bsize; + var dv0; + var dv1; + var dv2; + var dv3; + var dv4; + var dv5; + var ov1; + var ov2; + var ov3; + var ov4; + var ov5; + var sh; + var s0; + var s1; + var s2; + var s3; + var s4; + var s5; + var sv; + var ov; + var iv; + var i0; + var i1; + var i2; + var i3; + var i4; + var i5; + var j0; + var j1; + var j2; + var j3; + var j4; + var j5; + var N; + var x; + var v; + var o; + var k; + var f; + + // Note on variable naming convention: S#, dv#, i#, j# where # corresponds to the loop number, with `0` being the innermost loop... + + N = arrays.length; + x = arrays[ 0 ]; + + // Resolve the loop interchange order: + o = loopOrder( shape, stridesX ); + sh = o.sh; + sv = [ o.sx ]; + for ( k = 1; k < N; k++ ) { + sv.push( takeIndexed( arrays[k].strides, o.idx ) ); + } + // Determine the block size: + bsize = blockSize( x.dtype ); + + // Resolve a list of pointers to the first indexed elements in the respective ndarrays: + ov = offsets( arrays ); + + // Cache offset increments for the innermost loop... + dv0 = []; + for ( k = 0; k < N; k++ ) { + dv0.push( sv[k][0] ); + } + // Initialize loop variables... + ov1 = zeros( N ); + ov2 = zeros( N ); + ov3 = zeros( N ); + ov4 = zeros( N ); + ov5 = zeros( N ); + dv1 = zeros( N ); + dv2 = zeros( N ); + dv3 = zeros( N ); + dv4 = zeros( N ); + dv5 = zeros( N ); + iv = zeros( N ); + + // Shallow copy the list of views to an internal array so that we can update with reshaped views without impacting the original list of views: + v = copyIndexed( views ); + + // Iterate over blocks... + for ( j5 = sh[5]; j5 > 0; ) { + if ( j5 < bsize ) { + s5 = j5; + j5 = 0; + } else { + s5 = bsize; + j5 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov5[ k ] = ov[k] + ( j5*sv[k][5] ); + } + for ( j4 = sh[4]; j4 > 0; ) { + if ( j4 < bsize ) { + s4 = j4; + j4 = 0; + } else { + s4 = bsize; + j4 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov4[ k ] = ov5[k] + ( j4*sv[k][4] ); + dv5[ k ] = sv[k][5] - ( s4*sv[k][4] ); + } + for ( j3 = sh[3]; j3 > 0; ) { + if ( j3 < bsize ) { + s3 = j3; + j3 = 0; + } else { + s3 = bsize; + j3 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov3[ k ] = ov4[k] + ( j3*sv[k][3] ); + dv4[ k ] = sv[k][4] - ( s3*sv[k][3] ); + } + for ( j2 = sh[2]; j2 > 0; ) { + if ( j2 < bsize ) { + s2 = j2; + j2 = 0; + } else { + s2 = bsize; + j2 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov2[ k ] = ov3[k] + ( j2*sv[k][2] ); + dv3[ k ] = sv[k][3] - ( s2*sv[k][2] ); + } + for ( j1 = sh[1]; j1 > 0; ) { + if ( j1 < bsize ) { + s1 = j1; + j1 = 0; + } else { + s1 = bsize; + j1 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov1[ k ] = ov2[k] + ( j1*sv[k][1] ); + dv2[ k ] = sv[k][2] - ( s1*sv[k][1] ); + } + for ( j0 = sh[0]; j0 > 0; ) { + if ( j0 < bsize ) { + s0 = j0; + j0 = 0; + } else { + s0 = bsize; + j0 -= bsize; + } + // Compute index offsets and loop offset increments for the first ndarray elements in the current block... + for ( k = 0; k < N; k++ ) { + iv[ k ] = ov1[k] + ( j0*sv[k][0] ); + dv1[ k ] = sv[k][1] - ( s0*sv[k][0] ); + } + // Iterate over the loop dimensions... + for ( i5 = 0; i5 < s5; i5++ ) { + for ( i4 = 0; i4 < s4; i4++ ) { + for ( i3 = 0; i3 < s3; i3++ ) { + for ( i2 = 0; i2 < s2; i2++ ) { + for ( i1 = 0; i1 < s1; i1++ ) { + for ( i0 = 0; i0 < s0; i0++ ) { + setViewOffsets( views, iv ); + v[ 0 ] = strategyX.input( views[ 0 ] ); + f = wrap( ref, views[ 0 ], ibuf, ldims, [ j5+i5, j4+i4, j3+i3, j2+i2, j1+i1, j0+i0 ], cdims, clbk, thisArg ); + if ( hasOpts ) { + fcn( v, opts, f ); + } else { + fcn( v, f ); + } + strategyX.output( views[ 0 ] ); + incrementOffsets( iv, dv0 ); + } + incrementOffsets( iv, dv1 ); + } + incrementOffsets( iv, dv2 ); + } + incrementOffsets( iv, dv3 ); + } + incrementOffsets( iv, dv4 ); + } + incrementOffsets( iv, dv5 ); + } + } + } + } + } + } + } +} + + +// EXPORTS // + +module.exports = kernel; diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/7d.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/7d.js new file mode 100644 index 000000000000..797fa1e122f7 --- /dev/null +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/7d.js @@ -0,0 +1,371 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +/* eslint-disable max-len, max-params, max-depth */ + +'use strict'; + +// MODULES // + +var loopOrder = require( '@stdlib/ndarray/base/nullary-loop-interchange-order' ); +var blockSize = require( '@stdlib/ndarray/base/nullary-tiling-block-size' ); +var takeIndexed = require( '@stdlib/array/base/take-indexed' ); +var copyIndexed = require( '@stdlib/array/base/copy-indexed' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var incrementOffsets = require( '@stdlib/ndarray/base/kernels/utils/increment-offsets' ); +var setViewOffsets = require( '@stdlib/ndarray/base/set-descriptor-offsets' ); +var offsets = require( '@stdlib/ndarray/base/offsets' ); +var wrap = require( './callback_wrapper.js' ); + + +// MAIN // + +/** +* Applies a one-dimensional strided array function to an input ndarray according to a callback function via loop blocking. +* +* @private +* @param {Function} fcn - wrapper for a one-dimensional strided array function +* @param {Array} arrays - ndarray descriptors +* @param {ndarrayLike} ref - original input ndarray provided to a callback function +* @param {Array} views - initialized ndarray descriptors representing subarray views +* @param {Array} shape - loop dimensions +* @param {Array} stridesX - loop dimension strides for the input ndarray +* @param {Object} strategyX - strategy for marshaling data to and from an input ndarray view +* @param {NonNegativeIntegerArray} ibuf - workspace for storing iteration indices +* @param {NonNegativeIntegerArray} ldims - list of loop dimensions +* @param {NonNegativeIntegerArray} cdims - list of "core" dimensions +* @param {Options} opts - function options +* @param {boolean} hasOpts - boolean indicating whether to pass an options argument to a strided array function +* @param {Function} clbk - callback function +* @param {*} thisArg - callback execution context +* @returns {void} +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +* var zeros = require( '@stdlib/array/base/zeros' ); +* var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +* var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); +* +* function clbk( value, idx ) { +* return ( idx[ 6 ]*100 ) + ( idx[ 7 ]*10 ) + idx[ 8 ]; +* } +* +* // Create a data buffer: +* var xbuf = new Float64Array( 12 ); +* +* // Define the array shape: +* var xsh = [ 1, 1, 1, 1, 1, 1, 3, 2, 2 ]; +* +* // Define the array strides: +* var sx = [ 12, 12, 12, 12, 12, 12, 4, 2, 1 ]; +* +* // Define the index offset: +* var ox = 0; +* +* // Create an input ndarray descriptor: +* var x = { +* 'dtype': 'float64', +* 'data': xbuf, +* 'shape': xsh, +* 'strides': sx, +* 'offset': ox, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the starting index (inclusive): +* var start = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 0.0 ] ), +* 'shape': [ 1, 1, 1, 1, 1, 1, 3 ], +* 'strides': [ 0, 0, 0, 0, 0, 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the ending index (exclusive): +* var end = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 4.0 ] ), +* 'shape': [ 1, 1, 1, 1, 1, 1, 3 ], +* 'strides': [ 0, 0, 0, 0, 0, 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Initialize ndarray descriptors representing subarray views: +* var views = [ +* { +* 'dtype': x.dtype, +* 'data': x.data, +* 'shape': [ 2, 2 ], +* 'strides': [ 2, 1 ], +* 'offset': x.offset, +* 'order': x.order +* }, +* { +* 'dtype': start.dtype, +* 'data': start.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': start.offset, +* 'order': start.order +* }, +* { +* 'dtype': end.dtype, +* 'data': end.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': end.offset, +* 'order': end.order +* } +* ]; +* +* // Resolve an input strategy when iterating over subarray views: +* var strategyX = strategy( views[ 0 ] ); +* +* // Create a workspace array for storing iteration indices: +* var ibuf = zeros( xsh.length ); +* +* // Define the loop and core dimensions: +* var ldims = [ 0, 1, 2, 3, 4, 5, 6 ]; +* var cdims = [ 7, 8 ]; +* +* // Apply strided function: +* kernel( gfillBy, [ x, start, end ], x, views, [ 1, 1, 1, 1, 1, 1, 3 ], [ 12, 12, 12, 12, 12, 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); +* +* var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +* // returns [ [ [ [ [ [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] ] ] ] ] ] +*/ +function kernel( fcn, arrays, ref, views, shape, stridesX, strategyX, ibuf, ldims, cdims, opts, hasOpts, clbk, thisArg ) { + var bsize; + var dv0; + var dv1; + var dv2; + var dv3; + var dv4; + var dv5; + var dv6; + var ov1; + var ov2; + var ov3; + var ov4; + var ov5; + var ov6; + var sh; + var s0; + var s1; + var s2; + var s3; + var s4; + var s5; + var s6; + var sv; + var ov; + var iv; + var i0; + var i1; + var i2; + var i3; + var i4; + var i5; + var i6; + var j0; + var j1; + var j2; + var j3; + var j4; + var j5; + var j6; + var N; + var x; + var v; + var o; + var k; + var f; + + // Note on variable naming convention: S#, dv#, i#, j# where # corresponds to the loop number, with `0` being the innermost loop... + + N = arrays.length; + x = arrays[ 0 ]; + + // Resolve the loop interchange order: + o = loopOrder( shape, stridesX ); + sh = o.sh; + sv = [ o.sx ]; + for ( k = 1; k < N; k++ ) { + sv.push( takeIndexed( arrays[k].strides, o.idx ) ); + } + // Determine the block size: + bsize = blockSize( x.dtype ); + + // Resolve a list of pointers to the first indexed elements in the respective ndarrays: + ov = offsets( arrays ); + + // Cache offset increments for the innermost loop... + dv0 = []; + for ( k = 0; k < N; k++ ) { + dv0.push( sv[k][0] ); + } + // Initialize loop variables... + ov1 = zeros( N ); + ov2 = zeros( N ); + ov3 = zeros( N ); + ov4 = zeros( N ); + ov5 = zeros( N ); + ov6 = zeros( N ); + dv1 = zeros( N ); + dv2 = zeros( N ); + dv3 = zeros( N ); + dv4 = zeros( N ); + dv5 = zeros( N ); + dv6 = zeros( N ); + iv = zeros( N ); + + // Shallow copy the list of views to an internal array so that we can update with reshaped views without impacting the original list of views: + v = copyIndexed( views ); + + // Iterate over blocks... + for ( j6 = sh[6]; j6 > 0; ) { + if ( j6 < bsize ) { + s6 = j6; + j6 = 0; + } else { + s6 = bsize; + j6 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov6[ k ] = ov[k] + ( j6*sv[k][6] ); + } + for ( j5 = sh[5]; j5 > 0; ) { + if ( j5 < bsize ) { + s5 = j5; + j5 = 0; + } else { + s5 = bsize; + j5 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov5[ k ] = ov6[k] + ( j5*sv[k][5] ); + dv6[ k ] = sv[k][6] - ( s5*sv[k][5] ); + } + for ( j4 = sh[4]; j4 > 0; ) { + if ( j4 < bsize ) { + s4 = j4; + j4 = 0; + } else { + s4 = bsize; + j4 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov4[ k ] = ov5[k] + ( j4*sv[k][4] ); + dv5[ k ] = sv[k][5] - ( s4*sv[k][4] ); + } + for ( j3 = sh[3]; j3 > 0; ) { + if ( j3 < bsize ) { + s3 = j3; + j3 = 0; + } else { + s3 = bsize; + j3 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov3[ k ] = ov4[k] + ( j3*sv[k][3] ); + dv4[ k ] = sv[k][4] - ( s3*sv[k][3] ); + } + for ( j2 = sh[2]; j2 > 0; ) { + if ( j2 < bsize ) { + s2 = j2; + j2 = 0; + } else { + s2 = bsize; + j2 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov2[ k ] = ov3[k] + ( j2*sv[k][2] ); + dv3[ k ] = sv[k][3] - ( s2*sv[k][2] ); + } + for ( j1 = sh[1]; j1 > 0; ) { + if ( j1 < bsize ) { + s1 = j1; + j1 = 0; + } else { + s1 = bsize; + j1 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov1[ k ] = ov2[k] + ( j1*sv[k][1] ); + dv2[ k ] = sv[k][2] - ( s1*sv[k][1] ); + } + for ( j0 = sh[0]; j0 > 0; ) { + if ( j0 < bsize ) { + s0 = j0; + j0 = 0; + } else { + s0 = bsize; + j0 -= bsize; + } + // Compute index offsets and loop offset increments for the first ndarray elements in the current block... + for ( k = 0; k < N; k++ ) { + iv[ k ] = ov1[k] + ( j0*sv[k][0] ); + dv1[ k ] = sv[k][1] - ( s0*sv[k][0] ); + } + // Iterate over the loop dimensions... + for ( i6 = 0; i6 < s6; i6++ ) { + for ( i5 = 0; i5 < s5; i5++ ) { + for ( i4 = 0; i4 < s4; i4++ ) { + for ( i3 = 0; i3 < s3; i3++ ) { + for ( i2 = 0; i2 < s2; i2++ ) { + for ( i1 = 0; i1 < s1; i1++ ) { + for ( i0 = 0; i0 < s0; i0++ ) { + setViewOffsets( views, iv ); + v[ 0 ] = strategyX.input( views[ 0 ] ); + f = wrap( ref, views[ 0 ], ibuf, ldims, [ j6+i6, j5+i5, j4+i4, j3+i3, j2+i2, j1+i1, j0+i0 ], cdims, clbk, thisArg ); + if ( hasOpts ) { + fcn( v, opts, f ); + } else { + fcn( v, f ); + } + strategyX.output( views[ 0 ] ); + incrementOffsets( iv, dv0 ); + } + incrementOffsets( iv, dv1 ); + } + incrementOffsets( iv, dv2 ); + } + incrementOffsets( iv, dv3 ); + } + incrementOffsets( iv, dv4 ); + } + incrementOffsets( iv, dv5 ); + } + incrementOffsets( iv, dv6 ); + } + } + } + } + } + } + } + } +} + + +// EXPORTS // + +module.exports = kernel; diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/8d.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/8d.js new file mode 100644 index 000000000000..a8765c83dbe9 --- /dev/null +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/8d.js @@ -0,0 +1,394 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +/* eslint-disable max-len, max-params, max-depth */ + +'use strict'; + +// MODULES // + +var loopOrder = require( '@stdlib/ndarray/base/nullary-loop-interchange-order' ); +var blockSize = require( '@stdlib/ndarray/base/nullary-tiling-block-size' ); +var takeIndexed = require( '@stdlib/array/base/take-indexed' ); +var copyIndexed = require( '@stdlib/array/base/copy-indexed' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var incrementOffsets = require( '@stdlib/ndarray/base/kernels/utils/increment-offsets' ); +var setViewOffsets = require( '@stdlib/ndarray/base/set-descriptor-offsets' ); +var offsets = require( '@stdlib/ndarray/base/offsets' ); +var wrap = require( './callback_wrapper.js' ); + + +// MAIN // + +/** +* Applies a one-dimensional strided array function to an input ndarray according to a callback function via loop blocking. +* +* @private +* @param {Function} fcn - wrapper for a one-dimensional strided array function +* @param {Array} arrays - ndarray descriptors +* @param {ndarrayLike} ref - original input ndarray provided to a callback function +* @param {Array} views - initialized ndarray descriptors representing subarray views +* @param {Array} shape - loop dimensions +* @param {Array} stridesX - loop dimension strides for the input ndarray +* @param {Object} strategyX - strategy for marshaling data to and from an input ndarray view +* @param {NonNegativeIntegerArray} ibuf - workspace for storing iteration indices +* @param {NonNegativeIntegerArray} ldims - list of loop dimensions +* @param {NonNegativeIntegerArray} cdims - list of "core" dimensions +* @param {Options} opts - function options +* @param {boolean} hasOpts - boolean indicating whether to pass an options argument to a strided array function +* @param {Function} clbk - callback function +* @param {*} thisArg - callback execution context +* @returns {void} +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +* var zeros = require( '@stdlib/array/base/zeros' ); +* var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +* var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); +* +* function clbk( value, idx ) { +* return ( idx[ 7 ]*100 ) + ( idx[ 8 ]*10 ) + idx[ 9 ]; +* } +* +* // Create a data buffer: +* var xbuf = new Float64Array( 12 ); +* +* // Define the array shape: +* var xsh = [ 1, 1, 1, 1, 1, 1, 1, 3, 2, 2 ]; +* +* // Define the array strides: +* var sx = [ 12, 12, 12, 12, 12, 12, 12, 4, 2, 1 ]; +* +* // Define the index offset: +* var ox = 0; +* +* // Create an input ndarray descriptor: +* var x = { +* 'dtype': 'float64', +* 'data': xbuf, +* 'shape': xsh, +* 'strides': sx, +* 'offset': ox, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the starting index (inclusive): +* var start = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 0.0 ] ), +* 'shape': [ 1, 1, 1, 1, 1, 1, 1, 3 ], +* 'strides': [ 0, 0, 0, 0, 0, 0, 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the ending index (exclusive): +* var end = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 4.0 ] ), +* 'shape': [ 1, 1, 1, 1, 1, 1, 1, 3 ], +* 'strides': [ 0, 0, 0, 0, 0, 0, 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Initialize ndarray descriptors representing subarray views: +* var views = [ +* { +* 'dtype': x.dtype, +* 'data': x.data, +* 'shape': [ 2, 2 ], +* 'strides': [ 2, 1 ], +* 'offset': x.offset, +* 'order': x.order +* }, +* { +* 'dtype': start.dtype, +* 'data': start.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': start.offset, +* 'order': start.order +* }, +* { +* 'dtype': end.dtype, +* 'data': end.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': end.offset, +* 'order': end.order +* } +* ]; +* +* // Resolve an input strategy when iterating over subarray views: +* var strategyX = strategy( views[ 0 ] ); +* +* // Create a workspace array for storing iteration indices: +* var ibuf = zeros( xsh.length ); +* +* // Define the loop and core dimensions: +* var ldims = [ 0, 1, 2, 3, 4, 5, 6, 7 ]; +* var cdims = [ 8, 9 ]; +* +* // Apply strided function: +* kernel( gfillBy, [ x, start, end ], x, views, [ 1, 1, 1, 1, 1, 1, 1, 3 ], [ 12, 12, 12, 12, 12, 12, 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); +* +* var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +* // returns [ [ [ [ [ [ [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] ] ] ] ] ] ] +*/ +function kernel( fcn, arrays, ref, views, shape, stridesX, strategyX, ibuf, ldims, cdims, opts, hasOpts, clbk, thisArg ) { // eslint-disable-line max-statements, max-lines-per-function + var bsize; + var dv0; + var dv1; + var dv2; + var dv3; + var dv4; + var dv5; + var dv6; + var dv7; + var ov1; + var ov2; + var ov3; + var ov4; + var ov5; + var ov6; + var ov7; + var sh; + var s0; + var s1; + var s2; + var s3; + var s4; + var s5; + var s6; + var s7; + var sv; + var ov; + var iv; + var i0; + var i1; + var i2; + var i3; + var i4; + var i5; + var i6; + var i7; + var j0; + var j1; + var j2; + var j3; + var j4; + var j5; + var j6; + var j7; + var N; + var x; + var v; + var o; + var k; + var f; + + // Note on variable naming convention: S#, dv#, i#, j# where # corresponds to the loop number, with `0` being the innermost loop... + + N = arrays.length; + x = arrays[ 0 ]; + + // Resolve the loop interchange order: + o = loopOrder( shape, stridesX ); + sh = o.sh; + sv = [ o.sx ]; + for ( k = 1; k < N; k++ ) { + sv.push( takeIndexed( arrays[k].strides, o.idx ) ); + } + // Determine the block size: + bsize = blockSize( x.dtype ); + + // Resolve a list of pointers to the first indexed elements in the respective ndarrays: + ov = offsets( arrays ); + + // Cache offset increments for the innermost loop... + dv0 = []; + for ( k = 0; k < N; k++ ) { + dv0.push( sv[k][0] ); + } + // Initialize loop variables... + ov1 = zeros( N ); + ov2 = zeros( N ); + ov3 = zeros( N ); + ov4 = zeros( N ); + ov5 = zeros( N ); + ov6 = zeros( N ); + ov7 = zeros( N ); + dv1 = zeros( N ); + dv2 = zeros( N ); + dv3 = zeros( N ); + dv4 = zeros( N ); + dv5 = zeros( N ); + dv6 = zeros( N ); + dv7 = zeros( N ); + iv = zeros( N ); + + // Shallow copy the list of views to an internal array so that we can update with reshaped views without impacting the original list of views: + v = copyIndexed( views ); + + // Iterate over blocks... + for ( j7 = sh[7]; j7 > 0; ) { + if ( j7 < bsize ) { + s7 = j7; + j7 = 0; + } else { + s7 = bsize; + j7 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov7[ k ] = ov[k] + ( j7*sv[k][7] ); + } + for ( j6 = sh[6]; j6 > 0; ) { + if ( j6 < bsize ) { + s6 = j6; + j6 = 0; + } else { + s6 = bsize; + j6 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov6[ k ] = ov7[k] + ( j6*sv[k][6] ); + dv7[ k ] = sv[k][7] - ( s6*sv[k][6] ); + } + for ( j5 = sh[5]; j5 > 0; ) { + if ( j5 < bsize ) { + s5 = j5; + j5 = 0; + } else { + s5 = bsize; + j5 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov5[ k ] = ov6[k] + ( j5*sv[k][5] ); + dv6[ k ] = sv[k][6] - ( s5*sv[k][5] ); + } + for ( j4 = sh[4]; j4 > 0; ) { + if ( j4 < bsize ) { + s4 = j4; + j4 = 0; + } else { + s4 = bsize; + j4 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov4[ k ] = ov5[k] + ( j4*sv[k][4] ); + dv5[ k ] = sv[k][5] - ( s4*sv[k][4] ); + } + for ( j3 = sh[3]; j3 > 0; ) { + if ( j3 < bsize ) { + s3 = j3; + j3 = 0; + } else { + s3 = bsize; + j3 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov3[ k ] = ov4[k] + ( j3*sv[k][3] ); + dv4[ k ] = sv[k][4] - ( s3*sv[k][3] ); + } + for ( j2 = sh[2]; j2 > 0; ) { + if ( j2 < bsize ) { + s2 = j2; + j2 = 0; + } else { + s2 = bsize; + j2 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov2[ k ] = ov3[k] + ( j2*sv[k][2] ); + dv3[ k ] = sv[k][3] - ( s2*sv[k][2] ); + } + for ( j1 = sh[1]; j1 > 0; ) { + if ( j1 < bsize ) { + s1 = j1; + j1 = 0; + } else { + s1 = bsize; + j1 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov1[ k ] = ov2[k] + ( j1*sv[k][1] ); + dv2[ k ] = sv[k][2] - ( s1*sv[k][1] ); + } + for ( j0 = sh[0]; j0 > 0; ) { + if ( j0 < bsize ) { + s0 = j0; + j0 = 0; + } else { + s0 = bsize; + j0 -= bsize; + } + // Compute index offsets and loop offset increments for the first ndarray elements in the current block... + for ( k = 0; k < N; k++ ) { + iv[ k ] = ov1[k] + ( j0*sv[k][0] ); + dv1[ k ] = sv[k][1] - ( s0*sv[k][0] ); + } + // Iterate over the loop dimensions... + for ( i7 = 0; i7 < s7; i7++ ) { + for ( i6 = 0; i6 < s6; i6++ ) { + for ( i5 = 0; i5 < s5; i5++ ) { + for ( i4 = 0; i4 < s4; i4++ ) { + for ( i3 = 0; i3 < s3; i3++ ) { + for ( i2 = 0; i2 < s2; i2++ ) { + for ( i1 = 0; i1 < s1; i1++ ) { + for ( i0 = 0; i0 < s0; i0++ ) { + setViewOffsets( views, iv ); + v[ 0 ] = strategyX.input( views[ 0 ] ); + f = wrap( ref, views[ 0 ], ibuf, ldims, [ j7+i7, j6+i6, j5+i5, j4+i4, j3+i3, j2+i2, j1+i1, j0+i0 ], cdims, clbk, thisArg ); + if ( hasOpts ) { + fcn( v, opts, f ); + } else { + fcn( v, f ); + } + strategyX.output( views[ 0 ] ); + incrementOffsets( iv, dv0 ); + } + incrementOffsets( iv, dv1 ); + } + incrementOffsets( iv, dv2 ); + } + incrementOffsets( iv, dv3 ); + } + incrementOffsets( iv, dv4 ); + } + incrementOffsets( iv, dv5 ); + } + incrementOffsets( iv, dv6 ); + } + incrementOffsets( iv, dv7 ); + } + } + } + } + } + } + } + } + } +} + + +// EXPORTS // + +module.exports = kernel; diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/9d.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/9d.js new file mode 100644 index 000000000000..b8bf9c47b7db --- /dev/null +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/9d.js @@ -0,0 +1,417 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +/* eslint-disable max-len, max-params, max-depth */ + +'use strict'; + +// MODULES // + +var loopOrder = require( '@stdlib/ndarray/base/nullary-loop-interchange-order' ); +var blockSize = require( '@stdlib/ndarray/base/nullary-tiling-block-size' ); +var takeIndexed = require( '@stdlib/array/base/take-indexed' ); +var copyIndexed = require( '@stdlib/array/base/copy-indexed' ); +var zeros = require( '@stdlib/array/base/zeros' ); +var incrementOffsets = require( '@stdlib/ndarray/base/kernels/utils/increment-offsets' ); +var setViewOffsets = require( '@stdlib/ndarray/base/set-descriptor-offsets' ); +var offsets = require( '@stdlib/ndarray/base/offsets' ); +var wrap = require( './callback_wrapper.js' ); + + +// MAIN // + +/** +* Applies a one-dimensional strided array function to an input ndarray according to a callback function via loop blocking. +* +* @private +* @param {Function} fcn - wrapper for a one-dimensional strided array function +* @param {Array} arrays - ndarray descriptors +* @param {ndarrayLike} ref - original input ndarray provided to a callback function +* @param {Array} views - initialized ndarray descriptors representing subarray views +* @param {Array} shape - loop dimensions +* @param {Array} stridesX - loop dimension strides for the input ndarray +* @param {Object} strategyX - strategy for marshaling data to and from an input ndarray view +* @param {NonNegativeIntegerArray} ibuf - workspace for storing iteration indices +* @param {NonNegativeIntegerArray} ldims - list of loop dimensions +* @param {NonNegativeIntegerArray} cdims - list of "core" dimensions +* @param {Options} opts - function options +* @param {boolean} hasOpts - boolean indicating whether to pass an options argument to a strided array function +* @param {Function} clbk - callback function +* @param {*} thisArg - callback execution context +* @returns {void} +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); +* var zeros = require( '@stdlib/array/base/zeros' ); +* var gfillBy = require( '@stdlib/blas/ext/base/ndarray/gfill-by' ); +* var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' ); +* +* function clbk( value, idx ) { +* return ( idx[ 8 ]*100 ) + ( idx[ 9 ]*10 ) + idx[ 10 ]; +* } +* +* // Create a data buffer: +* var xbuf = new Float64Array( 12 ); +* +* // Define the array shape: +* var xsh = [ 1, 1, 1, 1, 1, 1, 1, 1, 3, 2, 2 ]; +* +* // Define the array strides: +* var sx = [ 12, 12, 12, 12, 12, 12, 12, 12, 4, 2, 1 ]; +* +* // Define the index offset: +* var ox = 0; +* +* // Create an input ndarray descriptor: +* var x = { +* 'dtype': 'float64', +* 'data': xbuf, +* 'shape': xsh, +* 'strides': sx, +* 'offset': ox, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the starting index (inclusive): +* var start = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 0.0 ] ), +* 'shape': [ 1, 1, 1, 1, 1, 1, 1, 1, 3 ], +* 'strides': [ 0, 0, 0, 0, 0, 0, 0, 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Create an ndarray descriptor for the ending index (exclusive): +* var end = { +* 'dtype': 'float64', +* 'data': new Float64Array( [ 4.0 ] ), +* 'shape': [ 1, 1, 1, 1, 1, 1, 1, 1, 3 ], +* 'strides': [ 0, 0, 0, 0, 0, 0, 0, 0, 0 ], +* 'offset': 0, +* 'order': 'row-major' +* }; +* +* // Initialize ndarray descriptors representing subarray views: +* var views = [ +* { +* 'dtype': x.dtype, +* 'data': x.data, +* 'shape': [ 2, 2 ], +* 'strides': [ 2, 1 ], +* 'offset': x.offset, +* 'order': x.order +* }, +* { +* 'dtype': start.dtype, +* 'data': start.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': start.offset, +* 'order': start.order +* }, +* { +* 'dtype': end.dtype, +* 'data': end.data, +* 'shape': [], +* 'strides': [ 0 ], +* 'offset': end.offset, +* 'order': end.order +* } +* ]; +* +* // Resolve an input strategy when iterating over subarray views: +* var strategyX = strategy( views[ 0 ] ); +* +* // Create a workspace array for storing iteration indices: +* var ibuf = zeros( xsh.length ); +* +* // Define the loop and core dimensions: +* var ldims = [ 0, 1, 2, 3, 4, 5, 6, 7, 8 ]; +* var cdims = [ 9, 10 ]; +* +* // Apply strided function: +* kernel( gfillBy, [ x, start, end ], x, views, [ 1, 1, 1, 1, 1, 1, 1, 1, 3 ], [ 12, 12, 12, 12, 12, 12, 12, 12, 4 ], strategyX, ibuf, ldims, cdims, {}, false, clbk, {} ); +* +* var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); +* // returns [ [ [ [ [ [ [ [ [ [ [ 0.0, 1.0 ], [ 10.0, 11.0 ] ], [ [ 100.0, 101.0 ], [ 110.0, 111.0 ] ], [ [ 200.0, 201.0 ], [ 210.0, 211.0 ] ] ] ] ] ] ] ] ] ] ] +*/ +function kernel( fcn, arrays, ref, views, shape, stridesX, strategyX, ibuf, ldims, cdims, opts, hasOpts, clbk, thisArg ) { // eslint-disable-line max-statements, max-lines-per-function + var bsize; + var dv0; + var dv1; + var dv2; + var dv3; + var dv4; + var dv5; + var dv6; + var dv7; + var dv8; + var ov1; + var ov2; + var ov3; + var ov4; + var ov5; + var ov6; + var ov7; + var ov8; + var sh; + var s0; + var s1; + var s2; + var s3; + var s4; + var s5; + var s6; + var s7; + var s8; + var sv; + var ov; + var iv; + var i0; + var i1; + var i2; + var i3; + var i4; + var i5; + var i6; + var i7; + var i8; + var j0; + var j1; + var j2; + var j3; + var j4; + var j5; + var j6; + var j7; + var j8; + var N; + var x; + var v; + var o; + var k; + var f; + + // Note on variable naming convention: S#, dv#, i#, j# where # corresponds to the loop number, with `0` being the innermost loop... + + N = arrays.length; + x = arrays[ 0 ]; + + // Resolve the loop interchange order: + o = loopOrder( shape, stridesX ); + sh = o.sh; + sv = [ o.sx ]; + for ( k = 1; k < N; k++ ) { + sv.push( takeIndexed( arrays[k].strides, o.idx ) ); + } + // Determine the block size: + bsize = blockSize( x.dtype ); + + // Resolve a list of pointers to the first indexed elements in the respective ndarrays: + ov = offsets( arrays ); + + // Cache offset increments for the innermost loop... + dv0 = []; + for ( k = 0; k < N; k++ ) { + dv0.push( sv[k][0] ); + } + // Initialize loop variables... + ov1 = zeros( N ); + ov2 = zeros( N ); + ov3 = zeros( N ); + ov4 = zeros( N ); + ov5 = zeros( N ); + ov6 = zeros( N ); + ov7 = zeros( N ); + ov8 = zeros( N ); + dv1 = zeros( N ); + dv2 = zeros( N ); + dv3 = zeros( N ); + dv4 = zeros( N ); + dv5 = zeros( N ); + dv6 = zeros( N ); + dv7 = zeros( N ); + dv8 = zeros( N ); + iv = zeros( N ); + + // Shallow copy the list of views to an internal array so that we can update with reshaped views without impacting the original list of views: + v = copyIndexed( views ); + + // Iterate over blocks... + for ( j8 = sh[8]; j8 > 0; ) { + if ( j8 < bsize ) { + s8 = j8; + j8 = 0; + } else { + s8 = bsize; + j8 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov8[ k ] = ov[k] + ( j8*sv[k][8] ); + } + for ( j7 = sh[7]; j7 > 0; ) { + if ( j7 < bsize ) { + s7 = j7; + j7 = 0; + } else { + s7 = bsize; + j7 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov7[ k ] = ov8[k] + ( j7*sv[k][7] ); + dv8[ k ] = sv[k][8] - ( s7*sv[k][7] ); + } + for ( j6 = sh[6]; j6 > 0; ) { + if ( j6 < bsize ) { + s6 = j6; + j6 = 0; + } else { + s6 = bsize; + j6 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov6[ k ] = ov7[k] + ( j6*sv[k][6] ); + dv7[ k ] = sv[k][7] - ( s6*sv[k][6] ); + } + for ( j5 = sh[5]; j5 > 0; ) { + if ( j5 < bsize ) { + s5 = j5; + j5 = 0; + } else { + s5 = bsize; + j5 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov5[ k ] = ov6[k] + ( j5*sv[k][5] ); + dv6[ k ] = sv[k][6] - ( s5*sv[k][5] ); + } + for ( j4 = sh[4]; j4 > 0; ) { + if ( j4 < bsize ) { + s4 = j4; + j4 = 0; + } else { + s4 = bsize; + j4 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov4[ k ] = ov5[k] + ( j4*sv[k][4] ); + dv5[ k ] = sv[k][5] - ( s4*sv[k][4] ); + } + for ( j3 = sh[3]; j3 > 0; ) { + if ( j3 < bsize ) { + s3 = j3; + j3 = 0; + } else { + s3 = bsize; + j3 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov3[ k ] = ov4[k] + ( j3*sv[k][3] ); + dv4[ k ] = sv[k][4] - ( s3*sv[k][3] ); + } + for ( j2 = sh[2]; j2 > 0; ) { + if ( j2 < bsize ) { + s2 = j2; + j2 = 0; + } else { + s2 = bsize; + j2 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov2[ k ] = ov3[k] + ( j2*sv[k][2] ); + dv3[ k ] = sv[k][3] - ( s2*sv[k][2] ); + } + for ( j1 = sh[1]; j1 > 0; ) { + if ( j1 < bsize ) { + s1 = j1; + j1 = 0; + } else { + s1 = bsize; + j1 -= bsize; + } + for ( k = 0; k < N; k++ ) { + ov1[ k ] = ov2[k] + ( j1*sv[k][1] ); + dv2[ k ] = sv[k][2] - ( s1*sv[k][1] ); + } + for ( j0 = sh[0]; j0 > 0; ) { + if ( j0 < bsize ) { + s0 = j0; + j0 = 0; + } else { + s0 = bsize; + j0 -= bsize; + } + // Compute index offsets and loop offset increments for the first ndarray elements in the current block... + for ( k = 0; k < N; k++ ) { + iv[ k ] = ov1[k] + ( j0*sv[k][0] ); + dv1[ k ] = sv[k][1] - ( s0*sv[k][0] ); + } + // Iterate over the loop dimensions... + for ( i8 = 0; i8 < s8; i8++ ) { + for ( i7 = 0; i7 < s7; i7++ ) { + for ( i6 = 0; i6 < s6; i6++ ) { + for ( i5 = 0; i5 < s5; i5++ ) { + for ( i4 = 0; i4 < s4; i4++ ) { + for ( i3 = 0; i3 < s3; i3++ ) { + for ( i2 = 0; i2 < s2; i2++ ) { + for ( i1 = 0; i1 < s1; i1++ ) { + for ( i0 = 0; i0 < s0; i0++ ) { + setViewOffsets( views, iv ); + v[ 0 ] = strategyX.input( views[ 0 ] ); + f = wrap( ref, views[ 0 ], ibuf, ldims, [ j8+i8, j7+i7, j6+i6, j5+i5, j4+i4, j3+i3, j2+i2, j1+i1, j0+i0 ], cdims, clbk, thisArg ); + if ( hasOpts ) { + fcn( v, opts, f ); + } else { + fcn( v, f ); + } + strategyX.output( views[ 0 ] ); + incrementOffsets( iv, dv0 ); + } + incrementOffsets( iv, dv1 ); + } + incrementOffsets( iv, dv2 ); + } + incrementOffsets( iv, dv3 ); + } + incrementOffsets( iv, dv4 ); + } + incrementOffsets( iv, dv5 ); + } + incrementOffsets( iv, dv6 ); + } + incrementOffsets( iv, dv7 ); + } + incrementOffsets( iv, dv8 ); + } + } + } + } + } + } + } + } + } + } +} + + +// EXPORTS // + +module.exports = kernel; diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/index.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/index.js index 9a737f43d174..e6eb199aa5c1 100644 --- a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/index.js +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/index.js @@ -130,6 +130,13 @@ var setReadOnly = require( '@stdlib/utils/define-nonenumerable-read-only-property' ); var kernel2d = require( './2d.js' ); var kernel3d = require( './3d.js' ); +var kernel4d = require( './4d.js' ); +var kernel5d = require( './5d.js' ); +var kernel6d = require( './6d.js' ); +var kernel7d = require( './7d.js' ); +var kernel8d = require( './8d.js' ); +var kernel9d = require( './9d.js' ); +var kernel10d = require( './10d.js' ); var main = require( './main.js' ); @@ -137,6 +144,13 @@ var main = require( './main.js' ); setReadOnly( main, 'kernel2d', kernel2d ); setReadOnly( main, 'kernel3d', kernel3d ); +setReadOnly( main, 'kernel4d', kernel4d ); +setReadOnly( main, 'kernel5d', kernel5d ); +setReadOnly( main, 'kernel6d', kernel6d ); +setReadOnly( main, 'kernel7d', kernel7d ); +setReadOnly( main, 'kernel8d', kernel8d ); +setReadOnly( main, 'kernel9d', kernel9d ); +setReadOnly( main, 'kernel10d', kernel10d ); // EXPORTS // diff --git a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/main.js b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/main.js index 90714b9c2fea..a08f4bfb3982 100644 --- a/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/main.js +++ b/lib/node_modules/@stdlib/ndarray/base/kernels/generic/nullary-strided1d-by/blocked/lib/main.js @@ -22,13 +22,27 @@ var kernel2d = require( './2d.js' ); var kernel3d = require( './3d.js' ); +var kernel4d = require( './4d.js' ); +var kernel5d = require( './5d.js' ); +var kernel6d = require( './6d.js' ); +var kernel7d = require( './7d.js' ); +var kernel8d = require( './8d.js' ); +var kernel9d = require( './9d.js' ); +var kernel10d = require( './10d.js' ); // VARIABLES // var KERNELS = [ kernel2d, - kernel3d + kernel3d, + kernel4d, + kernel5d, + kernel6d, + kernel7d, + kernel8d, + kernel9d, + kernel10d ]; var MAX_DIMS = KERNELS.length + 1;