test.ndarray.js

/**
* @license Apache-2.0
*
* Copyright (c) 2025 The Stdlib Authors.
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
*    http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/

/* eslint-disable max-lines, max-len */

'use strict';

// MODULES //

var tape = require( 'tape' );
var Float64Array = require( '@stdlib/array/float64' );
var dgebak = require( './../lib/ndarray.js' );


// TESTS //

tape( 'main export is a function', function test( t ) {
	t.ok( true, __filename );
	t.strictEqual( typeof dgebak, 'function', 'main export is a function' );
	t.end();
});

tape( 'the function has an arity of 13', function test( t ) {
	t.strictEqual( dgebak.length, 13, 'returns expected value' );
	t.end();
});

tape( 'the function throws an error if provided a first argument which is not a valid order', function test( t ) {
	var values;
	var scale;
	var i;
	var V;

	values = [
		'foo',
		'bar',
		'beep',
		'boop',
		-5,
		NaN,
		true,
		false,
		null,
		void 0,
		[],
		{},
		function noop() {}
	];

	V = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
	scale = new Float64Array( [ 1.0, 2.0, 3.0 ] );

	for ( i = 0; i < values.length; i++ ) {
		t.throws( badValue( values[ i ] ), TypeError, 'throws an error when provided ' + values[ i ] );
	}
	t.end();

	function badValue( value ) {
		return function badValue() {
			dgebak( value, 'scale', 'right', 3, 2, 0, 2, scale, V, 2 );
		};
	}
});

tape( 'the function throws an error if provided a second argument which is not a valid job', function test( t ) {
	var values;
	var scale;
	var i;
	var V;

	values = [
		'foo',
		'bar',
		'beep',
		'boop',
		-5,
		NaN,
		true,
		false,
		null,
		void 0,
		[],
		{},
		function noop() {}
	];

	V = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
	scale = new Float64Array( [ 1.0, 2.0, 3.0 ] );

	for ( i = 0; i < values.length; i++ ) {
		t.throws( badValue( values[ i ] ), TypeError, 'throws an error when provided ' + values[ i ] );
	}
	t.end();

	function badValue( value ) {
		return function badValue() {
			dgebak( value, 'right', 3, 2, 0, 2, scale, 1, 0, V, 2, 1, 0 );
		};
	}
});

tape( 'the function throws an error if provided a third argument which is not a valid side', function test( t ) {
	var values;
	var scale;
	var i;
	var V;

	values = [
		'foo',
		'bar',
		'beep',
		'boop',
		-5,
		NaN,
		true,
		false,
		null,
		void 0,
		[],
		{},
		function noop() {}
	];

	V = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
	scale = new Float64Array( [ 1.0, 2.0, 3.0 ] );

	for ( i = 0; i < values.length; i++ ) {
		t.throws( badValue( values[ i ] ), TypeError, 'throws an error when provided ' + values[ i ] );
	}
	t.end();

	function badValue( value ) {
		return function badValue() {
			dgebak( 'scale', value, 3, 2, 0, 2, scale, 1, 0, V, 2, 1, 0 );
		};
	}
});

tape( 'the function returns the array unchanged when N = 0', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 1.0, 2.0, 3.0 ] );
	V = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );

	expected = new Float64Array( V );

	actual = dgebak( 'scale', 'right', 0, 2, 0, 2, scale, 1, 0, V, 2, 1, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns the array unchanged when M = 0', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 1.0, 2.0, 3.0 ] );
	V = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );

	expected = new Float64Array( V );

	actual = dgebak( 'both', 'right', 3, 0, 0, 2, scale, 1, 0, V, 2, 1, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns the array unchanged when job = none (row-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 1.0, 2.0, 3.0 ] );
	V = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );

	expected = new Float64Array( V );

	actual = dgebak( 'none', 'right', 3, 2, 0, 2, scale, 1, 0, V, 2, 1, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns the array unchanged when job = none (column-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 1.0, 2.0, 3.0 ] );
	V = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );

	expected = new Float64Array( V );

	actual = dgebak( 'none', 'right', 3, 2, 0, 2, scale, 1, 0, V, 1, 3, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = permute (row-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 2.0, 0.0 ] );
	V = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 1.0, 2.0, 5.0, 6.0, 3.0, 4.0, 7.0, 8.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 2, 2, scale, 1, 0, V, 2, 1, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = permute (column-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 2.0, 0.0 ] );
	V = new Float64Array( [ 1.0, 3.0, 5.0, 7.0, 2.0, 4.0, 6.0, 8.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 1.0, 5.0, 3.0, 7.0, 2.0, 6.0, 4.0, 8.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 2, 2, scale, 1, 0, V, 1, 4, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = both (row-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 2.0, 0.0 ] );
	V = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 1.0, 2.0, 5.0, 6.0, 3.0, 4.0, 7.0, 8.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 2, 2, scale, 1, 0, V, 2, 1, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = both (column-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 2.0, 0.0 ] );
	V = new Float64Array( [ 1.0, 3.0, 5.0, 7.0, 2.0, 4.0, 6.0, 8.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 1.0, 5.0, 3.0, 7.0, 2.0, 6.0, 4.0, 8.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 2, 2, scale, 1, 0, V, 1, 4, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (right side) (row-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 10.0, 1.0, 0.1, 0.1 ] );
	V = new Float64Array( [ 10.0, 20.0, 3.0, 4.0, 0.5, 0.6, 0.04, 0.08 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 100.0, 200.0, 3.0, 4.0, 0.05, 0.06, 0.004, 0.008 ] );

	/*
		expected = [
			[ 100.0, 200.0 ],
			[ 3.0, 4.0 ],
			[ 0.05, 0.06 ],
			[ 0.004, 0.008 ]
		]
	*/

	actual = dgebak( 'scale', 'right', 4, 2, 0, 3, scale, 1, 0, V, 2, 1, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (right side) (column-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 10.0, 1.0, 0.1, 0.1 ] );
	V = new Float64Array( [ 10.0, 3.0, 0.5, 0.04, 20.0, 4.0, 0.6, 0.08 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 100.0, 3.0, 0.05, 0.004, 200.0, 4.0, 0.06, 0.008 ] );

	/*
		expected = [
			[ 100.0, 200.0 ],
			[ 3.0, 4.0 ],
			[ 0.05, 0.06 ],
			[ 0.004, 0.008 ]
		]
	*/

	actual = dgebak( 'scale', 'right', 4, 2, 0, 3, scale, 1, 0, V, 1, 4, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (left side) (row-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 10.0, 1.0, 0.1, 0.01 ] );
	V = new Float64Array( [ 10.0, 20.0, 3.0, 4.0, 0.5, 0.6, 0.04, 0.08 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 4.0, 8.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 4.0, 8.0 ]
		]
	*/

	actual = dgebak( 'scale', 'left', 4, 2, 0, 3, scale, 1, 0, V, 2, 1, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (left side) (column-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 10.0, 1.0, 0.1, 0.01 ] );
	V = new Float64Array( [ 10.0, 3.0, 0.5, 0.04, 20.0, 4.0, 0.6, 0.08 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 1.0, 3.0, 5.0, 4.0, 2.0, 4.0, 6.0, 8.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 4.0, 8.0 ]
		]
	*/

	actual = dgebak( 'scale', 'left', 4, 2, 0, 3, scale, 1, 0, V, 1, 4, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (right side) (row-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 1.0, 0.0 ] );
	V = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 7.0, 8.0, 5.0, 6.0, 3.0, 4.0, 1.0, 2.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 0, 1, scale, 1, 0, V, 2, 1, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (right side) (column-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 1.0, 0.0 ] );
	V = new Float64Array( [ 1.0, 3.0, 5.0, 7.0, 2.0, 4.0, 6.0, 8.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 7.0, 5.0, 3.0, 1.0, 8.0, 6.0, 4.0, 2.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 0, 1, scale, 1, 0, V, 1, 4, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (left side) (row-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 1.0, 0.0 ] );
	V = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 7.0, 8.0, 5.0, 6.0, 3.0, 4.0, 1.0, 2.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'left', 4, 2, 0, 1, scale, 1, 0, V, 2, 1, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (left side) (column-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 1.0, 0.0 ] );
	V = new Float64Array( [ 1.0, 3.0, 5.0, 7.0, 2.0, 4.0, 6.0, 8.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 7.0, 5.0, 3.0, 1.0, 8.0, 6.0, 4.0, 2.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'left', 4, 2, 0, 1, scale, 1, 0, V, 1, 4, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (right side) (row-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 10.0, 0.1, 0.0 ] );
	V = new Float64Array( [ 7.0, 8.0, 3.0, 4.0, 0.5, 0.6, 1.0, 2.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 7.0, 8.0, 30.0, 40.0, 0.05, 0.06, 1.0, 2.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 30.0, 40.0 ],
			[ 0.05, 0.06 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 1, 2, scale, 1, 0, V, 2, 1, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (right side) (column-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 10.0, 0.1, 0.0 ] );
	V = new Float64Array( [ 7.0, 3.0, 0.5, 1.0, 8.0, 4.0, 0.6, 2.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 7.0, 30.0, 0.05, 1.0, 8.0, 40.0, 0.06, 2.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 30.0, 40.0 ],
			[ 0.05, 0.06 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 1, 2, scale, 1, 0, V, 1, 4, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (left side) (row-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 10.0, 0.1, 0.0 ] );
	V = new Float64Array( [ 7.0, 8.0, 8.0, 4.0, 0.5, 0.6, 1.0, 2.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 8.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 7.0, 8.0, 0.8, 0.4, 5.0, 6.0, 1.0, 2.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 0.8, 0.4 ],
			[ 5.0, 6.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'left', 4, 2, 1, 2, scale, 1, 0, V, 2, 1, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (left side) (column-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 10.0, 0.1, 0.0 ] );
	V = new Float64Array( [ 7.0, 8.0, 0.5, 1.0, 8.0, 4.0, 0.6, 2.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 8.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 7.0, 0.8, 5.0, 1.0, 8.0, 0.4, 6.0, 2.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 0.8, 0.4 ],
			[ 5.0, 6.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'left', 4, 2, 1, 2, scale, 1, 0, V, 1, 4, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns the array unchanged when job = none (row-major) (offsets)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.0, 1.0, 2.0, 3.0 ] );
	V = new Float64Array( [ 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );

	expected = new Float64Array( V );

	actual = dgebak( 'none', 'right', 3, 2, 0, 2, scale, 1, 2, V, 2, 1, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns the array unchanged when job = none (column-major) (offsets)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.0, 1.0, 2.0, 3.0 ] );
	V = new Float64Array( [ 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );

	expected = new Float64Array( V );

	actual = dgebak( 'none', 'right', 3, 2, 0, 2, scale, 1, 2, V, 1, 3, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = permute (row-major) (offsets)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.0, 3.0, 2.0, 2.0, 0.0 ] );
	V = new Float64Array( [ 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 0.0, 1.0, 2.0, 5.0, 6.0, 3.0, 4.0, 7.0, 8.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 2, 2, scale, 1, 2, V, 2, 1, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = permute (column-major) (offsets)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.0, 3.0, 2.0, 2.0, 0.0 ] );
	V = new Float64Array( [ 0.0, 1.0, 3.0, 5.0, 7.0, 2.0, 4.0, 6.0, 8.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 0.0, 1.0, 5.0, 3.0, 7.0, 2.0, 6.0, 4.0, 8.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 2, 2, scale, 1, 2, V, 1, 4, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = both (row-major) (offsets)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.0, 3.0, 2.0, 2.0, 0.0 ] );
	V = new Float64Array( [ 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 0.0, 1.0, 2.0, 5.0, 6.0, 3.0, 4.0, 7.0, 8.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 2, 2, scale, 1, 2, V, 2, 1, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = both (column-major) (offsets)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.0, 3.0, 2.0, 2.0, 0.0 ] );
	V = new Float64Array( [ 0.0, 1.0, 3.0, 5.0, 7.0, 2.0, 4.0, 6.0, 8.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 0.0, 1.0, 5.0, 3.0, 7.0, 2.0, 6.0, 4.0, 8.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 2, 2, scale, 1, 2, V, 1, 4, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (right side) (row-major) (offsets)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.0, 10.0, 1.0, 0.1, 0.1 ] );
	V = new Float64Array( [ 0.0, 10.0, 20.0, 3.0, 4.0, 0.5, 0.6, 0.04, 0.08 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 0.0, 100.0, 200.0, 3.0, 4.0, 0.05, 0.06, 0.004, 0.008 ] );

	/*
		expected = [
			[ 100.0, 200.0 ],
			[ 3.0, 4.0 ],
			[ 0.05, 0.06 ],
			[ 0.004, 0.008 ]
		]
	*/

	actual = dgebak( 'scale', 'right', 4, 2, 0, 3, scale, 1, 2, V, 2, 1, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (right side) (column-major) (offsets)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.0, 10.0, 1.0, 0.1, 0.1, 0.0, 0.0 ] );
	V = new Float64Array( [ 0.0, 10.0, 3.0, 0.5, 0.04, 20.0, 4.0, 0.6, 0.08 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 0.0, 100.0, 3.0, 0.05, 0.004, 200.0, 4.0, 0.06, 0.008 ] );

	/*
		expected = [
			[ 100.0, 200.0 ],
			[ 3.0, 4.0 ],
			[ 0.05, 0.06 ],
			[ 0.004, 0.008 ]
		]
	*/

	actual = dgebak( 'scale', 'right', 4, 2, 0, 3, scale, 1, 2, V, 1, 4, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (left side) (row-major) (offsets)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.0, 10.0, 1.0, 0.1, 0.01, 0.0, 0.0 ] );
	V = new Float64Array( [ 0.0, 10.0, 20.0, 3.0, 4.0, 0.5, 0.6, 0.04, 0.08 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 4.0, 8.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 4.0, 8.0 ]
		]
	*/

	actual = dgebak( 'scale', 'left', 4, 2, 0, 3, scale, 1, 2, V, 2, 1, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (left side) (column-major) (offsets)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.0, 10.0, 1.0, 0.1, 0.01, 0.0, 0.0 ] );
	V = new Float64Array( [ 0.0, 10.0, 3.0, 0.5, 0.04, 20.0, 4.0, 0.6, 0.08 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 0.0, 1.0, 3.0, 5.0, 4.0, 2.0, 4.0, 6.0, 8.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 4.0, 8.0 ]
		]
	*/

	actual = dgebak( 'scale', 'left', 4, 2, 0, 3, scale, 1, 2, V, 1, 4, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (right side) (row-major) (offsets)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.0, 3.0, 2.0, 1.0, 0.0, 0.0, 0.0 ] );
	V = new Float64Array( [ 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 0.0, 7.0, 8.0, 5.0, 6.0, 3.0, 4.0, 1.0, 2.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 0, 1, scale, 1, 2, V, 2, 1, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (right side) (column-major) (offsets)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.0, 3.0, 2.0, 1.0, 0.0, 0.0, 0.0 ] );
	V = new Float64Array( [ 0.0, 1.0, 3.0, 5.0, 7.0, 2.0, 4.0, 6.0, 8.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 0.0, 7.0, 5.0, 3.0, 1.0, 8.0, 6.0, 4.0, 2.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 0, 1, scale, 1, 2, V, 1, 4, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (left side) (row-major) (offsets)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.0, 3.0, 2.0, 1.0, 0.0, 0.0, 0.0 ] );
	V = new Float64Array( [ 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 0.0, 7.0, 8.0, 5.0, 6.0, 3.0, 4.0, 1.0, 2.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'left', 4, 2, 0, 1, scale, 1, 2, V, 2, 1, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (left side) (column-major) (offsets)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.0, 3.0, 2.0, 1.0, 0.0 ] );
	V = new Float64Array( [ 0.0, 1.0, 3.0, 5.0, 7.0, 2.0, 4.0, 6.0, 8.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 0.0, 7.0, 5.0, 3.0, 1.0, 8.0, 6.0, 4.0, 2.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'left', 4, 2, 0, 1, scale, 1, 2, V, 1, 4, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (right side) (row-major) (offsets)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.0, 3.0, 10.0, 0.1, 0.0 ] );
	V = new Float64Array( [ 0.0, 7.0, 8.0, 3.0, 4.0, 0.5, 0.6, 1.0, 2.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 0.0, 7.0, 8.0, 30.0, 40.0, 0.05, 0.06, 1.0, 2.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 30.0, 40.0 ],
			[ 0.05, 0.06 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 1, 2, scale, 1, 2, V, 2, 1, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (right side) (column-major) (offsets)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.0, 3.0, 10.0, 0.1, 0.0 ] );
	V = new Float64Array( [ 0.0, 7.0, 3.0, 0.5, 1.0, 8.0, 4.0, 0.6, 2.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 0.0, 7.0, 30.0, 0.05, 1.0, 8.0, 40.0, 0.06, 2.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 30.0, 40.0 ],
			[ 0.05, 0.06 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 1, 2, scale, 1, 2, V, 1, 4, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (left side) (row-major) (offsets)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.0, 3.0, 10.0, 0.1, 0.0 ] );
	V = new Float64Array( [ 0.0, 7.0, 8.0, 8.0, 4.0, 0.5, 0.6, 1.0, 2.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 8.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 0.0, 7.0, 8.0, 0.8, 0.4, 5.0, 6.0, 1.0, 2.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 0.8, 0.4 ],
			[ 5.0, 6.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'left', 4, 2, 1, 2, scale, 1, 2, V, 2, 1, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (left side) (column-major) (offsets)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.0, 3.0, 10.0, 0.1, 0.0 ] );
	V = new Float64Array( [ 0.0, 7.0, 8.0, 0.5, 1.0, 8.0, 4.0, 0.6, 2.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 8.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 0.0, 7.0, 0.8, 5.0, 1.0, 8.0, 0.4, 6.0, 2.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 0.8, 0.4 ],
			[ 5.0, 6.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'left', 4, 2, 1, 2, scale, 1, 2, V, 1, 4, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns the array unchanged when job = none (row-major) (negative strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 1.0 ] );
	V = new Float64Array( [ 6.0, 5.0, 4.0, 3.0, 2.0, 1.0 ] );

	expected = new Float64Array( V );

	actual = dgebak( 'none', 'right', 3, 2, 0, 2, scale, -1, 2, V, -2, -1, 5 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns the array unchanged when job = none (column-major) (negative strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 1.0 ] );
	V = new Float64Array( [ 6.0, 5.0, 4.0, 3.0, 2.0, 1.0 ] );

	expected = new Float64Array( V );

	actual = dgebak( 'none', 'right', 3, 2, 0, 2, scale, -1, 2, V, -1, -3, 5 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = permute (row-major) (negative strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 2.0, 2.0, 3.0 ] );
	V = new Float64Array( [ 8.0, 7.0, 6.0, 5.0, 4.0, 3.0, 2.0, 1.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 8.0, 7.0, 4.0, 3.0, 6.0, 5.0, 2.0, 1.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 2, 2, scale, -1, 3, V, -2, -1, 7 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = permute (column-major) (negative strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 2.0, 2.0, 3.0 ] );
	V = new Float64Array( [ 8.0, 6.0, 4.0, 2.0, 7.0, 5.0, 3.0, 1.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 8.0, 4.0, 6.0, 2.0, 7.0, 3.0, 5.0, 1.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 2, 2, scale, -1, 3, V, -1, -4, 7 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = both (row-major) (negative strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 2.0, 2.0, 3.0 ] );
	V = new Float64Array( [ 8.0, 7.0, 6.0, 5.0, 4.0, 3.0, 2.0, 1.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 8.0, 7.0, 4.0, 3.0, 6.0, 5.0, 2.0, 1.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 2, 2, scale, -1, 3, V, -2, -1, 7 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = both (column-major) (negative strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 2.0, 2.0, 3.0 ] );
	V = new Float64Array( [ 8.0, 6.0, 4.0, 2.0, 7.0, 5.0, 3.0, 1.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 8.0, 4.0, 6.0, 2.0, 7.0, 3.0, 5.0, 1.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 2, 2, scale, -1, 3, V, -1, -4, 7 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (right side) (row-major) (negative strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.1, 0.1, 1.0, 10.0 ] );
	V = new Float64Array( [ 0.04, 0.08, 0.5, 0.6, 3.0, 4.0, 10.0, 20.0 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 0.004, 0.008, 0.05, 0.06, 3.0, 4.0, 100.0, 200.0 ] );

	/*
		expected = [
			[ 100.0, 200.0 ],
			[ 3.0, 4.0 ],
			[ 0.05, 0.06 ],
			[ 0.004, 0.008 ]
		]
	*/

	actual = dgebak( 'scale', 'right', 4, 2, 0, 3, scale, -1, 3, V, -2, -1, 7 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (right side) (column-major) (negative strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.1, 0.1, 1.0, 10.0 ] );
	V = new Float64Array( [ 0.04, 0.5, 3.0, 10.0, 0.08, 0.6, 4.0, 20.0 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 0.004, 0.05, 3.0, 100.0, 0.008, 0.06, 4.0, 200.0 ] );

	/*
		expected = [
			[ 100.0, 200.0 ],
			[ 3.0, 4.0 ],
			[ 0.05, 0.06 ],
			[ 0.004, 0.008 ]
		]
	*/

	actual = dgebak( 'scale', 'right', 4, 2, 0, 3, scale, -1, 3, V, -1, -4, 7 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (left side) (row-major) (negative strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.01, 0.1, 1.0, 10.0 ] );
	V = new Float64Array( [ 0.08, 0.04, 0.6, 0.5, 4.0, 3.0, 20.0, 10.0 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 8.0, 4.0, 6.0, 5.0, 4.0, 3.0, 2.0, 1.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 4.0, 8.0 ]
		]
	*/

	actual = dgebak( 'scale', 'left', 4, 2, 0, 3, scale, -1, 3, V, -2, -1, 7 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (left side) (column-major) (negative strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.01, 0.1, 1.0, 10.0 ] );
	V = new Float64Array( [ 10.0, 3.0, 0.5, 0.04, 20.0, 4.0, 0.6, 0.08 ] );
	V = new Float64Array( [ 0.08, 0.6, 4.0, 20.0, 0.04, 0.5, 3.0, 10.0 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 8.0, 6.0, 4.0, 2.0, 4.0, 5.0, 3.0, 1.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 4.0, 8.0 ]
		]
	*/

	actual = dgebak( 'scale', 'left', 4, 2, 0, 3, scale, -1, 3, V, -1, -4, 7 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (right side) (row-major) (negative strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 1.0, 0.0 ] );
	V = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 7.0, 8.0, 5.0, 6.0, 3.0, 4.0, 1.0, 2.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 0, 1, scale, 1, 0, V, 2, 1, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (right side) (column-major) (negative strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 1.0, 2.0, 3.0 ] );
	V = new Float64Array( [ 8.0, 6.0, 4.0, 2.0, 7.0, 5.0, 3.0, 1.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 2.0, 4.0, 6.0, 8.0, 1.0, 3.0, 5.0, 7.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 0, 1, scale, -1, 3, V, -1, -4, 7 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (left side) (row-major) (negative strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 1.0, 2.0, 3.0 ] );
	V = new Float64Array( [ 8.0, 7.0, 6.0, 5.0, 4.0, 3.0, 2.0, 1.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 2.0, 1.0, 4.0, 3.0, 6.0, 5.0, 8.0, 7.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'left', 4, 2, 0, 1, scale, -1, 3, V, -2, -1, 7 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (left side) (column-major) (negative strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 1.0, 2.0, 3.0 ] );
	V = new Float64Array( [ 8.0, 6.0, 4.0, 2.0, 7.0, 5.0, 3.0, 1.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 2.0, 4.0, 6.0, 8.0, 1.0, 3.0, 5.0, 7.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'left', 4, 2, 0, 1, scale, -1, 3, V, -1, -4, 7 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (right side) (row-major) (negative strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.1, 10.0, 3.0 ] );
	V = new Float64Array( [ 2.0, 1.0, 0.6, 0.5, 4.0, 3.0, 8.0, 7.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 2.0, 1.0, 0.06, 0.05, 40.0, 30.0, 8.0, 7.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 30.0, 40.0 ],
			[ 0.05, 0.06 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 1, 2, scale, -1, 3, V, -2, -1, 7 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (right side) (column-major) (negative strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.1, 10.0, 3.0 ] );
	V = new Float64Array( [ 2.0, 0.6, 4.0, 8.0, 1.0, 0.5, 3.0, 7.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 2.0, 0.06, 40.0, 8.0, 1.0, 0.05, 30.0, 7.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 30.0, 40.0 ],
			[ 0.05, 0.06 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 1, 2, scale, -1, 3, V, -1, -4, 7 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (left side) (row-major)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.1, 10.0, 3.0 ] );
	V = new Float64Array( [ 2.0, 1.0, 0.6, 0.5, 4.0, 8.0, 8.0, 7.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 8.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 2.0, 1.0, 6.0, 5.0, 0.4, 0.8, 8.0, 7.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 0.8, 0.4 ],
			[ 5.0, 6.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'left', 4, 2, 1, 2, scale, -1, 3, V, -2, -1, 7 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (left side) (column-major) (negative strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 0.0, 0.1, 10.0, 3.0 ] );
	V = new Float64Array( [ 2.0, 0.6, 4.0, 8.0, 1.0, 0.5, 8.0, 7.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 8.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 2.0, 6.0, 0.4, 8.0, 1.0, 5.0, 0.8, 7.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 0.8, 0.4 ],
			[ 5.0, 6.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'left', 4, 2, 1, 2, scale, -1, 3, V, -1, -4, 7 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = permute (row-major) (mixed strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 2.0, 0.0 ] );
	V = new Float64Array( [ 2.0, 1.0, 4.0, 3.0, 6.0, 5.0, 8.0, 7.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 2.0, 1.0, 6.0, 5.0, 4.0, 3.0, 8.0, 7.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 2, 2, scale, 1, 0, V, 2, -1, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = permute (column-major) (mixed strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 2.0, 0.0 ] );
	V = new Float64Array( [ 2.0, 4.0, 6.0, 8.0, 1.0, 3.0, 5.0, 7.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 2.0, 6.0, 4.0, 8.0, 1.0, 5.0, 3.0, 7.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 2, 2, scale, 1, 0, V, 1, -4, 4 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = both (row-major) (mixed strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 2.0, 0.0 ] );
	V = new Float64Array( [ 2.0, 1.0, 4.0, 3.0, 6.0, 5.0, 8.0, 7.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 2.0, 1.0, 6.0, 5.0, 4.0, 3.0, 8.0, 7.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 2, 2, scale, 1, 0, V, 2, -1, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = both (column-major) (mixed strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 2.0, 0.0 ] );
	V = new Float64Array( [ 2.0, 4.0, 6.0, 8.0, 1.0, 3.0, 5.0, 7.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 2.0, 6.0, 4.0, 8.0, 1.0, 5.0, 3.0, 7.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 2, 2, scale, 1, 0, V, 1, -4, 4 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (right side) (row-major) (mixed strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 10.0, 1.0, 0.1, 0.1 ] );
	V = new Float64Array( [ 20.0, 10.0, 4.0, 3.0, 0.6, 0.5, 0.08, 0.04 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 200.0, 100.0, 4.0, 3.0, 0.06, 0.05, 0.008, 0.004 ] );

	/*
		expected = [
			[ 100.0, 200.0 ],
			[ 3.0, 4.0 ],
			[ 0.05, 0.06 ],
			[ 0.004, 0.008 ]
		]
	*/

	actual = dgebak( 'scale', 'right', 4, 2, 0, 3, scale, 1, 0, V, 2, -1, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (right side) (column-major) (mixed strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 10.0, 1.0, 0.1, 0.1 ] );
	V = new Float64Array( [ 20.0, 4.0, 0.6, 0.08, 10.0, 3.0, 0.5, 0.04 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 200.0, 4.0, 0.06, 0.008, 100.0, 3.0, 0.05, 0.004 ] );

	/*
		expected = [
			[ 100.0, 200.0 ],
			[ 3.0, 4.0 ],
			[ 0.05, 0.06 ],
			[ 0.004, 0.008 ]
		]
	*/

	actual = dgebak( 'scale', 'right', 4, 2, 0, 3, scale, 1, 0, V, 1, -4, 4 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (left side) (row-major) (mixed strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 10.0, 1.0, 0.1, 0.01 ] );
	V = new Float64Array( [ 20.0, 10.0, 4.0, 3.0, 0.6, 0.5, 0.08, 0.04 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 2.0, 1.0, 4.0, 3.0, 6.0, 5.0, 8.0, 4.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 4.0, 8.0 ]
		]
	*/

	actual = dgebak( 'scale', 'left', 4, 2, 0, 3, scale, 1, 0, V, 2, -1, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (left side) (column-major) (mixed strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 10.0, 1.0, 0.1, 0.01 ] );
	V = new Float64Array( [ 20.0, 4.0, 0.6, 0.08, 10.0, 3.0, 0.5, 0.04 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 2.0, 4.0, 6.0, 8.0, 1.0, 3.0, 5.0, 4.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 4.0, 8.0 ]
		]
	*/

	actual = dgebak( 'scale', 'left', 4, 2, 0, 3, scale, 1, 0, V, 1, -4, 4 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (right side) (row-major) (mixed strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 1.0, 0.0 ] );
	V = new Float64Array( [ 2.0, 1.0, 4.0, 3.0, 6.0, 5.0, 8.0, 7.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 8.0, 7.0, 6.0, 5.0, 4.0, 3.0, 2.0, 1.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 0, 1, scale, 1, 0, V, 2, -1, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (right side) (column-major) (mixed strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 1.0, 0.0 ] );
	V = new Float64Array( [ 2.0, 4.0, 6.0, 8.0, 1.0, 3.0, 5.0, 7.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 8.0, 6.0, 4.0, 2.0, 7.0, 5.0, 3.0, 1.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 0, 1, scale, 1, 0, V, 1, -4, 4 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (left side) (row-major) (mixed strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 1.0, 0.0 ] );
	V = new Float64Array( [ 2.0, 1.0, 4.0, 3.0, 6.0, 5.0, 8.0, 7.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 8.0, 7.0, 6.0, 5.0, 4.0, 3.0, 2.0, 1.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'left', 4, 2, 0, 1, scale, 1, 0, V, 2, -1, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (left side) (column-major) (mixed strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 2.0, 1.0, 0.0 ] );
	V = new Float64Array( [ 2.0, 4.0, 6.0, 8.0, 1.0, 3.0, 5.0, 7.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 8.0, 6.0, 4.0, 2.0, 7.0, 5.0, 3.0, 1.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'left', 4, 2, 0, 1, scale, 1, 0, V, 1, -4, 4 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (right side) (row-major) (mixed strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 10.0, 0.1, 0.0 ] );
	V = new Float64Array( [ 8.0, 7.0, 4.0, 3.0, 0.6, 0.5, 2.0, 1.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 8.0, 7.0, 40.0, 30.0, 0.06, 0.05, 2.0, 1.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 30.0, 40.0 ],
			[ 0.05, 0.06 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 1, 2, scale, 1, 0, V, 2, -1, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (right side) (column-major) (mixed strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 10.0, 0.1, 0.0 ] );
	V = new Float64Array( [ 8.0, 4.0, 0.6, 2.0, 7.0, 3.0, 0.5, 1.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 8.0, 40.0, 0.06, 2.0, 7.0, 30.0, 0.05, 1.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 30.0, 40.0 ],
			[ 0.05, 0.06 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 1, 2, scale, 1, 0, V, 1, -4, 4 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (left side) (row-major) (mixed strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 10.0, 0.1, 0.0 ] );
	V = new Float64Array( [ 8.0, 7.0, 4.0, 8.0, 0.6, 0.5, 2.0, 1.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 8.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 8.0, 7.0, 0.4, 0.8, 6.0, 5.0, 2.0, 1.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 0.8, 0.4 ],
			[ 5.0, 6.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'left', 4, 2, 1, 2, scale, 1, 0, V, 2, -1, 1 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (left side) (column-major) (mixed strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 10.0, 0.1, 0.0 ] );
	V = new Float64Array( [ 8.0, 4.0, 0.6, 2.0, 7.0, 8.0, 0.5, 1.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 8.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 8.0, 0.4, 6.0, 2.0, 7.0, 0.8, 5.0, 1.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 0.8, 0.4 ],
			[ 5.0, 6.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'left', 4, 2, 1, 2, scale, 1, 0, V, 1, -4, 4 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = permute (row-major) (large strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 0.0, 2.0, 0.0, 2.0, 0.0, 0.0, 0.0 ] );
	V = new Float64Array( [ 1.0, 0.0, 2.0, 0.0, 3.0, 0.0, 4.0, 0.0, 5.0, 0.0, 6.0, 0.0, 7.0, 0.0, 8.0, 0.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 1.0, 0.0, 2.0, 0.0, 5.0, 0.0, 6.0, 0.0, 3.0, 0.0, 4.0, 0.0, 7.0, 0.0, 8.0, 0.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 2, 2, scale, 2, 0, V, 4, 2, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = permute (column-major) (large strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 0.0, 2.0, 0.0, 2.0, 0.0, 0.0, 0.0 ] );
	V = new Float64Array( [ 1.0, 0.0, 3.0, 0.0, 5.0, 0.0, 7.0, 0.0, 2.0, 0.0, 4.0, 0.0, 6.0, 0.0, 8.0, 0.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 1.0, 0.0, 5.0, 0.0, 3.0, 0.0, 7.0, 0.0, 2.0, 0.0, 6.0, 0.0, 4.0, 0.0, 8.0, 0.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 2, 2, scale, 2, 0, V, 2, 8, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = both (row-major) (mixed strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 0.0, 2.0, 0.0, 2.0, 0.0, 0.0, 0.0 ] );
	V = new Float64Array( [ 1.0, 0.0, 2.0, 0.0, 3.0, 0.0, 4.0, 0.0, 5.0, 0.0, 6.0, 0.0, 7.0, 0.0, 8.0, 0.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 1.0, 0.0, 2.0, 0.0, 5.0, 0.0, 6.0, 0.0, 3.0, 0.0, 4.0, 0.0, 7.0, 0.0, 8.0, 0.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 2, 2, scale, 2, 0, V, 4, 2, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output when ilo = ihi and job = both (column-major) (large strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 0.0, 2.0, 0.0, 2.0, 0.0, 0.0, 0.0 ] );
	V = new Float64Array( [ 1.0, 0.0, 3.0, 0.0, 5.0, 0.0, 7.0, 0.0, 2.0, 0.0, 4.0, 0.0, 6.0, 0.0, 8.0, 0.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 1.0, 0.0, 5.0, 0.0, 3.0, 0.0, 7.0, 0.0, 2.0, 0.0, 6.0, 0.0, 4.0, 0.0, 8.0, 0.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 2, 2, scale, 2, 0, V, 2, 8, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (right side) (row-major) (large strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 10.0, 0.0, 1.0, 0.0, 0.1, 0.0, 0.1, 0.0 ] );
	V = new Float64Array( [ 10.0, 0.0, 20.0, 0.0, 3.0, 0.0, 4.0, 0.0, 0.5, 0.0, 0.6, 0.0, 0.04, 0.0, 0.08, 0.0 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 100.0, 0.0, 200.0, 0.0, 3.0, 0.0, 4.0, 0.0, 0.05, 0.0, 0.06, 0.0, 0.004, 0.0, 0.008, 0.0 ] );

	/*
		expected = [
			[ 100.0, 200.0 ],
			[ 3.0, 4.0 ],
			[ 0.05, 0.06 ],
			[ 0.004, 0.008 ]
		]
	*/

	actual = dgebak( 'scale', 'right', 4, 2, 0, 3, scale, 2, 0, V, 4, 2, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (right side) (column-major) (large strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 10.0, 0.0, 1.0, 0.0, 0.1, 0.0, 0.1, 0.0 ] );
	V = new Float64Array( [ 10.0, 0.0, 3.0, 0.0, 0.5, 0.0, 0.04, 0.0, 20.0, 0.0, 4.0, 0.0, 0.6, 0.0, 0.08, 0.0 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 100.0, 0.0, 3.0, 0.0, 0.05, 0.0, 0.004, 0.0, 200.0, 0.0, 4.0, 0.0, 0.06, 0.0, 0.008, 0.0 ] );

	/*
		expected = [
			[ 100.0, 200.0 ],
			[ 3.0, 4.0 ],
			[ 0.05, 0.06 ],
			[ 0.004, 0.008 ]
		]
	*/

	actual = dgebak( 'scale', 'right', 4, 2, 0, 3, scale, 2, 0, V, 2, 8, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (left side) (row-major) (large strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 10.0, 0.0, 1.0, 0.0, 0.1, 0.0, 0.01, 0.0 ] );
	V = new Float64Array( [ 10.0, 0.0, 20.0, 0.0, 3.0, 0.0, 4.0, 0.0, 0.5, 0.0, 0.6, 0.0, 0.04, 0.0, 0.08, 0.0 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 1.0, 0.0, 2.0, 0.0, 3.0, 0.0, 4.0, 0.0, 5.0, 0.0, 6.0, 0.0, 4.0, 0.0, 8.0, 0.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 4.0, 8.0 ]
		]
	*/

	actual = dgebak( 'scale', 'left', 4, 2, 0, 3, scale, 2, 0, V, 4, 2, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = scale (left side) (column-major) (large strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 10.0, 0.0, 1.0, 0.0, 0.1, 0.0, 0.01, 0.0 ] );
	V = new Float64Array( [ 10.0, 0.0, 3.0, 0.0, 0.5, 0.0, 0.04, 0.0, 20.0, 0.0, 4.0, 0.0, 0.6, 0.0, 0.08, 0.0 ] );

	/*
		V = [
			[ 10.0, 20.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 0.04, 0.08 ]
		]
	*/

	expected = new Float64Array( [ 1.0, 0.0, 3.0, 0.0, 5.0, 0.0, 4.0, 0.0, 2.0, 0.0, 4.0, 0.0, 6.0, 0.0, 8.0, 0.0 ] );

	/*
		expected = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 4.0, 8.0 ]
		]
	*/

	actual = dgebak( 'scale', 'left', 4, 2, 0, 3, scale, 2, 0, V, 2, 8, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (right side) (row-major) (large strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 0.0, 2.0, 0.0, 1.0, 0.0, 0.0, 0.0 ] );
	V = new Float64Array( [ 1.0, 0.0, 2.0, 0.0, 3.0, 0.0, 4.0, 0.0, 5.0, 0.0, 6.0, 0.0, 7.0, 0.0, 8.0, 0.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 7.0, 0.0, 8.0, 0.0, 5.0, 0.0, 6.0, 0.0, 3.0, 0.0, 4.0, 0.0, 1.0, 0.0, 2.0, 0.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 0, 1, scale, 2, 0, V, 4, 2, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (right side) (column-major) (large strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 0.0, 2.0, 0.0, 1.0, 0.0, 0.0, 0.0 ] );
	V = new Float64Array( [ 1.0, 0.0, 3.0, 0.0, 5.0, 0.0, 7.0, 0.0, 2.0, 0.0, 4.0, 0.0, 6.0, 0.0, 8.0, 0.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 7.0, 0.0, 5.0, 0.0, 3.0, 0.0, 1.0, 0.0, 8.0, 0.0, 6.0, 0.0, 4.0, 0.0, 2.0, 0.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'right', 4, 2, 0, 1, scale, 2, 0, V, 2, 8, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (left side) (row-major) (large strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 0.0, 2.0, 0.0, 1.0, 0.0, 0.0, 0.0 ] );
	V = new Float64Array( [ 1.0, 0.0, 2.0, 0.0, 3.0, 0.0, 4.0, 0.0, 5.0, 0.0, 6.0, 0.0, 7.0, 0.0, 8.0, 0.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 7.0, 0.0, 8.0, 0.0, 5.0, 0.0, 6.0, 0.0, 3.0, 0.0, 4.0, 0.0, 1.0, 0.0, 2.0, 0.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'left', 4, 2, 0, 1, scale, 2, 0, V, 4, 2, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = permute (left side) (column-major) (large strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 0.0, 2.0, 0.0, 1.0, 0.0, 0.0, 0.0 ] );
	V = new Float64Array( [ 1.0, 0.0, 3.0, 0.0, 5.0, 0.0, 7.0, 0.0, 2.0, 0.0, 4.0, 0.0, 6.0, 0.0, 8.0, 0.0 ] );

	/*
		V = [
			[ 1.0, 2.0 ],
			[ 3.0, 4.0 ],
			[ 5.0, 6.0 ],
			[ 7.0, 8.0 ]
		]
	*/

	expected = new Float64Array( [ 7.0, 0.0, 5.0, 0.0, 3.0, 0.0, 1.0, 0.0, 8.0, 0.0, 6.0, 0.0, 4.0, 0.0, 2.0, 0.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 5.0, 6.0 ],
			[ 3.0, 4.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'permute', 'left', 4, 2, 0, 1, scale, 2, 0, V, 2, 8, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (right side) (row-major) (large strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 0.0, 10.0, 0.0, 0.1, 0.0, 0.0, 0.0 ] );
	V = new Float64Array( [ 7.0, 0.0, 8.0, 0.0, 3.0, 0.0, 4.0, 0.0, 0.5, 0.0, 0.6, 0.0, 1.0, 0.0, 2.0, 0.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 7.0, 0.0, 8.0, 0.0, 30.0, 0.0, 40.0, 0.0, 0.05, 0.0, 0.06, 0.0, 1.0, 0.0, 2.0, 0.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 30.0, 40.0 ],
			[ 0.05, 0.06 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 1, 2, scale, 2, 0, V, 4, 2, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (right side) (column-major) (large strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 0.0, 10.0, 0.0, 0.1, 0.0, 0.0, 0.0 ] );
	V = new Float64Array( [ 7.0, 0.0, 3.0, 0.0, 0.5, 0.0, 1.0, 0.0, 8.0, 0.0, 4.0, 0.0, 0.6, 0.0, 2.0, 0.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 3.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 7.0, 0.0, 30.0, 0.0, 0.05, 0.0, 1.0, 0.0, 8.0, 0.0, 40.0, 0.0, 0.06, 0.0, 2.0, 0.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 30.0, 40.0 ],
			[ 0.05, 0.06 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'right', 4, 2, 1, 2, scale, 2, 0, V, 2, 8, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (left side) (row-major) (large strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 0.0, 10.0, 0.0, 0.1, 0.0, 0.0, 0.0 ] );
	V = new Float64Array( [ 7.0, 0.0, 8.0, 0.0, 8.0, 0.0, 4.0, 0.0, 0.5, 0.0, 0.6, 0.0, 1.0, 0.0, 2.0, 0.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 8.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 7.0, 0.0, 8.0, 0.0, 0.8, 0.0, 0.4, 0.0, 5.0, 0.0, 6.0, 0.0, 1.0, 0.0, 2.0, 0.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 0.8, 0.4 ],
			[ 5.0, 6.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'left', 4, 2, 1, 2, scale, 2, 0, V, 4, 2, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});

tape( 'the function returns expected output for job = both (left side) (column-major) (large strides)', function test( t ) {
	var expected;
	var actual;
	var scale;
	var V;

	scale = new Float64Array( [ 3.0, 0.0, 10.0, 0.0, 0.1, 0.0, 0.0, 0.0 ] );
	V = new Float64Array( [ 7.0, 0.0, 8.0, 0.0, 0.5, 0.0, 1.0, 0.0, 8.0, 0.0, 4.0, 0.0, 0.6, 0.0, 2.0, 0.0 ] );

	/*
		V = [
			[ 7.0, 8.0 ],
			[ 8.0, 4.0 ],
			[ 0.5, 0.6 ],
			[ 1.0, 2.0 ]
		]
	*/

	expected = new Float64Array( [ 7.0, 0.0, 0.8, 0.0, 5.0, 0.0, 1.0, 0.0, 8.0, 0.0, 0.4, 0.0, 6.0, 0.0, 2.0, 0.0 ] );

	/*
		expected = [
			[ 7.0, 8.0 ],
			[ 0.8, 0.4 ],
			[ 5.0, 6.0 ],
			[ 1.0, 2.0 ]
		]
	*/

	actual = dgebak( 'both', 'left', 4, 2, 1, 2, scale, 2, 0, V, 2, 8, 0 );

	t.deepEqual( actual, expected, 'returns expected value' );
	t.end();
});