diff --git a/lib/node_modules/@stdlib/ndarray/base/nullary/test/test.8d.js b/lib/node_modules/@stdlib/ndarray/base/nullary/test/test.8d.js
index 0f51650441d..779263cb6e6 100644
--- a/lib/node_modules/@stdlib/ndarray/base/nullary/test/test.8d.js
+++ b/lib/node_modules/@stdlib/ndarray/base/nullary/test/test.8d.js
@@ -32,6 +32,9 @@ var ndarray = require( '@stdlib/ndarray/ctor' );
 var shape2strides = require( '@stdlib/ndarray/base/shape2strides' );
 var strides2offset = require( '@stdlib/ndarray/base/strides2offset' );
 var numel = require( '@stdlib/ndarray/base/numel' );
+var dfill = require( '@stdlib/blas/ext/base/dfill' );
+var gfill = require( '@stdlib/blas/ext/base/gfill' );
+var blockSize = require( '@stdlib/ndarray/base/nullary-tiling-block-size' );
 var nullary = require( './../lib' );
 
 
@@ -107,6 +110,25 @@ tape( 'the function applies a nullary callback to each indexed element of an 8-d
 	t.end();
 });
 
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, empty)', function test( t ) {
+	var expected;
+	var x;
+
+	x = ndarray( 'float64', zeros( 4, 'float64' ), [ 0, 0, 0, 0, 0, 0, 0, 0 ], [ 1, 1, 1, 1, 1, 1, 1, 1 ], 0, 'row-major' );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array([
+		0.0,
+		0.0,
+		0.0,
+		0.0
+	]);
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
 tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, contiguous)', function test( t ) {
 	var expected;
 	var ord;
@@ -141,6 +163,40 @@ tape( 'the function applies a nullary callback to each indexed element of an 8-d
 	t.end();
 });
 
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, contiguous, negative strides)', function test( t ) {
+	var expected;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+	sh = [ 1, 1, 1, 1, 2, 2, 1, 2 ];
+	st = [ -8, -8, -8, -8, -4, -2, -2, -1 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array([
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0
+	]);
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
 tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, same sign strides)', function test( t ) {
 	var expected;
 	var ord;
@@ -205,27 +261,1242 @@ tape( 'the function applies a nullary callback to each indexed element of an 8-d
 	expected = new Float64Array([
 		10.0,
 		10.0,
-		0.0,
-		0.0,
+		0.0,
+		0.0,
+		10.0,
+		10.0,
+		0.0,
+		0.0,
+		10.0,
+		10.0,
+		0.0,
+		0.0,
+		10.0,
+		10.0,
+		0.0,
+		0.0
+	]);
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, large arrays)', function test( t ) {
+	var expected;
+	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, 2, 1, 1, 2, 1, 2 ];
+	st = [ 32, -16, 8, -8, -8, 4, 4, 2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array( x.length*2 );
+	dfill( x.length, 10.0, expected, st[ 7 ] );
+
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, large arrays)', function test( t ) {
+	var expected;
+	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, 2, 1, 1, 2, 1, 2 ];
+	st = [ -bsize*32, -16, 8, -8, -8, 4, -4, 2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array( x.length*2 );
+	dfill( x.length, 10.0, expected, st[ 7 ] );
+
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, large arrays)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, bsize*2, 2, 1, 2, 1, 2 ];
+	st = [ -bsize*32, -bsize*32, 16, -8, -8, 4, -4, 2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array( x.length*2 );
+	dfill( x.length, 10.0, expected, st[ 7 ] );
+
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, large arrays)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 1, bsize*2, 2, 2, 1, 2 ];
+	st = [ -bsize*32, -bsize*32, bsize*32, -16, -8, 4, -4, 2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array( x.length*2 );
+	dfill( x.length, 10.0, expected, st[ 7 ] );
+
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, large arrays)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 2, 1, bsize*2, 2, 1, 2 ];
+	st = [ -bsize*32, -bsize*32, bsize*16, -bsize*16, -8, 4, -4, 2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array( x.length*2 );
+	dfill( x.length, 10.0, expected, st[ 7 ] );
+
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, large arrays)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 2, 1, 1, bsize*2, 2, 2 ];
+	st = [ -bsize*32, -bsize*32, bsize*16, -bsize*16, -bsize*16, 8, -4, 2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array( x.length*2 );
+	dfill( x.length, 10.0, expected, st[ 7 ] );
+
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, large arrays)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 2, 1, 1, 1, bsize*2, 2 ];
+	st = [ -bsize*16, -bsize*16, bsize*8, -bsize*8, -bsize*8, bsize*8, -4, 2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array( x.length*2 );
+	dfill( x.length, 10.0, expected, st[ 7 ] );
+
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, large arrays)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 2, 1, 1, 2, 1, bsize*2 ];
+	st = [
+		-bsize*16,
+		-bsize*16,
+		bsize*8,
+		-bsize*8,
+		-bsize*8,
+		bsize*4,
+		-bsize*4,
+		2
+	];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array( x.length*2 );
+	dfill( x.length, 10.0, expected, st[ 7 ] );
+
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, contiguous, accessors)', function test( t ) {
+	var expected;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 1, 1, 1, 1, 2, 2, 1, 2 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
+
+	expected = new Complex128Array([
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0
+	]);
+	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, contiguous, accessors, negative strides)', function test( t ) {
+	var expected;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 1, 1, 1, 1, 2, 2, 1, 2 ];
+	st = [ -8, -8, -8, -8, -4, -2, -2, -1 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
+
+	expected = new Complex128Array([
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0
+	]);
+	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, same sign strides, accessors)', function test( t ) {
+	var expected;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 1, 1, 1, 1, 2, 2, 1, 2 ];
+	st = [ 16, 16, 16, 16, 8, 4, 2, 1 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( 16, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
+
+	expected = new Complex128Array([
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		0.0,
+		0.0,
+		0.0,
+		0.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		0.0,
+		0.0,
+		0.0,
+		0.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		0.0,
+		0.0,
+		0.0,
+		0.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		0.0,
+		0.0,
+		0.0,
+		0.0
+	]);
+	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, mixed sign strides, accessors)', function test( t ) {
+	var expected;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+	sh = [ 1, 1, 1, 1, 2, 2, 1, 2 ];
+	st = [ 16, 16, -16, 16, -8, -4, -2, 1 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( 16, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
+
+	expected = new Complex128Array([
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		0.0,
+		0.0,
+		0.0,
+		0.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		0.0,
+		0.0,
+		0.0,
+		0.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		0.0,
+		0.0,
+		0.0,
+		0.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		0.0,
+		0.0,
+		0.0,
+		0.0
+	]);
+	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, large arrays, accessors)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ bsize*2, 2, 2, 1, 1, 2, 1, 2 ];
+	st = [ 32, 16, 8, -8, -8, 4, 4, 2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
+
+	expected = new Complex128Array( x.length*2 );
+	gfill( x.length, new Complex128( 10.0, 10.0 ), expected, st[ 7 ] );
+
+	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, large arrays, accessors)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, bsize*2, 2, 1, 1, 2, 1, 2 ];
+	st = [ -bsize*32, -16, 8, -8, -8, 4, -4, 2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
+
+	expected = new Complex128Array( x.length*2 );
+	gfill( x.length, new Complex128( 10.0, 10.0 ), expected, st[ 7 ] );
+
+	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, large arrays, accessors)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, bsize*2, 2, 1, 2, 1, 2 ];
+	st = [ -bsize*32, -bsize*32, 16, -8, -8, 4, -4, 2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
+
+	expected = new Complex128Array( x.length*2 );
+	gfill( x.length, new Complex128( 10.0, 10.0 ), expected, st[ 7 ] );
+
+	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, large arrays, accessors)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 1, bsize*2, 2, 2, 1, 2 ];
+	st = [ -bsize*32, -bsize*32, bsize*32, -16, -8, 4, -4, 2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
+
+	expected = new Complex128Array( x.length*2 );
+	gfill( x.length, new Complex128( 10.0, 10.0 ), expected, st[ 7 ] );
+
+	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, large arrays, accessors)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 2, 1, bsize*2, 2, 1, 2 ];
+	st = [ -bsize*32, -bsize*32, bsize*16, -bsize*16, -8, 4, -4, 2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
+
+	expected = new Complex128Array( x.length*2 );
+	gfill( x.length, new Complex128( 10.0, 10.0 ), expected, st[ 7 ] );
+
+	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, large arrays, accessors)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 2, 1, 1, bsize*2, 2, 2 ];
+	st = [ -bsize*32, -bsize*32, bsize*16, -bsize*16, -bsize*16, 8, -4, 2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
+
+	expected = new Complex128Array( x.length*2 );
+	gfill( x.length, new Complex128( 10.0, 10.0 ), expected, st[ 7 ] );
+
+	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, large arrays, accessors)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 2, 1, 1, 1, bsize*2, 2 ];
+	st = [ -bsize*16, -bsize*16, bsize*8, -bsize*8, -bsize*8, bsize*8, -4, 2 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
+
+	expected = new Complex128Array( x.length*2 );
+	gfill( x.length, new Complex128( 10.0, 10.0 ), expected, st[ 7 ] );
+
+	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, large arrays, accessors)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'row-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 2, 1, 1, 2, 1, bsize*2 ];
+	st = [
+		-bsize*16,
+		-bsize*16,
+		bsize*8,
+		-bsize*8,
+		-bsize*8,
+		bsize*4,
+		-bsize*4,
+		2
+	];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
+
+	expected = new Complex128Array( x.length*2 );
+	gfill( x.length, new Complex128( 10.0, 10.0 ), expected, st[ 7 ] );
+
+	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, singleton dimensions)', function test( t ) {
+	var expected;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-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 );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array([
+		10.0,
+		10.0,
+		10.0,
+		10.0
+	]);
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, singleton dimensions, accessors)', function test( t ) {
+	var expected;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-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 );
+
+	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
+
+	expected = new Complex128Array([
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0
+	]);
+	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, empty)', function test( t ) {
+	var expected;
+	var x;
+
+	x = ndarray( 'float64', zeros( 4, 'float64' ), [ 0, 0, 0, 0, 0, 0, 0, 0 ], [ 1, 1, 1, 1, 1, 1, 1, 1 ], 0, 'column-major' );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array([
+		0.0,
+		0.0,
+		0.0,
+		0.0
+	]);
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, contiguous)', function test( t ) {
+	var expected;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 1, 1, 1, 2, 2, 1, 2, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array([
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0
+	]);
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, contiguous, negative strides)', function test( t ) {
+	var expected;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 1, 1, 1, 2, 2, 1, 2, 1 ];
+	st = [ -1, -1, -1, -1, -2, -4, -4, -8 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array([
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0
+	]);
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape('the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, same sign strides)', function test( t ) {
+	var expected;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 1, 1, 1, 1, 2, 2, 1, 2 ];
+	st = [ 1, 1, 1, 1, 1, 4, 8, 8 ];
+	o = strides2offset(sh, st);
+
+	x = ndarray( dt, zeros( 16, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array( [
+		10.0,
+		10.0,
+		0.0,
+		0.0,
+		10.0,
+		10.0,
+		0.0,
+		0.0,
+		10.0,
+		10.0,
+		0.0,
+		0.0,
+		10.0,
+		10.0,
+		0.0,
+		0.0
+	]);
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape('the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, mixed sign strides)', function test( t ) {
+	var expected;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+	sh = [ 1, 1, 1, 1, 2, 2, 1, 2 ];
+	st = [ 1, 1, 1, 1, -1, -4, 8, -8 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( 16, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array( [
+		10.0,
+		10.0,
+		0.0,
+		0.0,
+		10.0,
+		10.0,
+		0.0,
+		0.0,
+		10.0,
+		10.0,
+		0.0,
+		0.0,
+		10.0,
+		10.0,
+		0.0,
+		0.0
+	]);
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, large arrays)', function test( t ) {
+	var expected;
+	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, 1, 2, 1, 1, 2, 1, 2 ];
+	st = [ 2, bsize*4, bsize*4, bsize*8, bsize*8, bsize*8, bsize*16, bsize*16 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array( x.length*2 );
+	dfill( x.length, 10.0, expected, st[ 0 ] );
+
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, large arrays)', function test( t ) {
+	var expected;
+	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, 1, 1, 1, 2, 1, 2 ];
+	st = [ 2, -4, bsize*8, bsize*8, bsize*8, bsize*8, bsize*16, bsize*16 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array( x.length*2 );
+	dfill( x.length, 10.0, expected, st[ 0 ] );
+
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, large arrays)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, 2, bsize*2, 1, 1, 2, 1, 2 ];
+	st = [ 2, 4, 8, bsize*16, bsize*16, bsize*16, bsize*32, bsize*32 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array( x.length*2 );
+	dfill( x.length, 10.0, expected, st[ 0 ] );
+
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, large arrays)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 2, bsize*2, 1, 2, 1, 2 ];
+	st = [ 2, 4, 4, 8, bsize*16, bsize*16, bsize*32, bsize*32 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array( x.length*2 );
+	dfill( x.length, 10.0, expected, st[ 0 ] );
+
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, large arrays)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 2, 2, bsize*2, 1, 1, 2 ];
+	st = [ 2, 4, 4, 8, 16, bsize*32, bsize*32, bsize*32 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array( x.length*2 );
+	dfill( x.length, 10.0, expected, st[ 0 ] );
+
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, large arrays)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 2, 1, 2, bsize*2, 1, 2 ];
+	st = [ 2, 4, 4, 8, 8, 16, bsize*32, bsize*32 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array( x.length*2 );
+	dfill( x.length, 10.0, expected, st[ 0 ] );
+
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, large arrays)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 2, 1, 1, 2, bsize*2, 2 ];
+	st = [ 2, 4, 4, 8, 8, 8, 16, bsize*32 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array( x.length*2 );
+	dfill( x.length, 10.0, expected, st[ 0 ] );
+
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, large arrays)', function test( t ) {
+	var expected;
+	var bsize;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'float64';
+	ord = 'column-major';
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 2, 1, 1, 2, 2, bsize*2 ];
+	st = [ 2, 4, 4, 8, 8, 8, 16, 32 ];
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( 10.0 ) );
+
+	expected = new Float64Array( x.length*2 );
+	dfill( x.length, 10.0, expected, st[ 0 ] );
+
+	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, contiguous, accessors)', function test( t ) {
+	var expected;
+	var ord;
+	var sh;
+	var st;
+	var dt;
+	var o;
+	var x;
+
+	dt = 'complex128';
+	ord = 'column-major';
+	sh = [ 1, 1, 1, 2, 2, 1, 2, 1 ];
+	st = shape2strides( sh, ord );
+	o = strides2offset( sh, st );
+
+	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+
+	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
+
+	expected = new Complex128Array([
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
+		10.0,
 		10.0,
 		10.0,
-		0.0,
-		0.0,
 		10.0,
 		10.0,
-		0.0,
-		0.0,
 		10.0,
 		10.0,
-		0.0,
-		0.0
+		10.0,
+		10.0,
+		10.0,
+		10.0
 	]);
-	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, contiguous, accessors)', function test( t ) {
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, contiguous, accessors, negative strides)', function test( t ) {
 	var expected;
 	var ord;
 	var sh;
@@ -235,9 +1506,9 @@ tape( 'the function applies a nullary callback to each indexed element of an 8-d
 	var x;
 
 	dt = 'complex128';
-	ord = 'row-major';
-	sh = [ 1, 1, 1, 1, 2, 2, 1, 2 ];
-	st = shape2strides( sh, ord );
+	ord = 'column-major';
+	sh = [ 1, 1, 1, 2, 2, 1, 2, 1 ];
+	st = [ -1, -1, -1, -1, -2, -4, -4, -8 ];
 	o = strides2offset( sh, st );
 
 	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
@@ -267,7 +1538,7 @@ tape( 'the function applies a nullary callback to each indexed element of an 8-d
 	t.end();
 });
 
-tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, same sign strides, accessors)', function test( t ) {
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, same sign strides, accessors)', function test( t ) {
 	var expected;
 	var ord;
 	var sh;
@@ -277,16 +1548,16 @@ tape( 'the function applies a nullary callback to each indexed element of an 8-d
 	var x;
 
 	dt = 'complex128';
-	ord = 'row-major';
+	ord = 'column-major';
 	sh = [ 1, 1, 1, 1, 2, 2, 1, 2 ];
-	st = [ 16, 16, 16, 16, 8, 4, 2, 1 ];
+	st = [ 1, 1, 1, 1, 1, 4, 8, 8 ];
 	o = strides2offset( sh, st );
 
 	x = ndarray( dt, zeros( 16, dt ), sh, st, o, ord );
 
 	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
 
-	expected = new Complex128Array([
+	expected = new Complex128Array( [
 		10.0,
 		10.0,
 		10.0,
@@ -325,7 +1596,7 @@ tape( 'the function applies a nullary callback to each indexed element of an 8-d
 	t.end();
 });
 
-tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (row-major, non-contiguous, mixed sign strides, accessors)', function test( t ) {
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, mixed sign strides, accessors)', function test( t ) {
 	var expected;
 	var ord;
 	var sh;
@@ -335,9 +1606,9 @@ tape( 'the function applies a nullary callback to each indexed element of an 8-d
 	var x;
 
 	dt = 'complex128';
-	ord = 'row-major';
+	ord = 'column-major';
 	sh = [ 1, 1, 1, 1, 2, 2, 1, 2 ];
-	st = [ 16, 16, -16, 16, -8, -4, -2, 1 ];
+	st = [ -1, 1, -1, 1, 1, -4, -8, 8 ];
 	o = strides2offset( sh, st );
 
 	x = ndarray( dt, zeros( 16, dt ), sh, st, o, ord );
@@ -383,8 +1654,9 @@ tape( 'the function applies a nullary callback to each indexed element of an 8-d
 	t.end();
 });
 
-tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, singleton dimensions)', function test( t ) {
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, large arrays, accessors)', function test( t ) {
 	var expected;
+	var bsize;
 	var ord;
 	var sh;
 	var st;
@@ -392,29 +1664,29 @@ tape( 'the function applies a nullary callback to each indexed element of an 8-d
 	var o;
 	var x;
 
-	dt = 'float64';
+	dt = 'complex128';
 	ord = 'column-major';
-	sh = [ 4, 1, 1, 1, 1, 1, 1, 1 ];
-	st = shape2strides( sh, ord );
+
+	bsize = blockSize( dt );
+	sh = [ bsize*2, 1, 2, 1, 1, 2, 1, 2 ];
+	st = [ 2, bsize*4, bsize*4, bsize*8, bsize*8, bsize*8, bsize*16, bsize*16 ];
 	o = strides2offset( sh, st );
 
-	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
 
-	nullary( [ x ], constantFunction( 10.0 ) );
+	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
 
-	expected = new Float64Array([
-		10.0,
-		10.0,
-		10.0,
-		10.0
-	]);
-	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+	expected = new Complex128Array( x.length*2 );
+	gfill( x.length, new Complex128( 10.0, 10.0 ), expected, st[ 0 ] );
+
+	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, singleton dimensions, accessors)', function test( t ) {
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, large arrays, accessors)', function test( t ) {
 	var expected;
+	var bsize;
 	var ord;
 	var sh;
 	var st;
@@ -424,31 +1696,27 @@ tape( 'the function applies a nullary callback to each indexed element of an 8-d
 
 	dt = 'complex128';
 	ord = 'column-major';
-	sh = [ 4, 1, 1, 1, 1, 1, 1, 1 ];
-	st = shape2strides( sh, ord );
+
+	bsize = blockSize( dt );
+	sh = [ 2, bsize*2, 1, 1, 1, 2, 1, 2 ];
+	st = [ 2, -4, bsize*8, bsize*8, bsize*8, bsize*8, bsize*16, bsize*16 ];
 	o = strides2offset( sh, st );
 
-	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
 
 	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
 
-	expected = new Complex128Array([
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0
-	]);
+	expected = new Complex128Array( x.length*2 );
+	gfill( x.length, new Complex128( 10.0, 10.0 ), expected, st[ 0 ] );
+
 	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, contiguous)', function test( t ) {
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, large arrays, accessors)', function test( t ) {
 	var expected;
+	var bsize;
 	var ord;
 	var sh;
 	var st;
@@ -456,33 +1724,29 @@ tape( 'the function applies a nullary callback to each indexed element of an 8-d
 	var o;
 	var x;
 
-	dt = 'float64';
+	dt = 'complex128';
 	ord = 'column-major';
-	sh = [ 1, 1, 1, 2, 2, 1, 2, 1 ];
-	st = shape2strides( sh, ord );
+
+	bsize = blockSize( dt );
+	sh = [ 2, 2, bsize*2, 1, 1, 2, 1, 2 ];
+	st = [ 2, 4, 8, bsize*16, bsize*16, bsize*16, bsize*32, bsize*32 ];
 	o = strides2offset( sh, st );
 
-	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
 
-	nullary( [ x ], constantFunction( 10.0 ) );
+	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
 
-	expected = new Float64Array([
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0
-	]);
-	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+	expected = new Complex128Array( x.length*2 );
+	gfill( x.length, new Complex128( 10.0, 10.0 ), expected, st[ 0 ] );
+
+	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
 
 	t.end();
 });
 
-tape('the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, same sign strides)', function test( t ) {
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, large arrays, accessors)', function test( t ) {
 	var expected;
+	var bsize;
 	var ord;
 	var sh;
 	var st;
@@ -490,41 +1754,29 @@ tape('the function applies a nullary callback to each indexed element of an 8-di
 	var o;
 	var x;
 
-	dt = 'float64';
+	dt = 'complex128';
 	ord = 'column-major';
-	sh = [ 1, 1, 1, 1, 2, 2, 1, 2 ];
-	st = [ 1, 1, 1, 1, 1, 4, 8, 8 ];
-	o = strides2offset(sh, st);
 
-	x = ndarray( dt, zeros( 16, dt ), sh, st, o, ord );
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 2, bsize*2, 1, 2, 1, 2 ];
+	st = [ 2, 4, 4, 8, bsize*16, bsize*16, bsize*32, bsize*32 ];
+	o = strides2offset( sh, st );
 
-	nullary( [ x ], constantFunction( 10.0 ) );
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
 
-	expected = new Float64Array( [
-		10.0,
-		10.0,
-		0.0,
-		0.0,
-		10.0,
-		10.0,
-		0.0,
-		0.0,
-		10.0,
-		10.0,
-		0.0,
-		0.0,
-		10.0,
-		10.0,
-		0.0,
-		0.0
-	]);
-	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
+
+	expected = new Complex128Array( x.length*2 );
+	gfill( x.length, new Complex128( 10.0, 10.0 ), expected, st[ 0 ] );
+
+	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
 
 	t.end();
 });
 
-tape('the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, mixed sign strides)', function test( t ) {
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, large arrays, accessors)', function test( t ) {
 	var expected;
+	var bsize;
 	var ord;
 	var sh;
 	var st;
@@ -532,41 +1784,29 @@ tape('the function applies a nullary callback to each indexed element of an 8-di
 	var o;
 	var x;
 
-	dt = 'float64';
+	dt = 'complex128';
 	ord = 'column-major';
-	sh = [ 1, 1, 1, 1, 2, 2, 1, 2 ];
-	st = [ 1, 1, 1, 1, -1, -4, 8, -8 ];
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 2, 2, bsize*2, 1, 1, 2 ];
+	st = [ 2, 4, 4, 8, 16, bsize*32, bsize*32, bsize*32 ];
 	o = strides2offset( sh, st );
 
-	x = ndarray( dt, zeros( 16, dt ), sh, st, o, ord );
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
 
-	nullary( [ x ], constantFunction( 10.0 ) );
+	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
 
-	expected = new Float64Array( [
-		10.0,
-		10.0,
-		0.0,
-		0.0,
-		10.0,
-		10.0,
-		0.0,
-		0.0,
-		10.0,
-		10.0,
-		0.0,
-		0.0,
-		10.0,
-		10.0,
-		0.0,
-		0.0
-	]);
-	t.strictEqual( isSameFloat64Array( x.data, expected ), true, 'returns expected value' );
+	expected = new Complex128Array( x.length*2 );
+	gfill( x.length, new Complex128( 10.0, 10.0 ), expected, st[ 0 ] );
+
+	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, contiguous, accessors)', function test( t ) {
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, large arrays, accessors)', function test( t ) {
 	var expected;
+	var bsize;
 	var ord;
 	var sh;
 	var st;
@@ -576,39 +1816,27 @@ tape( 'the function applies a nullary callback to each indexed element of an 8-d
 
 	dt = 'complex128';
 	ord = 'column-major';
-	sh = [ 1, 1, 1, 2, 2, 1, 2, 1 ];
-	st = shape2strides( sh, ord );
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 2, 1, 2, bsize*2, 1, 2 ];
+	st = [ 2, 4, 4, 8, 8, 16, bsize*32, bsize*32 ];
 	o = strides2offset( sh, st );
 
-	x = ndarray( dt, zeros( numel( sh ), dt ), sh, st, o, ord );
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
 
 	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
 
-	expected = new Complex128Array([
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0
-	]);
+	expected = new Complex128Array( x.length*2 );
+	gfill( x.length, new Complex128( 10.0, 10.0 ), expected, st[ 0 ] );
+
 	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, same sign strides, accessors)', function test( t ) {
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, large arrays, accessors)', function test( t ) {
 	var expected;
+	var bsize;
 	var ord;
 	var sh;
 	var st;
@@ -618,55 +1846,27 @@ tape( 'the function applies a nullary callback to each indexed element of an 8-d
 
 	dt = 'complex128';
 	ord = 'column-major';
-	sh = [ 1, 1, 1, 1, 2, 2, 1, 2 ];
-	st = [ 1, 1, 1, 1, 1, 4, 8, 8 ];
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 2, 1, 1, 2, bsize*2, 2 ];
+	st = [ 2, 4, 4, 8, 8, 8, 16, bsize*32 ];
 	o = strides2offset( sh, st );
 
-	x = ndarray( dt, zeros( 16, dt ), sh, st, o, ord );
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
 
 	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
 
-	expected = new Complex128Array( [
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		0.0,
-		0.0,
-		0.0,
-		0.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		0.0,
-		0.0,
-		0.0,
-		0.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		0.0,
-		0.0,
-		0.0,
-		0.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		0.0,
-		0.0,
-		0.0,
-		0.0
-	]);
+	expected = new Complex128Array( x.length*2 );
+	gfill( x.length, new Complex128( 10.0, 10.0 ), expected, st[ 0 ] );
+
 	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, mixed sign strides, accessors)', function test( t ) {
+tape( 'the function applies a nullary callback to each indexed element of an 8-dimensional ndarray (column-major, non-contiguous, large arrays, accessors)', function test( t ) {
 	var expected;
+	var bsize;
 	var ord;
 	var sh;
 	var st;
@@ -676,48 +1876,19 @@ tape( 'the function applies a nullary callback to each indexed element of an 8-d
 
 	dt = 'complex128';
 	ord = 'column-major';
-	sh = [ 1, 1, 1, 1, 2, 2, 1, 2 ];
-	st = [ -1, 1, -1, 1, 1, -4, -8, 8 ];
+
+	bsize = blockSize( dt );
+	sh = [ 2, 1, 2, 1, 1, 2, 2, bsize*2 ];
+	st = [ 2, 4, 4, 8, 8, 8, 16, 32 ];
 	o = strides2offset( sh, st );
 
-	x = ndarray( dt, zeros( 16, dt ), sh, st, o, ord );
+	x = ndarray( dt, zeros( numel( sh )*2, dt ), sh, st, o, ord );
 
 	nullary( [ x ], constantFunction( new Complex128( 10.0, 10.0 ) ) );
 
-	expected = new Complex128Array([
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		0.0,
-		0.0,
-		0.0,
-		0.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		0.0,
-		0.0,
-		0.0,
-		0.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		0.0,
-		0.0,
-		0.0,
-		0.0,
-		10.0,
-		10.0,
-		10.0,
-		10.0,
-		0.0,
-		0.0,
-		0.0,
-		0.0
-	]);
+	expected = new Complex128Array( x.length*2 );
+	gfill( x.length, new Complex128( 10.0, 10.0 ), expected, st[ 0 ] );
+
 	t.strictEqual( isSameComplex128Array( x.data, expected ), true, 'returns expected value' );
 
 	t.end();