diff --git a/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.0d.js b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.0d.js
new file mode 100644
index 000000000000..d198042a8869
--- /dev/null
+++ b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.0d.js
@@ -0,0 +1,111 @@
+/**
+* @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 Float64Array = require( '@stdlib/array/float64' );
+var toAccessorArray = require( '@stdlib/array/base/to-accessor-array' );
+var Complex128 = require( '@stdlib/complex/float64/ctor' );
+var real = require( '@stdlib/complex/float64/real' );
+var imag = require( '@stdlib/complex/float64/imag' );
+var scalar2ndarray = require( '@stdlib/ndarray/from-scalar' );
+var ndarray = require( '@stdlib/ndarray/ctor' );
+var accumulateUnary = require( './../lib' );
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+	t.ok( true, __filename );
+	t.strictEqual( typeof accumulateUnary, 'function', 'main export is a function' );
+	t.end();
+});
+
+tape( 'the function performs accumulation on a 0-dimensional ndarray', function test( t ) {
+	var actual;
+	var x;
+
+	x = scalar2ndarray( 0.0, {
+		'dtype': 'float64'
+	});
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = scalar2ndarray( 5.0, {
+		'dtype': 'float64'
+	});
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 5.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 0-dimensional ndarray (accessors)', function test( t ) {
+	var actual;
+	var x;
+
+	x = ndarray( 'generic', toAccessorArray( new Float64Array( [ 0.0 ] ) ), [], [ 0 ], 0, 'row-major' );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( 'generic', toAccessorArray( new Float64Array( [ 5.0 ] ) ), [], [ 0 ], 0, 'row-major' );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 5.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 0-dimensional ndarray (complex)', function test( t ) {
+	var actual;
+	var x;
+
+	x = scalar2ndarray( new Complex128( 0.0, 0.0 ), {
+		'dtype': 'complex128'
+	});
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = scalar2ndarray( new Complex128( 3.0, 4.0 ), {
+		'dtype': 'complex128'
+	});
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 7.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
diff --git a/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.10d.js b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.10d.js
new file mode 100644
index 000000000000..736f3d8dc2ee
--- /dev/null
+++ b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.10d.js
@@ -0,0 +1,243 @@
+/**
+* @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 zeros = require( '@stdlib/array/zeros' );
+var ones = require( '@stdlib/array/ones' );
+var real = require( '@stdlib/complex/float64/real' );
+var imag = require( '@stdlib/complex/float64/imag' );
+var Complex128Array = require( '@stdlib/array/complex128' );
+var numel = require( '@stdlib/ndarray/base/numel' );
+var shape2strides = require( '@stdlib/ndarray/base/shape2strides' );
+var strides2offset = require( '@stdlib/ndarray/base/strides2offset' );
+var ndarray = require( '@stdlib/ndarray/ctor' );
+var accumulateUnary = require( './../lib' );
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+	t.ok( true, __filename );
+	t.strictEqual( typeof accumulateUnary, 'function', 'main export is a function' );
+	t.end();
+});
+
+tape( 'the function performs accumulation on a 10-dimensional ndarray (row-major, singleton dimensions)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 4, 1, 1, 1, 1, 1, 1, 1, 1, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 10-dimensional ndarray (row-major, singleton dimensions, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 4, 1, 1, 1, 1, 1, 1, 1, 1, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 10-dimensional ndarray (row-major, contiguous)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 1024.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 10-dimensional ndarray (row-major, contiguous, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 2048.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 10-dimensional ndarray (column-major, contiguous)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 1024.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 10-dimensional ndarray (column-major, contiguous, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-major';
+	sh = [ 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 2048.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
diff --git a/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.1d.js b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.1d.js
new file mode 100644
index 000000000000..2edd1b8689cd
--- /dev/null
+++ b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.1d.js
@@ -0,0 +1,84 @@
+/**
+* @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 ones = require( '@stdlib/array/ones' );
+var zeros = require( '@stdlib/array/zeros' );
+var real = require( '@stdlib/complex/float64/real' );
+var imag = require( '@stdlib/complex/float64/imag' );
+var Complex128Array = require( '@stdlib/array/complex128' );
+var ndarray = require( '@stdlib/ndarray/ctor' );
+var accumulateUnary = require( './../lib' );
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+	t.ok( true, __filename );
+	t.strictEqual( typeof accumulateUnary, 'function', 'main export is a function' );
+	t.end();
+});
+
+tape( 'the function performs accumulation on a 1-dimensional ndarray', function test( t ) {
+	var actual;
+	var x;
+
+	x = ndarray( 'float64', zeros( 8, 'float64' ), [ 4 ], [ 2 ], 1, 'row-major' );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( 'float64', ones( 8, 'float64' ), [ 4 ], [ 2 ], 1, 'row-major' );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 1-dimensional ndarray (accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var x;
+
+	xbuf = zeros( 6*2, 'float64' );
+	x = ndarray( 'complex128', new Complex128Array( xbuf ), [ 4 ], [ 1 ], 1, 'row-major' );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( 6*2, 'float64' );
+	x = ndarray( 'complex128', new Complex128Array( xbuf ), [ 4 ], [ 1 ], 1, 'row-major' );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
diff --git a/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.2d.js b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.2d.js
new file mode 100644
index 000000000000..4478a85701d9
--- /dev/null
+++ b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.2d.js
@@ -0,0 +1,1005 @@
+/**
+* @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 zeros = require( '@stdlib/array/zeros' );
+var ones = require( '@stdlib/array/ones' );
+var real = require( '@stdlib/complex/float64/real' );
+var imag = require( '@stdlib/complex/float64/imag' );
+var Complex128Array = require( '@stdlib/array/complex128' );
+var numel = require( '@stdlib/ndarray/base/numel' );
+var shape2strides = require( '@stdlib/ndarray/base/shape2strides' );
+var strides2offset = require( '@stdlib/ndarray/base/strides2offset' );
+var ndarray = require( '@stdlib/ndarray/ctor' );
+var blockSize = require( '@stdlib/ndarray/base/nullary-tiling-block-size' );
+var accumulateUnary = require( './../lib' );
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+	t.ok( true, __filename );
+	t.strictEqual( typeof accumulateUnary, 'function', 'main export is a function' );
+	t.end();
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (row-major, singleton dimensions)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 4, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (row-major, singleton dimensions, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 4, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (row-major, contiguous)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (row-major, contiguous, negative strides)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 2, 2 ];
+	st = [ -2, -1 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (row-major, non-contiguous, same sign strides)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 2, 2 ];
+	st = [ 4, 2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( 8, dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( 8, dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (row-major, non-contiguous, mixed sign strides)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 2, 2 ];
+	st = [ 4, -2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( 8, dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( 8, dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (row-major, non-contiguous, large arrays)', function test( t ) {
+	var actual;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ bsize*2, 2 ];
+	st = [ -4, 2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh )*2, dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, bsize*4, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (row-major, non-contiguous, large arrays)', function test( t ) {
+	var actual;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, bsize*2 ];
+	st = [ bsize*4, -2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh )*2, dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, bsize*4, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (row-major, contiguous, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (row-major, contiguous, negative strides, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 2, 2 ];
+	st = [ -2, -1 ];
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (row-major, non-contiguous, same sign strides, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 2, 2 ];
+	st = [ 4, 2 ];
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( 16, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( 16, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (row-major, non-contiguous, mixed sign strides, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 2, 2 ];
+	st = [ 4, -2 ];
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( 16, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( 16, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (row-major, non-contiguous, large arrays, accessors)', function test( t ) {
+	var actual;
+	var bsize;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ bsize*2, 2 ];
+	st = [ -4, 2 ];
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*4, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*4, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, bsize*8, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (row-major, non-contiguous, large arrays, accessors)', function test( t ) {
+	var actual;
+	var bsize;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, bsize*2 ];
+	st = [ bsize*4, -2 ];
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*4, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*4, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, bsize*8, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (column-major, singleton dimensions)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 4, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (column-major, singleton dimensions, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-major';
+	sh = [ 4, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (column-major, contiguous)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (column-major, contiguous, negative strides)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 2, 2 ];
+	st = [ -1, -2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (column-major, non-contiguous, same sign strides)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 2, 2 ];
+	st = [ 2, 4 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( 8, dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( 8, dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (column-major, non-contiguous, mixed sign strides)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 2, 2 ];
+	st = [ -2, 4 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( 8, dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( 8, dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (column-major, non-contiguous, large arrays)', function test( t ) {
+	var actual;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+
+	bsize = blockSize( dt );
+	sh = [ bsize*2, 2 ];
+	st = [ -2, bsize*4 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh )*2, dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, bsize*4, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (column-major, non-contiguous, large arrays)', function test( t ) {
+	var actual;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, bsize*2 ];
+	st = [ bsize*4, -2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh )*2, dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, bsize*4, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (column-major, contiguous, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-major';
+	sh = [ 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (column-major, contiguous, negative strides, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-major';
+	sh = [ 2, 2 ];
+	st = [ -1, -2 ];
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (column-major, non-contiguous, same sign strides, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-major';
+	sh = [ 2, 2 ];
+	st = [ 2, 4 ];
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( 16, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( 16, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (column-major, non-contiguous, mixed sign strides, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-major';
+	sh = [ 2, 2 ];
+	st = [ -2, 4 ];
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( 16, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( 16, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (column-major, non-contiguous, large arrays, accessors)', function test( t ) {
+	var actual;
+	var bsize;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-major';
+
+	bsize = blockSize( dt );
+	sh = [ bsize*2, 2 ];
+	st = [ -2, bsize*4 ];
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*4, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*4, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, bsize*8, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 2-dimensional ndarray (column-major, non-contiguous, large arrays, accessors)', function test( t ) {
+	var actual;
+	var bsize;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, bsize*2 ];
+	st = [ bsize*4, -2 ];
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*4, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*4, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, bsize*8, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
diff --git a/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.3d.js b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.3d.js
new file mode 100644
index 000000000000..32d478bb3d85
--- /dev/null
+++ b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.3d.js
@@ -0,0 +1,346 @@
+/**
+* @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 zeros = require( '@stdlib/array/zeros' );
+var ones = require( '@stdlib/array/ones' );
+var real = require( '@stdlib/complex/float64/real' );
+var imag = require( '@stdlib/complex/float64/imag' );
+var Complex128Array = require( '@stdlib/array/complex128' );
+var numel = require( '@stdlib/ndarray/base/numel' );
+var shape2strides = require( '@stdlib/ndarray/base/shape2strides' );
+var strides2offset = require( '@stdlib/ndarray/base/strides2offset' );
+var ndarray = require( '@stdlib/ndarray/ctor' );
+var blockSize = require( '@stdlib/ndarray/base/nullary-tiling-block-size' );
+var accumulateUnary = require( './../lib' );
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+	t.ok( true, __filename );
+	t.strictEqual( typeof accumulateUnary, 'function', 'main export is a function' );
+	t.end();
+});
+
+tape( 'the function performs accumulation on a 3-dimensional ndarray (row-major, singleton dimensions)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 4, 1, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 3-dimensional ndarray (row-major, singleton dimensions, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 4, 1, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 3-dimensional ndarray (row-major, contiguous)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 3-dimensional ndarray (row-major, non-contiguous, large arrays)', function test( t ) {
+	var actual;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ bsize, 2, 2 ];
+	st = [ 8, 4, 1 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh )*2, dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, bsize*4, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 3-dimensional ndarray (row-major, contiguous, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 16.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 3-dimensional ndarray (column-major, singleton dimensions)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 4, 1, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 3-dimensional ndarray (column-major, singleton dimensions, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-major';
+	sh = [ 4, 1, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 3-dimensional ndarray (column-major, contiguous)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 3-dimensional ndarray (column-major, contiguous, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-major';
+	sh = [ 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 16.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
diff --git a/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.4d.js b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.4d.js
new file mode 100644
index 000000000000..0c6946933c41
--- /dev/null
+++ b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.4d.js
@@ -0,0 +1,243 @@
+/**
+* @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 zeros = require( '@stdlib/array/zeros' );
+var ones = require( '@stdlib/array/ones' );
+var real = require( '@stdlib/complex/float64/real' );
+var imag = require( '@stdlib/complex/float64/imag' );
+var Complex128Array = require( '@stdlib/array/complex128' );
+var numel = require( '@stdlib/ndarray/base/numel' );
+var shape2strides = require( '@stdlib/ndarray/base/shape2strides' );
+var strides2offset = require( '@stdlib/ndarray/base/strides2offset' );
+var ndarray = require( '@stdlib/ndarray/ctor' );
+var accumulateUnary = require( './../lib' );
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+	t.ok( true, __filename );
+	t.strictEqual( typeof accumulateUnary, 'function', 'main export is a function' );
+	t.end();
+});
+
+tape( 'the function performs accumulation on a 4-dimensional ndarray (row-major, singleton dimensions)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 4, 1, 1, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 4-dimensional ndarray (row-major, singleton dimensions, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 4, 1, 1, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 4-dimensional ndarray (row-major, contiguous)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 16.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 4-dimensional ndarray (row-major, contiguous, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 32.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 4-dimensional ndarray (column-major, contiguous)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 16.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 4-dimensional ndarray (column-major, contiguous, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-major';
+	sh = [ 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 32.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
diff --git a/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.5d.js b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.5d.js
new file mode 100644
index 000000000000..04d7f2284d8b
--- /dev/null
+++ b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.5d.js
@@ -0,0 +1,243 @@
+/**
+* @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 zeros = require( '@stdlib/array/zeros' );
+var ones = require( '@stdlib/array/ones' );
+var real = require( '@stdlib/complex/float64/real' );
+var imag = require( '@stdlib/complex/float64/imag' );
+var Complex128Array = require( '@stdlib/array/complex128' );
+var numel = require( '@stdlib/ndarray/base/numel' );
+var shape2strides = require( '@stdlib/ndarray/base/shape2strides' );
+var strides2offset = require( '@stdlib/ndarray/base/strides2offset' );
+var ndarray = require( '@stdlib/ndarray/ctor' );
+var accumulateUnary = require( './../lib' );
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+	t.ok( true, __filename );
+	t.strictEqual( typeof accumulateUnary, 'function', 'main export is a function' );
+	t.end();
+});
+
+tape( 'the function performs accumulation on a 5-dimensional ndarray (row-major, singleton dimensions)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 4, 1, 1, 1, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 5-dimensional ndarray (row-major, singleton dimensions, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 4, 1, 1, 1, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 5-dimensional ndarray (row-major, contiguous)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 32.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 5-dimensional ndarray (row-major, contiguous, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 64.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 5-dimensional ndarray (column-major, contiguous)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 32.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 5-dimensional ndarray (column-major, contiguous, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-major';
+	sh = [ 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 64.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
diff --git a/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.6d.js b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.6d.js
new file mode 100644
index 000000000000..87afdf46a481
--- /dev/null
+++ b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.6d.js
@@ -0,0 +1,243 @@
+/**
+* @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 zeros = require( '@stdlib/array/zeros' );
+var ones = require( '@stdlib/array/ones' );
+var real = require( '@stdlib/complex/float64/real' );
+var imag = require( '@stdlib/complex/float64/imag' );
+var Complex128Array = require( '@stdlib/array/complex128' );
+var numel = require( '@stdlib/ndarray/base/numel' );
+var shape2strides = require( '@stdlib/ndarray/base/shape2strides' );
+var strides2offset = require( '@stdlib/ndarray/base/strides2offset' );
+var ndarray = require( '@stdlib/ndarray/ctor' );
+var accumulateUnary = require( './../lib' );
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+	t.ok( true, __filename );
+	t.strictEqual( typeof accumulateUnary, 'function', 'main export is a function' );
+	t.end();
+});
+
+tape( 'the function performs accumulation on a 6-dimensional ndarray (row-major, singleton dimensions)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 4, 1, 1, 1, 1, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 6-dimensional ndarray (row-major, singleton dimensions, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 4, 1, 1, 1, 1, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 6-dimensional ndarray (row-major, contiguous)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 64.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 6-dimensional ndarray (row-major, contiguous, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 128.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 6-dimensional ndarray (column-major, contiguous)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 64.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 6-dimensional ndarray (column-major, contiguous, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-major';
+	sh = [ 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 128.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
diff --git a/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.7d.js b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.7d.js
new file mode 100644
index 000000000000..c29a3ba8f480
--- /dev/null
+++ b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.7d.js
@@ -0,0 +1,243 @@
+/**
+* @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 zeros = require( '@stdlib/array/zeros' );
+var ones = require( '@stdlib/array/ones' );
+var real = require( '@stdlib/complex/float64/real' );
+var imag = require( '@stdlib/complex/float64/imag' );
+var Complex128Array = require( '@stdlib/array/complex128' );
+var numel = require( '@stdlib/ndarray/base/numel' );
+var shape2strides = require( '@stdlib/ndarray/base/shape2strides' );
+var strides2offset = require( '@stdlib/ndarray/base/strides2offset' );
+var ndarray = require( '@stdlib/ndarray/ctor' );
+var accumulateUnary = require( './../lib' );
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+	t.ok( true, __filename );
+	t.strictEqual( typeof accumulateUnary, 'function', 'main export is a function' );
+	t.end();
+});
+
+tape( 'the function performs accumulation on a 7-dimensional ndarray (row-major, singleton dimensions)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 4, 1, 1, 1, 1, 1, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 7-dimensional ndarray (row-major, singleton dimensions, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 4, 1, 1, 1, 1, 1, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 7-dimensional ndarray (row-major, contiguous)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 2, 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 128.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 7-dimensional ndarray (row-major, contiguous, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 2, 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 256.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 7-dimensional ndarray (column-major, contiguous)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 2, 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 128.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 7-dimensional ndarray (column-major, contiguous, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-major';
+	sh = [ 2, 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 256.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
diff --git a/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.8d.js b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.8d.js
new file mode 100644
index 000000000000..69cb7474b8d5
--- /dev/null
+++ b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.8d.js
@@ -0,0 +1,243 @@
+/**
+* @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 zeros = require( '@stdlib/array/zeros' );
+var ones = require( '@stdlib/array/ones' );
+var real = require( '@stdlib/complex/float64/real' );
+var imag = require( '@stdlib/complex/float64/imag' );
+var Complex128Array = require( '@stdlib/array/complex128' );
+var numel = require( '@stdlib/ndarray/base/numel' );
+var shape2strides = require( '@stdlib/ndarray/base/shape2strides' );
+var strides2offset = require( '@stdlib/ndarray/base/strides2offset' );
+var ndarray = require( '@stdlib/ndarray/ctor' );
+var accumulateUnary = require( './../lib' );
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+	t.ok( true, __filename );
+	t.strictEqual( typeof accumulateUnary, 'function', 'main export is a function' );
+	t.end();
+});
+
+tape( 'the function performs accumulation on a 8-dimensional ndarray (row-major, singleton dimensions)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 4, 1, 1, 1, 1, 1, 1, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 8-dimensional ndarray (row-major, singleton dimensions, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 4, 1, 1, 1, 1, 1, 1, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 8-dimensional ndarray (row-major, contiguous)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 2, 2, 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 256.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 8-dimensional ndarray (row-major, contiguous, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 2, 2, 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 512.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 8-dimensional ndarray (column-major, contiguous)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 2, 2, 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 256.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 8-dimensional ndarray (column-major, contiguous, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-major';
+	sh = [ 2, 2, 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 512.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
diff --git a/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.9d.js b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.9d.js
new file mode 100644
index 000000000000..dd7f8d30b1cf
--- /dev/null
+++ b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.9d.js
@@ -0,0 +1,243 @@
+/**
+* @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 zeros = require( '@stdlib/array/zeros' );
+var ones = require( '@stdlib/array/ones' );
+var real = require( '@stdlib/complex/float64/real' );
+var imag = require( '@stdlib/complex/float64/imag' );
+var Complex128Array = require( '@stdlib/array/complex128' );
+var numel = require( '@stdlib/ndarray/base/numel' );
+var shape2strides = require( '@stdlib/ndarray/base/shape2strides' );
+var strides2offset = require( '@stdlib/ndarray/base/strides2offset' );
+var ndarray = require( '@stdlib/ndarray/ctor' );
+var accumulateUnary = require( './../lib' );
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+	t.ok( true, __filename );
+	t.strictEqual( typeof accumulateUnary, 'function', 'main export is a function' );
+	t.end();
+});
+
+tape( 'the function performs accumulation on a 9-dimensional ndarray (row-major, singleton dimensions)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 4, 1, 1, 1, 1, 1, 1, 1, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 4.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 9-dimensional ndarray (row-major, singleton dimensions, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 4, 1, 1, 1, 1, 1, 1, 1, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 9-dimensional ndarray (row-major, contiguous)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 2, 2, 2, 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 512.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 9-dimensional ndarray (row-major, contiguous, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 2, 2, 2, 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 1024.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on a 9-dimensional ndarray (column-major, contiguous)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 2, 2, 2, 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 512.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on a 9-dimensional ndarray (column-major, contiguous, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-major';
+	sh = [ 2, 2, 2, 2, 2, 2, 2, 2, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 1024.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
diff --git a/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.js b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.js
index 3fea05a50627..6e1508a997db 100644
--- a/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.js
+++ b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.js
@@ -1,7 +1,7 @@
 /**
 * @license Apache-2.0
 *
-* Copyright (c) 2025 The Stdlib Authors.
+* 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.
@@ -31,5 +31,3 @@ tape( 'main export is a function', function test( t ) {
 	t.strictEqual( typeof accumulateUnary, 'function', 'main export is a function' );
 	t.end();
 });
-
-// FIXME: add tests
diff --git a/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.nd.js b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.nd.js
new file mode 100644
index 000000000000..90cfecde9397
--- /dev/null
+++ b/lib/node_modules/@stdlib/ndarray/base/unary-accumulate/test/test.nd.js
@@ -0,0 +1,176 @@
+/**
+* @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 zeros = require( '@stdlib/array/zeros' );
+var ones = require( '@stdlib/array/ones' );
+var real = require( '@stdlib/complex/float64/real' );
+var imag = require( '@stdlib/complex/float64/imag' );
+var Complex128Array = require( '@stdlib/array/complex128' );
+var numel = require( '@stdlib/ndarray/base/numel' );
+var shape2strides = require( '@stdlib/ndarray/base/shape2strides' );
+var strides2offset = require( '@stdlib/ndarray/base/strides2offset' );
+var ndarray = require( '@stdlib/ndarray/ctor' );
+var accumulateUnary = require( './../lib' );
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+	t.ok( true, __filename );
+	t.strictEqual( typeof accumulateUnary, 'function', 'main export is a function' );
+	t.end();
+});
+
+tape( 'the function performs accumulation on an ndarray with more than 10 dimensions (row-major)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on an ndarray with more than 10 dimensions (row-major, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 16.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});
+
+tape( 'the function performs accumulation on an ndarray with more than 10 dimensions (column-major)', function test( t ) {
+	var actual;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	x = ndarray( dt, ones( numel( sh ), dt ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 8.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + v;
+	}
+});
+
+tape( 'the function performs accumulation on an ndarray with more than 10 dimensions (column-major, accessors)', function test( t ) {
+	var actual;
+	var xbuf;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-major';
+	sh = [ 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	xbuf = zeros( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 0.0, 'returns expected value' );
+
+	xbuf = ones( numel( sh )*2, 'float64' );
+	x = ndarray( dt, new Complex128Array( xbuf ), sh, st, o, ord );
+
+	actual = accumulateUnary( [ x ], 0.0, clbk );
+	t.strictEqual( actual, 16.0, 'returns expected value' );
+
+	t.end();
+
+	function clbk( acc, v ) {
+		return acc + real( v ) + imag( v );
+	}
+});