From 0d71234301ca239a989d749741055c4f52622692 Mon Sep 17 00:00:00 2001 From: aman-095 Date: Sun, 25 Aug 2024 06:39:54 +0530 Subject: [PATCH 1/6] refactor: update implementation for sspmv --- .../@stdlib/blas/base/sspmv/lib/base.js | 139 ++++++++++++++++++ .../@stdlib/blas/base/sspmv/lib/index.js | 4 +- .../@stdlib/blas/base/sspmv/lib/ndarray.js | 106 +++---------- .../@stdlib/blas/base/sspmv/lib/sspmv.js | 107 ++------------ 4 files changed, 175 insertions(+), 181 deletions(-) create mode 100644 lib/node_modules/@stdlib/blas/base/sspmv/lib/base.js diff --git a/lib/node_modules/@stdlib/blas/base/sspmv/lib/base.js b/lib/node_modules/@stdlib/blas/base/sspmv/lib/base.js new file mode 100644 index 000000000000..2ef5c879a9dd --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/sspmv/lib/base.js @@ -0,0 +1,139 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2024 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var sfill = require( '@stdlib/blas/ext/base/sfill' ).ndarray; +var sscal = require( '@stdlib/blas/base/sscal' ).ndarray; +var f32 = require( '@stdlib/number/float64/base/to-float32' ); + + +// MAIN // + +/** +* 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 +* @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 {Float32Array} AP - packed form of a symmetric matrix `A` +* @param {integer} strideAP - `AP` stride length +* @param {NonNegativeInteger} offsetAP - starting index for `AP` +* @param {Float32Array} x - first input array +* @param {integer} strideX - `x` stride length +* @param {NonNegativeInteger} offsetX - starting `x` index +* @param {number} beta - scalar constant +* @param {Float32Array} y - second input array +* @param {integer} strideY - `y` stride length +* @param {NonNegativeInteger} offsetY - starting `y` index +* @returns {Float32Array} `y` +* +* @example +* var Float32Array = require( '@stdlib/array/float32' ); +* +* var AP = new Float32Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); +* var x = new Float32Array( [ 1.0, 1.0, 1.0 ] ); +* var y = new Float32Array( [ 1.0, 1.0, 1.0 ] ); +* +* sspmv( 'column-major', 'lower', 3, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); +* // y => [ ~7.0, ~12.0, ~15.0 ] +*/ +function sspmv( 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 ox; + var oy; + + if ( beta !== 1.0 ) { + if ( beta === 0.0 ) { + sfill( N, beta, y, strideY, offsetY ); + } else { + sscal( N, beta, y, strideY, offsetY ); + } + } + if ( alpha === 0.0 ) { + return y; + } + isrm = ( order === 'row-major' ); + kk = offsetAP; + ox = offsetX; + oy = offsetY; + if ( ( !isrm && uplo === 'upper' ) || ( isrm && uplo === 'lower' ) ) { + ix1 = ox; + iy1 = oy; + for ( i1 = 0; i1 < N; i1++ ) { + tmp1 = f32( 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[ iy1 ] += f32( f32( tmp1 * AP[ iap ] ) + f32( alpha * tmp2 ) ); + ix1 += strideX; + iy1 += strideY; + kk += ( i1 + 1 ) * strideAP; + } + return y; + } + ix1 = ox; + iy1 = oy; + for ( i1 = 0; i1 < N; i1++ ) { + tmp1 = f32( alpha * x[ ix1 ] ); + tmp2 = 0.0; + iap = kk; + y[ iy1 ] += f32( tmp1 * AP[ iap ] ); + 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[ iy1 ] += f32( alpha * tmp2 ); + ix1 += strideX; + iy1 += strideY; + kk += ( N - i1 ) * strideAP; + } + return y; +} + + +// EXPORTS // + +module.exports = sspmv; diff --git a/lib/node_modules/@stdlib/blas/base/sspmv/lib/index.js b/lib/node_modules/@stdlib/blas/base/sspmv/lib/index.js index 7313620df248..5755a8df251e 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/lib/index.js +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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/sspmv * @@ -42,7 +42,7 @@ * var x = new Float32Array( [ 1.0, 1.0, 1.0 ] ); * var y = new Float32Array( [ 1.0, 1.0, 1.0 ] ); * -* sspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, x, 1, 0, 1.0, y, 1, 0 ); +* sspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); * // y => [ ~7.0, ~12.0, ~15.0 ] */ diff --git a/lib/node_modules/@stdlib/blas/base/sspmv/lib/ndarray.js b/lib/node_modules/@stdlib/blas/base/sspmv/lib/ndarray.js index f27999fabd92..9975cb022f1b 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/lib/ndarray.js +++ b/lib/node_modules/@stdlib/blas/base/sspmv/lib/ndarray.js @@ -20,23 +20,24 @@ // MODULES // -var sfill = require( '@stdlib/blas/ext/base/sfill' ).ndarray; -var sscal = require( '@stdlib/blas/base/sscal' ).ndarray; -var f32 = require( '@stdlib/number/float64/base/to-float32' ); var isLayout = require( '@stdlib/blas/base/assert/is-layout' ); var isMatrixTriangle = require( '@stdlib/blas/base/assert/is-matrix-triangle' ); +var format = require( '@stdlib/string/format' ); +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 {Float32Array} AP - packed form of a symmetric matrix `A` +* @param {integer} strideAP - `AP` stride length +* @param {NonNegativeInteger} offsetAP - starting index for `AP` * @param {Float32Array} x - first input array * @param {integer} strideX - `x` stride length * @param {NonNegativeInteger} offsetX - starting `x` index @@ -47,111 +48,44 @@ var isMatrixTriangle = require( '@stdlib/blas/base/assert/is-matrix-triangle' ); * @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} seventh argument must be non-zero -* @throws {RangeError} eleventh 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 {Float32Array} `y` * * @example * var Float32Array = require( '@stdlib/array/float32' ); * -* var A = new Float32Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); +* var AP = new Float32Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); * var x = new Float32Array( [ 1.0, 1.0, 1.0 ] ); * var y = new Float32Array( [ 1.0, 1.0, 1.0 ] ); * -* sspmv( 'column-major', 'lower', 3, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); +* sspmv( 'column-major', 'lower', 3, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); * // y => [ ~7.0, ~12.0, ~15.0 ] */ -function sspmv( order, uplo, N, alpha, AP, 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; - var kk; - var kx; - var ky; - var j; - var k; - +function sspmv( 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( 'invalid argument. First argument must be a valid order. Value: `%s`.', order ); + throw new TypeError( format( 'invalid argument. First argument must be a valid order. Value: `%s`.', order ) ); } if ( !isMatrixTriangle( uplo ) ) { - throw new TypeError( 'invalid argument. Second argument must specify whether the lower or upper triangular matrix is supplied. Value: `%s`.', uplo ); + throw new TypeError( format( 'invalid argument. Second argument must specify whether the lower or upper triangular matrix is supplied. Value: `%s`.', uplo ) ); } if ( N < 0 ) { - throw new RangeError( 'invalid argument. Third argument must be a nonnegative integer. Value: `%d`.', N ); + 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( 'invalid argument. Seventh 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( 'invalid argument. Eleventh 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; } - // Form: y = beta*y - if ( beta !== 1.0 ) { - if ( beta === 0.0 ) { - sfill( N, 0.0, y, strideY, offsetY ); - } else { - sscal( N, beta, y, strideY, offsetY ); - } - } - if ( alpha === 0.0 ) { - return y; - } - // Form: y = alpha*A*x + y - kx = offsetX; - ky = offsetY; - kk = 0; - if ( - ( order === 'row-major' && uplo === 'upper' ) || - ( order === 'column-major' && uplo === 'lower' ) - ) { - jx = kx; - jy = ky; - for ( j = 0; j < N; j++ ) { - temp1 = f32( alpha * x[ jx ] ); - temp2 = 0.0; - y[ jy ] += f32( temp1 * AP[ kk ] ); - ix = jx; - iy = jy; - for ( k = kk + 1; k < kk + N - j; k++ ) { - ix += strideX; - iy += strideY; - y[ iy ] += f32( temp1 * AP[ k ] ); - temp2 += f32( AP[ k ] * x[ ix ] ); - } - y[ jy ] += f32( alpha * temp2 ); - jx += strideX; - jy += strideY; - kk += N - j; - } - return y; - } - // ( order === 'row-major' && uplo === 'lower') || ( order === 'column-major' && uplo === 'upper' ) - jx = kx; - jy = ky; - for ( j = 0; j < N; j++ ) { - temp1 = f32( alpha * x[ jx ] ); - temp2 = 0.0; - ix = kx; - iy = ky; - for ( k = kk; k < kk + j; k++ ) { - y[ iy ] += f32( temp1 * AP[ k ] ); - temp2 += f32( AP[ k ] * x[ ix ] ); - ix += strideX; - iy += strideY; - } - y[ jy ] += f32( f32( temp1 * AP[ kk + j ] ) + f32( alpha * temp2 ) ); - jx += strideX; - jy += strideY; - kk += j + 1; - } - return y; + return base( order, uplo, N, alpha, AP, strideAP, offsetAP, x, strideX, offsetX, beta, y, strideY, offsetY ); } diff --git a/lib/node_modules/@stdlib/blas/base/sspmv/lib/sspmv.js b/lib/node_modules/@stdlib/blas/base/sspmv/lib/sspmv.js index 1982ab3929b1..204b53f7d0c3 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/lib/sspmv.js +++ b/lib/node_modules/@stdlib/blas/base/sspmv/lib/sspmv.js @@ -20,17 +20,17 @@ // MODULES // -var sfill = require( '@stdlib/blas/ext/base/sfill' ); -var sscal = require( '@stdlib/blas/base/sscal' ); -var f32 = require( '@stdlib/number/float64/base/to-float32' ); +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 format = require( '@stdlib/string/format' ); +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,109 +60,30 @@ var isMatrixTriangle = require( '@stdlib/blas/base/assert/is-matrix-triangle' ); * // y => [ ~7.0, ~12.0, ~15.0 ] */ function sspmv( order, uplo, N, alpha, AP, x, strideX, beta, y, strideY ) { - var temp1; - var temp2; - var ix; - var iy; - var jx; - var jy; - var kk; - var kx; - var ky; - var sy; - var j; - var k; + var ox; + var oy; if ( !isLayout( order ) ) { - throw new TypeError( 'invalid argument. First argument must be a valid order. Value: `%s`.', order ); + throw new TypeError( format( 'invalid argument. First argument must be a valid order. Value: `%s`.', order ) ); } if ( !isMatrixTriangle( uplo ) ) { - throw new TypeError( 'invalid argument. Second argument must specify whether the lower or upper triangular matrix is supplied. Value: `%s`.', uplo ); + throw new TypeError( format( 'invalid argument. Second argument must specify whether the lower or upper triangular matrix is supplied. Value: `%s`.', uplo ) ); } if ( N < 0 ) { - throw new RangeError( 'invalid argument. Third argument must be a nonnegative integer. Value: `%d`.', N ); + throw new RangeError( format( 'invalid argument. Third argument must be a nonnegative integer. Value: `%d`.', N ) ); } if ( strideX === 0 ) { - throw new RangeError( 'invalid argument. Seventh argument must be non-zero. Value: `%d`.', strideX ); + throw new RangeError( format( 'invalid argument. Seventh argument must be non-zero. Value: `%d`.', strideX ) ); } if ( strideY === 0 ) { - throw new RangeError( 'invalid argument. Tenth argument must be non-zero. Value: `%d`.', strideY ); + throw new RangeError( format( 'invalid argument. Tenth argument must be non-zero. Value: `%d`.', strideY ) ); } if ( N === 0 || ( alpha === 0.0 && beta === 1.0 ) ) { return y; } - // Form: y = beta*y - sy = strideY; - if ( beta !== 1.0 ) { - if ( beta === 0.0 ) { - sfill( N, 0.0, y, strideY ); - } else { - if ( strideY < 0 ) { - sy = -sy; - } - sscal( N, beta, y, sy ); - } - } - if ( alpha === 0.0 ) { - return y; - } - if ( strideX > 0 ) { - kx = 0; - } else { - kx = ( 1 - N ) * strideX; - } - if ( strideY > 0 ) { - ky = 0; - } else { - ky = ( 1 - N ) * strideY; - } - // Form: y = alpha*A*x + y - kk = 0; - if ( - ( order === 'row-major' && uplo === 'upper' ) || - ( order === 'column-major' && uplo === 'lower' ) - ) { - jx = kx; - jy = ky; - for ( j = 0; j < N; j++ ) { - temp1 = f32( alpha * x[ jx ] ); - temp2 = 0.0; - y[ jy ] += f32( temp1 * AP[ kk ] ); - ix = jx; - iy = jy; - for ( k = kk + 1; k < kk + N - j; k++ ) { - ix += strideX; - iy += strideY; - y[ iy ] += f32( temp1 * AP[ k ] ); - temp2 += f32( AP[ k ] * x[ ix ] ); - } - y[ jy ] += f32( alpha * temp2 ); - jx += strideX; - jy += strideY; - kk += N - j; - } - return y; - } - // ( order === 'row-major' && uplo === 'lower') || ( order === 'column-major' && uplo === 'upper' ) - jx = kx; - jy = ky; - for ( j = 0; j < N; j++ ) { - temp1 = f32( alpha * x[ jx ] ); - temp2 = 0.0; - ix = kx; - iy = ky; - for ( k = kk; k < kk + j; k++ ) { - y[ iy ] += f32( temp1 * AP[ k ] ); - temp2 += f32( AP[ k ] * x[ ix ] ); - ix += strideX; - iy += strideY; - } - y[ jy ] += f32( f32( temp1 * AP[ kk + j ] ) + f32( alpha * temp2 ) ); - jx += strideX; - jy += strideY; - kk += j + 1; - } - 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 ); } From 6812b08f6f45a29c46de5732958103530c919fe1 Mon Sep 17 00:00:00 2001 From: aman-095 Date: Sun, 25 Aug 2024 06:42:11 +0530 Subject: [PATCH 2/6] docs: update comments --- lib/node_modules/@stdlib/blas/base/sspmv/lib/ndarray.js | 2 +- lib/node_modules/@stdlib/blas/base/sspmv/lib/sspmv.js | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/base/sspmv/lib/ndarray.js b/lib/node_modules/@stdlib/blas/base/sspmv/lib/ndarray.js index 9975cb022f1b..7a70a875e598 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/lib/ndarray.js +++ b/lib/node_modules/@stdlib/blas/base/sspmv/lib/ndarray.js @@ -85,7 +85,7 @@ function sspmv( order, uplo, N, alpha, AP, strideAP, offsetAP, x, strideX, offse if ( N === 0 || ( alpha === 0.0 && beta === 1.0 ) ) { return y; } - return base( order, uplo, N, alpha, AP, strideAP, offsetAP, x, strideX, offsetX, beta, y, strideY, offsetY ); + 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/sspmv/lib/sspmv.js b/lib/node_modules/@stdlib/blas/base/sspmv/lib/sspmv.js index 204b53f7d0c3..caed43203ca3 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/lib/sspmv.js +++ b/lib/node_modules/@stdlib/blas/base/sspmv/lib/sspmv.js @@ -83,7 +83,7 @@ function sspmv( order, uplo, N, alpha, AP, x, strideX, beta, y, strideY ) { } ox = stride2offset( N, strideX ); oy = stride2offset( N, strideY ); - return base( order, uplo, N, alpha, AP, 1, 0, x, strideX, ox, beta, y, strideY, oy ); + return base( order, uplo, N, alpha, AP, 1, 0, x, strideX, ox, beta, y, strideY, oy ); // eslint-disable-line max-len } From ab9cf40763885fcf0939c77af285e68eec5f65e9 Mon Sep 17 00:00:00 2001 From: aman-095 Date: Mon, 26 Aug 2024 10:34:29 +0530 Subject: [PATCH 3/6] bench: update benchmarks to match updated implementation --- .../@stdlib/blas/base/sspmv/benchmark/benchmark.js | 12 ++++++------ .../blas/base/sspmv/benchmark/benchmark.ndarray.js | 12 ++++++------ 2 files changed, 12 insertions(+), 12 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/base/sspmv/benchmark/benchmark.js b/lib/node_modules/@stdlib/blas/base/sspmv/benchmark/benchmark.js index fee324e074ab..1ae7437c5b8e 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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 = sspmv( 'row-major', 'upper', len, 1.0, AP, x, 1, 1.0, y, 1 ); + z = sspmv( 'row-major', 'upper', N, 1.0, AP, x, 1, 1.0, y, 1 ); if ( isnanf( z[ i%z.length ] ) ) { b.fail( 'should not return NaN' ); } diff --git a/lib/node_modules/@stdlib/blas/base/sspmv/benchmark/benchmark.ndarray.js b/lib/node_modules/@stdlib/blas/base/sspmv/benchmark/benchmark.ndarray.js index ed179effd94b..ee6e01e99ae8 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/benchmark/benchmark.ndarray.js +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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 = sspmv( 'row-major', 'upper', len, 1.0, AP, x, 1, 0, 1.0, y, 1, 0 ); + z = sspmv( 'row-major', 'upper', N, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); if ( isnanf( z[ i%z.length ] ) ) { b.fail( 'should not return NaN' ); } From 647e7b02d602f42907a345eef82cfd2a5129f2ea Mon Sep 17 00:00:00 2001 From: aman-095 Date: Mon, 26 Aug 2024 10:39:17 +0530 Subject: [PATCH 4/6] docs: update descriptions --- .../@stdlib/blas/base/sspmv/README.md | 17 +++++++++++------ .../@stdlib/blas/base/sspmv/examples/index.js | 7 +++++-- .../@stdlib/blas/base/sspmv/package.json | 2 +- 3 files changed, 17 insertions(+), 9 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/base/sspmv/README.md b/lib/node_modules/@stdlib/blas/base/sspmv/README.md index d302ac8e5542..a38addab8fb3 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/README.md +++ b/lib/node_modules/@stdlib/blas/base/sspmv/README.md @@ -20,7 +20,7 @@ limitations under the License. # sspmv -> 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`.
@@ -91,7 +91,7 @@ sspmv( 'row-major', 'upper', 2, 1.0, AP, x1, -1, 1.0, y1, -1 ); // y0 => [ 0.0, 6.0, 4.0 ] ``` -#### sspmv.ndarray( order, uplo, N, α, AP, x, sx, ox, β, y, sy, oy ) +#### sspmv.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,12 +102,14 @@ var AP = new Float32Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); var x = new Float32Array( [ 1.0, 1.0, 1.0 ] ); var y = new Float32Array( [ 1.0, 1.0, 1.0 ] ); -sspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, x, 1, 0, 1.0, y, 1, 0 ); +sspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // y => [ 7.0, 12.0, 15.0 ] ``` The function has the following additional parameters: +- **sap**: `AP` stride length. +- **oap**: starting index for `AP`. - **ox**: starting index for `x`. - **oy**: starting index for `y`. @@ -120,7 +122,7 @@ var AP = new Float32Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); var x = new Float32Array( [ 1.0, 1.0, 1.0 ] ); var y = new Float32Array( [ 1.0, 1.0, 1.0 ] ); -sspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, x, 1, 0, 1.0, y, -1, 2 ); +sspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, -1, 2 ); // y => [ 15.0, 12.0, 7.0 ] ``` @@ -152,13 +154,16 @@ var opts = { 'dtype': 'float32' }; -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 ); -sspmv.ndarray( 'row-major', 'upper', N, 1.0, AP, x, 1, 0, 1.0, y, 1, 0 ); +sspmv( 'row-major', 'upper', N, 1.0, AP, x, 1, 1.0, y, 1 ); +console.log( y ); + +sspmv.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/sspmv/examples/index.js b/lib/node_modules/@stdlib/blas/base/sspmv/examples/index.js index 23391078e6b8..977c095c94a6 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/examples/index.js +++ b/lib/node_modules/@stdlib/blas/base/sspmv/examples/index.js @@ -25,11 +25,14 @@ var opts = { 'dtype': 'float32' }; -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 ); -sspmv.ndarray( 'row-major', 'upper', N, 1.0, AP, x, 1, 0, 1.0, y, 1, 0 ); +sspmv( 'row-major', 'upper', N, 1.0, AP, x, 1, 1.0, y, 1 ); +console.log( y ); + +sspmv.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/sspmv/package.json b/lib/node_modules/@stdlib/blas/base/sspmv/package.json index 9bcc532b0903..d255d18f619f 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/package.json +++ b/lib/node_modules/@stdlib/blas/base/sspmv/package.json @@ -1,7 +1,7 @@ { "name": "@stdlib/blas/base/sspmv", "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", From bad1e5b4c6cfbc4d408646bf60a159ba2efb8a14 Mon Sep 17 00:00:00 2001 From: aman-095 Date: Wed, 28 Aug 2024 12:56:37 +0530 Subject: [PATCH 5/6] test: update and add test cases --- .../@stdlib/blas/base/sspmv/lib/base.js | 2 +- .../column_major_complex_access_pattern.json | 17 + ...ajor_xoyt.json => column_major_lower.json} | 4 +- .../sspmv/test/fixtures/column_major_oap.json | 17 + .../sspmv/test/fixtures/column_major_ox.json | 17 + .../sspmv/test/fixtures/column_major_oy.json | 17 + .../sspmv/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 | 2 + .../test/fixtures/column_major_xnyp.json | 2 + .../test/fixtures/column_major_xpyn.json | 2 + .../test/fixtures/column_major_xpyp.json | 2 + .../row_major_complex_access_pattern.json | 17 + ...w_major_xoyt.json => row_major_lower.json} | 8 +- .../sspmv/test/fixtures/row_major_oap.json | 17 + .../sspmv/test/fixtures/row_major_ox.json | 17 + .../sspmv/test/fixtures/row_major_oy.json | 17 + .../sspmv/test/fixtures/row_major_sap.json | 17 + .../sspmv/test/fixtures/row_major_sapn.json | 17 + .../sspmv/test/fixtures/row_major_upper.json | 17 + .../sspmv/test/fixtures/row_major_xnyn.json | 2 + .../sspmv/test/fixtures/row_major_xnyp.json | 2 + .../sspmv/test/fixtures/row_major_xpyn.json | 2 + .../sspmv/test/fixtures/row_major_xpyp.json | 2 + .../blas/base/sspmv/test/test.ndarray.js | 733 +++++++++++++----- .../blas/base/sspmv/test/test.sspmv.js | 405 ++++++---- 27 files changed, 1060 insertions(+), 346 deletions(-) create mode 100644 lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_complex_access_pattern.json rename lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/{column_major_xoyt.json => column_major_lower.json} (90%) create mode 100644 lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_oap.json create mode 100644 lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_ox.json create mode 100644 lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_oy.json create mode 100644 lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_sap.json create mode 100644 lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_sapn.json create mode 100644 lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_upper.json create mode 100644 lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_complex_access_pattern.json rename lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/{row_major_xoyt.json => row_major_lower.json} (72%) create mode 100644 lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_oap.json create mode 100644 lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_ox.json create mode 100644 lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_oy.json create mode 100644 lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_sap.json create mode 100644 lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_sapn.json create mode 100644 lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_upper.json diff --git a/lib/node_modules/@stdlib/blas/base/sspmv/lib/base.js b/lib/node_modules/@stdlib/blas/base/sspmv/lib/base.js index 2ef5c879a9dd..43b4883bef9c 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/lib/base.js +++ b/lib/node_modules/@stdlib/blas/base/sspmv/lib/base.js @@ -114,8 +114,8 @@ function sspmv( order, uplo, N, alpha, AP, strideAP, offsetAP, x, strideX, offse for ( i1 = 0; i1 < N; i1++ ) { tmp1 = f32( alpha * x[ ix1 ] ); tmp2 = 0.0; + y[ iy1 ] += f32( tmp1 * AP[ kk ] ); iap = kk; - y[ iy1 ] += f32( tmp1 * AP[ iap ] ); ix0 = ix1; iy0 = iy1; for ( i0 = i1+1; i0 < N; i0++ ) { diff --git a/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_complex_access_pattern.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_complex_access_pattern.json new file mode 100644 index 000000000000..8fd5802c13c7 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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/sspmv/test/fixtures/column_major_xoyt.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_lower.json similarity index 90% rename from lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_xoyt.json rename to lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_lower.json index 7ce888f02a01..37c7121ef6dd 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_xoyt.json +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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, @@ -9,6 +9,8 @@ "strideY": 2, "offsetY": 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/sspmv/test/fixtures/column_major_oap.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_oap.json new file mode 100644 index 000000000000..b54652dadca0 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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/sspmv/test/fixtures/column_major_ox.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_ox.json new file mode 100644 index 000000000000..2dc0f16614b4 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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/sspmv/test/fixtures/column_major_oy.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_oy.json new file mode 100644 index 000000000000..7a90138b116a --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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/sspmv/test/fixtures/column_major_sap.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_sap.json new file mode 100644 index 000000000000..4ba7aad930ad --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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/sspmv/test/fixtures/column_major_sapn.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_sapn.json new file mode 100644 index 000000000000..410bec99c5e8 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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/sspmv/test/fixtures/column_major_upper.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_upper.json new file mode 100644 index 000000000000..f2965d50f6c1 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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/sspmv/test/fixtures/column_major_xnyn.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_xnyn.json index dd535606e057..07cc43b13d85 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_xnyn.json +++ b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_xnyn.json @@ -9,6 +9,8 @@ "strideY": -1, "offsetY": 2, "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/sspmv/test/fixtures/column_major_xnyp.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_xnyp.json index 88b876d94788..1807dbf1b0be 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_xnyp.json +++ b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_xnyp.json @@ -9,6 +9,8 @@ "strideY": 1, "offsetY": 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/sspmv/test/fixtures/column_major_xpyn.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_xpyn.json index 8fedbd443972..4e08b89290a8 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_xpyn.json +++ b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_xpyn.json @@ -9,6 +9,8 @@ "strideY": -1, "offsetY": 2, "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/sspmv/test/fixtures/column_major_xpyp.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_xpyp.json index 0e57578aba8e..09d16df7ab73 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_xpyp.json +++ b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/column_major_xpyp.json @@ -9,6 +9,8 @@ "strideY": 1, "offsetY": 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/sspmv/test/fixtures/row_major_complex_access_pattern.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_complex_access_pattern.json new file mode 100644 index 000000000000..88b89e2cc53c --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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/sspmv/test/fixtures/row_major_xoyt.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_lower.json similarity index 72% rename from lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_xoyt.json rename to lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_lower.json index f77e7da3427f..a393e4279a66 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_xoyt.json +++ b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_lower.json @@ -1,15 +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, "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/sspmv/test/fixtures/row_major_oap.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_oap.json new file mode 100644 index 000000000000..fb34b7e18d1b --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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/sspmv/test/fixtures/row_major_ox.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_ox.json new file mode 100644 index 000000000000..b39813d836c9 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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/sspmv/test/fixtures/row_major_oy.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_oy.json new file mode 100644 index 000000000000..660a96d15dd4 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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/sspmv/test/fixtures/row_major_sap.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_sap.json new file mode 100644 index 000000000000..b827f7fa9a43 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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/sspmv/test/fixtures/row_major_sapn.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_sapn.json new file mode 100644 index 000000000000..b73e985d8388 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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/sspmv/test/fixtures/row_major_upper.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_upper.json new file mode 100644 index 000000000000..d9fc2e5f7d05 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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/sspmv/test/fixtures/row_major_xnyn.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_xnyn.json index 4763913667f8..d9aacf91386f 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_xnyn.json +++ b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_xnyn.json @@ -9,6 +9,8 @@ "strideY": -1, "offsetY": 2, "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/sspmv/test/fixtures/row_major_xnyp.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_xnyp.json index fc5df6d7ce93..ca868017cb2c 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_xnyp.json +++ b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_xnyp.json @@ -9,6 +9,8 @@ "strideY": 1, "offsetY": 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/sspmv/test/fixtures/row_major_xpyn.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_xpyn.json index 82387497ca33..66696f9f69e4 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_xpyn.json +++ b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_xpyn.json @@ -9,6 +9,8 @@ "strideY": -1, "offsetY": 2, "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/sspmv/test/fixtures/row_major_xpyp.json b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_xpyp.json index d4deab60a3fd..4878a217089b 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_xpyp.json +++ b/lib/node_modules/@stdlib/blas/base/sspmv/test/fixtures/row_major_xpyp.json @@ -9,6 +9,8 @@ "strideY": 1, "offsetY": 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/sspmv/test/test.ndarray.js b/lib/node_modules/@stdlib/blas/base/sspmv/test/test.ndarray.js index c3c1d111cc4c..fa9c0e5ad89e 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/test/test.ndarray.js +++ b/lib/node_modules/@stdlib/blas/base/sspmv/test/test.ndarray.js @@ -24,54 +24,37 @@ var tape = require( 'tape' ); var Float32Array = require( '@stdlib/array/float32' ); -var EPS = require( '@stdlib/constants/float32/eps' ); -var abs = require( '@stdlib/math/base/special/abs' ); var ones = require( '@stdlib/array/ones' ); var sspmv = 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 12', function test( t ) { - t.strictEqual( sspmv.length, 12, 'returns expected value' ); +tape( 'the function has an arity of 14', function test( t ) { + t.strictEqual( sspmv.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() { - sspmv( value, rxpyp.uplo, rxpyp.N, rxpyp.alpha, new Float32Array( rxpyp.AP ), new Float32Array( rxpyp.x ), rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, new Float32Array( rxpyp.y ), rxpyp.strideY, rxpyp.offsetY ); + sspmv( value, data.uplo, data.N, data.alpha, new Float32Array( data.AP ), data.strideAP, data.offsetAP, new Float32Array( data.x ), data.strideX, data.offsetX, data.beta, new Float32Array( 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() { - sspmv( rxpyp.order, value, rxpyp.N, rxpyp.alpha, new Float32Array( rxpyp.AP ), new Float32Array( rxpyp.x ), rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, new Float32Array( rxpyp.y ), rxpyp.strideY, rxpyp.offsetY ); + sspmv( data.order, value, data.N, data.alpha, new Float32Array( data.AP ), data.strideAP, data.offsetAP, new Float32Array( data.x ), data.strideX, data.offsetX, data.beta, new Float32Array( 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,41 @@ tape( 'the function throws an error if provided an invalid third argument', func function badValue( value ) { return function badValue() { - sspmv( rxpyp.order, rxpyp.uplo, value, rxpyp.alpha, new Float32Array( rxpyp.AP ), new Float32Array( rxpyp.x ), rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, new Float32Array( rxpyp.y ), rxpyp.strideY, rxpyp.offsetY ); + sspmv( data.order, data.uplo, value, data.alpha, new Float32Array( data.AP ), data.strideAP, data.offsetAP, new Float32Array( data.x ), data.strideX, data.offsetX, data.beta, new Float32Array( data.y ), data.strideY, data.offsetY ); + }; + } +}); + +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 + ]; + + 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() { + sspmv( data.order, data.uplo, data.N, data.alpha, new Float32Array( data.AP ), value, data.offsetAP, new Float32Array( data.x ), data.strideX, data.offsetX, data.beta, new Float32Array( data.y ), data.strideY, data.offsetY ); }; } }); -tape( 'the function throws an error if provided an invalid seventh 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 ]; @@ -170,15 +188,18 @@ tape( 'the function throws an error if provided an invalid seventh argument', fu function badValue( value ) { return function badValue() { - sspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, new Float32Array( rxpyp.AP ), new Float32Array( rxpyp.x ), value, rxpyp.offsetX, rxpyp.beta, new Float32Array( rxpyp.y ), rxpyp.strideY, rxpyp.offsetY ); + sspmv( data.order, data.uplo, data.N, data.alpha, new Float32Array( data.AP ), data.strideAP, data.offsetAP, new Float32Array( data.x ), value, data.offsetX, data.beta, new Float32Array( data.y ), data.strideY, data.offsetY ); }; } }); -tape( 'the function throws an error if provided an invalid eleventh argument', 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 ]; @@ -190,375 +211,739 @@ tape( 'the function throws an error if provided an invalid eleventh argument', f function badValue( value ) { return function badValue() { - sspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, new Float32Array( rxpyp.AP ), new Float32Array( rxpyp.x ), rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, new Float32Array( rxpyp.y ), value, rxpyp.offsetY ); + sspmv( data.order, data.uplo, data.N, data.alpha, new Float32Array( data.AP ), data.strideAP, data.offsetAP, new Float32Array( data.x ), data.strideX, data.offsetX, data.beta, new Float32Array( data.y ), value, 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 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 Float32Array( rxpyp.AP ); - x = new Float32Array( rxpyp.x ); - y = new Float32Array( rxpyp.y ); + data = rl; + + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - expected = new Float32Array( rxpyp.y_out ); + expected = new Float32Array( data.y_out ); - out = sspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, ap, x, rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, y, rxpyp.strideY, rxpyp.offsetY ); + out = sspmv( 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 Float32Array( cxpyp.AP ); - x = new Float32Array( cxpyp.x ); - y = new Float32Array( cxpyp.y ); + data = cl; - expected = new Float32Array( cxpyp.y_out ); + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - out = sspmv( cxpyp.order, cxpyp.uplo, cxpyp.N, cxpyp.alpha, ap, x, cxpyp.strideX, cxpyp.offsetX, cxpyp.beta, y, cxpyp.strideY, cxpyp.offsetY ); + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 Float32Array( rxoyt.AP ); - x = new Float32Array( rxoyt.x ); - y = new Float32Array( rxoyt.y ); + data = ru; + + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - expected = new Float32Array( rxoyt.y_out ); + expected = new Float32Array( data.y_out ); - out = sspmv( rxoyt.order, rxoyt.uplo, rxoyt.N, rxoyt.alpha, ap, x, rxoyt.strideX, rxoyt.offsetX, rxoyt.beta, y, rxoyt.strideY, rxoyt.offsetY ); + out = sspmv( 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 Float32Array( cxoyt.AP ); - x = new Float32Array( cxoyt.x ); - y = new Float32Array( cxoyt.y ); + data = cu; + + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - expected = new Float32Array( cxoyt.y_out ); + expected = new Float32Array( data.y_out ); - out = sspmv( cxoyt.order, cxoyt.uplo, cxoyt.N, cxoyt.alpha, ap, x, cxoyt.strideX, cxoyt.offsetX, cxoyt.beta, y, cxoyt.strideY, cxoyt.offsetY ); + out = sspmv( 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 Float32Array( rxpyn.AP ); - x = new Float32Array( rxpyn.x ); - y = new Float32Array( rxpyn.y ); + data = rl; - expected = new Float32Array( rxpyn.y_out ); + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - out = sspmv( rxpyn.order, rxpyn.uplo, rxpyn.N, rxpyn.alpha, ap, x, rxpyn.strideX, rxpyn.offsetX, rxpyn.beta, y, rxpyn.strideY, rxpyn.offsetY ); + out = sspmv( 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 Float32Array( cxpyn.AP ); - x = new Float32Array( cxpyn.x ); - y = new Float32Array( cxpyn.y ); + data = rl; + + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); + + expected = new Float32Array( data.y ); - expected = new Float32Array( cxpyn.y_out ); + out = sspmv( 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' ); + t.deepEqual( y, expected, 'returns expected value' ); - out = sspmv( cxpyn.order, cxpyn.uplo, cxpyn.N, cxpyn.alpha, ap, x, cxpyn.strideX, cxpyn.offsetX, cxpyn.beta, y, cxpyn.strideY, cxpyn.offsetY ); + out = sspmv( 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' ); - 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 Float32Array( rxnyp.AP ); - x = new Float32Array( rxnyp.x ); - y = new Float32Array( rxnyp.y ); + data = cl; + + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - expected = new Float32Array( rxnyp.y_out ); + expected = new Float32Array( data.y ); + + out = sspmv( 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' ); + t.deepEqual( y, expected, 'returns expected value' ); - out = sspmv( rxnyp.order, rxnyp.uplo, rxnyp.N, rxnyp.alpha, ap, x, rxnyp.strideX, rxnyp.offsetX, rxnyp.beta, y, rxnyp.strideY, rxnyp.offsetY ); + out = sspmv( 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' ); - 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 (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 Float32Array( cxnyp.AP ); - x = new Float32Array( cxnyp.x ); - y = new Float32Array( cxnyp.y ); + ap = ones( 6, 'float32' ); + x = ones( 3, 'float32' ); + y = ones( 3, 'float32' ); - expected = new Float32Array( cxnyp.y_out ); + expected = new Float32Array( [ 5.0, 5.0, 5.0 ] ); - out = sspmv( cxnyp.order, cxnyp.uplo, cxnyp.N, cxnyp.alpha, ap, x, cxnyp.strideX, cxnyp.offsetX, cxnyp.beta, y, cxnyp.strideY, cxnyp.offsetY ); + out = sspmv( 'row-major', 'upper', 3, 0.0, ap, 1, 0, x, 1, 0, 5.0, y, 1, 0 ); t.strictEqual( out, y, 'returns expected value' ); - isApprox( t, y, expected, 2.0 ); + t.deepEqual( y, expected, 'returns expected value' ); + + expected = new Float32Array( [ 0.0, 0.0, 0.0 ] ); + + out = sspmv( 'row-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( '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 Float32Array( rxpyp.AP ); - x = new Float32Array( rxpyp.x ); - y = new Float32Array( rxpyp.y ); + data = rl; - out = sspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, ap, 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 Float32Array( [ 5.0, 5.0, 5.0 ] ); + + out = sspmv( '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 Float32Array( [ 0.0, 0.0, 0.0 ] ); + + out = sspmv( '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 Float32Array( rxpyp.AP ); - x = new Float32Array( rxpyp.x ); - y = new Float32Array( rxpyp.y ); + ap = ones( 6, 'float32' ); + x = ones( 3, 'float32' ); + y = ones( 3, 'float32' ); - expected = new Float32Array( rxpyp.y ); + expected = new Float32Array( [ 5.0, 5.0, 5.0 ] ); - out = sspmv( rxpyp.order, rxpyp.uplo, 0, rxpyp.alpha, ap, x, rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, y, rxpyp.strideY, rxpyp.offsetY ); + out = sspmv( '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 = sspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, 0.0, ap, x, rxpyp.strideX, rxpyp.offsetX, 1.0, y, rxpyp.strideY, rxpyp.offsetY ); + expected = new Float32Array( [ 0.0, 0.0, 0.0 ] ); + + out = sspmv( '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 Float32Array( cxpyp.AP ); - x = new Float32Array( cxpyp.x ); - y = new Float32Array( cxpyp.y ); + ap = ones( 6, 'float32' ); + x = ones( 3, 'float32' ); + y = ones( 3, 'float32' ); - expected = new Float32Array( cxpyp.y ); + expected = new Float32Array( [ 5.0, 5.0, 5.0 ] ); - out = sspmv( cxpyp.order, cxpyp.uplo, 0, cxpyp.alpha, ap, x, cxpyp.strideX, cxpyp.offsetX, cxpyp.beta, y, cxpyp.strideY, cxpyp.offsetY ); + out = sspmv( '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 = sspmv( cxpyp.order, cxpyp.uplo, cxpyp.N, 0.0, ap, x, cxpyp.strideX, cxpyp.offsetX, 1.0, y, cxpyp.strideY, cxpyp.offsetY ); + expected = new Float32Array( [ 0.0, 0.0, 0.0 ] ); + + out = sspmv( '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, 'float32' ); - x = ones( 3, 'float32' ); - y = ones( 3, 'float32' ); + data = rxpyp; - expected = new Float32Array( [ 5.0, 5.0, 5.0 ] ); + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - out = sspmv( 'row-major', 'upper', 3, 0.0, ap, x, 1, 0, 5.0, y, 1, 0 ); + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 Float32Array( [ 0.0, 0.0, 0.0 ] ); + t.end(); +}); + +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 Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - out = sspmv( 'row-major', 'upper', 3, 0.0, ap, x, 1, 0, 0.0, y, -1, 2 ); + expected = new Float32Array( data.y_out ); + + out = sspmv( 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, 'float32' ); - x = ones( 3, 'float32' ); - y = ones( 3, 'float32' ); + data = rxpyn; - expected = new Float32Array( [ 5.0, 5.0, 5.0 ] ); + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - out = sspmv( 'row-major', 'lower', 3, 0.0, ap, x, 1, 0, 5.0, y, -1, 2 ); + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 Float32Array( [ 0.0, 0.0, 0.0 ] ); + t.end(); +}); + +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 Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - out = sspmv( 'row-major', 'lower', 3, 0.0, ap, x, 1, 0, 0.0, y, 1, 0 ); + expected = new Float32Array( data.y_out ); + + out = sspmv( 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, 'float32' ); - x = ones( 3, 'float32' ); - y = ones( 3, 'float32' ); + data = rxnyp; - expected = new Float32Array( [ 5.0, 5.0, 5.0 ] ); + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - out = sspmv( 'column-major', 'upper', 3, 0.0, ap, x, 1, 0, 5.0, y, 1, 0 ); + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 Float32Array( [ 0.0, 0.0, 0.0 ] ); + t.end(); +}); - out = sspmv( 'column-major', 'upper', 3, 0.0, ap, x, 1, 0, 0.0, y, -1, 2 ); +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 Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); + + expected = new Float32Array( data.y_out ); + + out = sspmv( 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, 'float32' ); - x = ones( 3, 'float32' ); - y = ones( 3, 'float32' ); + data = rxnyn; - expected = new Float32Array( [ 5.0, 5.0, 5.0 ] ); + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); + + expected = new Float32Array( data.y_out ); - out = sspmv( 'column-major', 'lower', 3, 0.0, ap, x, 1, 0, 5.0, y, -1, 2 ); + out = sspmv( 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 Float32Array( [ 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; + + data = cxnyn; + + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); + + expected = new Float32Array( data.y_out ); - out = sspmv( 'column-major', 'lower', 3, 0.0, ap, x, 1, 0, 0.0, y, 1, 0 ); + out = sspmv( 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 Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); + + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); + + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); + + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); + + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); + + expected = new Float32Array( data.y_out ); + + out = sspmv( 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` (column-major)', function test( t ) { + var expected; + var data; + var out; + var ap; + var x; + var y; + + data = csap; + + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); + + expected = new Float32Array( data.y_out ); + + out = sspmv( 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` (row-major)', function test( t ) { + var expected; + var data; + var out; + var ap; + var x; + var y; + + data = rsapn; + + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); + + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 negative stride for `AP` (column-major)', function test( t ) { var expected; + var data; var out; var ap; var x; var y; - ap = new Float32Array( rxnyn.AP ); - x = new Float32Array( rxnyn.x ); - y = new Float32Array( rxnyn.y ); + data = csapn; - expected = new Float32Array( rxnyn.y_out ); + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - out = sspmv( rxnyn.order, rxnyn.uplo, rxnyn.N, rxnyn.alpha, ap, x, rxnyn.strideX, rxnyn.offsetX, rxnyn.beta, y, rxnyn.strideY, rxnyn.offsetY ); + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 an offset for `AP` (row-major)', function test( t ) { var expected; + var data; var out; var ap; var x; var y; - ap = new Float32Array( cxnyn.AP ); - x = new Float32Array( cxnyn.x ); - y = new Float32Array( cxnyn.y ); + data = roap; + + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); + + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); + + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); + + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - expected = new Float32Array( cxnyn.y_out ); + expected = new Float32Array( data.y_out ); - out = sspmv( cxnyn.order, cxnyn.uplo, cxnyn.N, cxnyn.alpha, ap, x, cxnyn.strideX, cxnyn.offsetX, cxnyn.beta, y, cxnyn.strideY, cxnyn.offsetY ); + out = sspmv( 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(); }); diff --git a/lib/node_modules/@stdlib/blas/base/sspmv/test/test.sspmv.js b/lib/node_modules/@stdlib/blas/base/sspmv/test/test.sspmv.js index 605e1a626986..b619fb269833 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/test/test.sspmv.js +++ b/lib/node_modules/@stdlib/blas/base/sspmv/test/test.sspmv.js @@ -24,56 +24,27 @@ var tape = require( 'tape' ); var Float32Array = require( '@stdlib/array/float32' ); -var EPS = require( '@stdlib/constants/float32/eps' ); -var abs = require( '@stdlib/math/base/special/abs' ); var ones = require( '@stdlib/array/ones' ); var sspmv = require( './../lib/sspmv.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() { - sspmv( value, rxpyp.uplo, rxpyp.N, rxpyp.alpha, new Float32Array( rxpyp.AP ), new Float32Array( rxpyp.x ), rxpyp.strideX, rxpyp.beta, new Float32Array( rxpyp.y ), rxpyp.strideY ); + sspmv( value, data.uplo, data.N, data.alpha, new Float32Array( data.AP ), new Float32Array( data.x ), data.strideX, data.beta, new Float32Array( 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() { - sspmv( rxpyp.order, value, rxpyp.N, rxpyp.alpha, new Float32Array( rxpyp.AP ), new Float32Array( rxpyp.x ), rxpyp.strideX, rxpyp.beta, new Float32Array( rxpyp.y ), rxpyp.strideY ); + sspmv( data.order, value, data.N, data.alpha, new Float32Array( data.AP ), new Float32Array( data.x ), data.strideX, data.beta, new Float32Array( 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() { - sspmv( rxpyp.order, rxpyp.uplo, value, rxpyp.alpha, new Float32Array( rxpyp.AP ), new Float32Array( rxpyp.x ), rxpyp.strideX, rxpyp.beta, new Float32Array( rxpyp.y ), rxpyp.strideY ); + sspmv( data.order, data.uplo, value, data.alpha, new Float32Array( data.AP ), new Float32Array( data.x ), data.strideX, data.beta, new Float32Array( 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() { - sspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, new Float32Array( rxpyp.AP ), new Float32Array( rxpyp.x ), value, rxpyp.beta, new Float32Array( rxpyp.y ), rxpyp.strideY ); + sspmv( data.order, data.uplo, data.N, data.alpha, new Float32Array( data.AP ), new Float32Array( data.x ), value, data.beta, new Float32Array( 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,236 +176,255 @@ tape( 'the function throws an error if provided an invalid tenth argument', func function badValue( value ) { return function badValue() { - sspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, new Float32Array( rxpyp.AP ), new Float32Array( rxpyp.x ), rxpyp.strideX, rxpyp.beta, new Float32Array( rxpyp.y ), value ); + sspmv( data.order, data.uplo, data.N, data.alpha, new Float32Array( data.AP ), new Float32Array( data.x ), data.strideX, data.beta, new Float32Array( 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 Float32Array( rxpyp.AP ); - x = new Float32Array( rxpyp.x ); - y = new Float32Array( rxpyp.y ); + data = rl; + + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - expected = new Float32Array( rxpyp.y_out ); + expected = new Float32Array( data.y_out ); - out = sspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, ap, x, rxpyp.strideX, rxpyp.beta, y, rxpyp.strideY ); + out = sspmv( 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 Float32Array( cxpyp.AP ); - x = new Float32Array( cxpyp.x ); - y = new Float32Array( cxpyp.y ); + data = cl; - expected = new Float32Array( cxpyp.y_out ); + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - out = sspmv( cxpyp.order, cxpyp.uplo, cxpyp.N, cxpyp.alpha, ap, x, cxpyp.strideX, cxpyp.beta, y, cxpyp.strideY ); + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 Float32Array( rxoyt.AP ); - x = new Float32Array( rxoyt.x ); - y = new Float32Array( rxoyt.y ); + data = ru; + + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - expected = new Float32Array( rxoyt.y_out ); + expected = new Float32Array( data.y_out ); - out = sspmv( rxoyt.order, rxoyt.uplo, rxoyt.N, rxoyt.alpha, ap, x, rxoyt.strideX, rxoyt.beta, y, rxoyt.strideY ); + out = sspmv( 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 Float32Array( cxoyt.AP ); - x = new Float32Array( cxoyt.x ); - y = new Float32Array( cxoyt.y ); + data = cu; - expected = new Float32Array( cxoyt.y_out ); + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - out = sspmv( cxoyt.order, cxoyt.uplo, cxoyt.N, cxoyt.alpha, ap, x, cxoyt.strideX, cxoyt.beta, y, cxoyt.strideY ); + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 Float32Array( rxpyn.AP ); - x = new Float32Array( rxpyn.x ); - y = new Float32Array( rxpyn.y ); + data = rl; - expected = new Float32Array( rxpyn.y_out ); + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - out = sspmv( rxpyn.order, rxpyn.uplo, rxpyn.N, rxpyn.alpha, ap, x, rxpyn.strideX, rxpyn.beta, y, rxpyn.strideY ); + out = sspmv( 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 Float32Array( cxpyn.AP ); - x = new Float32Array( cxpyn.x ); - y = new Float32Array( cxpyn.y ); + data = rl; + + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - expected = new Float32Array( cxpyn.y_out ); + expected = new Float32Array( data.y ); - out = sspmv( cxpyn.order, cxpyn.uplo, cxpyn.N, cxpyn.alpha, ap, x, cxpyn.strideX, cxpyn.beta, y, cxpyn.strideY ); + out = sspmv( 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 = sspmv( 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 (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 Float32Array( rxnyp.AP ); - x = new Float32Array( rxnyp.x ); - y = new Float32Array( rxnyp.y ); + data = cl; + + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - expected = new Float32Array( rxnyp.y_out ); + expected = new Float32Array( data.y ); - out = sspmv( rxnyp.order, rxnyp.uplo, rxnyp.N, rxnyp.alpha, ap, x, rxnyp.strideX, rxnyp.beta, y, rxnyp.strideY ); + out = sspmv( 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 = sspmv( 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 Float32Array( cxnyp.AP ); - x = new Float32Array( cxnyp.x ); - y = new Float32Array( cxnyp.y ); + ap = ones( 6, 'float32' ); + x = ones( 3, 'float32' ); + y = ones( 3, 'float32' ); - expected = new Float32Array( cxnyp.y_out ); + expected = new Float32Array( [ 5.0, 5.0, 5.0 ] ); - out = sspmv( cxnyp.order, cxnyp.uplo, cxnyp.N, cxnyp.alpha, ap, x, cxnyp.strideX, cxnyp.beta, y, cxnyp.strideY ); + out = sspmv( '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 Float32Array( rxpyp.AP ); - x = new Float32Array( rxpyp.x ); - y = new Float32Array( rxpyp.y ); + expected = new Float32Array( [ 0.0, 0.0, 0.0 ] ); - out = sspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, ap, x, rxpyp.strideX, rxpyp.beta, y, rxpyp.strideY ); + out = sspmv( '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 Float32Array( rxpyp.AP ); - x = new Float32Array( rxpyp.x ); - y = new Float32Array( rxpyp.y ); + ap = ones( 6, 'float32' ); + x = ones( 3, 'float32' ); + y = ones( 3, 'float32' ); - expected = new Float32Array( rxpyp.y ); + expected = new Float32Array( [ 5.0, 5.0, 5.0 ] ); - out = sspmv( rxpyp.order, rxpyp.uplo, 0, rxpyp.alpha, ap, x, rxpyp.strideX, rxpyp.beta, y, rxpyp.strideY ); + out = sspmv( '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 = sspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, 0.0, ap, x, rxpyp.strideX, 1.0, y, rxpyp.strideY ); + expected = new Float32Array( [ 0.0, 0.0, 0.0 ] ); + + out = sspmv( '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 Float32Array( cxpyp.AP ); - x = new Float32Array( cxpyp.x ); - y = new Float32Array( cxpyp.y ); + ap = ones( 6, 'float32' ); + x = ones( 3, 'float32' ); + y = ones( 3, 'float32' ); - expected = new Float32Array( cxpyp.y ); + expected = new Float32Array( [ 5.0, 5.0, 5.0 ] ); - out = sspmv( cxpyp.order, cxpyp.uplo, 0, cxpyp.alpha, ap, x, cxpyp.strideX, cxpyp.beta, y, cxpyp.strideY ); + out = sspmv( '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 = sspmv( cxpyp.order, cxpyp.uplo, cxpyp.N, 0.0, ap, x, cxpyp.strideX, 1.0, y, cxpyp.strideY ); + expected = new Float32Array( [ 0.0, 0.0, 0.0 ] ); + + out = sspmv( '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; @@ -432,133 +437,199 @@ tape( 'when `α` is zero, the function scales the second input vector (row-major expected = new Float32Array( [ 5.0, 5.0, 5.0 ] ); - out = sspmv( 'row-major', 'upper', 3, 0.0, ap, x, 1, 5.0, y, 1 ); + out = sspmv( '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 Float32Array( [ 0.0, 0.0, 0.0 ] ); - out = sspmv( 'row-major', 'upper', 3, 0.0, ap, x, 1, 0.0, y, -1 ); + out = sspmv( '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, 'float32' ); - x = ones( 3, 'float32' ); - y = ones( 3, 'float32' ); + data = rxpyp; - expected = new Float32Array( [ 5.0, 5.0, 5.0 ] ); + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - out = sspmv( 'row-major', 'lower', 3, 0.0, ap, x, 1, 5.0, y, -1 ); + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 Float32Array( [ 0.0, 0.0, 0.0 ] ); + t.end(); +}); - out = sspmv( '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 Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); + + expected = new Float32Array( data.y_out ); + + out = sspmv( 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, 'float32' ); - x = ones( 3, 'float32' ); - y = ones( 3, 'float32' ); + data = rxpyn; - expected = new Float32Array( [ 5.0, 5.0, 5.0 ] ); + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - out = sspmv( 'column-major', 'upper', 3, 0.0, ap, x, 1, 5.0, y, 1 ); + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 Float32Array( [ 0.0, 0.0, 0.0 ] ); + t.end(); +}); - out = sspmv( '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 Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); + + expected = new Float32Array( data.y_out ); + + out = sspmv( 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, 'float32' ); - x = ones( 3, 'float32' ); - y = ones( 3, 'float32' ); + data = rxnyp; - expected = new Float32Array( [ 5.0, 5.0, 5.0 ] ); + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - out = sspmv( 'column-major', 'lower', 3, 0.0, ap, x, 1, 5.0, y, -1 ); + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 Float32Array( [ 0.0, 0.0, 0.0 ] ); + t.end(); +}); - out = sspmv( '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 Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); + + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 Float32Array( rxnyn.AP ); - x = new Float32Array( rxnyn.x ); - y = new Float32Array( rxnyn.y ); + data = rxnyn; - expected = new Float32Array( rxnyn.y_out ); + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - out = sspmv( rxnyn.order, rxnyn.uplo, rxnyn.N, rxnyn.alpha, ap, x, rxnyn.strideX, rxnyn.beta, y, rxnyn.strideY ); + expected = new Float32Array( data.y_out ); + + out = sspmv( 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 Float32Array( cxnyn.AP ); - x = new Float32Array( cxnyn.x ); - y = new Float32Array( cxnyn.y ); + data = cxnyn; + + ap = new Float32Array( data.AP ); + x = new Float32Array( data.x ); + y = new Float32Array( data.y ); - expected = new Float32Array( cxnyn.y_out ); + expected = new Float32Array( data.y_out ); - out = sspmv( cxnyn.order, cxnyn.uplo, cxnyn.N, cxnyn.alpha, ap, x, cxnyn.strideX, cxnyn.beta, y, cxnyn.strideY ); + out = sspmv( 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(); }); From fa1403a10c8ee02d11f1f59efcb06c010d71b0ae Mon Sep 17 00:00:00 2001 From: aman-095 Date: Wed, 28 Aug 2024 13:16:52 +0530 Subject: [PATCH 6/6] docs: update descriptions --- .../@stdlib/blas/base/sspmv/docs/repl.txt | 14 +- .../blas/base/sspmv/docs/types/index.d.ts | 18 +- .../blas/base/sspmv/docs/types/test.ts | 258 ++++++++++-------- 3 files changed, 166 insertions(+), 124 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/base/sspmv/docs/repl.txt b/lib/node_modules/@stdlib/blas/base/sspmv/docs/repl.txt index 424383dd11f1..a5af0798faf5 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/docs/repl.txt +++ b/lib/node_modules/@stdlib/blas/base/sspmv/docs/repl.txt @@ -72,12 +72,12 @@ > AP = new {{alias:@stdlib/array/float32}}( [ 1.0, 2.0, 3.0 ] ); > var x1 = new {{alias:@stdlib/array/float32}}( x0.buffer, x0.BYTES_PER_ELEMENT*0 ); > var y1 = new {{alias:@stdlib/array/float32}}( y0.buffer, y0.BYTES_PER_ELEMENT*0 ); - > {{alias}}( 'row-major', 'upper', 2, 1.0, AP, x1, -1, 1.0, y1, -1 ) + > {{alias}}( 'row-major', 'upper', 2, 1.0, AP, x1, -1, 1.0, y1, -1 ); > y0 [ ~6.0, ~4.0 ] -{{alias}}.ndarray( ord, uplo, N, α, AP, 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,6 +106,12 @@ AP: Float32Array Matrix in packed form. + sap: integer + Index increment for `AP`. + + oap: integer + Starting index for `AP`. + x: Float32Array Input vector `x`. @@ -140,14 +146,14 @@ > var AP = new {{alias:@stdlib/array/float32}}( [ 1.0, 2.0, 3.0, 4.0 ] ); > var ord = 'row-major'; > var uplo = 'upper'; - > {{alias}}.ndarray( ord, uplo, 2, 1.0, AP, 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 ) [ ~4.0, ~6.0 ] // Advanced indexing: > x = new {{alias:@stdlib/array/float32}}( [ 1.0, 1.0 ] ); > y = new {{alias:@stdlib/array/float32}}( [ 1.0, 1.0 ] ); > AP = new {{alias:@stdlib/array/float32}}( [ 1.0, 2.0, 3.0, 4.0 ] ); - > {{alias}}.ndarray( ord, uplo, 2, 1.0, AP, 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 ) [ ~6.0, ~4.0 ] See Also diff --git a/lib/node_modules/@stdlib/blas/base/sspmv/docs/types/index.d.ts b/lib/node_modules/@stdlib/blas/base/sspmv/docs/types/index.d.ts index cf08d1ccba8d..4de375e92aa7 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/docs/types/index.d.ts +++ b/lib/node_modules/@stdlib/blas/base/sspmv/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,20 +54,22 @@ interface Routine { ( order: Layout, uplo: MatrixTriangle, N: number, alpha: number, AP: Float32Array, x: Float32Array, strideX: number, beta: number, y: Float32Array, strideY: number ): Float32Array; /** - * 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 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 @@ -77,14 +79,14 @@ interface Routine { * var x = new Float32Array( [ 1.0, 1.0, 1.0 ] ); * var y = new Float32Array( [ 1.0, 1.0, 1.0 ] ); * - * sspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, x, 1, 0, 1.0, y, 1, 0 ); + * sspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); * // y => [ 7.0, 12.0, 15.0 ] */ - ndarray( order: Layout, uplo: MatrixTriangle, N: number, alpha: number, AP: Float32Array, x: Float32Array, strideX: number, offsetX: number, beta: number, y: Float32Array, strideY: number, offsetY: number ): Float32Array; + ndarray( order: Layout, uplo: MatrixTriangle, N: number, alpha: number, AP: Float32Array, strideAP: number, offsetAP: number, x: Float32Array, strideX: number, offsetX: number, beta: number, y: Float32Array, strideY: number, offsetY: number ): Float32Array; } /** -* 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 @@ -115,7 +117,7 @@ interface Routine { * var x = new Float32Array( [ 1.0, 1.0, 1.0 ] ); * var y = new Float32Array( [ 1.0, 1.0, 1.0 ] ); * -* sspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, x, 1, 0, 1.0, y, -1, 2 ); +* sspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 1, 0, x, 1, 0, 1.0, y, -1, 2 ); * // y => [ 15.0, 12.0, 7.0 ] */ declare var sspmv: Routine; diff --git a/lib/node_modules/@stdlib/blas/base/sspmv/docs/types/test.ts b/lib/node_modules/@stdlib/blas/base/sspmv/docs/types/test.ts index 930e5f394280..f55d492cdb5c 100644 --- a/lib/node_modules/@stdlib/blas/base/sspmv/docs/types/test.ts +++ b/lib/node_modules/@stdlib/blas/base/sspmv/docs/types/test.ts @@ -215,7 +215,7 @@ import sspmv = require( './index' ); const y = new Float32Array( 10 ); const A = new Float32Array( 55 ); - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectType Float32Array + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectType Float32Array } // The compiler throws an error if the function is provided a first argument which is not a string... @@ -224,14 +224,14 @@ import sspmv = require( './index' ); const y = new Float32Array( 10 ); const A = new Float32Array( 55 ); - sspmv.ndarray( 10, 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( true, 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( false, 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( null, 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( void 0, 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( [], 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( {}, 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( ( x: number ): number => x, 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 10, 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( true, 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( false, 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( null, 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( void 0, 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( [], 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( {}, 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.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 sspmv = require( './index' ); const y = new Float32Array( 10 ); const A = new Float32Array( 55 ); - sspmv.ndarray( 'row-major', 10, 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', true, 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', false, 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', null, 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', void 0, 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', [ '1' ], 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', {}, 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', ( x: number ): number => x, 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 10, 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', true, 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', false, 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', null, 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', void 0, 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', [ '1' ], 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', {}, 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.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 sspmv = require( './index' ); const y = new Float32Array( 10 ); const A = new Float32Array( 55 ); - sspmv.ndarray( 'row-major', 'upper', '10', 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', true, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', false, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', null, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', void 0, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', [], 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', {}, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', ( x: number ): number => x, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', '10', 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', true, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', false, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', null, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', void 0, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', [], 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', {}, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.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 sspmv = require( './index' ); const y = new Float32Array( 10 ); const A = new Float32Array( 55 ); - sspmv.ndarray( 'row-major', 'upper', 10, '10', A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, true, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, false, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, null, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, void 0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, [], A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, {}, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, ( x: number ): number => x, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, '10', A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, true, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, false, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, null, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, void 0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, [], A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, {}, A, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.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 Float32Array... @@ -287,31 +287,31 @@ import sspmv = require( './index' ); const x = new Float32Array( 10 ); const y = new Float32Array( 10 ); - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, 10, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, '10', x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, true, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, false, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, null, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, void 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, [ '1' ], x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, {}, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, ( x: number ): number => x, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, 10, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, '10', 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, true, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, false, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, null, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, void 0, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, [ '1' ], 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, {}, 1, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.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 Float32Array... +// The compiler throws an error if the function is provided a sixth argument which is not a number... { + const x = new Float32Array( 10 ); const y = new Float32Array( 10 ); const A = new Float32Array( 55 ); - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 10, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, '10', 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, true, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, false, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, null, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, void 0, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, [ '1' ], 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, {}, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, ( x: number ): number => x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, '10', 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, true, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, false, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, null, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, void 0, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, [], 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, {}, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.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 number... @@ -320,30 +320,30 @@ import sspmv = require( './index' ); const y = new Float32Array( 10 ); const A = new Float32Array( 55 ); - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, '10', 0, 1.0, y, 1 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, true, 0, 1.0, y, 1 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, false, 0, 1.0, y, 1 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, null, 0, 1.0, y, 1 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, void 0, 0, 1.0, y, 1 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, [], 0, 1.0, y, 1 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, {}, 0, 1.0, y, 1 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, ( x: number ): number => x, 0, 1.0, y, 1 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, '10', x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, true, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, false, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, null, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, void 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, [], x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, {}, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.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 Float32Array... { - const x = new Float32Array( 10 ); const y = new Float32Array( 10 ); const A = new Float32Array( 55 ); - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, '10', 1.0, y, 1 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, true, 1.0, y, 1 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, false, 1.0, y, 1 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, null, 1.0, y, 1 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, void 0, 1.0, y, 1 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, [], 1.0, y, 1 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, {}, 1.0, y, 1 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, ( x: number ): number => x, 1.0, y, 1 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, 10, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, '10', 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, true, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, false, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, null, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, void 0, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, [ '1' ], 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, {}, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.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,30 +352,30 @@ import sspmv = require( './index' ); const y = new Float32Array( 10 ); const A = new Float32Array( 55 ); - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, '10', y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, true, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, false, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, null, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, void 0, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, [], y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, {}, y, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, ( x: number ): number => x, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, '10', 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, true, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, false, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, null, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, void 0, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, [], 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, {}, 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.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 Float32Array... +// The compiler throws an error if the function is provided a tenth argument which is not a number... { const x = new Float32Array( 10 ); + const y = new Float32Array( 10 ); const A = new Float32Array( 55 ); - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, 10, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, '10', 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, true, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, false, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, null, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, void 0, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, [ '1' ], 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, {}, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, ( x: number ): number => x, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, '10', 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, true, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, false, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, null, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, void 0, 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, [], 1.0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, {}, 1.0, y, 1, 0 ); // $ExpectError + sspmv.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 number... @@ -384,30 +384,62 @@ import sspmv = require( './index' ); const y = new Float32Array( 10 ); const A = new Float32Array( 55 ); - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, '10', 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, true, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, false, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, null, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, void 0, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, [], 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, {}, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, ( x: number ): number => x, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, '10', y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, true, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, false, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, null, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, void 0, y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, [], y, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, {}, y, 1, 0 ); // $ExpectError + sspmv.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 Float32Array... +{ + const x = new Float32Array( 10 ); + const A = new Float32Array( 55 ); + + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, 10, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, '10', 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, true, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, false, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, null, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, void 0, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, [ '1' ], 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, {}, 1, 0 ); // $ExpectError + sspmv.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... +{ + const x = new Float32Array( 10 ); + const y = new Float32Array( 10 ); + const A = new Float32Array( 55 ); + + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, '10', 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, true, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, false, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, null, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, void 0, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, [], 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, {}, 0 ); // $ExpectError + sspmv.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 twelfth argument which is not a number... +// The compiler throws an error if the function is provided a fourteenth argument which is not a number... { const x = new Float32Array( 10 ); const y = new Float32Array( 10 ); const A = new Float32Array( 55 ); - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, '10' ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, true ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, false ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, null ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, void 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, [] ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, {} ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, ( x: number ): number => x ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, '10' ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, true ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, false ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, null ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, void 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, [] ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, {} ); // $ExpectError + sspmv.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... @@ -422,11 +454,13 @@ import sspmv = require( './index' ); sspmv.ndarray( 'row-major', 'upper', 10 ); // $ExpectError sspmv.ndarray( 'row-major', 'upper', 10, 1.0 ); // $ExpectError sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1 ); // $ExpectError - sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0, 10 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1 ); // $ExpectError + sspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 1, 0, x, 1, 0, 1.0, y, 1, 0, 10 ); // $ExpectError }