refactor: add support for ancillary ndarray arguments having trailing dimensions

---
type: pre_commit_static_analysis_report
description: Results of running static analysis checks when committing changes.
report:
  - task: lint_filenames
    status: passed
  - task: lint_editorconfig
    status: passed
  - task: lint_markdown
    status: passed
  - task: lint_package_json
    status: na
  - task: lint_repl_help
    status: passed
  - task: lint_javascript_src
    status: passed
  - task: lint_javascript_cli
    status: na
  - task: lint_javascript_examples
    status: na
  - task: lint_javascript_tests
    status: na
  - task: lint_javascript_benchmarks
    status: na
  - task: lint_python
    status: na
  - task: lint_r
    status: na
  - task: lint_c_src
    status: na
  - task: lint_c_examples
    status: na
  - task: lint_c_benchmarks
    status: na
  - task: lint_c_tests_fixtures
    status: na
  - task: lint_shell
    status: na
  - task: lint_typescript_declarations
    status: passed
  - task: lint_typescript_tests
    status: na
  - task: lint_license_headers
    status: passed
---

Author: kgryte

lib/node_modules/@stdlib/ndarray/base/binary-reduce-strided1d/README.md

@@ -128,7 +128,11 @@ Each provided ndarray should be an object with the following properties:
 
 ## Notes
 
--   The output ndarray and any additional ndarray arguments are expected to have the same dimensions as the non-reduced dimensions of the input ndarrays. When calling the reduction function, any additional ndarray arguments are provided as zero-dimensional ndarray-like objects.
+-   The output ndarray is expected to have the same dimensions as the non-reduced dimensions of the input ndarrays.
+
+-   Any additional ndarray arguments are expected to have the same leading dimensions as the non-reduced dimensions of the input ndarrays.
+
+-   When calling the reduction function, any additional ndarray arguments are provided as k-dimensional subarrays, where `k = M - N` with `M` being the number of dimensions in an ndarray argument and `N` being the number of non-reduced dimensions in the input ndarrays. For example, if an input ndarrays have three dimensions, the number of reduced dimensions is two, and an additional ndarray argument has one dimension, thus matching the number of non-reduced dimensions in the input ndarrays, the reduction function is provided a zero-dimensional subarray as an additional ndarray argument. In the same scenario but where an additional ndarray argument has two dimensions, thus exceeding the number of non-reduced dimensions in the input ndarrays, the reduction function is provided a one-dimensional subarray as an additional ndarray argument.
 
 -   The reduction function is expected to have the following signature:
 
@@ -138,7 +142,7 @@ Each provided ndarray should be an object with the following properties:
 
     where
 
-    -   **arrays**: array containing a one-dimensional subarray for each input ndarray and any additional ndarray arguments as zero-dimensional ndarrays.
+    -   **arrays**: array containing a one-dimensional subarray for each input ndarray and any additional ndarray arguments as subarrays.
     -   **options**: function options (_optional_).
 
 -   For very high-dimensional ndarrays which are non-contiguous, one should consider copying the underlying data to contiguous memory before performing a reduction in order to achieve better performance.

lib/node_modules/@stdlib/ndarray/base/binary-reduce-strided1d/docs/repl.txt

@@ -14,10 +14,24 @@
     - order: specifies whether an ndarray is row-major (C-style) or column-major
     (Fortran-style).
 
-    The output ndarray and any additional ndarray arguments are expected to have
-    the same dimensions as the non-reduced dimensions of the input ndarray. When
-    calling the reduction function, any additional ndarray arguments are
-    provided as zero-dimensional ndarray-like objects.
+    The output ndarray is expected to have the same dimensions as the non-
+    reduced dimensions of the input ndarrays.
+
+    Any additional ndarray arguments are expected to have the same leading
+    dimensions as the non-reduced dimensions of the input ndarrays.
+
+    When calling the reduction function, any additional ndarray arguments are
+    provided as k-dimensional subarrays, where `k = M - N` with `M` being the
+    number of dimensions in an ndarray argument and `N` being the number of non-
+    reduced dimensions in the input ndarrays. For example, if an input ndarrays
+    have three dimensions, the number of reduced dimensions is two, and an
+    additional ndarray argument has one dimension, thus matching the number of
+    non-reduced dimensions in the input ndarrays, the reduction function is
+    provided a zero-dimensional subarray as an additional ndarray argument. In
+    the same scenario but where an additional ndarray argument has two
+    dimensions, thus exceeding the number of non-reduced dimensions in the input
+    ndarrays, the reduction function is provided a one-dimensional subarray as
+    an additional ndarray argument.
 
     Parameters
     ----------
@@ -31,8 +45,7 @@
         where
 
         - arrays: array containing a one-dimensional subarrays for each input
-        ndarray and any additional ndarray arguments as zero-dimensional
-        ndarrays.
+        ndarray and any additional ndarray arguments as subarrays.
         - options: function options.
 
     arrays: ArrayLikeObject<ndarray>

lib/node_modules/@stdlib/ndarray/base/binary-reduce-strided1d/lib/initialize_array_views.js

@@ -18,10 +18,15 @@
 
 'use strict';
 
+// MODULES //
+
+var slice = require( '@stdlib/array/base/slice' );
+
+
 // MAIN //
 
 /**
-* Initialize ndarray-like objects for representing zero-dimensional sub-array views of ancillary ndarray arguments.
+* Initialize ndarray-like objects for representing sub-array views of ancillary ndarray arguments.
 *
 * ## Notes
 *
@@ -30,20 +35,25 @@
 *
 * @private
 * @param {ArrayLikeObject<Object>} arrays - list of ndarray-like objects
+* @param {NonNegativeInteger} k - number of non-reduced dimensions
 * @param {Array<Object>} out - output array
 * @returns {Array<Object>} output array
 */
-function initializeViews( arrays, out ) {
+function initializeViews( arrays, k, out ) {
+	var sh;
+	var N;
 	var v;
 	var i;
 
 	for ( i = 3; i < arrays.length; i++ ) {
 		v = arrays[ i ];
+		sh = v.shape;
+		N = sh.length;
 		out.push({
 			'dtype': v.dtype,
 			'data': v.data,
-			'shape': [],
-			'strides': [ 0 ],
+			'shape': slice( sh, k, N ),
+			'strides': ( N === k ) ? [ 0 ] : slice( v.strides, k, N ),
 			'offset': v.offset,
 			'order': v.order
 		});

lib/node_modules/@stdlib/ndarray/base/binary-reduce-strided1d/lib/main.js

@@ -418,7 +418,6 @@ function binaryReduceStrided1d( fcn, arrays, dims, options ) { // eslint-disable
 			throw new Error( format( 'invalid argument. Input arrays must have the same shape. First array shape: [%s]. Second array shape: [%s].', join( shx, ',' ), join( shy, ',' ) ) );
 		}
 	}
-
 	// Verify that we've been provided a list of unique dimension indices...
 	M = dims.length;
 	d = normalizeIndices( dims, ndims-1 );
@@ -433,12 +432,19 @@ function binaryReduceStrided1d( fcn, arrays, dims, options ) { // eslint-disable
 	if ( M > ndims ) {
 		throw new RangeError( format( 'invalid argument. Number of specified dimensions cannot exceed the number of dimensions in the input array. Number of dimensions: %d. Value: [%s].', ndims, join( dims, ',' ) ) );
 	}
-	// Verify that provided ndarrays have the expected number of dimensions...
+	// Compute the number of non-reduced dimensions:
 	K = ndims - M;
+
+	// Verify that the output ndarray has the expected number of dimensions...
 	if ( z.shape.length !== K ) {
-		throw new Error( format( 'invalid argument. Output array must have the same number of non-reduced dimensions as input arrays. Input array shape: [%s]. Number of non-reduced dimensions: %d. Output array shape: [%s].', join( shx, ',' ), K, join( z.shape, ',' ) ) );
+		throw new Error( format( 'invalid argument. Arrays which are not being reduced must have the same number of non-reduced dimensions. Input array shape: [%s]. Number of non-reduced dimensions: %d. Array shape: [%s] (index: %d).', join( shx, ',' ), K, join( z.shape, ',' ), 2 ) );
+	}
+	// Verify that any ancillary ndarrays have at least the number of non-reduced dimensions...
+	for ( i = 3; i < N; i++ ) {
+		if ( arr[ i ].shape.length < K ) {
+			throw new Error( format( 'invalid argument. Arrays which are not being reduced must have the same number of non-reduced dimensions. Input array shape: [%s]. Number of non-reduced dimensions: %d. Array shape: [%s] (index: %d).', join( shx, ',' ), K, join( arr[ i ].shape, ',' ), i ) );
+		}
 	}
-
 	// Resolve the non-reduced ("loop") dimensions and associated strides:
 	ldims = indicesComplement( shx.length, d );
 	tmpx = takeIndexed2( shx, x.strides, ldims );
@@ -495,7 +501,7 @@ function binaryReduceStrided1d( fcn, arrays, dims, options ) { // eslint-disable
 			'order': y.order
 		}
 	];
-	initializeViews( arr, views );
+	initializeViews( arr, K, views );
 
 	// Determine the strategy for reshaping sub-array views of the input arrays prior to performing a reduction:
 	strategyX = reshapeStrategy( views[ 0 ] );
@@ -515,8 +521,6 @@ function binaryReduceStrided1d( fcn, arrays, dims, options ) { // eslint-disable
 		}
 		return BINARY[ K ]( fcn, arr, views, slx, sly, strategyX, strategyY, opts ); // eslint-disable-line max-len
 	}
-	sz = z.strides;
-
 	// Determine whether the loop dimensions have only **one** non-singleton dimension (e.g., shape=[10,1,1,1]) so that we can treat loop iteration as being equivalent to one-dimensional iteration...
 	if ( ns === K-1 ) {
 		// Get the index of the non-singleton dimension...
@@ -536,6 +540,7 @@ function binaryReduceStrided1d( fcn, arrays, dims, options ) { // eslint-disable
 		}
 		return BINARY[ 1 ]( fcn, arr, views, slx, sly, strategyX, strategyY, opts ); // eslint-disable-line max-len
 	}
+	sz = z.strides;
 	iox = iterationOrder( slx ); // +/-1
 	ioy = iterationOrder( sly ); // +/-1
 	ioz = iterationOrder( sz ); // +/-1