From 3942c89cc466c548eb6b16f8eec8cc22a3b6281f Mon Sep 17 00:00:00 2001
From: aman-095 <bhansali.1@iitj.ac.in>
Date: Wed, 28 Aug 2024 13:36:36 +0530
Subject: [PATCH] refactor: update implementation, descriptions, examples and
 add tests

---
 .../@stdlib/blas/base/dspmv/README.md         |  20 +-
 .../blas/base/dspmv/benchmark/benchmark.js    |  12 +-
 .../base/dspmv/benchmark/benchmark.ndarray.js |  12 +-
 .../@stdlib/blas/base/dspmv/docs/repl.txt     |  11 +-
 .../blas/base/dspmv/docs/types/index.d.ts     |  19 +-
 .../blas/base/dspmv/docs/types/test.ts        | 261 ++++---
 .../@stdlib/blas/base/dspmv/examples/index.js |   7 +-
 .../@stdlib/blas/base/dspmv/lib/base.js       | 125 ++-
 .../@stdlib/blas/base/dspmv/lib/dspmv.js      |  17 +-
 .../@stdlib/blas/base/dspmv/lib/index.js      |   4 +-
 .../@stdlib/blas/base/dspmv/lib/ndarray.js    |  26 +-
 .../@stdlib/blas/base/dspmv/package.json      |   2 +-
 .../column_major_complex_access_pattern.json  |  17 +
 ...ajor_xoyt.json => column_major_lower.json} |   5 +-
 .../dspmv/test/fixtures/column_major_oap.json |  17 +
 .../dspmv/test/fixtures/column_major_ox.json  |  17 +
 .../dspmv/test/fixtures/column_major_oy.json  |  17 +
 .../dspmv/test/fixtures/column_major_sap.json |  17 +
 .../test/fixtures/column_major_sapn.json      |  17 +
 .../test/fixtures/column_major_upper.json     |  17 +
 .../test/fixtures/column_major_xnyn.json      |   5 +-
 .../test/fixtures/column_major_xnyp.json      |   3 +-
 .../test/fixtures/column_major_xpyn.json      |   3 +-
 .../test/fixtures/column_major_xpyp.json      |   3 +-
 .../row_major_complex_access_pattern.json     |  17 +
 ...w_major_xoyt.json => row_major_lower.json} |   9 +-
 .../dspmv/test/fixtures/row_major_oap.json    |  17 +
 .../dspmv/test/fixtures/row_major_ox.json     |  17 +
 .../dspmv/test/fixtures/row_major_oy.json     |  17 +
 .../dspmv/test/fixtures/row_major_sap.json    |  17 +
 .../dspmv/test/fixtures/row_major_sapn.json   |  17 +
 .../dspmv/test/fixtures/row_major_upper.json  |  17 +
 .../dspmv/test/fixtures/row_major_xnyn.json   |   5 +-
 .../dspmv/test/fixtures/row_major_xnyp.json   |   3 +-
 .../dspmv/test/fixtures/row_major_xpyn.json   |   3 +-
 .../dspmv/test/fixtures/row_major_xpyp.json   |   3 +-
 .../blas/base/dspmv/test/test.dspmv.js        | 411 ++++++----
 .../blas/base/dspmv/test/test.ndarray.js      | 733 +++++++++++++-----
 38 files changed, 1341 insertions(+), 599 deletions(-)
 create mode 100644 lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_complex_access_pattern.json
 rename lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/{column_major_xoyt.json => column_major_lower.json} (95%)
 create mode 100644 lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_oap.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_ox.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_oy.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_sap.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_sapn.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_upper.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_complex_access_pattern.json
 rename lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/{row_major_xoyt.json => row_major_lower.json} (77%)
 create mode 100644 lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_oap.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_ox.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_oy.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_sap.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_sapn.json
 create mode 100644 lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_upper.json

diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/README.md b/lib/node_modules/@stdlib/blas/base/dspmv/README.md
index 2b2aa972558..0149404af19 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/README.md
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/README.md
@@ -20,7 +20,7 @@ limitations under the License.
 
 # dspmv
 
-> Perform the matrix-vector operation `y = α*A*x + β*y` where `α` and `β` are scalars, `x` and `y` are `N` element vectors and, `A` is an `N` by `N` symmetric matrix supplied in packed form.
+> Perform the matrix-vector operation `y = α*A*x + β*y`.
 
 <section class = "usage">
 
@@ -91,7 +91,7 @@ dspmv( 'row-major', 'upper', 2, 1.0, AP, x1, -1, 1.0, y1, -1 );
 // y0 => <Float64Array>[ 0.0, 6.0, 4.0 ]
 ```
 
-#### dspmv.ndarray( order, uplo, N, α, AP, oa, x, sx, ox, β, y, sy, oy )
+#### dspmv.ndarray( order, uplo, N, α, AP, sap, oap x, sx, ox, β, y, sy, oy )
 
 Performs the matrix-vector operation `y = α*A*x + β*y` using alternative indexing semantics and where `α` and `β` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form `AP`.
 
@@ -102,13 +102,14 @@ var AP = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
 var x = new Float64Array( [ 1.0, 1.0, 1.0 ] );
 var y = new Float64Array( [ 1.0, 1.0, 1.0 ] );
 
-dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 0, x, 1, 0, 1.0, y, 1, 0 );
+dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, 1, 0 );
 // y => <Float64Array>[ 7.0, 12.0, 15.0 ]
 ```
 
 The function has the following additional parameters:
 
--   **oa**: starting index for `AP`.
+-   **sap**: `AP` stride length.
+-   **oap**: starting index for `AP`.
 -   **ox**: starting index for `x`.
 -   **oy**: starting index for `y`.
 
@@ -117,11 +118,11 @@ While [`typed array`][mdn-typed-array] views mandate a view offset based on the
 ```javascript
 var Float64Array = require( '@stdlib/array/float64' );
 
-var AP = new Float64Array( [ 0.0, 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
+var AP = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
 var x = new Float64Array( [ 1.0, 1.0, 1.0 ] );
 var y = new Float64Array( [ 1.0, 1.0, 1.0 ] );
 
-dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 2, x, 1, 0, 1.0, y, -1, 2 );
+dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, -1, 2 );
 // y => <Float64Array>[ 15.0, 12.0, 7.0 ]
 ```
 
@@ -153,13 +154,16 @@ var opts = {
     'dtype': 'float64'
 };
 
-var N = 3;
+var N = 5;
 var AP = discreteUniform( N * ( N + 1 ) / 2, -10, 10, opts );
 
 var x = discreteUniform( N, -10, 10, opts );
 var y = discreteUniform( N, -10, 10, opts );
 
-dspmv.ndarray( 'row-major', 'upper', N, 1.0, AP, 0, x, 1, 0, 1.0, y, 1, 0 );
+dspmv( 'row-major', 'upper', N, 1.0, AP, x, 1, 1.0, y, 1 );
+console.log( y );
+
+dspmv.ndarray( 'row-major', 'upper', N, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, 1, 0 );
 console.log( y );
 ```
 
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/benchmark/benchmark.js b/lib/node_modules/@stdlib/blas/base/dspmv/benchmark/benchmark.js
index 7a27e9cbf7a..4b115eda079 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/benchmark/benchmark.js
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/benchmark/benchmark.js
@@ -42,13 +42,13 @@ var options = {
 * Creates a benchmark function.
 *
 * @private
-* @param {PositiveInteger} len - array length
+* @param {PositiveInteger} N - array length
 * @returns {Function} benchmark function
 */
-function createBenchmark( len ) {
-	var AP = uniform( len * ( len + 1 ) / 2, -10.0, 10.0, options );
-	var x = uniform( len, -10.0, 10.0, options );
-	var y = uniform( len, -10.0, 10.0, options );
+function createBenchmark( N ) {
+	var AP = uniform( N*(N+1)/2, -10.0, 10.0, options );
+	var x = uniform( N, -10.0, 10.0, options );
+	var y = uniform( N, -10.0, 10.0, options );
 	return benchmark;
 
 	/**
@@ -63,7 +63,7 @@ function createBenchmark( len ) {
 
 		b.tic();
 		for ( i = 0; i < b.iterations; i++ ) {
-			z = dspmv( 'row-major', 'upper', len, 1.0, AP, x, 1, 1.0, y, 1 );
+			z = dspmv( 'row-major', 'upper', N, 1.0, AP, x, 1, 1.0, y, 1 );
 			if ( isnan( z[ i%z.length ] ) ) {
 				b.fail( 'should not return NaN' );
 			}
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/benchmark/benchmark.ndarray.js b/lib/node_modules/@stdlib/blas/base/dspmv/benchmark/benchmark.ndarray.js
index a2877463004..ace2bed1eff 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/benchmark/benchmark.ndarray.js
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/benchmark/benchmark.ndarray.js
@@ -42,13 +42,13 @@ var options = {
 * Creates a benchmark function.
 *
 * @private
-* @param {PositiveInteger} len - array length
+* @param {PositiveInteger} N - array length
 * @returns {Function} benchmark function
 */
-function createBenchmark( len ) {
-	var AP = uniform( len * ( len + 1 ) / 2, -10.0, 10.0, options );
-	var x = uniform( len, -10.0, 10.0, options );
-	var y = uniform( len, -10.0, 10.0, options );
+function createBenchmark( N ) {
+	var AP = uniform( N * ( N + 1 ) / 2, -10.0, 10.0, options );
+	var x = uniform( N, -10.0, 10.0, options );
+	var y = uniform( N, -10.0, 10.0, options );
 	return benchmark;
 
 	/**
@@ -63,7 +63,7 @@ function createBenchmark( len ) {
 
 		b.tic();
 		for ( i = 0; i < b.iterations; i++ ) {
-			z = dspmv( 'row-major', 'upper', len, 1.0, AP, 0, x, 1, 0, 1.0, y, 1, 0 );
+			z = dspmv( 'row-major', 'upper', N, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, 1, 0 );
 			if ( isnan( z[ i%z.length ] ) ) {
 				b.fail( 'should not return NaN' );
 			}
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/docs/repl.txt b/lib/node_modules/@stdlib/blas/base/dspmv/docs/repl.txt
index 5340eeefa55..568959135d1 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/docs/repl.txt
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/docs/repl.txt
@@ -77,7 +77,7 @@
     <Float64Array>[ ~6.0, ~4.0 ]
 
 
-{{alias}}.ndarray( ord, uplo, N, α, AP, oa, x, sx, ox, β, y, sy, oy )
+{{alias}}.ndarray( ord, uplo, N, α, AP, sap, oap, x, sx, ox, β, y, sy, oy )
     Performs the matrix-vector operation `y = α*A*x + β*y` using alternative
     indexing semantics and where `α` and `β` are scalars, `x` and `y` are `N`
     element vectors, and `A` is an `N` by `N` symmetric matrix supplied in
@@ -106,7 +106,10 @@
     AP: Float64Array
         Matrix in packed form.
 
-    oa: integer
+    sap: integer
+        Index increment for `AP`.
+
+    oap: integer
         Starting index for `AP`.
 
     x: Float64Array
@@ -143,14 +146,14 @@
     > var AP = new {{alias:@stdlib/array/float64}}( [ 1.0, 2.0, 3.0, 4.0 ] );
     > var ord = 'row-major';
     > var uplo = 'upper';
-    > {{alias}}.ndarray( ord, uplo, 2, 1.0, AP, 0, x, 1, 0, 1.0, y, 1, 0 )
+    > {{alias}}.ndarray( ord, uplo, 2, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, 1, 0 )
     <Float64Array>[ ~4.0, ~6.0 ]
 
     // Advanced indexing:
     > x = new {{alias:@stdlib/array/float64}}( [ 1.0, 1.0 ] );
     > y = new {{alias:@stdlib/array/float64}}( [ 1.0, 1.0 ] );
     > AP = new {{alias:@stdlib/array/float64}}( [ 1.0, 2.0, 3.0, 4.0 ] );
-    > {{alias}}.ndarray( ord, uplo, 2, 1.0, AP, 0, x, -1, 1, 1.0, y, -1, 1 )
+    > {{alias}}.ndarray( ord, uplo, 2, 1.0, AP, 1, 0, x, -1, 1, 1.0, y, -1, 1 )
     <Float64Array>[ ~6.0, ~4.0 ]
 
     See Also
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/docs/types/index.d.ts b/lib/node_modules/@stdlib/blas/base/dspmv/docs/types/index.d.ts
index e755ce04ec5..93fac01e31d 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/docs/types/index.d.ts
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/docs/types/index.d.ts
@@ -27,7 +27,7 @@ import { Layout, MatrixTriangle } from '@stdlib/types/blas';
 */
 interface Routine {
 	/**
-	* Performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form.
+	* Performs the matrix-vector operation `y = α*A*x + β*y` where `α` and `β` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form.
 	*
 	* @param order - storage layout
 	* @param uplo - specifies whether the upper or lower triangular part of the symmetric matrix `A` is being supplied
@@ -54,21 +54,22 @@ interface Routine {
 	( order: Layout, uplo: MatrixTriangle, N: number, alpha: number, AP: Float64Array, x: Float64Array, strideX: number, beta: number, y: Float64Array, strideY: number ): Float64Array;
 
 	/**
-	* Performs the matrix-vector operation `y = alpha*A*x + beta*y` using alternative indexing semantics and where `alpha` and `beta` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form.
+	* Performs the matrix-vector operation `y = α*A*x + β*y`, using alternative indexing semantics and where `α` and `β` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form.
 	*
 	* @param order - storage layout
 	* @param uplo - specifies whether the upper or lower triangular part of the symmetric matrix `A` is being supplied
 	* @param N - number of columns in the matrix `A`
 	* @param alpha - scalar constant
 	* @param AP - packed form of a symmetric matrix `A`
-	* @param offsetAP - starting `AP` index
+	* @param sap - `AP` stride length
+	* @param oap - starting index for `AP`
 	* @param x - first input array
 	* @param strideX - `x` stride length
-	* @param offsetX - starting `x` index
+	* @param offsetX - starting index for `x`
 	* @param beta - scalar constant
 	* @param y - second input array
 	* @param strideY - `y` stride length
-	* @param offsetY - starting `y` index
+	* @param offsetY - starting index for `y`
 	* @returns `y`
 	*
 	* @example
@@ -78,14 +79,14 @@ interface Routine {
 	* var x = new Float64Array( [ 1.0, 1.0, 1.0 ] );
 	* var y = new Float64Array( [ 1.0, 1.0, 1.0 ] );
 	*
-	* dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 0, x, 1, 0, 1.0, y, 1, 0 );
+	* dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, 1, 0 );
 	* // y => <Float64Array>[ 7.0, 12.0, 15.0 ]
 	*/
-	ndarray( order: Layout, uplo: MatrixTriangle, N: number, alpha: number, AP: Float64Array, offsetAP: number, x: Float64Array, strideX: number, offsetX: number, beta: number, y: Float64Array, strideY: number, offsetY: number ): Float64Array;
+	ndarray( order: Layout, uplo: MatrixTriangle, N: number, alpha: number, AP: Float64Array, strideAP: number, offsetAP: number, x: Float64Array, strideX: number, offsetX: number, beta: number, y: Float64Array, strideY: number, offsetY: number ): Float64Array;
 }
 
 /**
-* Performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form.
+* Performs the matrix-vector operation `y = α*A*x + β*y` where `α` and `β` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form.
 *
 * @param order - storage layout
 * @param uplo - specifies whether the upper or lower triangular part of the symmetric matrix `A` is being supplied
@@ -116,7 +117,7 @@ interface Routine {
 * var x = new Float64Array( [ 1.0, 1.0, 1.0 ] );
 * var y = new Float64Array( [ 1.0, 1.0, 1.0 ] );
 *
-* dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 0, x, 1, 0, 1.0, y, -1, 2 );
+* dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, -1, 2 );
 * // y => <Float64Array>[ 15.0, 12.0, 7.0 ]
 */
 declare var dspmv: Routine;
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/docs/types/test.ts b/lib/node_modules/@stdlib/blas/base/dspmv/docs/types/test.ts
index d717efff965..535cd59006f 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/docs/types/test.ts
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/docs/types/test.ts
@@ -215,7 +215,7 @@ import dspmv = require( './index' );
 	const y = new Float64Array( 10 );
 	const A = new Float64Array( 55 );
 
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectType Float64Array
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectType Float64Array
 }
 
 // The compiler throws an error if the function is provided a first argument which is not a string...
@@ -224,14 +224,14 @@ import dspmv = require( './index' );
 	const y = new Float64Array( 10 );
 	const A = new Float64Array( 55 );
 
-	dspmv.ndarray( 10, 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( true, 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( false, 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( null, 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( void 0, 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( [], 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( {}, 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( ( x: number ): number => x, 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 10, 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( true, 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( false, 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( null, 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( void 0, 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( [], 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( {}, 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( ( x: number ): number => x, 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
 }
 
 // The compiler throws an error if the function is provided a second argument which is not a string...
@@ -240,14 +240,14 @@ import dspmv = require( './index' );
 	const y = new Float64Array( 10 );
 	const A = new Float64Array( 55 );
 
-	dspmv.ndarray( 'row-major', 10, 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', true, 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', false, 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', null, 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', void 0, 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', [ '1' ], 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', {}, 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', ( x: number ): number => x, 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 10, 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', true, 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', false, 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', null, 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', void 0, 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', [ '1' ], 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', {}, 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', ( x: number ): number => x, 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
 }
 
 // The compiler throws an error if the function is provided a third argument which is not a number...
@@ -256,14 +256,14 @@ import dspmv = require( './index' );
 	const y = new Float64Array( 10 );
 	const A = new Float64Array( 55 );
 
-	dspmv.ndarray( 'row-major', 'upper', '10', 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', true, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', false, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', null, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', void 0, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', [], 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', {}, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', ( x: number ): number => x, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', '10', 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', true, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', false, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', null, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', void 0, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', [], 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', {}, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', ( x: number ): number => x, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
 }
 
 // The compiler throws an error if the function is provided a fourth argument which is not a number...
@@ -272,14 +272,14 @@ import dspmv = require( './index' );
 	const y = new Float64Array( 10 );
 	const A = new Float64Array( 55 );
 
-	dspmv.ndarray( 'row-major', 'upper', 10, '10', A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, true, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, false, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, null, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, void 0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, [], A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, {}, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, ( x: number ): number => x, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, '10', A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, true, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, false, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, null, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, void 0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, [], A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, {}, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, ( x: number ): number => x, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
 }
 
 // The compiler throws an error if the function is provided a fifth argument which is not a Float64Array...
@@ -287,15 +287,15 @@ import dspmv = require( './index' );
 	const x = new Float64Array( 10 );
 	const y = new Float64Array( 10 );
 
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, 10, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, '10', 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, true, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, false, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, null, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, void 0, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, [ '1' ], 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, {}, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, ( x: number ): number => x, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, 10, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, '10', 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, true, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, false, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, null, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, void 0, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, [ '1' ], 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, {}, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, ( x: number ): number => x, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
 }
 
 // The compiler throws an error if the function is provided a sixth argument which is not a number...
@@ -304,46 +304,46 @@ import dspmv = require( './index' );
 	const y = new Float64Array( 10 );
 	const A = new Float64Array( 55 );
 
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, '10', x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, true, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, false, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, null, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, void 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, [ '1' ], x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, {}, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, ( x: number ): number => x, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, '10', 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, true, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, false, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, null, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, void 0, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, [], 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, {}, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, ( x: number ): number => x, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
 }
 
-// The compiler throws an error if the function is provided a seventh argument which is not a Float64Array...
+// The compiler throws an error if the function is provided a seventh argument which is not a number...
 {
+	const x = new Float64Array( 10 );
 	const y = new Float64Array( 10 );
 	const A = new Float64Array( 55 );
 
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, 10, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, '10', 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, true, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, false, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, null, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, void 0, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, [ '1' ], 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, {}, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, ( x: number ): number => x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, '10', x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, true, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, false, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, null, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, void 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, [], x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, {}, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, ( x: number ): number => x, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
 }
 
-// The compiler throws an error if the function is provided an eighth argument which is not a number...
+// The compiler throws an error if the function is provided an eighth argument which is not a Float64Array...
 {
-	const x = new Float64Array( 10 );
 	const y = new Float64Array( 10 );
 	const A = new Float64Array( 55 );
 
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, '10', 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, true, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, false, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, null, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, void 0, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, [], 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, {}, 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, ( x: number ): number => x, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, 10, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, '10', 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, true, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, false, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, null, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, void 0, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, [ '1' ], 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, {}, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, ( x: number ): number => x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError
 }
 
 // The compiler throws an error if the function is provided a ninth argument which is not a number...
@@ -352,14 +352,14 @@ import dspmv = require( './index' );
 	const y = new Float64Array( 10 );
 	const A = new Float64Array( 55 );
 
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, '10', 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, true, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, false, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, null, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, void 0, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, [], 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, {}, 1.0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, ( x: number ): number => x, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, '10', 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, true, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, false, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, null, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, void 0, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, [], 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, {}, 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, ( x: number ): number => x, 0, 1.0, y, 1, 0 ); // $ExpectError
 }
 
 // The compiler throws an error if the function is provided a tenth argument which is not a number...
@@ -368,46 +368,46 @@ import dspmv = require( './index' );
 	const y = new Float64Array( 10 );
 	const A = new Float64Array( 55 );
 
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, '10', y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, true, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, false, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, null, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, void 0, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, [], y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, {}, y, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, ( x: number ): number => x, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, '10', 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, true, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, false, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, null, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, void 0, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, [], 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, {}, 1.0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, ( x: number ): number => x, 1.0, y, 1, 0 ); // $ExpectError
 }
 
-// The compiler throws an error if the function is provided an eleventh argument which is not a Float64Array...
+// The compiler throws an error if the function is provided an eleventh argument which is not a number...
 {
 	const x = new Float64Array( 10 );
+	const y = new Float64Array( 10 );
 	const A = new Float64Array( 55 );
 
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, 10, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, '10', 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, true, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, false, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, null, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, void 0, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, [ '1' ], 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, {}, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, ( x: number ): number => x, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, '10', y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, true, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, false, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, null, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, void 0, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, [], y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, {}, y, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, ( x: number ): number => x, y, 1, 0 ); // $ExpectError
 }
 
-// The compiler throws an error if the function is provided a twelfth argument which is not a number...
+// The compiler throws an error if the function is provided a twelfth argument which is not a Float64Array...
 {
 	const x = new Float64Array( 10 );
-	const y = new Float64Array( 10 );
 	const A = new Float64Array( 55 );
 
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, '10', 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, true, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, false, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, null, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, void 0, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, [], 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, {}, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, ( x: number ): number => x, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, 10, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, '10', 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, true, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, false, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, null, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, void 0, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, [ '1' ], 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, {}, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, ( x: number ): number => x, 1, 0 ); // $ExpectError
 }
 
 // The compiler throws an error if the function is provided a thirteenth argument which is not a number...
@@ -416,14 +416,30 @@ import dspmv = require( './index' );
 	const y = new Float64Array( 10 );
 	const A = new Float64Array( 55 );
 
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, '10' ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, true ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, false ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, null ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, void 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, [] ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, {} ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, ( x: number ): number => x ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, '10', 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, true, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, false, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, null, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, void 0, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, [], 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, {}, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, ( x: number ): number => x, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a fourteenth argument which is not a number...
+{
+	const x = new Float64Array( 10 );
+	const y = new Float64Array( 10 );
+	const A = new Float64Array( 55 );
+
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, '10' ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, true ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, false ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, null ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, void 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, [] ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, {} ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, ( x: number ): number => x ); // $ExpectError
 }
 
 // The compiler throws an error if the `ndarray` method is provided an unsupported number of arguments...
@@ -438,12 +454,13 @@ import dspmv = require( './index' );
 	dspmv.ndarray( 'row-major', 'upper', 10 ); // $ExpectError
 	dspmv.ndarray( 'row-major', 'upper', 10, 1.0 ); // $ExpectError
 	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1 ); // $ExpectError
-	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0, 10 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1 ); // $ExpectError
+	dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0, 10 ); // $ExpectError
 }
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/examples/index.js b/lib/node_modules/@stdlib/blas/base/dspmv/examples/index.js
index 96199e8b634..48e06d6aee8 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/examples/index.js
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/examples/index.js
@@ -25,11 +25,14 @@ var opts = {
 	'dtype': 'float64'
 };
 
-var N = 3;
+var N = 5;
 var AP = discreteUniform( N * ( N + 1 ) / 2, -10, 10, opts );
 
 var x = discreteUniform( N, -10, 10, opts );
 var y = discreteUniform( N, -10, 10, opts );
 
-dspmv.ndarray( 'row-major', 'upper', N, 1.0, AP, 0, x, 1, 0, 1.0, y, 1, 0 );
+dspmv( 'row-major', 'upper', N, 1.0, AP, x, 1, 1.0, y, 1 );
+console.log( y );
+
+dspmv.ndarray( 'row-major', 'upper', N, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, 1, 0 );
 console.log( y );
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/lib/base.js b/lib/node_modules/@stdlib/blas/base/dspmv/lib/base.js
index 5b07549d6f0..97680b79dd2 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/lib/base.js
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/lib/base.js
@@ -20,14 +20,14 @@
 
 // MODULES //
 
-var dfill = require( '@stdlib/blas/ext/base/dfill' ).ndarray;
-var dscal = require( '@stdlib/blas/base/dscal' ).ndarray;
+var sfill = require( '@stdlib/blas/ext/base/sfill' ).ndarray;
+var sscal = require( '@stdlib/blas/base/sscal' ).ndarray;
 
 
 // MAIN //
 
 /**
-* Performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form.
+* Performs the matrix-vector operation `y = α*A*x + β*y` where `α` and `β` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form.
 *
 * @private
 * @param {string} order - storage layout
@@ -35,7 +35,8 @@ var dscal = require( '@stdlib/blas/base/dscal' ).ndarray;
 * @param {NonNegativeInteger} N - number of elements along each dimension of `A`
 * @param {number} alpha - scalar constant
 * @param {Float64Array} AP - packed form of a symmetric matrix `A`
-* @param {NonNegativeInteger} offsetAP - starting `AP` index
+* @param {integer} strideAP - `AP` stride length
+* @param {NonNegativeInteger} offsetAP - starting index for `AP`
 * @param {Float64Array} x - first input array
 * @param {integer} strideX - `x` stride length
 * @param {NonNegativeInteger} offsetX - starting `x` index
@@ -52,83 +53,81 @@ var dscal = require( '@stdlib/blas/base/dscal' ).ndarray;
 * var x = new Float64Array( [ 1.0, 1.0, 1.0 ] );
 * var y = new Float64Array( [ 1.0, 1.0, 1.0 ] );
 *
-* dspmv( 'column-major', 'lower', 3, 1.0, AP, 0, x, 1, 0, 1.0, y, 1, 0 );
+* dspmv( 'column-major', 'lower', 3, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, 1, 0 );
 * // y => <Float64Array>[ ~7.0, ~12.0, ~15.0 ]
 */
-function dspmv( order, uplo, N, alpha, AP, offsetAP, x, strideX, offsetX, beta, y, strideY, offsetY ) { // eslint-disable-line max-params, max-len
-	var temp1;
-	var temp2;
-	var ix;
-	var iy;
-	var jx;
-	var jy;
+function dspmv( order, uplo, N, alpha, AP, strideAP, offsetAP, x, strideX, offsetX, beta, y, strideY, offsetY ) { // eslint-disable-line max-params, max-len
+	var isrm;
+	var tmp1;
+	var tmp2;
+	var iap;
+	var ix0;
+	var ix1;
+	var iy0;
+	var iy1;
+	var i1;
+	var i0;
 	var kk;
-	var kx;
-	var ky;
-	var j;
-	var k;
+	var ox;
+	var oy;
 
-	if ( N === 0 || ( alpha === 0.0 && beta === 1.0 ) ) {
-		return y;
-	}
-	// Form: y = beta*y
 	if ( beta !== 1.0 ) {
 		if ( beta === 0.0 ) {
-			dfill( N, 0.0, y, strideY, offsetY );
+			sfill( N, beta, y, strideY, offsetY );
 		} else {
-			dscal( N, beta, y, strideY, offsetY );
+			sscal( N, beta, y, strideY, offsetY );
 		}
 	}
 	if ( alpha === 0.0 ) {
 		return y;
 	}
-	// Form: y = alpha*A*x + y
-	kx = offsetX;
-	ky = offsetY;
+	isrm = ( order === 'row-major' );
 	kk = offsetAP;
-	if (
-		( order === 'row-major' && uplo === 'upper' ) ||
-		( order === 'column-major' && uplo === 'lower' )
-	) {
-		jx = kx;
-		jy = ky;
-		for ( j = 0; j < N; j++ ) {
-			temp1 = alpha * x[ jx ];
-			temp2 = 0.0;
-			y[ jy ] += temp1 * AP[ kk ];
-			ix = jx;
-			iy = jy;
-			for ( k = kk + 1; k < kk + N - j; k++ ) {
-				ix += strideX;
-				iy += strideY;
-				y[ iy ] += temp1 * AP[ k ];
-				temp2 += AP[ k ] * x[ ix ];
+	ox = offsetX;
+	oy = offsetY;
+	if ( ( !isrm && uplo === 'upper' ) || ( isrm && uplo === 'lower' ) ) {
+		ix1 = ox;
+		iy1 = oy;
+		for ( i1 = 0; i1 < N; i1++ ) {
+			tmp1 = alpha * x[ ix1 ];
+			tmp2 = 0.0;
+			ix0 = ox;
+			iy0 = oy;
+			iap = kk;
+			for ( i0 = 0; i0 < i1; i0++ ) {
+				y[ iy0 ] += tmp1 * AP[ iap ];
+				tmp2 += AP[ iap ] * x[ ix0 ];
+				ix0 += strideX;
+				iy0 += strideY;
+				iap += strideAP;
 			}
-			y[ jy ] += alpha * temp2;
-			jx += strideX;
-			jy += strideY;
-			kk += N - j;
+			y[ iy1 ] += ( tmp1 * AP[ iap ] ) + ( alpha * tmp2 );
+			ix1 += strideX;
+			iy1 += strideY;
+			kk += ( i1 + 1 ) * strideAP;
 		}
 		return y;
 	}
-	// ( order === 'row-major' && uplo === 'lower') || ( order === 'column-major' && uplo === 'upper' )
-	jx = kx;
-	jy = ky;
-	for ( j = 0; j < N; j++ ) {
-		temp1 = alpha * x[ jx ];
-		temp2 = 0.0;
-		ix = kx;
-		iy = ky;
-		for ( k = kk; k < kk + j; k++ ) {
-			y[ iy ] += temp1 * AP[ k ];
-			temp2 += AP[ k ] * x[ ix ];
-			ix += strideX;
-			iy += strideY;
+	ix1 = ox;
+	iy1 = oy;
+	for ( i1 = 0; i1 < N; i1++ ) {
+		tmp1 = alpha * x[ ix1 ];
+		tmp2 = 0.0;
+		y[ iy1 ] += tmp1 * AP[ kk ];
+		iap = kk;
+		ix0 = ix1;
+		iy0 = iy1;
+		for ( i0 = i1+1; i0 < N; i0++ ) {
+			ix0 += strideX;
+			iy0 += strideY;
+			iap += strideAP;
+			y[ iy0 ] += tmp1 * AP[ iap ];
+			tmp2 += AP[ iap ] * x[ ix0 ];
 		}
-		y[ jy ] += ( temp1 * AP[ kk + j ] ) + ( alpha * temp2 );
-		jx += strideX;
-		jy += strideY;
-		kk += j + 1;
+		y[ iy1 ] += alpha * tmp2;
+		ix1 += strideX;
+		iy1 += strideY;
+		kk += ( N - i1 ) * strideAP;
 	}
 	return y;
 }
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/lib/dspmv.js b/lib/node_modules/@stdlib/blas/base/dspmv/lib/dspmv.js
index 32b865e6331..e9803e9b4c0 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/lib/dspmv.js
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/lib/dspmv.js
@@ -20,9 +20,9 @@
 
 // MODULES //
 
+var stride2offset = require( '@stdlib/strided/base/stride2offset' );
 var isLayout = require( '@stdlib/blas/base/assert/is-layout' );
 var isMatrixTriangle = require( '@stdlib/blas/base/assert/is-matrix-triangle' );
-var stride2offset = require( '@stdlib/strided/base/stride2offset' );
 var format = require( '@stdlib/string/format' );
 var base = require( './base.js' );
 
@@ -30,7 +30,7 @@ var base = require( './base.js' );
 // MAIN //
 
 /**
-* Performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form.
+* Performs the matrix-vector operation `y = α*A*x + β*y` where `α` and `β` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form.
 *
 * @param {string} order - storage layout
 * @param {string} uplo - specifies whether the upper or lower triangular part of the symmetric matrix `A` is supplied
@@ -60,8 +60,8 @@ var base = require( './base.js' );
 * // y => <Float64Array>[ ~7.0, ~12.0, ~15.0 ]
 */
 function dspmv( order, uplo, N, alpha, AP, x, strideX, beta, y, strideY ) {
-	var offsetX;
-	var offsetY;
+	var ox;
+	var oy;
 
 	if ( !isLayout( order ) ) {
 		throw new TypeError( format( 'invalid argument. First argument must be a valid order. Value: `%s`.', order ) );
@@ -78,9 +78,12 @@ function dspmv( order, uplo, N, alpha, AP, x, strideX, beta, y, strideY ) {
 	if ( strideY === 0 ) {
 		throw new RangeError( format( 'invalid argument. Tenth argument must be non-zero. Value: `%d`.', strideY ) );
 	}
-	offsetX = stride2offset( N, strideX );
-	offsetY = stride2offset( N, strideY );
-	return base( order, uplo, N, alpha, AP, 0, x, strideX, offsetX, beta, y, strideY, offsetY ); // eslint-disable-line max-len
+	if ( N === 0 || ( alpha === 0.0 && beta === 1.0 ) ) {
+		return y;
+	}
+	ox = stride2offset( N, strideX );
+	oy = stride2offset( N, strideY );
+	return base( order, uplo, N, alpha, AP, 1, 0, x, strideX, ox, beta, y, strideY, oy ); // eslint-disable-line max-len
 }
 
 
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/lib/index.js b/lib/node_modules/@stdlib/blas/base/dspmv/lib/index.js
index c18159951df..caf718651d2 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/lib/index.js
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/lib/index.js
@@ -19,7 +19,7 @@
 'use strict';
 
 /**
-* BLAS level 2 routine to perform the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form.
+* BLAS level 2 routine to perform the matrix-vector operation `y = α*A*x + β*y` where `α` and `β` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form.
 *
 * @module @stdlib/blas/base/dspmv
 *
@@ -42,7 +42,7 @@
 * var x = new Float64Array( [ 1.0, 1.0, 1.0 ] );
 * var y = new Float64Array( [ 1.0, 1.0, 1.0 ] );
 *
-* dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 0, x, 1, 0, 1.0, y, 1, 0 );
+* dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, 1, 0 );
 * // y => <Float64Array>[ ~7.0, ~12.0, ~15.0 ]
 */
 
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/lib/ndarray.js b/lib/node_modules/@stdlib/blas/base/dspmv/lib/ndarray.js
index 2c5a1479e79..1d12b90dcd3 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/lib/ndarray.js
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/lib/ndarray.js
@@ -29,14 +29,15 @@ var base = require( './base.js' );
 // MAIN //
 
 /**
-* Performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form.
+* Performs the matrix-vector operation `y = α*A*x + β*y` where `α` and `β` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form.
 *
 * @param {string} order - storage layout
 * @param {string} uplo - specifies whether the upper or lower triangular part of the symmetric matrix `A` is supplied
 * @param {NonNegativeInteger} N - number of elements along each dimension of `A`
 * @param {number} alpha - scalar constant
 * @param {Float64Array} AP - packed form of a symmetric matrix `A`
-* @param {NonNegativeInteger} offsetAP - starting `AP` index
+* @param {integer} strideAP - `AP` stride length
+* @param {NonNegativeInteger} offsetAP - starting index for `AP`
 * @param {Float64Array} x - first input array
 * @param {integer} strideX - `x` stride length
 * @param {NonNegativeInteger} offsetX - starting `x` index
@@ -47,8 +48,9 @@ var base = require( './base.js' );
 * @throws {TypeError} first argument must be a valid order
 * @throws {TypeError} second argument must specify whether the lower or upper triangular matrix is supplied
 * @throws {RangeError} third argument must be a nonnegative integer
-* @throws {RangeError} eighth argument must be non-zero
-* @throws {RangeError} twelfth argument must be non-zero
+* @throws {RangeError} sixth argument must be non-zero
+* @throws {RangeError} ninth argument must be non-zero
+* @throws {RangeError} thirteenth argument must be non-zero
 * @returns {Float64Array} `y`
 *
 * @example
@@ -58,10 +60,10 @@ var base = require( './base.js' );
 * var x = new Float64Array( [ 1.0, 1.0, 1.0 ] );
 * var y = new Float64Array( [ 1.0, 1.0, 1.0 ] );
 *
-* dspmv( 'column-major', 'lower', 3, 1.0, AP, 0, x, 1, 0, 1.0, y, 1, 0 );
+* dspmv( 'column-major', 'lower', 3, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, 1, 0 );
 * // y => <Float64Array>[ ~7.0, ~12.0, ~15.0 ]
 */
-function dspmv( order, uplo, N, alpha, AP, offsetAP, x, strideX, offsetX, beta, y, strideY, offsetY ) { // eslint-disable-line max-params, max-len
+function dspmv( order, uplo, N, alpha, AP, strideAP, offsetAP, x, strideX, offsetX, beta, y, strideY, offsetY ) { // eslint-disable-line max-params, max-len
 	if ( !isLayout( order ) ) {
 		throw new TypeError( format( 'invalid argument. First argument must be a valid order. Value: `%s`.', order ) );
 	}
@@ -71,13 +73,19 @@ function dspmv( order, uplo, N, alpha, AP, offsetAP, x, strideX, offsetX, beta,
 	if ( N < 0 ) {
 		throw new RangeError( format( 'invalid argument. Third argument must be a nonnegative integer. Value: `%d`.', N ) );
 	}
+	if ( strideAP === 0 ) {
+		throw new RangeError( format( 'invalid argument. Sixth argument must be non-zero. Value: `%d`.', strideAP ) );
+	}
 	if ( strideX === 0 ) {
-		throw new RangeError( format( 'invalid argument. Eighth argument must be non-zero. Value: `%d`.', strideX ) );
+		throw new RangeError( format( 'invalid argument. Ninth argument must be non-zero. Value: `%d`.', strideX ) );
 	}
 	if ( strideY === 0 ) {
-		throw new RangeError( format( 'invalid argument. Twelfth argument must be non-zero. Value: `%d`.', strideY ) );
+		throw new RangeError( format( 'invalid argument. Thirteenth argument must be non-zero. Value: `%d`.', strideY ) );
+	}
+	if ( N === 0 || ( alpha === 0.0 && beta === 1.0 ) ) {
+		return y;
 	}
-	return base( order, uplo, N, alpha, AP, offsetAP, x, strideX, offsetX, beta, y, strideY, offsetY ); // eslint-disable-line max-len
+	return base( order, uplo, N, alpha, AP, strideAP, offsetAP, x, strideX, offsetX, beta, y, strideY, offsetY ); // eslint-disable-line max-len
 }
 
 
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/package.json b/lib/node_modules/@stdlib/blas/base/dspmv/package.json
index a3796aeb02c..2226291ced4 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/package.json
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/package.json
@@ -1,7 +1,7 @@
 {
   "name": "@stdlib/blas/base/dspmv",
   "version": "0.0.0",
-  "description": "Perform the matrix-vector operation `y = α*A*x + β*y` where `α` and `β` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form.",
+  "description": "Perform the matrix-vector operation `y = α*A*x + β*y`.",
   "license": "Apache-2.0",
   "author": {
     "name": "The Stdlib Authors",
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_complex_access_pattern.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_complex_access_pattern.json
new file mode 100644
index 00000000000..8fd5802c13c
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_complex_access_pattern.json
@@ -0,0 +1,17 @@
+{
+  "order": "column-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "beta": 1.0,
+  "strideX": -1,
+  "offsetX": 2,
+  "strideY": -2,
+  "offsetY": 4,
+  "AP": [ 6.0, 5.0, 3.0, 4.0, 2.0, 1.0 ],
+  "strideAP": -1,
+  "offsetAP": 5,
+  "x": [ 3.0, 2.0, 1.0 ],
+  "y": [ 1.0, 0.0, 1.0, 0.0, 1.0, 0.0 ],
+  "y_out": [ 32.0, 0.0, 26.0, 0.0, 15.0, 0.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xoyt.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_lower.json
similarity index 95%
rename from lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xoyt.json
rename to lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_lower.json
index b05e6fe81a6..37c7121ef6d 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xoyt.json
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_lower.json
@@ -1,6 +1,6 @@
 {
-  "uplo": "lower",
   "order": "column-major",
+  "uplo": "lower",
   "N": 3,
   "alpha": 1.0,
   "beta": 1.0,
@@ -8,8 +8,9 @@
   "offsetX": 0,
   "strideY": 2,
   "offsetY": 0,
-  "offsetAP": 0,
   "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ],
+  "strideAP": 1,
+  "offsetAP": 0,
   "x": [ 1.0, 1.0, 1.0 ],
   "y": [ 1.0, 0.0, 1.0, 0.0, 1.0, 0.0 ],
   "y_out": [ 7.0, 0.0, 12.0, 0.0, 15.0, 0.0 ]
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_oap.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_oap.json
new file mode 100644
index 00000000000..b54652dadca
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_oap.json
@@ -0,0 +1,17 @@
+{
+  "order": "column-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "beta": 1.0,
+  "strideX": 1,
+  "offsetX": 0,
+  "strideY": 2,
+  "offsetY": 0,
+  "AP": [ 0.0, 1.0, 0.0, 2.0, 0.0, 4.0, 0.0, 3.0, 0.0, 5.0, 0.0, 6.0 ],
+  "strideAP": 2,
+  "offsetAP": 1,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "y": [ 1.0, 0.0, 1.0, 0.0, 1.0, 0.0 ],
+  "y_out": [ 15.0, 0.0, 26.0, 0.0, 32.0, 0.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_ox.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_ox.json
new file mode 100644
index 00000000000..2dc0f16614b
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_ox.json
@@ -0,0 +1,17 @@
+{
+  "order": "column-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "beta": 1.0,
+  "strideX": 1,
+  "offsetX": 2,
+  "strideY": 2,
+  "offsetY": 0,
+  "AP": [ 1.0, 2.0, 4.0, 3.0, 5.0, 6.0 ],
+  "strideAP": 1,
+  "offsetAP": 0,
+  "x": [ 0.0, 0.0, 1.0, 2.0, 3.0 ],
+  "y": [ 1.0, 0.0, 1.0, 0.0, 1.0, 0.0 ],
+  "y_out": [ 15.0, 0.0, 26.0, 0.0, 32.0, 0.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_oy.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_oy.json
new file mode 100644
index 00000000000..7a90138b116
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_oy.json
@@ -0,0 +1,17 @@
+{
+  "order": "column-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "beta": 1.0,
+  "strideX": 1,
+  "offsetX": 0,
+  "strideY": 2,
+  "offsetY": 1,
+  "AP": [ 1.0, 2.0, 4.0, 3.0, 5.0, 6.0 ],
+  "strideAP": 1,
+  "offsetAP": 0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "y": [ 0.0, 1.0, 0.0, 1.0, 0.0, 1.0, 0.0 ],
+  "y_out": [ 0.0, 15.0, 0.0, 26.0, 0.0, 32.0, 0.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_sap.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_sap.json
new file mode 100644
index 00000000000..4ba7aad930a
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_sap.json
@@ -0,0 +1,17 @@
+{
+  "order": "column-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "beta": 1.0,
+  "strideX": 1,
+  "offsetX": 0,
+  "strideY": 2,
+  "offsetY": 0,
+  "AP": [ 1.0, 0.0, 2.0, 0.0, 4.0, 0.0, 3.0, 0.0, 5.0, 0.0, 6.0 ],
+  "strideAP": 2,
+  "offsetAP": 0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "y": [ 1.0, 0.0, 1.0, 0.0, 1.0, 0.0 ],
+  "y_out": [ 15.0, 0.0, 26.0, 0.0, 32.0, 0.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_sapn.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_sapn.json
new file mode 100644
index 00000000000..410bec99c5e
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_sapn.json
@@ -0,0 +1,17 @@
+{
+  "order": "column-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "beta": 1.0,
+  "strideX": 1,
+  "offsetX": 0,
+  "strideY": 2,
+  "offsetY": 0,
+  "AP": [ 6.0, 5.0, 3.0, 4.0, 2.0, 1.0 ],
+  "strideAP": -1,
+  "offsetAP": 5,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "y": [ 1.0, 0.0, 1.0, 0.0, 1.0, 0.0 ],
+  "y_out": [ 15.0, 0.0, 26.0, 0.0, 32.0, 0.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_upper.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_upper.json
new file mode 100644
index 00000000000..f2965d50f6c
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_upper.json
@@ -0,0 +1,17 @@
+{
+  "order": "column-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "beta": 1.0,
+  "strideX": 1,
+  "offsetX": 0,
+  "strideY": 2,
+  "offsetY": 0,
+  "AP": [ 1.0, 2.0, 4.0, 3.0, 5.0, 6.0 ],
+  "strideAP": 1,
+  "offsetAP": 0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "y": [ 1.0, 0.0, 1.0, 0.0, 1.0, 0.0 ],
+  "y_out": [ 15.0, 0.0, 26.0, 0.0, 32.0, 0.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xnyn.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xnyn.json
index 06b6850679f..07cc43b13d8 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xnyn.json
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xnyn.json
@@ -8,8 +8,9 @@
   "offsetX": 2,
   "strideY": -1,
   "offsetY": 2,
-  "offsetAP": 2,
-  "AP": [ 0.0, 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ],
+  "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ],
+  "strideAP": 1,
+  "offsetAP": 0,
   "x": [ 1.0, 1.0, 1.0 ],
   "y": [ 1.0, 1.0, 1.0 ],
   "y_out": [ 16.0, 11.0, 8.0 ]
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xnyp.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xnyp.json
index c2d83166074..1807dbf1b0b 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xnyp.json
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xnyp.json
@@ -8,8 +8,9 @@
   "offsetX": 2,
   "strideY": 1,
   "offsetY": 0,
-  "offsetAP": 0,
   "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ],
+  "strideAP": 1,
+  "offsetAP": 0,
   "x": [ 1.0, 1.0, 1.0 ],
   "y": [ 1.0, 1.0, 1.0 ],
   "y_out": [ 8.0, 11.0, 16.0 ]
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xpyn.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xpyn.json
index f4aad6d8b84..4e08b89290a 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xpyn.json
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xpyn.json
@@ -8,8 +8,9 @@
   "offsetX": 0,
   "strideY": -1,
   "offsetY": 2,
-  "offsetAP": 0,
   "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ],
+  "strideAP": 1,
+  "offsetAP": 0,
   "x": [ 1.0, 1.0, 1.0 ],
   "y": [ 1.0, 1.0, 1.0 ],
   "y_out": [ 16.0, 11.0, 8.0 ]
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xpyp.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xpyp.json
index b3123f8e2f1..09d16df7ab7 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xpyp.json
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xpyp.json
@@ -8,8 +8,9 @@
   "offsetX": 0,
   "strideY": 1,
   "offsetY": 0,
-  "offsetAP": 0,
   "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ],
+  "strideAP": 1,
+  "offsetAP": 0,
   "x": [ 1.0, 1.0, 1.0 ],
   "y": [ 1.0, 1.0, 1.0 ],
   "y_out": [ 8.0, 11.0, 16.0 ]
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_complex_access_pattern.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_complex_access_pattern.json
new file mode 100644
index 00000000000..88b89e2cc53
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_complex_access_pattern.json
@@ -0,0 +1,17 @@
+{
+  "order": "row-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "beta": 1.0,
+  "strideX": -1,
+  "offsetX": 2,
+  "strideY": -1,
+  "offsetY": 2,
+  "AP": [ 6.0, 5.0, 4.0, 3.0, 2.0, 1.0 ],
+  "strideAP": -1,
+  "offsetAP": 5,
+  "x": [ 3.0, 2.0, 1.0 ],
+  "y": [ 1.0, 1.0, 1.0 ],
+  "y_out": [ 32.0, 26.0, 15.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xoyt.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_lower.json
similarity index 77%
rename from lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xoyt.json
rename to lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_lower.json
index 211518d861c..a393e4279a6 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xoyt.json
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_lower.json
@@ -1,16 +1,17 @@
 {
-  "uplo": "lower",
   "order": "row-major",
+  "uplo": "lower",
   "N": 3,
   "alpha": 1.0,
-  "beta": 1.0,
+  "beta": 2.0,
   "strideX": 1,
   "offsetX": 0,
   "strideY": 2,
   "offsetY": 0,
-  "offsetAP": 0,
   "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ],
+  "strideAP": 1,
+  "offsetAP": 0,
   "x": [ 1.0, 1.0, 1.0 ],
   "y": [ 1.0, 0.0, 1.0, 0.0, 1.0, 0.0 ],
-  "y_out": [ 8.0, 0.0, 11.0, 0.0, 16.0, 0.0 ]
+  "y_out": [ 9.0, 0.0, 12.0, 0.0, 17.0, 0.0 ]
 }
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_oap.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_oap.json
new file mode 100644
index 00000000000..fb34b7e18d1
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_oap.json
@@ -0,0 +1,17 @@
+{
+  "order": "row-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "beta": 1.0,
+  "strideX": 1,
+  "offsetX": 0,
+  "strideY": 1,
+  "offsetY": 0,
+  "AP": [ 0.0, 1.0, 0.0, 2.0, 0.0, 3.0, 0.0, 4.0, 0.0, 5.0, 0.0, 6.0 ],
+  "strideAP": 2,
+  "offsetAP": 1,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "y": [ 1.0, 1.0, 1.0 ],
+  "y_out": [ 15.0, 26.0, 32.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_ox.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_ox.json
new file mode 100644
index 00000000000..b39813d836c
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_ox.json
@@ -0,0 +1,17 @@
+{
+  "order": "row-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "beta": 1.0,
+  "strideX": 1,
+  "offsetX": 2,
+  "strideY": 1,
+  "offsetY": 0,
+  "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ],
+  "strideAP": 1,
+  "offsetAP": 0,
+  "x": [ 0.0, 0.0, 1.0, 2.0, 3.0 ],
+  "y": [ 1.0, 1.0, 1.0 ],
+  "y_out": [ 15.0, 26.0, 32.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_oy.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_oy.json
new file mode 100644
index 00000000000..660a96d15dd
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_oy.json
@@ -0,0 +1,17 @@
+{
+  "order": "row-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "beta": 1.0,
+  "strideX": 1,
+  "offsetX": 0,
+  "strideY": 1,
+  "offsetY": 2,
+  "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ],
+  "strideAP": 1,
+  "offsetAP": 0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "y": [ 0.0, 0.0, 1.0, 1.0, 1.0 ],
+  "y_out": [ 0.0, 0.0, 15.0, 26.0, 32.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_sap.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_sap.json
new file mode 100644
index 00000000000..b827f7fa9a4
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_sap.json
@@ -0,0 +1,17 @@
+{
+  "order": "row-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "beta": 1.0,
+  "strideX": 1,
+  "offsetX": 0,
+  "strideY": 1,
+  "offsetY": 0,
+  "AP": [ 1.0, 0.0, 2.0, 0.0, 3.0, 0.0, 4.0, 0.0, 5.0, 0.0, 6.0 ],
+  "strideAP": 2,
+  "offsetAP": 0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "y": [ 1.0, 1.0, 1.0 ],
+  "y_out": [ 15.0, 26.0, 32.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_sapn.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_sapn.json
new file mode 100644
index 00000000000..b73e985d838
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_sapn.json
@@ -0,0 +1,17 @@
+{
+  "order": "row-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "beta": 1.0,
+  "strideX": 1,
+  "offsetX": 0,
+  "strideY": 1,
+  "offsetY": 0,
+  "AP": [ 6.0, 5.0, 4.0, 3.0, 2.0, 1.0 ],
+  "strideAP": -1,
+  "offsetAP": 5,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "y": [ 1.0, 1.0, 1.0 ],
+  "y_out": [ 15.0, 26.0, 32.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_upper.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_upper.json
new file mode 100644
index 00000000000..d9fc2e5f7d0
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_upper.json
@@ -0,0 +1,17 @@
+{
+  "order": "row-major",
+  "uplo": "upper",
+  "N": 3,
+  "alpha": 1.0,
+  "beta": 1.0,
+  "strideX": 1,
+  "offsetX": 0,
+  "strideY": 1,
+  "offsetY": 0,
+  "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ],
+  "strideAP": 1,
+  "offsetAP": 0,
+  "x": [ 1.0, 2.0, 3.0 ],
+  "y": [ 1.0, 1.0, 1.0 ],
+  "y_out": [ 15.0, 26.0, 32.0 ]
+}
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xnyn.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xnyn.json
index f1631104ae5..d9aacf91386 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xnyn.json
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xnyn.json
@@ -8,8 +8,9 @@
   "offsetX": 2,
   "strideY": -1,
   "offsetY": 2,
-  "offsetAP": 2,
-  "AP": [ 0.0, 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ],
+  "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ],
+  "strideAP": 1,
+  "offsetAP": 0,
   "x": [ 1.0, 1.0, 1.0 ],
   "y": [ 1.0, 1.0, 1.0 ],
   "y_out": [ 15.0, 12.0, 7.0 ]
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xnyp.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xnyp.json
index 475ee4d4919..ca868017cb2 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xnyp.json
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xnyp.json
@@ -8,8 +8,9 @@
   "offsetX": 2,
   "strideY": 1,
   "offsetY": 0,
-  "offsetAP": 0,
   "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ],
+  "strideAP": 1,
+  "offsetAP": 0,
   "x": [ 1.0, 1.0, 1.0 ],
   "y": [ 1.0, 1.0, 1.0 ],
   "y_out": [ 7.0, 12.0, 15.0 ]
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xpyn.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xpyn.json
index 2a6711742d6..66696f9f69e 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xpyn.json
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xpyn.json
@@ -8,8 +8,9 @@
   "offsetX": 0,
   "strideY": -1,
   "offsetY": 2,
-  "offsetAP": 0,
   "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ],
+  "strideAP": 1,
+  "offsetAP": 0,
   "x": [ 1.0, 1.0, 1.0 ],
   "y": [ 1.0, 1.0, 1.0 ],
   "y_out": [ 15.0, 12.0, 7.0 ]
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xpyp.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xpyp.json
index ddee2b15624..4878a217089 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xpyp.json
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xpyp.json
@@ -8,8 +8,9 @@
   "offsetX": 0,
   "strideY": 1,
   "offsetY": 0,
-  "offsetAP": 0,
   "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ],
+  "strideAP": 1,
+  "offsetAP": 0,
   "x": [ 1.0, 1.0, 1.0 ],
   "y": [ 1.0, 1.0, 1.0 ],
   "y_out": [ 7.0, 12.0, 15.0 ]
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/test.dspmv.js b/lib/node_modules/@stdlib/blas/base/dspmv/test/test.dspmv.js
index 7328a8656fe..42a63e61897 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/test/test.dspmv.js
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/test.dspmv.js
@@ -24,56 +24,27 @@
 
 var tape = require( 'tape' );
 var Float64Array = require( '@stdlib/array/float64' );
-var EPS = require( '@stdlib/constants/float64/eps' );
-var abs = require( '@stdlib/math/base/special/abs' );
 var ones = require( '@stdlib/array/ones' );
 var dspmv = require( './../lib/dspmv.js' );
 
 
 // FIXTURES //
 
-var rxoyt = require( './fixtures/row_major_xoyt.json' );
+var rl = require( './fixtures/row_major_lower.json' );
+var ru = require( './fixtures/row_major_upper.json' );
 var rxpyp = require( './fixtures/row_major_xpyp.json' );
 var rxnyp = require( './fixtures/row_major_xnyp.json' );
 var rxpyn = require( './fixtures/row_major_xpyn.json' );
 var rxnyn = require( './fixtures/row_major_xnyn.json' );
 
-var cxoyt = require( './fixtures/column_major_xoyt.json' );
+var cl = require( './fixtures/column_major_lower.json' );
+var cu = require( './fixtures/column_major_upper.json' );
 var cxpyp = require( './fixtures/column_major_xpyp.json' );
 var cxnyp = require( './fixtures/column_major_xnyp.json' );
 var cxpyn = require( './fixtures/column_major_xpyn.json' );
 var cxnyn = require( './fixtures/column_major_xnyn.json' );
 
 
-// FUNCTIONS //
-
-/**
-* Tests for element-wise approximate equality.
-*
-* @private
-* @param {Object} t - test object
-* @param {Collection} actual - actual values
-* @param {Collection} expected - expected values
-* @param {number} rtol - relative tolerance
-*/
-function isApprox( t, actual, expected, rtol ) {
-	var delta;
-	var tol;
-	var i;
-
-	t.strictEqual( actual.length, expected.length, 'returns expected value' );
-	for ( i = 0; i < expected.length; i++ ) {
-		if ( actual[ i ] === expected[ i ] ) {
-			t.strictEqual( actual[ i ], expected[ i ], 'returns expected value' );
-		} else {
-			delta = abs( actual[ i ] - expected[ i ] );
-			tol = rtol * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. actual: '+actual[ i ]+'. expected: '+expected[ i ]+'. delta: '+delta+'. tol: '+tol+'.' );
-		}
-	}
-}
-
-
 // TESTS //
 
 tape( 'main export is a function', function test( t ) {
@@ -89,8 +60,11 @@ tape( 'the function has an arity of 10', function test( t ) {
 
 tape( 'the function throws an error if provided an invalid first argument', function test( t ) {
 	var values;
+	var data;
 	var i;
 
+	data = rl;
+
 	values = [
 		'foo',
 		'bar',
@@ -105,15 +79,18 @@ tape( 'the function throws an error if provided an invalid first argument', func
 
 	function badValue( value ) {
 		return function badValue() {
-			dspmv( value, rxpyp.uplo, rxpyp.N, rxpyp.alpha, new Float64Array( rxpyp.AP ), new Float64Array( rxpyp.x ), rxpyp.strideX, rxpyp.beta, new Float64Array( rxpyp.y ), rxpyp.strideY );
+			dspmv( value, data.uplo, data.N, data.alpha, new Float64Array( data.AP ), new Float64Array( data.x ), data.strideX, data.beta, new Float64Array( data.y ), data.strideY );
 		};
 	}
 });
 
 tape( 'the function throws an error if provided an invalid second argument', function test( t ) {
 	var values;
+	var data;
 	var i;
 
+	data = rl;
+
 	values = [
 		'foo',
 		'bar',
@@ -128,15 +105,18 @@ tape( 'the function throws an error if provided an invalid second argument', fun
 
 	function badValue( value ) {
 		return function badValue() {
-			dspmv( rxpyp.order, value, rxpyp.N, rxpyp.alpha, new Float64Array( rxpyp.AP ), new Float64Array( rxpyp.x ), rxpyp.strideX, rxpyp.beta, new Float64Array( rxpyp.y ), rxpyp.strideY );
+			dspmv( data.order, value, data.N, data.alpha, new Float64Array( data.AP ), new Float64Array( data.x ), data.strideX, data.beta, new Float64Array( data.y ), data.strideY );
 		};
 	}
 });
 
 tape( 'the function throws an error if provided an invalid third argument', function test( t ) {
 	var values;
+	var data;
 	var i;
 
+	data = rl;
+
 	values = [
 		-1,
 		-2,
@@ -150,15 +130,18 @@ tape( 'the function throws an error if provided an invalid third argument', func
 
 	function badValue( value ) {
 		return function badValue() {
-			dspmv( rxpyp.order, rxpyp.uplo, value, rxpyp.alpha, new Float64Array( rxpyp.AP ), new Float64Array( rxpyp.x ), rxpyp.strideX, rxpyp.beta, new Float64Array( rxpyp.y ), rxpyp.strideY );
+			dspmv( data.order, data.uplo, value, data.alpha, new Float64Array( data.AP ), new Float64Array( data.x ), data.strideX, data.beta, new Float64Array( data.y ), data.strideY );
 		};
 	}
 });
 
 tape( 'the function throws an error if provided an invalid seventh argument', function test( t ) {
 	var values;
+	var data;
 	var i;
 
+	data = rl;
+
 	values = [
 		0
 	];
@@ -170,15 +153,18 @@ tape( 'the function throws an error if provided an invalid seventh argument', fu
 
 	function badValue( value ) {
 		return function badValue() {
-			dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, new Float64Array( rxpyp.AP ), new Float64Array( rxpyp.x ), value, rxpyp.beta, new Float64Array( rxpyp.y ), rxpyp.strideY );
+			dspmv( data.order, data.uplo, data.N, data.alpha, new Float64Array( data.AP ), new Float64Array( data.x ), value, data.beta, new Float64Array( data.y ), data.strideY );
 		};
 	}
 });
 
 tape( 'the function throws an error if provided an invalid tenth argument', function test( t ) {
 	var values;
+	var data;
 	var i;
 
+	data = rl;
+
 	values = [
 		0
 	];
@@ -190,375 +176,460 @@ tape( 'the function throws an error if provided an invalid tenth argument', func
 
 	function badValue( value ) {
 		return function badValue() {
-			dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, new Float64Array( rxpyp.AP ), new Float64Array( rxpyp.x ), rxpyp.strideX, rxpyp.beta, new Float64Array( rxpyp.y ), value );
+			dspmv( data.order, data.uplo, data.N, data.alpha, new Float64Array( data.AP ), new Float64Array( data.x ), data.strideX, data.beta, new Float64Array( data.y ), value );
 		};
 	}
 });
 
-tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form (row-major, sx=1, sy=1)', function test( t ) {
+tape( 'the function performs the matrix-vector operation `y = α*A*x + β*y (row-major, lower)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( rxpyp.AP );
-	x = new Float64Array( rxpyp.x );
-	y = new Float64Array( rxpyp.y );
+	data = rl;
 
-	expected = new Float64Array( rxpyp.y_out );
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	out = dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, ap, x, rxpyp.strideX, rxpyp.beta, y, rxpyp.strideY );
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, x, data.strideX, data.beta, y, data.strideY );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form (column-major, sx=1, sy=1)', function test( t ) {
+tape( 'the function performs the matrix-vector operation `y = α*A*x + β*y (column-major, lower)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( cxpyp.AP );
-	x = new Float64Array( cxpyp.x );
-	y = new Float64Array( cxpyp.y );
+	data = cl;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	expected = new Float64Array( cxpyp.y_out );
+	expected = new Float64Array( data.y_out );
 
-	out = dspmv( cxpyp.order, cxpyp.uplo, cxpyp.N, cxpyp.alpha, ap, x, cxpyp.strideX, cxpyp.beta, y, cxpyp.strideY );
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, x, data.strideX, data.beta, y, data.strideY );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form (row-major, sx=1, sy=2)', function test( t ) {
+tape( 'the function performs the matrix-vector operation `y = α*A*x + β*y (row-major, upper)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( rxoyt.AP );
-	x = new Float64Array( rxoyt.x );
-	y = new Float64Array( rxoyt.y );
+	data = ru;
 
-	expected = new Float64Array( rxoyt.y_out );
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	out = dspmv( rxoyt.order, rxoyt.uplo, rxoyt.N, rxoyt.alpha, ap, x, rxoyt.strideX, rxoyt.beta, y, rxoyt.strideY );
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, x, data.strideX, data.beta, y, data.strideY );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form (column-major, sx=1, sy=2)', function test( t ) {
+tape( 'the function performs the matrix-vector operation `y = α*A*x + β*y (column-major, upper)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( cxoyt.AP );
-	x = new Float64Array( cxoyt.x );
-	y = new Float64Array( cxoyt.y );
+	data = cu;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	expected = new Float64Array( cxoyt.y_out );
+	expected = new Float64Array( data.y_out );
 
-	out = dspmv( cxoyt.order, cxoyt.uplo, cxoyt.N, cxoyt.alpha, ap, x, cxoyt.strideX, cxoyt.beta, y, cxoyt.strideY );
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, x, data.strideX, data.beta, y, data.strideY );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form (row-major, sx=1, sy=-1)', function test( t ) {
-	var expected;
+tape( 'the function returns a reference to the second input vector', function test( t ) {
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( rxpyn.AP );
-	x = new Float64Array( rxpyn.x );
-	y = new Float64Array( rxpyn.y );
+	data = rl;
 
-	expected = new Float64Array( rxpyn.y_out );
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	out = dspmv( rxpyn.order, rxpyn.uplo, rxpyn.N, rxpyn.alpha, ap, x, rxpyn.strideX, rxpyn.beta, y, rxpyn.strideY );
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, x, data.strideX, data.beta, y, data.strideY );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
 
 	t.end();
 });
 
-tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form (column-major, sx=1, sy=-1)', function test( t ) {
+tape( 'if `N` is zero or the scalar constants are zero and unity, respectively, the function returns the second input vector unchanged (row-major)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( cxpyn.AP );
-	x = new Float64Array( cxpyn.x );
-	y = new Float64Array( cxpyn.y );
+	data = rl;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	expected = new Float64Array( cxpyn.y_out );
+	expected = new Float64Array( data.y );
+
+	out = dspmv( data.order, data.uplo, 0, data.alpha, ap, x, data.strideX, data.beta, y, data.strideY );
+	t.strictEqual( out, y, 'returns expected value' );
+	t.deepEqual( y, expected, 'returns expected value' );
 
-	out = dspmv( cxpyn.order, cxpyn.uplo, cxpyn.N, cxpyn.alpha, ap, x, cxpyn.strideX, cxpyn.beta, y, cxpyn.strideY );
+	out = dspmv( data.order, data.uplo, data.N, 0.0, ap, x, data.strideX, 1.0, y, data.strideY );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
+	t.deepEqual( y, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form (row-major, sx=-1, sy=1)', function test( t ) {
+tape( 'if `N` is zero or the scalar constants are zero and unity, respectively, the function returns the second input vector unchanged (column-major)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( rxnyp.AP );
-	x = new Float64Array( rxnyp.x );
-	y = new Float64Array( rxnyp.y );
+	data = cl;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	expected = new Float64Array( rxnyp.y_out );
+	expected = new Float64Array( data.y );
 
-	out = dspmv( rxnyp.order, rxnyp.uplo, rxnyp.N, rxnyp.alpha, ap, x, rxnyp.strideX, rxnyp.beta, y, rxnyp.strideY );
+	out = dspmv( data.order, data.uplo, 0, data.alpha, ap, x, data.strideX, data.beta, y, data.strideY );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
+	t.deepEqual( y, expected, 'returns expected value' );
+
+	out = dspmv( data.order, data.uplo, data.N, 0.0, ap, x, data.strideX, 1.0, y, data.strideY );
+	t.strictEqual( out, y, 'returns expected value' );
+	t.deepEqual( y, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `N` element vectors, and `A` is an `N` by `N` symmetric matrix supplied in packed form (column-major, sx=-1, sy=1)', function test( t ) {
+tape( 'when `α` is zero, the function scales the second input vector (row-major, upper)', function test( t ) {
 	var expected;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( cxnyp.AP );
-	x = new Float64Array( cxnyp.x );
-	y = new Float64Array( cxnyp.y );
+	ap = ones( 6, 'float32' );
+	x = ones( 3, 'float32' );
+	y = ones( 3, 'float32' );
 
-	expected = new Float64Array( cxnyp.y_out );
+	expected = new Float64Array( [ 5.0, 5.0, 5.0 ] );
 
-	out = dspmv( cxnyp.order, cxnyp.uplo, cxnyp.N, cxnyp.alpha, ap, x, cxnyp.strideX, cxnyp.beta, y, cxnyp.strideY );
+	out = dspmv( 'row-major', 'upper', 3, 0.0, ap, x, 1, 5.0, y, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
-
-	t.end();
-});
-
-tape( 'the function returns a reference to the second input vector', function test( t ) {
-	var out;
-	var ap;
-	var x;
-	var y;
+	t.deepEqual( y, expected, 'returns expected value' );
 
-	ap = new Float64Array( rxpyp.AP );
-	x = new Float64Array( rxpyp.x );
-	y = new Float64Array( rxpyp.y );
+	expected = new Float64Array( [ 0.0, 0.0, 0.0 ] );
 
-	out = dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, ap, x, rxpyp.strideX, rxpyp.beta, y, rxpyp.strideY );
+	out = dspmv( 'row-major', 'upper', 3, 0.0, ap, x, 1, 0.0, y, -1 );
 	t.strictEqual( out, y, 'returns expected value' );
+	t.deepEqual( y, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'if `N` is zero or the scalar constants are zero and unity, respectively, the function returns the second input vector unchanged (row-major)', function test( t ) {
+tape( 'when `α` is zero, the function scales the second input vector (row-major, lower)', function test( t ) {
 	var expected;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( rxpyp.AP );
-	x = new Float64Array( rxpyp.x );
-	y = new Float64Array( rxpyp.y );
+	ap = ones( 6, 'float32' );
+	x = ones( 3, 'float32' );
+	y = ones( 3, 'float32' );
 
-	expected = new Float64Array( rxpyp.y );
+	expected = new Float64Array( [ 5.0, 5.0, 5.0 ] );
 
-	out = dspmv( rxpyp.order, rxpyp.uplo, 0, rxpyp.alpha, ap, x, rxpyp.strideX, rxpyp.beta, y, rxpyp.strideY );
+	out = dspmv( 'row-major', 'lower', 3, 0.0, ap, x, 1, 5.0, y, -1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.deepEqual( y, expected, 'returns expected value' );
 
-	out = dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, 0.0, ap, x, rxpyp.strideX, 1.0, y, rxpyp.strideY );
+	expected = new Float64Array( [ 0.0, 0.0, 0.0 ] );
+
+	out = dspmv( 'row-major', 'lower', 3, 0.0, ap, x, 1, 0.0, y, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.deepEqual( y, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'if `N` is zero or the scalar constants are zero and unity, respectively, the function returns the second input vector unchanged (column-major)', function test( t ) {
+tape( 'when `α` is zero, the function scales the second input vector (column-major, upper)', function test( t ) {
 	var expected;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( cxpyp.AP );
-	x = new Float64Array( cxpyp.x );
-	y = new Float64Array( cxpyp.y );
+	ap = ones( 6, 'float32' );
+	x = ones( 3, 'float32' );
+	y = ones( 3, 'float32' );
 
-	expected = new Float64Array( cxpyp.y );
+	expected = new Float64Array( [ 5.0, 5.0, 5.0 ] );
 
-	out = dspmv( cxpyp.order, cxpyp.uplo, 0, cxpyp.alpha, ap, x, cxpyp.strideX, cxpyp.beta, y, cxpyp.strideY );
+	out = dspmv( 'column-major', 'upper', 3, 0.0, ap, x, 1, 5.0, y, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.deepEqual( y, expected, 'returns expected value' );
 
-	out = dspmv( cxpyp.order, cxpyp.uplo, cxpyp.N, 0.0, ap, x, cxpyp.strideX, 1.0, y, cxpyp.strideY );
+	expected = new Float64Array( [ 0.0, 0.0, 0.0 ] );
+
+	out = dspmv( 'column-major', 'upper', 3, 0.0, ap, x, 1, 0.0, y, -1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.deepEqual( y, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'when `α` is zero, the function scales the second input vector (row-major, upper)', function test( t ) {
+tape( 'when `α` is zero, the function scales the second input vector (column-major, lower)', function test( t ) {
 	var expected;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = ones( 6, 'float64' );
-	x = ones( 3, 'float64' );
-	y = ones( 3, 'float64' );
+	ap = ones( 6, 'float32' );
+	x = ones( 3, 'float32' );
+	y = ones( 3, 'float32' );
 
 	expected = new Float64Array( [ 5.0, 5.0, 5.0 ] );
 
-	out = dspmv( 'row-major', 'upper', 3, 0.0, ap, x, 1, 5.0, y, 1 );
+	out = dspmv( 'column-major', 'lower', 3, 0.0, ap, x, 1, 5.0, y, -1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.deepEqual( y, expected, 'returns expected value' );
 
 	expected = new Float64Array( [ 0.0, 0.0, 0.0 ] );
 
-	out = dspmv( 'row-major', 'upper', 3, 0.0, ap, x, 1, 0.0, y, -1 );
+	out = dspmv( 'column-major', 'lower', 3, 0.0, ap, x, 1, 0.0, y, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.deepEqual( y, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'when `α` is zero, the function scales the second input vector (row-major, lower)', function test( t ) {
+tape( 'the function supports specifying `x` and `y` strides (row-major)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = ones( 6, 'float64' );
-	x = ones( 3, 'float64' );
-	y = ones( 3, 'float64' );
+	data = rxpyp;
 
-	expected = new Float64Array( [ 5.0, 5.0, 5.0 ] );
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	out = dspmv( 'row-major', 'lower', 3, 0.0, ap, x, 1, 5.0, y, -1 );
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, x, data.strideX, data.beta, y, data.strideY );
 	t.strictEqual( out, y, 'returns expected value' );
-	t.deepEqual( y, expected, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
 
-	expected = new Float64Array( [ 0.0, 0.0, 0.0 ] );
+	t.end();
+});
 
-	out = dspmv( 'row-major', 'lower', 3, 0.0, ap, x, 1, 0.0, y, 1 );
+tape( 'the function supports specifying `x` and `y` strides (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var ap;
+	var x;
+	var y;
+
+	data = cxpyp;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
+
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, x, data.strideX, data.beta, y, data.strideY );
 	t.strictEqual( out, y, 'returns expected value' );
-	t.deepEqual( y, expected, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'when `α` is zero, the function scales the second input vector (column-major, upper)', function test( t ) {
+tape( 'the function supports specifying a negative `y` stride (row-major)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = ones( 6, 'float64' );
-	x = ones( 3, 'float64' );
-	y = ones( 3, 'float64' );
+	data = rxpyn;
 
-	expected = new Float64Array( [ 5.0, 5.0, 5.0 ] );
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	out = dspmv( 'column-major', 'upper', 3, 0.0, ap, x, 1, 5.0, y, 1 );
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, x, data.strideX, data.beta, y, data.strideY );
 	t.strictEqual( out, y, 'returns expected value' );
-	t.deepEqual( y, expected, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
 
-	expected = new Float64Array( [ 0.0, 0.0, 0.0 ] );
+	t.end();
+});
 
-	out = dspmv( 'column-major', 'upper', 3, 0.0, ap, x, 1, 0.0, y, -1 );
+tape( 'the function supports specifying a negative `y` stride (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var ap;
+	var x;
+	var y;
+
+	data = cxpyn;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
+
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, x, data.strideX, data.beta, y, data.strideY );
 	t.strictEqual( out, y, 'returns expected value' );
-	t.deepEqual( y, expected, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'when `α` is zero, the function scales the second input vector (column-major, lower)', function test( t ) {
+tape( 'the function supports specifying a negative `x` stride (row-major)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = ones( 6, 'float64' );
-	x = ones( 3, 'float64' );
-	y = ones( 3, 'float64' );
+	data = rxnyp;
 
-	expected = new Float64Array( [ 5.0, 5.0, 5.0 ] );
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	out = dspmv( 'column-major', 'lower', 3, 0.0, ap, x, 1, 5.0, y, -1 );
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, x, data.strideX, data.beta, y, data.strideY );
 	t.strictEqual( out, y, 'returns expected value' );
-	t.deepEqual( y, expected, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
 
-	expected = new Float64Array( [ 0.0, 0.0, 0.0 ] );
+	t.end();
+});
 
-	out = dspmv( 'column-major', 'lower', 3, 0.0, ap, x, 1, 0.0, y, 1 );
+tape( 'the function supports specifying a negative `x` stride (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var ap;
+	var x;
+	var y;
+
+	data = cxnyp;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
+
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, x, data.strideX, data.beta, y, data.strideY );
 	t.strictEqual( out, y, 'returns expected value' );
-	t.deepEqual( y, expected, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });
 
 tape( 'the function supports complex access patterns (row-major)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( rxnyn.AP.slice( rxnyn.offsetAP ) );
-	x = new Float64Array( rxnyn.x );
-	y = new Float64Array( rxnyn.y );
+	data = rxnyn;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	expected = new Float64Array( rxnyn.y_out );
+	expected = new Float64Array( data.y_out );
 
-	out = dspmv( rxnyn.order, rxnyn.uplo, rxnyn.N, rxnyn.alpha, ap, x, rxnyn.strideX, rxnyn.beta, y, rxnyn.strideY );
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, x, data.strideX, data.beta, y, data.strideY );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });
 
 tape( 'the function supports complex access patterns (column-major)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( cxnyn.AP.slice( cxnyn.offsetAP ) );
-	x = new Float64Array( cxnyn.x );
-	y = new Float64Array( cxnyn.y );
+	data = cxnyn;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	expected = new Float64Array( cxnyn.y_out );
+	expected = new Float64Array( data.y_out );
 
-	out = dspmv( cxnyn.order, cxnyn.uplo, cxnyn.N, cxnyn.alpha, ap, x, cxnyn.strideX, cxnyn.beta, y, cxnyn.strideY );
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, x, data.strideX, data.beta, y, data.strideY );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/test.ndarray.js b/lib/node_modules/@stdlib/blas/base/dspmv/test/test.ndarray.js
index e02d81ec309..42fdf0c2769 100644
--- a/lib/node_modules/@stdlib/blas/base/dspmv/test/test.ndarray.js
+++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/test.ndarray.js
@@ -24,54 +24,37 @@
 
 var tape = require( 'tape' );
 var Float64Array = require( '@stdlib/array/float64' );
-var EPS = require( '@stdlib/constants/float64/eps' );
-var abs = require( '@stdlib/math/base/special/abs' );
 var ones = require( '@stdlib/array/ones' );
 var dspmv = require( './../lib/ndarray.js' );
 
 
 // FIXTURES //
 
-var rxoyt = require( './fixtures/row_major_xoyt.json' );
+var rl = require( './fixtures/row_major_lower.json' );
+var ru = require( './fixtures/row_major_upper.json' );
 var rxpyp = require( './fixtures/row_major_xpyp.json' );
 var rxnyp = require( './fixtures/row_major_xnyp.json' );
 var rxpyn = require( './fixtures/row_major_xpyn.json' );
 var rxnyn = require( './fixtures/row_major_xnyn.json' );
-
-var cxoyt = require( './fixtures/column_major_xoyt.json' );
+var rsap = require( './fixtures/row_major_sap.json' );
+var rsapn = require( './fixtures/row_major_sapn.json' );
+var roap = require( './fixtures/row_major_oap.json' );
+var rox = require( './fixtures/row_major_ox.json' );
+var roy = require( './fixtures/row_major_oy.json' );
+var rcap = require( './fixtures/row_major_complex_access_pattern.json' );
+
+var cl = require( './fixtures/column_major_lower.json' );
+var cu = require( './fixtures/column_major_upper.json' );
 var cxpyp = require( './fixtures/column_major_xpyp.json' );
 var cxnyp = require( './fixtures/column_major_xnyp.json' );
 var cxpyn = require( './fixtures/column_major_xpyn.json' );
 var cxnyn = require( './fixtures/column_major_xnyn.json' );
-
-
-// FUNCTIONS //
-
-/**
-* Tests for element-wise approximate equality.
-*
-* @private
-* @param {Object} t - test object
-* @param {Collection} actual - actual values
-* @param {Collection} expected - expected values
-* @param {number} rtol - relative tolerance
-*/
-function isApprox( t, actual, expected, rtol ) {
-	var delta;
-	var tol;
-	var i;
-
-	t.strictEqual( actual.length, expected.length, 'returns expected value' );
-	for ( i = 0; i < expected.length; i++ ) {
-		if ( actual[ i ] === expected[ i ] ) {
-			t.strictEqual( actual[ i ], expected[ i ], 'returns expected value' );
-		} else {
-			delta = abs( actual[ i ] - expected[ i ] );
-			tol = rtol * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. actual: '+actual[ i ]+'. expected: '+expected[ i ]+'. delta: '+delta+'. tol: '+tol+'.' );
-		}
-	}
-}
+var csap = require( './fixtures/column_major_sap.json' );
+var csapn = require( './fixtures/column_major_sapn.json' );
+var coap = require( './fixtures/column_major_oap.json' );
+var cox = require( './fixtures/column_major_ox.json' );
+var coy = require( './fixtures/column_major_oy.json' );
+var ccap = require( './fixtures/column_major_complex_access_pattern.json' );
 
 
 // TESTS //
@@ -82,15 +65,18 @@ tape( 'main export is a function', function test( t ) {
 	t.end();
 });
 
-tape( 'the function has an arity of 13', function test( t ) {
-	t.strictEqual( dspmv.length, 13, 'returns expected value' );
+tape( 'the function has an arity of 14', function test( t ) {
+	t.strictEqual( dspmv.length, 14, 'returns expected value' );
 	t.end();
 });
 
 tape( 'the function throws an error if provided an invalid first argument', function test( t ) {
 	var values;
+	var data;
 	var i;
 
+	data = rl;
+
 	values = [
 		'foo',
 		'bar',
@@ -105,15 +91,18 @@ tape( 'the function throws an error if provided an invalid first argument', func
 
 	function badValue( value ) {
 		return function badValue() {
-			dspmv( value, rxpyp.uplo, rxpyp.N, rxpyp.alpha, new Float64Array( rxpyp.AP ), rxpyp.offsetAP, new Float64Array( rxpyp.x ), rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, new Float64Array( rxpyp.y ), rxpyp.strideY, rxpyp.offsetY );
+			dspmv( value, data.uplo, data.N, data.alpha, new Float64Array( data.AP ), data.strideAP, data.offsetAP, new Float64Array( data.x ), data.strideX, data.offsetX, data.beta, new Float64Array( data.y ), data.strideY, data.offsetY );
 		};
 	}
 });
 
 tape( 'the function throws an error if provided an invalid second argument', function test( t ) {
 	var values;
+	var data;
 	var i;
 
+	data = rl;
+
 	values = [
 		'foo',
 		'bar',
@@ -128,15 +117,18 @@ tape( 'the function throws an error if provided an invalid second argument', fun
 
 	function badValue( value ) {
 		return function badValue() {
-			dspmv( rxpyp.order, value, rxpyp.N, rxpyp.alpha, new Float64Array( rxpyp.AP ), rxpyp.offsetAP, new Float64Array( rxpyp.x ), rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, new Float64Array( rxpyp.y ), rxpyp.strideY, rxpyp.offsetY );
+			dspmv( data.order, value, data.N, data.alpha, new Float64Array( data.AP ), data.strideAP, data.offsetAP, new Float64Array( data.x ), data.strideX, data.offsetX, data.beta, new Float64Array( data.y ), data.strideY, data.offsetY );
 		};
 	}
 });
 
 tape( 'the function throws an error if provided an invalid third argument', function test( t ) {
 	var values;
+	var data;
 	var i;
 
+	data = rl;
+
 	values = [
 		-1,
 		-2,
@@ -150,15 +142,18 @@ tape( 'the function throws an error if provided an invalid third argument', func
 
 	function badValue( value ) {
 		return function badValue() {
-			dspmv( rxpyp.order, rxpyp.uplo, value, rxpyp.alpha, new Float64Array( rxpyp.AP ), rxpyp.offsetAP, new Float64Array( rxpyp.x ), rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, new Float64Array( rxpyp.y ), rxpyp.strideY, rxpyp.offsetY );
+			dspmv( data.order, data.uplo, value, data.alpha, new Float64Array( data.AP ), data.strideAP, data.offsetAP, new Float64Array( data.x ), data.strideX, data.offsetX, data.beta, new Float64Array( data.y ), data.strideY, data.offsetY );
 		};
 	}
 });
 
-tape( 'the function throws an error if provided an invalid eighth argument', function test( t ) {
+tape( 'the function throws an error if provided an invalid sixth argument', function test( t ) {
 	var values;
+	var data;
 	var i;
 
+	data = rl;
+
 	values = [
 		0
 	];
@@ -170,15 +165,18 @@ tape( 'the function throws an error if provided an invalid eighth argument', fun
 
 	function badValue( value ) {
 		return function badValue() {
-			dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, new Float64Array( rxpyp.AP ), rxpyp.offsetAP, new Float64Array( rxpyp.x ), value, rxpyp.offsetX, rxpyp.beta, new Float64Array( rxpyp.y ), rxpyp.strideY, rxpyp.offsetY );
+			dspmv( data.order, data.uplo, data.N, data.alpha, new Float64Array( data.AP ), value, data.offsetAP, new Float64Array( data.x ), data.strideX, data.offsetX, data.beta, new Float64Array( data.y ), data.strideY, data.offsetY );
 		};
 	}
 });
 
-tape( 'the function throws an error if provided an invalid twelfth argument', function test( t ) {
+tape( 'the function throws an error if provided an invalid ninth argument', function test( t ) {
 	var values;
+	var data;
 	var i;
 
+	data = rl;
+
 	values = [
 		0
 	];
@@ -190,375 +188,762 @@ tape( 'the function throws an error if provided an invalid twelfth argument', fu
 
 	function badValue( value ) {
 		return function badValue() {
-			dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, new Float64Array( rxpyp.AP ), rxpyp.offsetAP, new Float64Array( rxpyp.x ), rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, new Float64Array( rxpyp.y ), value, rxpyp.offsetY );
+			dspmv( data.order, data.uplo, data.N, data.alpha, new Float64Array( data.AP ), data.strideAP, data.offsetAP, new Float64Array( data.x ), value, data.offsetX, data.beta, new Float64Array( data.y ), data.strideY, data.offsetY );
 		};
 	}
 });
 
-tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `n` element vectors, and `A` is an `n` by `n` symmetric matrix supplied in packed form (row-major, sx=1, sy=1)', function test( t ) {
+tape( 'the function throws an error if provided an invalid thirteenth argument', function test( t ) {
+	var values;
+	var data;
+	var i;
+
+	data = rl;
+
+	values = [
+		0
+	];
+
+	for ( i = 0; i < values.length; i++ ) {
+		t.throws( badValue( values[ i ] ), RangeError, 'throws an error when provided ' + values[ i ] );
+	}
+	t.end();
+
+	function badValue( value ) {
+		return function badValue() {
+			dspmv( data.order, data.uplo, data.N, data.alpha, new Float64Array( data.AP ), data.strideAP, data.offsetAP, new Float64Array( data.x ), data.strideX, data.offsetX, data.beta, new Float64Array( data.y ), value, data.offsetY );
+		};
+	}
+});
+
+tape( 'the function performs the matrix-vector operation `y = α*A*x + β*y (row-major, lower)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( rxpyp.AP );
-	x = new Float64Array( rxpyp.x );
-	y = new Float64Array( rxpyp.y );
+	data = rl;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	expected = new Float64Array( rxpyp.y_out );
+	expected = new Float64Array( data.y_out );
 
-	out = dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, ap, rxpyp.offsetAP, x, rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, y, rxpyp.strideY, rxpyp.offsetY );
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `n` element vectors, and `A` is an `n` by `n` symmetric matrix supplied in packed form (column-major, sx=1, sy=1)', function test( t ) {
+tape( 'the function performs the matrix-vector operation `y = α*A*x + β*y (column-major, lower)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( cxpyp.AP );
-	x = new Float64Array( cxpyp.x );
-	y = new Float64Array( cxpyp.y );
+	data = cl;
 
-	expected = new Float64Array( cxpyp.y_out );
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	out = dspmv( cxpyp.order, cxpyp.uplo, cxpyp.N, cxpyp.alpha, ap, cxpyp.offsetAP, x, cxpyp.strideX, cxpyp.offsetX, cxpyp.beta, y, cxpyp.strideY, cxpyp.offsetY );
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `n` element vectors, and `A` is an `n` by `n` symmetric matrix supplied in packed form (row-major, sx=1, sy=2)', function test( t ) {
+tape( 'the function performs the matrix-vector operation `y = α*A*x + β*y (row-major, upper)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( rxoyt.AP );
-	x = new Float64Array( rxoyt.x );
-	y = new Float64Array( rxoyt.y );
+	data = ru;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	expected = new Float64Array( rxoyt.y_out );
+	expected = new Float64Array( data.y_out );
 
-	out = dspmv( rxoyt.order, rxoyt.uplo, rxoyt.N, rxoyt.alpha, ap, rxoyt.offsetAP, x, rxoyt.strideX, rxoyt.offsetX, rxoyt.beta, y, rxoyt.strideY, rxoyt.offsetY );
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `n` element vectors, and `A` is an `n` by `n` symmetric matrix supplied in packed form (column-major, sx=1, sy=2)', function test( t ) {
+tape( 'the function performs the matrix-vector operation `y = α*A*x + β*y (column-major, upper)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( cxoyt.AP );
-	x = new Float64Array( cxoyt.x );
-	y = new Float64Array( cxoyt.y );
+	data = cu;
 
-	expected = new Float64Array( cxoyt.y_out );
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	out = dspmv( cxoyt.order, cxoyt.uplo, cxoyt.N, cxoyt.alpha, ap, cxoyt.offsetAP, x, cxoyt.strideX, cxoyt.offsetX, cxoyt.beta, y, cxoyt.strideY, cxoyt.offsetY );
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `n` element vectors, and `A` is an `n` by `n` symmetric matrix supplied in packed form (row-major, sx=1, sy=-1)', function test( t ) {
-	var expected;
+tape( 'the function returns a reference to the second input vector', function test( t ) {
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( rxpyn.AP );
-	x = new Float64Array( rxpyn.x );
-	y = new Float64Array( rxpyn.y );
+	data = rl;
 
-	expected = new Float64Array( rxpyn.y_out );
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	out = dspmv( rxpyn.order, rxpyn.uplo, rxpyn.N, rxpyn.alpha, ap, rxpyn.offsetAP, x, rxpyn.strideX, rxpyn.offsetX, rxpyn.beta, y, rxpyn.strideY, rxpyn.offsetY );
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
 
 	t.end();
 });
 
-tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `n` element vectors, and `A` is an `n` by `n` symmetric matrix supplied in packed form (column-major, sx=1, sy=-1)', function test( t ) {
+tape( 'if `N` is zero or the scalar constants are zero and unity, respectively, the function returns the second input vector unchanged (row-major)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( cxpyn.AP );
-	x = new Float64Array( cxpyn.x );
-	y = new Float64Array( cxpyn.y );
+	data = rl;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	expected = new Float64Array( cxpyn.y_out );
+	expected = new Float64Array( data.y );
 
-	out = dspmv( cxpyn.order, cxpyn.uplo, cxpyn.N, cxpyn.alpha, ap, cxpyn.offsetAP, x, cxpyn.strideX, cxpyn.offsetX, cxpyn.beta, y, cxpyn.strideY, cxpyn.offsetY );
+	out = dspmv( data.order, data.uplo, 0, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
+	t.deepEqual( y, expected, 'returns expected value' );
+
+	out = dspmv( data.order, data.uplo, data.N, 0.0, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, 1.0, y, data.strideY, data.offsetY );
+	t.strictEqual( out, y, 'returns expected value' );
+	t.deepEqual( y, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `n` element vectors, and `A` is an `n` by `n` symmetric matrix supplied in packed form (row-major, sx=-1, sy=1)', function test( t ) {
+tape( 'if `N` is zero or the scalar constants are zero and unity, respectively, the function returns the second input vector unchanged (column-major)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( rxnyp.AP );
-	x = new Float64Array( rxnyp.x );
-	y = new Float64Array( rxnyp.y );
+	data = cl;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	expected = new Float64Array( rxnyp.y_out );
+	expected = new Float64Array( data.y );
 
-	out = dspmv( rxnyp.order, rxnyp.uplo, rxnyp.N, rxnyp.alpha, ap, rxnyp.offsetAP, x, rxnyp.strideX, rxnyp.offsetX, rxnyp.beta, y, rxnyp.strideY, rxnyp.offsetY );
+	out = dspmv( data.order, data.uplo, 0, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
+	t.deepEqual( y, expected, 'returns expected value' );
+
+	out = dspmv( data.order, data.uplo, data.N, 0.0, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, 1.0, y, data.strideY, data.offsetY );
+	t.strictEqual( out, y, 'returns expected value' );
+	t.deepEqual( y, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y` where `alpha` and `beta` are scalars, `x` and `y` are `n` element vectors, and `A` is an `n` by `n` symmetric matrix supplied in packed form (column-major, sx=-1, sy=1)', function test( t ) {
+tape( 'when `α` is zero, the function scales the second input vector (row-major, upper)', function test( t ) {
 	var expected;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( cxnyp.AP );
-	x = new Float64Array( cxnyp.x );
-	y = new Float64Array( cxnyp.y );
+	ap = ones( 6, 'float32' );
+	x = ones( 3, 'float32' );
+	y = ones( 3, 'float32' );
+
+	expected = new Float64Array( [ 5.0, 5.0, 5.0 ] );
+
+	out = dspmv( 'row-major', 'upper', 3, 0.0, ap, 1, 0, x, 1, 0, 5.0, y, 1, 0 );
+	t.strictEqual( out, y, 'returns expected value' );
+	t.deepEqual( y, expected, 'returns expected value' );
 
-	expected = new Float64Array( cxnyp.y_out );
+	expected = new Float64Array( [ 0.0, 0.0, 0.0 ] );
 
-	out = dspmv( cxnyp.order, cxnyp.uplo, cxnyp.N, cxnyp.alpha, ap, cxnyp.offsetAP, x, cxnyp.strideX, cxnyp.offsetX, cxnyp.beta, y, cxnyp.strideY, cxnyp.offsetY );
+	out = dspmv( 'row-major', 'upper', 3, 0.0, ap, 1, 0, x, 1, 0, 0.0, y, -1, 2 );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
+	t.deepEqual( y, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function returns a reference to the second input vector', function test( t ) {
+tape( 'when `α` is zero, the function scales the second input vector (row-major, lower)', function test( t ) {
+	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( rxpyp.AP );
-	x = new Float64Array( rxpyp.x );
-	y = new Float64Array( rxpyp.y );
+	data = rl;
 
-	out = dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, ap, rxpyp.offsetAP, x, rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, y, rxpyp.strideY, rxpyp.offsetY );
-	t.strictEqual( out, y, 'returns expected value' );
+	ap = ones( 6, 'float32' );
+	x = ones( 3, 'float32' );
+	y = ones( 3, 'float32' );
+
+	expected = new Float64Array( [ 5.0, 5.0, 5.0 ] );
+
+	out = dspmv( 'row-major', 'lower', 3, 0.0, ap, 1, 0, x, 1, 0, 5.0, y, -1, 2 );
+	t.strictEqual( out, y, 'returns expected value', data.offsetY );
+	t.deepEqual( y, expected, 'returns expected value' );
+
+	expected = new Float64Array( [ 0.0, 0.0, 0.0 ] );
+
+	out = dspmv( 'row-major', 'lower', 3, 0.0, ap, 1, 0, x, 1, 0, 0.0, y, 1, 0 );
+	t.strictEqual( out, y, 'returns expected value', data.offsetY );
+	t.deepEqual( y, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'if `N` is zero or the scalar constants are zero and unity, respectively, the function returns the second input vector unchanged (row-major)', function test( t ) {
+tape( 'when `α` is zero, the function scales the second input vector (column-major, upper)', function test( t ) {
 	var expected;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( rxpyp.AP );
-	x = new Float64Array( rxpyp.x );
-	y = new Float64Array( rxpyp.y );
+	ap = ones( 6, 'float32' );
+	x = ones( 3, 'float32' );
+	y = ones( 3, 'float32' );
 
-	expected = new Float64Array( rxpyp.y );
+	expected = new Float64Array( [ 5.0, 5.0, 5.0 ] );
 
-	out = dspmv( rxpyp.order, rxpyp.uplo, 0, rxpyp.alpha, ap, rxpyp.offsetAP, x, rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, y, rxpyp.strideY, rxpyp.offsetY );
+	out = dspmv( 'column-major', 'upper', 3, 0.0, ap, 1, 0, x, 1, 0, 5.0, y, 1, 0 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.deepEqual( y, expected, 'returns expected value' );
 
-	out = dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, 0.0, ap, rxpyp.offsetAP, x, rxpyp.strideX, rxpyp.offsetX, 1.0, y, rxpyp.strideY, rxpyp.offsetY );
+	expected = new Float64Array( [ 0.0, 0.0, 0.0 ] );
+
+	out = dspmv( 'column-major', 'upper', 3, 0.0, ap, 1, 0, x, 1, 0, 0.0, y, -1, 2 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.deepEqual( y, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'if `N` is zero or the scalar constants are zero and unity, respectively, the function returns the second input vector unchanged (column-major)', function test( t ) {
+tape( 'when `α` is zero, the function scales the second input vector (column-major, lower)', function test( t ) {
 	var expected;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( cxpyp.AP );
-	x = new Float64Array( cxpyp.x );
-	y = new Float64Array( cxpyp.y );
+	ap = ones( 6, 'float32' );
+	x = ones( 3, 'float32' );
+	y = ones( 3, 'float32' );
 
-	expected = new Float64Array( cxpyp.y );
+	expected = new Float64Array( [ 5.0, 5.0, 5.0 ] );
 
-	out = dspmv( cxpyp.order, cxpyp.uplo, 0, cxpyp.alpha, ap, cxpyp.offsetAP, x, cxpyp.strideX, cxpyp.offsetX, cxpyp.beta, y, cxpyp.strideY, cxpyp.offsetY );
+	out = dspmv( 'column-major', 'lower', 3, 0.0, ap, 1, 0, x, 1, 0, 5.0, y, -1, 2 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.deepEqual( y, expected, 'returns expected value' );
 
-	out = dspmv( cxpyp.order, cxpyp.uplo, cxpyp.N, 0.0, ap, cxpyp.offsetAP, x, cxpyp.strideX, cxpyp.offsetX, 1.0, y, cxpyp.strideY, cxpyp.offsetY );
+	expected = new Float64Array( [ 0.0, 0.0, 0.0 ] );
+
+	out = dspmv( 'column-major', 'lower', 3, 0.0, ap, 1, 0, x, 1, 0, 0.0, y, 1, 0 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.deepEqual( y, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'when `α` is zero, the function scales the second input vector (row-major, upper)', function test( t ) {
+tape( 'the function supports specifying `x` and `y` strides (row-major)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = ones( 6, 'float64' );
-	x = ones( 3, 'float64' );
-	y = ones( 3, 'float64' );
+	data = rxpyp;
 
-	expected = new Float64Array( [ 5.0, 5.0, 5.0 ] );
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
+
+	expected = new Float64Array( data.y_out );
 
-	out = dspmv( 'row-major', 'upper', 3, 0.0, ap, 0, x, 1, 0, 5.0, y, 1, 0 );
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
 	t.strictEqual( out, y, 'returns expected value' );
-	t.deepEqual( y, expected, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
 
-	expected = new Float64Array( [ 0.0, 0.0, 0.0 ] );
+	t.end();
+});
 
-	out = dspmv( 'row-major', 'upper', 3, 0.0, ap, 0, x, 1, 0, 0.0, y, -1, 2 );
+tape( 'the function supports specifying `x` and `y` strides (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var ap;
+	var x;
+	var y;
+
+	data = cxpyp;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
+
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
 	t.strictEqual( out, y, 'returns expected value' );
-	t.deepEqual( y, expected, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'when `α` is zero, the function scales the second input vector (row-major, lower)', function test( t ) {
+tape( 'the function supports specifying a negative `y` stride (row-major)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = ones( 6, 'float64' );
-	x = ones( 3, 'float64' );
-	y = ones( 3, 'float64' );
+	data = rxpyn;
 
-	expected = new Float64Array( [ 5.0, 5.0, 5.0 ] );
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	out = dspmv( 'row-major', 'lower', 3, 0.0, ap, 0, x, 1, 0, 5.0, y, -1, 2 );
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
 	t.strictEqual( out, y, 'returns expected value' );
-	t.deepEqual( y, expected, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
 
-	expected = new Float64Array( [ 0.0, 0.0, 0.0 ] );
+	t.end();
+});
 
-	out = dspmv( 'row-major', 'lower', 3, 0.0, ap, 0, x, 1, 0, 0.0, y, 1, 0 );
+tape( 'the function supports specifying a negative `y` stride (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var ap;
+	var x;
+	var y;
+
+	data = cxpyn;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
+
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
 	t.strictEqual( out, y, 'returns expected value' );
-	t.deepEqual( y, expected, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'when `α` is zero, the function scales the second input vector (column-major, upper)', function test( t ) {
+tape( 'the function supports specifying a negative `x` stride (row-major)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = ones( 6, 'float64' );
-	x = ones( 3, 'float64' );
-	y = ones( 3, 'float64' );
+	data = rxnyp;
 
-	expected = new Float64Array( [ 5.0, 5.0, 5.0 ] );
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
+
+	expected = new Float64Array( data.y_out );
 
-	out = dspmv( 'column-major', 'upper', 3, 0.0, ap, 0, x, 1, 0, 5.0, y, 1, 0 );
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
 	t.strictEqual( out, y, 'returns expected value' );
-	t.deepEqual( y, expected, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
 
-	expected = new Float64Array( [ 0.0, 0.0, 0.0 ] );
+	t.end();
+});
+
+tape( 'the function supports specifying a negative `x` stride (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var ap;
+	var x;
+	var y;
+
+	data = cxnyp;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	out = dspmv( 'column-major', 'upper', 3, 0.0, ap, 0, x, 1, 0, 0.0, y, -1, 2 );
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
 	t.strictEqual( out, y, 'returns expected value' );
-	t.deepEqual( y, expected, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'when `α` is zero, the function scales the second input vector (column-major, lower)', function test( t ) {
+tape( 'the function supports specifying negative strides for `x` and `y` (row-major)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = ones( 6, 'float64' );
-	x = ones( 3, 'float64' );
-	y = ones( 3, 'float64' );
+	data = rxnyn;
 
-	expected = new Float64Array( [ 5.0, 5.0, 5.0 ] );
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	out = dspmv( 'column-major', 'lower', 3, 0.0, ap, 0, x, 1, 0, 5.0, y, -1, 2 );
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
 	t.strictEqual( out, y, 'returns expected value' );
-	t.deepEqual( y, expected, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
 
-	expected = new Float64Array( [ 0.0, 0.0, 0.0 ] );
+	t.end();
+});
+
+tape( 'the function supports specifying negative strides for `x` and `y` (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var ap;
+	var x;
+	var y;
 
-	out = dspmv( 'column-major', 'lower', 3, 0.0, ap, 0, x, 1, 0, 0.0, y, 1, 0 );
+	data = cxnyn;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
+
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
 	t.strictEqual( out, y, 'returns expected value' );
-	t.deepEqual( y, expected, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying an `x` offset (row-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var ap;
+	var x;
+	var y;
+
+	data = rox;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
+
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
+	t.strictEqual( out, y, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying an `x` offset (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var ap;
+	var x;
+	var y;
+
+	data = cox;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
+
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
+	t.strictEqual( out, y, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a `y` offset (row-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var ap;
+	var x;
+	var y;
+
+	data = roy;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
+
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
+	t.strictEqual( out, y, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a `y` offset (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var ap;
+	var x;
+	var y;
+
+	data = coy;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
+
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
+	t.strictEqual( out, y, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a stride for `AP` (row-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var ap;
+	var x;
+	var y;
+
+	data = rsap;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
+
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
+	t.strictEqual( out, y, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function supports complex access patterns (row-major)', function test( t ) {
+tape( 'the function supports specifying a stride for `AP` (column-major)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( rxnyn.AP );
-	x = new Float64Array( rxnyn.x );
-	y = new Float64Array( rxnyn.y );
+	data = csap;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	expected = new Float64Array( rxnyn.y_out );
+	expected = new Float64Array( data.y_out );
 
-	out = dspmv( rxnyn.order, rxnyn.uplo, rxnyn.N, rxnyn.alpha, ap, rxnyn.offsetAP, x, rxnyn.strideX, rxnyn.offsetX, rxnyn.beta, y, rxnyn.strideY, rxnyn.offsetY );
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });
 
-tape( 'the function supports complex access patterns (column-major)', function test( t ) {
+tape( 'the function supports specifying a negative stride for `AP` (row-major)', function test( t ) {
 	var expected;
+	var data;
 	var out;
 	var ap;
 	var x;
 	var y;
 
-	ap = new Float64Array( cxnyn.AP );
-	x = new Float64Array( cxnyn.x );
-	y = new Float64Array( cxnyn.y );
+	data = rsapn;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
+
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
+	t.strictEqual( out, y, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying a negative stride for `AP` (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var ap;
+	var x;
+	var y;
+
+	data = csapn;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
+
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
+	t.strictEqual( out, y, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying an offset for `AP` (row-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var ap;
+	var x;
+	var y;
+
+	data = roap;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
+
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
+	t.strictEqual( out, y, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying an offset for `AP` (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var ap;
+	var x;
+	var y;
+
+	data = coap;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
+
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
+	t.strictEqual( out, y, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying complex access patterns (row-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var ap;
+	var x;
+	var y;
+
+	data = rcap;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
+
+	expected = new Float64Array( data.y_out );
+
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
+	t.strictEqual( out, y, 'returns expected value' );
+	t.deepEqual( out, expected, 'returns expected value' );
+
+	t.end();
+});
+
+tape( 'the function supports specifying complex access patterns (column-major)', function test( t ) {
+	var expected;
+	var data;
+	var out;
+	var ap;
+	var x;
+	var y;
+
+	data = ccap;
+
+	ap = new Float64Array( data.AP );
+	x = new Float64Array( data.x );
+	y = new Float64Array( data.y );
 
-	expected = new Float64Array( cxnyn.y_out );
+	expected = new Float64Array( data.y_out );
 
-	out = dspmv( cxnyn.order, cxnyn.uplo, cxnyn.N, cxnyn.alpha, ap, cxnyn.offsetAP, x, cxnyn.strideX, cxnyn.offsetX, cxnyn.beta, y, cxnyn.strideY, cxnyn.offsetY );
+	out = dspmv( data.order, data.uplo, data.N, data.alpha, ap, data.strideAP, data.offsetAP, x, data.strideX, data.offsetX, data.beta, y, data.strideY, data.offsetY );
 	t.strictEqual( out, y, 'returns expected value' );
-	isApprox( t, y, expected, 2.0 );
+	t.deepEqual( out, expected, 'returns expected value' );
 
 	t.end();
 });