From bc2355981ffa28082d8e23273887501bd3ba674a Mon Sep 17 00:00:00 2001 From: Athan Reines Date: Mon, 1 Jul 2024 00:59:37 -0700 Subject: [PATCH] feat!: add support for specifying the index offset for `AP` BREAKING CHANGE: add offset parameter to `ndarray` method To migrate, users should set the `offsetAP` parameter. For most cases, this parameter will be zero, but supporting this parameter allows users to specify alternative starting indices, such as needed when working with ndarray views. --- .../@stdlib/blas/base/dspmv/README.md | 11 +- .../base/dspmv/benchmark/benchmark.ndarray.js | 2 +- .../@stdlib/blas/base/dspmv/docs/repl.txt | 9 +- .../blas/base/dspmv/docs/types/index.d.ts | 7 +- .../blas/base/dspmv/docs/types/test.ts | 243 ++++++++++-------- .../@stdlib/blas/base/dspmv/examples/index.js | 2 +- .../@stdlib/blas/base/dspmv/lib/base.js | 9 +- .../@stdlib/blas/base/dspmv/lib/dspmv.js | 2 +- .../@stdlib/blas/base/dspmv/lib/index.js | 2 +- .../@stdlib/blas/base/dspmv/lib/ndarray.js | 17 +- .../test/fixtures/column_major_xnyn.json | 3 +- .../test/fixtures/column_major_xnyp.json | 1 + .../test/fixtures/column_major_xoyt.json | 1 + .../test/fixtures/column_major_xpyn.json | 1 + .../test/fixtures/column_major_xpyp.json | 1 + .../dspmv/test/fixtures/row_major_xnyn.json | 3 +- .../dspmv/test/fixtures/row_major_xnyp.json | 1 + .../dspmv/test/fixtures/row_major_xoyt.json | 1 + .../dspmv/test/fixtures/row_major_xpyn.json | 1 + .../dspmv/test/fixtures/row_major_xpyp.json | 1 + .../blas/base/dspmv/test/test.dspmv.js | 4 +- .../blas/base/dspmv/test/test.ndarray.js | 64 ++--- 22 files changed, 210 insertions(+), 176 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/README.md b/lib/node_modules/@stdlib/blas/base/dspmv/README.md index 7d8149aa8e9..f4986fba5b5 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/README.md +++ b/lib/node_modules/@stdlib/blas/base/dspmv/README.md @@ -91,7 +91,7 @@ dspmv( 'row-major', 'upper', 2, 1.0, AP, x1, -1, 1.0, y1, -1 ); // y0 => [ 6.0, 4.0 ] ``` -#### dspmv.ndarray( order, uplo, N, α, AP, x, sx, ox, β, y, sy, oy ) +#### dspmv.ndarray( order, uplo, N, α, AP, oa, 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,13 @@ var AP = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); var x = new Float64Array( [ 1.0, 1.0, 1.0 ] ); var y = new Float64Array( [ 1.0, 1.0, 1.0 ] ); -dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, x, 1, 0, 1.0, y, 1, 0 ); +dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 0, x, 1, 0, 1.0, y, 1, 0 ); // y => [ 7.0, 12.0, 15.0 ] ``` The function has the following additional parameters: +- **oa**: starting index for `AP`. - **ox**: starting index for `x`. - **oy**: starting index for `y`. @@ -116,11 +117,11 @@ While [`typed array`][mdn-typed-array] views mandate a view offset based on the ```javascript var Float64Array = require( '@stdlib/array/float64' ); -var AP = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); +var AP = new Float64Array( [ 0.0, 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); var x = new Float64Array( [ 1.0, 1.0, 1.0 ] ); var y = new Float64Array( [ 1.0, 1.0, 1.0 ] ); -dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, x, 1, 0, 1.0, y, -1, 2 ); +dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 2, x, 1, 0, 1.0, y, -1, 2 ); // y => [ 15.0, 12.0, 7.0 ] ``` @@ -158,7 +159,7 @@ var AP = discreteUniform( N * ( N + 1 ) / 2, -10, 10, opts ); var x = discreteUniform( N, -10, 10, opts ); var y = discreteUniform( N, -10, 10, opts ); -dspmv.ndarray( 'row-major', 'upper', N, 1.0, AP, x, 1, 0, 1.0, y, 1, 0 ); +dspmv.ndarray( 'row-major', 'upper', N, 1.0, AP, 0, x, 1, 0, 1.0, y, 1, 0 ); console.log( y ); ``` diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/benchmark/benchmark.ndarray.js b/lib/node_modules/@stdlib/blas/base/dspmv/benchmark/benchmark.ndarray.js index 10433a29453..a2877463004 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/benchmark/benchmark.ndarray.js +++ b/lib/node_modules/@stdlib/blas/base/dspmv/benchmark/benchmark.ndarray.js @@ -63,7 +63,7 @@ function createBenchmark( len ) { b.tic(); for ( i = 0; i < b.iterations; i++ ) { - z = dspmv( 'row-major', 'upper', len, 1.0, AP, x, 1, 0, 1.0, y, 1, 0 ); + z = dspmv( 'row-major', 'upper', len, 1.0, AP, 0, x, 1, 0, 1.0, y, 1, 0 ); if ( isnan( z[ i%z.length ] ) ) { b.fail( 'should not return NaN' ); } diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/docs/repl.txt b/lib/node_modules/@stdlib/blas/base/dspmv/docs/repl.txt index 9cd3209ff83..5340eeefa55 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/docs/repl.txt +++ b/lib/node_modules/@stdlib/blas/base/dspmv/docs/repl.txt @@ -77,7 +77,7 @@ [ ~6.0, ~4.0 ] -{{alias}}.ndarray( ord, uplo, N, α, AP, x, sx, ox, β, y, sy, oy ) +{{alias}}.ndarray( ord, uplo, N, α, AP, oa, 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,9 @@ AP: Float64Array Matrix in packed form. + oa: integer + Starting index for `AP`. + x: Float64Array Input vector `x`. @@ -140,14 +143,14 @@ > var AP = new {{alias:@stdlib/array/float64}}( [ 1.0, 2.0, 3.0, 4.0 ] ); > var ord = 'row-major'; > var uplo = 'upper'; - > {{alias}}.ndarray( ord, uplo, 2, 1.0, AP, x, 1, 0, 1.0, y, 1, 0 ) + > {{alias}}.ndarray( ord, uplo, 2, 1.0, AP, 0, x, 1, 0, 1.0, y, 1, 0 ) [ ~4.0, ~6.0 ] // Advanced indexing: > x = new {{alias:@stdlib/array/float64}}( [ 1.0, 1.0 ] ); > y = new {{alias:@stdlib/array/float64}}( [ 1.0, 1.0 ] ); > AP = new {{alias:@stdlib/array/float64}}( [ 1.0, 2.0, 3.0, 4.0 ] ); - > {{alias}}.ndarray( ord, uplo, 2, 1.0, AP, x, -1, 1, 1.0, y, -1, 1 ) + > {{alias}}.ndarray( ord, uplo, 2, 1.0, AP, 0, x, -1, 1, 1.0, y, -1, 1 ) [ ~6.0, ~4.0 ] See Also diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/docs/types/index.d.ts b/lib/node_modules/@stdlib/blas/base/dspmv/docs/types/index.d.ts index 37627586c98..e755ce04ec5 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/docs/types/index.d.ts +++ b/lib/node_modules/@stdlib/blas/base/dspmv/docs/types/index.d.ts @@ -61,6 +61,7 @@ interface Routine { * @param N - number of columns in the matrix `A` * @param alpha - scalar constant * @param AP - packed form of a symmetric matrix `A` + * @param offsetAP - starting `AP` index * @param x - first input array * @param strideX - `x` stride length * @param offsetX - starting `x` index @@ -77,10 +78,10 @@ interface Routine { * var x = new Float64Array( [ 1.0, 1.0, 1.0 ] ); * var y = new Float64Array( [ 1.0, 1.0, 1.0 ] ); * - * dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, x, 1, 0, 1.0, y, 1, 0 ); + * dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 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: Float64Array, x: Float64Array, strideX: number, offsetX: number, beta: number, y: Float64Array, strideY: number, offsetY: number ): Float64Array; + ndarray( order: Layout, uplo: MatrixTriangle, N: number, alpha: number, AP: Float64Array, offsetAP: number, x: Float64Array, strideX: number, offsetX: number, beta: number, y: Float64Array, strideY: number, offsetY: number ): Float64Array; } /** @@ -115,7 +116,7 @@ interface Routine { * var x = new Float64Array( [ 1.0, 1.0, 1.0 ] ); * var y = new Float64Array( [ 1.0, 1.0, 1.0 ] ); * -* dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, x, 1, 0, 1.0, y, -1, 2 ); +* dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 0, x, 1, 0, 1.0, y, -1, 2 ); * // y => [ 15.0, 12.0, 7.0 ] */ declare var dspmv: Routine; diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/docs/types/test.ts b/lib/node_modules/@stdlib/blas/base/dspmv/docs/types/test.ts index 1aefcc3f2df..d717efff965 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/docs/types/test.ts +++ b/lib/node_modules/@stdlib/blas/base/dspmv/docs/types/test.ts @@ -215,7 +215,7 @@ import dspmv = require( './index' ); const y = new Float64Array( 10 ); const A = new Float64Array( 55 ); - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectType Float64Array + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectType Float64Array } // The compiler throws an error if the function is provided a first argument which is not a string... @@ -224,14 +224,14 @@ import dspmv = require( './index' ); const y = new Float64Array( 10 ); const A = new Float64Array( 55 ); - dspmv.ndarray( 10, 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( true, 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( false, 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( null, 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( void 0, 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( [], 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( {}, 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( ( x: number ): number => x, 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 10, 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( true, 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( false, 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( null, 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( void 0, 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( [], 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( {}, 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( ( x: number ): number => x, 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError } // The compiler throws an error if the function is provided a second argument which is not a string... @@ -240,14 +240,14 @@ import dspmv = require( './index' ); const y = new Float64Array( 10 ); const A = new Float64Array( 55 ); - dspmv.ndarray( 'row-major', 10, 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', true, 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', false, 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', null, 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', void 0, 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', [ '1' ], 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', {}, 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', ( x: number ): number => x, 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 10, 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', true, 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', false, 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', null, 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', void 0, 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', [ '1' ], 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', {}, 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', ( x: number ): number => x, 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError } // The compiler throws an error if the function is provided a third argument which is not a number... @@ -256,14 +256,14 @@ import dspmv = require( './index' ); const y = new Float64Array( 10 ); const A = new Float64Array( 55 ); - dspmv.ndarray( 'row-major', 'upper', '10', 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', true, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', false, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', null, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', void 0, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', [], 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', {}, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', ( x: number ): number => x, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', '10', 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', true, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', false, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', null, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', void 0, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', [], 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', {}, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', ( x: number ): number => x, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError } // The compiler throws an error if the function is provided a fourth argument which is not a number... @@ -272,14 +272,14 @@ import dspmv = require( './index' ); const y = new Float64Array( 10 ); const A = new Float64Array( 55 ); - dspmv.ndarray( 'row-major', 'upper', 10, '10', A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, true, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, false, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, null, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, void 0, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, [], A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, {}, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, ( x: number ): number => x, A, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, '10', A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, true, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, false, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, null, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, void 0, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, [], A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, {}, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, ( x: number ): number => x, A, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError } // The compiler throws an error if the function is provided a fifth argument which is not a Float64Array... @@ -287,47 +287,47 @@ import dspmv = require( './index' ); const x = new Float64Array( 10 ); const y = new Float64Array( 10 ); - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, 10, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, '10', x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, true, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, false, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, null, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, void 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, [ '1' ], x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, {}, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, ( x: number ): number => x, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, 10, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, '10', 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, true, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, false, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, null, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, void 0, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, [ '1' ], 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, {}, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, ( x: number ): number => x, 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError } -// The compiler throws an error if the function is provided a sixth argument which is not a Float64Array... +// The compiler throws an error if the function is provided a sixth argument which is not a number... { + const x = new Float64Array( 10 ); const y = new Float64Array( 10 ); const A = new Float64Array( 55 ); - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 10, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, '10', 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, true, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, false, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, null, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, void 0, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, [ '1' ], 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, {}, 1, 0, 1.0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, ( x: number ): number => x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, '10', x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, true, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, false, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, null, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, void 0, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, [ '1' ], x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, {}, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, ( x: number ): number => x, x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError } -// The compiler throws an error if the function is provided a seventh argument which is not a number... +// The compiler throws an error if the function is provided a seventh argument which is not a Float64Array... { - const x = new Float64Array( 10 ); const y = new Float64Array( 10 ); const A = new Float64Array( 55 ); - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, '10', 0, 1.0, y, 1 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, true, 0, 1.0, y, 1 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, false, 0, 1.0, y, 1 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, null, 0, 1.0, y, 1 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, void 0, 0, 1.0, y, 1 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, [], 0, 1.0, y, 1 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, {}, 0, 1.0, y, 1 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, ( x: number ): number => x, 0, 1.0, y, 1 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, 10, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, '10', 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, true, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, false, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, null, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, void 0, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, [ '1' ], 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, {}, 1, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, ( x: number ): number => x, 1, 0, 1.0, y, 1, 0 ); // $ExpectError } // The compiler throws an error if the function is provided an eighth argument which is not a number... @@ -336,14 +336,14 @@ import dspmv = require( './index' ); const y = new Float64Array( 10 ); const A = new Float64Array( 55 ); - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, '10', 1.0, y, 1 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, true, 1.0, y, 1 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, false, 1.0, y, 1 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, null, 1.0, y, 1 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, void 0, 1.0, y, 1 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, [], 1.0, y, 1 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, {}, 1.0, y, 1 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, ( x: number ): number => x, 1.0, y, 1 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, '10', 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, true, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, false, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, null, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, void 0, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, [], 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, {}, 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, ( x: number ): number => x, 0, 1.0, y, 1, 0 ); // $ExpectError } // The compiler throws an error if the function is provided a ninth argument which is not a number... @@ -352,46 +352,46 @@ import dspmv = require( './index' ); const y = new Float64Array( 10 ); const A = new Float64Array( 55 ); - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, '10', y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, true, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, false, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, null, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, void 0, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, [], y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, {}, y, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, ( x: number ): number => x, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, '10', 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, true, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, false, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, null, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, void 0, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, [], 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, {}, 1.0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, ( x: number ): number => x, 1.0, y, 1, 0 ); // $ExpectError } -// The compiler throws an error if the function is provided a tenth argument which is not a Float64Array... +// The compiler throws an error if the function is provided a tenth argument which is not a number... { const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); const A = new Float64Array( 55 ); - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, 10, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, '10', 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, true, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, false, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, null, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, void 0, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, [ '1' ], 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, {}, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, ( x: number ): number => x, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, '10', y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, true, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, false, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, null, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, void 0, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, [], y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, {}, y, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, ( x: number ): number => x, y, 1, 0 ); // $ExpectError } -// The compiler throws an error if the function is provided an eleventh argument which is not a number... +// The compiler throws an error if the function is provided an eleventh argument which is not a Float64Array... { const x = new Float64Array( 10 ); - const y = new Float64Array( 10 ); const A = new Float64Array( 55 ); - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, '10', 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, true, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, false, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, null, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, void 0, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, [], 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, {}, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, ( x: number ): number => x, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, 10, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, '10', 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, true, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, false, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, null, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, void 0, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, [ '1' ], 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, {}, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, ( x: number ): number => x, 1, 0 ); // $ExpectError } // The compiler throws an error if the function is provided a twelfth argument which is not a number... @@ -400,14 +400,30 @@ import dspmv = require( './index' ); const y = new Float64Array( 10 ); const A = new Float64Array( 55 ); - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, '10' ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, true ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, false ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, null ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, void 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, [] ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, {} ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, ( x: number ): number => x ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, '10', 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, true, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, false, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, null, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, void 0, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, [], 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, {}, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, ( x: number ): number => x, 0 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a thirteenth argument which is not a number... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + const A = new Float64Array( 55 ); + + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, '10' ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, true ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, false ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, null ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, void 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, [] ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, {} ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, ( x: number ): number => x ); // $ExpectError } // The compiler throws an error if the `ndarray` method is provided an unsupported number of arguments... @@ -422,11 +438,12 @@ import dspmv = require( './index' ); dspmv.ndarray( 'row-major', 'upper', 10 ); // $ExpectError dspmv.ndarray( 'row-major', 'upper', 10, 1.0 ); // $ExpectError dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1 ); // $ExpectError - dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, x, 1, 0, 1.0, y, 1, 0, 10 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1 ); // $ExpectError + dspmv.ndarray( 'row-major', 'upper', 10, 1.0, A, 0, x, 1, 0, 1.0, y, 1, 0, 10 ); // $ExpectError } diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/examples/index.js b/lib/node_modules/@stdlib/blas/base/dspmv/examples/index.js index 6ea87c25168..96199e8b634 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/examples/index.js +++ b/lib/node_modules/@stdlib/blas/base/dspmv/examples/index.js @@ -31,5 +31,5 @@ var AP = discreteUniform( N * ( N + 1 ) / 2, -10, 10, opts ); var x = discreteUniform( N, -10, 10, opts ); var y = discreteUniform( N, -10, 10, opts ); -dspmv.ndarray( 'row-major', 'upper', N, 1.0, AP, x, 1, 0, 1.0, y, 1, 0 ); +dspmv.ndarray( 'row-major', 'upper', N, 1.0, AP, 0, x, 1, 0, 1.0, y, 1, 0 ); console.log( y ); diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/lib/base.js b/lib/node_modules/@stdlib/blas/base/dspmv/lib/base.js index 37d7c39b464..5b07549d6f0 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/lib/base.js +++ b/lib/node_modules/@stdlib/blas/base/dspmv/lib/base.js @@ -35,6 +35,7 @@ var dscal = require( '@stdlib/blas/base/dscal' ).ndarray; * @param {NonNegativeInteger} N - number of elements along each dimension of `A` * @param {number} alpha - scalar constant * @param {Float64Array} AP - packed form of a symmetric matrix `A` +* @param {NonNegativeInteger} offsetAP - starting `AP` index * @param {Float64Array} x - first input array * @param {integer} strideX - `x` stride length * @param {NonNegativeInteger} offsetX - starting `x` index @@ -47,14 +48,14 @@ var dscal = require( '@stdlib/blas/base/dscal' ).ndarray; * @example * var Float64Array = require( '@stdlib/array/float64' ); * -* var A = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); +* var AP = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); * var x = new Float64Array( [ 1.0, 1.0, 1.0 ] ); * var y = new Float64Array( [ 1.0, 1.0, 1.0 ] ); * -* dspmv( 'column-major', 'lower', 3, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); +* dspmv( 'column-major', 'lower', 3, 1.0, AP, 0, x, 1, 0, 1.0, y, 1, 0 ); * // y => [ ~7.0, ~12.0, ~15.0 ] */ -function dspmv( order, uplo, N, alpha, AP, x, strideX, offsetX, beta, y, strideY, offsetY ) { // eslint-disable-line max-params, max-len +function dspmv( order, uplo, N, alpha, AP, offsetAP, x, strideX, offsetX, beta, y, strideY, offsetY ) { // eslint-disable-line max-params, max-len var temp1; var temp2; var ix; @@ -84,7 +85,7 @@ function dspmv( order, uplo, N, alpha, AP, x, strideX, offsetX, beta, y, strideY // Form: y = alpha*A*x + y kx = offsetX; ky = offsetY; - kk = 0; + kk = offsetAP; if ( ( order === 'row-major' && uplo === 'upper' ) || ( order === 'column-major' && uplo === 'lower' ) diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/lib/dspmv.js b/lib/node_modules/@stdlib/blas/base/dspmv/lib/dspmv.js index 7ddfe3c3648..32b865e6331 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/lib/dspmv.js +++ b/lib/node_modules/@stdlib/blas/base/dspmv/lib/dspmv.js @@ -80,7 +80,7 @@ function dspmv( order, uplo, N, alpha, AP, x, strideX, beta, y, strideY ) { } offsetX = stride2offset( N, strideX ); offsetY = stride2offset( N, strideY ); - return base( order, uplo, N, alpha, AP, x, strideX, offsetX, beta, y, strideY, offsetY ); // eslint-disable-line max-len + return base( order, uplo, N, alpha, AP, 0, x, strideX, offsetX, beta, y, strideY, offsetY ); // eslint-disable-line max-len } diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/lib/index.js b/lib/node_modules/@stdlib/blas/base/dspmv/lib/index.js index 98196e09833..c18159951df 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/lib/index.js +++ b/lib/node_modules/@stdlib/blas/base/dspmv/lib/index.js @@ -42,7 +42,7 @@ * var x = new Float64Array( [ 1.0, 1.0, 1.0 ] ); * var y = new Float64Array( [ 1.0, 1.0, 1.0 ] ); * -* dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, x, 1, 0, 1.0, y, 1, 0 ); +* dspmv.ndarray( 'column-major', 'lower', 3, 1.0, AP, 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/dspmv/lib/ndarray.js b/lib/node_modules/@stdlib/blas/base/dspmv/lib/ndarray.js index 15965075d42..2c5a1479e79 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/lib/ndarray.js +++ b/lib/node_modules/@stdlib/blas/base/dspmv/lib/ndarray.js @@ -36,6 +36,7 @@ var base = require( './base.js' ); * @param {NonNegativeInteger} N - number of elements along each dimension of `A` * @param {number} alpha - scalar constant * @param {Float64Array} AP - packed form of a symmetric matrix `A` +* @param {NonNegativeInteger} offsetAP - starting `AP` index * @param {Float64Array} x - first input array * @param {integer} strideX - `x` stride length * @param {NonNegativeInteger} offsetX - starting `x` index @@ -46,21 +47,21 @@ var base = require( './base.js' ); * @throws {TypeError} first argument must be a valid order * @throws {TypeError} second argument must specify whether the lower or upper triangular matrix is supplied * @throws {RangeError} third argument must be a nonnegative integer -* @throws {RangeError} seventh argument must be non-zero -* @throws {RangeError} eleventh argument must be non-zero +* @throws {RangeError} eighth argument must be non-zero +* @throws {RangeError} twelfth argument must be non-zero * @returns {Float64Array} `y` * * @example * var Float64Array = require( '@stdlib/array/float64' ); * -* var A = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); +* var AP = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); * var x = new Float64Array( [ 1.0, 1.0, 1.0 ] ); * var y = new Float64Array( [ 1.0, 1.0, 1.0 ] ); * -* dspmv( 'column-major', 'lower', 3, 1.0, A, x, 1, 0, 1.0, y, 1, 0 ); +* dspmv( 'column-major', 'lower', 3, 1.0, AP, 0, x, 1, 0, 1.0, y, 1, 0 ); * // y => [ ~7.0, ~12.0, ~15.0 ] */ -function dspmv( order, uplo, N, alpha, AP, x, strideX, offsetX, beta, y, strideY, offsetY ) { // eslint-disable-line max-params, max-len +function dspmv( order, uplo, N, alpha, AP, offsetAP, x, strideX, offsetX, beta, y, strideY, offsetY ) { // eslint-disable-line max-params, max-len if ( !isLayout( order ) ) { throw new TypeError( format( 'invalid argument. First argument must be a valid order. Value: `%s`.', order ) ); } @@ -71,12 +72,12 @@ function dspmv( order, uplo, N, alpha, AP, x, strideX, offsetX, beta, y, strideY throw new RangeError( format( 'invalid argument. Third argument must be a nonnegative integer. Value: `%d`.', N ) ); } if ( strideX === 0 ) { - throw new RangeError( format( 'invalid argument. Seventh argument must be non-zero. Value: `%d`.', strideX ) ); + throw new RangeError( format( 'invalid argument. Eighth argument must be non-zero. Value: `%d`.', strideX ) ); } if ( strideY === 0 ) { - throw new RangeError( format( 'invalid argument. Eleventh argument must be non-zero. Value: `%d`.', strideY ) ); + throw new RangeError( format( 'invalid argument. Twelfth argument must be non-zero. Value: `%d`.', strideY ) ); } - return base( order, uplo, N, alpha, AP, x, strideX, offsetX, beta, y, strideY, offsetY ); // eslint-disable-line max-len + return base( order, uplo, N, alpha, AP, offsetAP, x, strideX, offsetX, beta, y, strideY, offsetY ); // eslint-disable-line max-len } diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xnyn.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xnyn.json index dd535606e05..06b6850679f 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xnyn.json +++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xnyn.json @@ -8,7 +8,8 @@ "offsetX": 2, "strideY": -1, "offsetY": 2, - "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ], + "offsetAP": 2, + "AP": [ 0.0, 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ], "x": [ 1.0, 1.0, 1.0 ], "y": [ 1.0, 1.0, 1.0 ], "y_out": [ 16.0, 11.0, 8.0 ] diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xnyp.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xnyp.json index 88b876d9478..c2d83166074 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xnyp.json +++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xnyp.json @@ -8,6 +8,7 @@ "offsetX": 2, "strideY": 1, "offsetY": 0, + "offsetAP": 0, "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ], "x": [ 1.0, 1.0, 1.0 ], "y": [ 1.0, 1.0, 1.0 ], diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xoyt.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xoyt.json index 7ce888f02a0..b05e6fe81a6 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xoyt.json +++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xoyt.json @@ -8,6 +8,7 @@ "offsetX": 0, "strideY": 2, "offsetY": 0, + "offsetAP": 0, "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ], "x": [ 1.0, 1.0, 1.0 ], "y": [ 1.0, 0.0, 1.0, 0.0, 1.0, 0.0 ], diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xpyn.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xpyn.json index 8fedbd44397..f4aad6d8b84 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xpyn.json +++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xpyn.json @@ -8,6 +8,7 @@ "offsetX": 0, "strideY": -1, "offsetY": 2, + "offsetAP": 0, "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ], "x": [ 1.0, 1.0, 1.0 ], "y": [ 1.0, 1.0, 1.0 ], diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xpyp.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xpyp.json index 0e57578aba8..b3123f8e2f1 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xpyp.json +++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/column_major_xpyp.json @@ -8,6 +8,7 @@ "offsetX": 0, "strideY": 1, "offsetY": 0, + "offsetAP": 0, "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ], "x": [ 1.0, 1.0, 1.0 ], "y": [ 1.0, 1.0, 1.0 ], diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xnyn.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xnyn.json index 4763913667f..f1631104ae5 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xnyn.json +++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xnyn.json @@ -8,7 +8,8 @@ "offsetX": 2, "strideY": -1, "offsetY": 2, - "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ], + "offsetAP": 2, + "AP": [ 0.0, 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ], "x": [ 1.0, 1.0, 1.0 ], "y": [ 1.0, 1.0, 1.0 ], "y_out": [ 15.0, 12.0, 7.0 ] diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xnyp.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xnyp.json index fc5df6d7ce9..475ee4d4919 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xnyp.json +++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xnyp.json @@ -8,6 +8,7 @@ "offsetX": 2, "strideY": 1, "offsetY": 0, + "offsetAP": 0, "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ], "x": [ 1.0, 1.0, 1.0 ], "y": [ 1.0, 1.0, 1.0 ], diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xoyt.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xoyt.json index f77e7da3427..211518d861c 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xoyt.json +++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xoyt.json @@ -8,6 +8,7 @@ "offsetX": 0, "strideY": 2, "offsetY": 0, + "offsetAP": 0, "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ], "x": [ 1.0, 1.0, 1.0 ], "y": [ 1.0, 0.0, 1.0, 0.0, 1.0, 0.0 ], diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xpyn.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xpyn.json index 82387497ca3..2a6711742d6 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xpyn.json +++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xpyn.json @@ -8,6 +8,7 @@ "offsetX": 0, "strideY": -1, "offsetY": 2, + "offsetAP": 0, "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ], "x": [ 1.0, 1.0, 1.0 ], "y": [ 1.0, 1.0, 1.0 ], diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xpyp.json b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xpyp.json index d4deab60a3f..ddee2b15624 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xpyp.json +++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/fixtures/row_major_xpyp.json @@ -8,6 +8,7 @@ "offsetX": 0, "strideY": 1, "offsetY": 0, + "offsetAP": 0, "AP": [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ], "x": [ 1.0, 1.0, 1.0 ], "y": [ 1.0, 1.0, 1.0 ], diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/test.dspmv.js b/lib/node_modules/@stdlib/blas/base/dspmv/test/test.dspmv.js index 22282eee40c..7328a8656fe 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/test/test.dspmv.js +++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/test.dspmv.js @@ -530,7 +530,7 @@ tape( 'the function supports complex access patterns (row-major)', function test var x; var y; - ap = new Float64Array( rxnyn.AP ); + ap = new Float64Array( rxnyn.AP.slice( rxnyn.offsetAP ) ); x = new Float64Array( rxnyn.x ); y = new Float64Array( rxnyn.y ); @@ -550,7 +550,7 @@ tape( 'the function supports complex access patterns (column-major)', function t var x; var y; - ap = new Float64Array( cxnyn.AP ); + ap = new Float64Array( cxnyn.AP.slice( cxnyn.offsetAP ) ); x = new Float64Array( cxnyn.x ); y = new Float64Array( cxnyn.y ); diff --git a/lib/node_modules/@stdlib/blas/base/dspmv/test/test.ndarray.js b/lib/node_modules/@stdlib/blas/base/dspmv/test/test.ndarray.js index 2f446b3d536..e02d81ec309 100644 --- a/lib/node_modules/@stdlib/blas/base/dspmv/test/test.ndarray.js +++ b/lib/node_modules/@stdlib/blas/base/dspmv/test/test.ndarray.js @@ -82,8 +82,8 @@ tape( 'main export is a function', function test( t ) { t.end(); }); -tape( 'the function has an arity of 12', function test( t ) { - t.strictEqual( dspmv.length, 12, 'returns expected value' ); +tape( 'the function has an arity of 13', function test( t ) { + t.strictEqual( dspmv.length, 13, 'returns expected value' ); t.end(); }); @@ -105,7 +105,7 @@ tape( 'the function throws an error if provided an invalid first argument', func function badValue( value ) { return function badValue() { - dspmv( value, rxpyp.uplo, rxpyp.N, rxpyp.alpha, new Float64Array( rxpyp.AP ), new Float64Array( rxpyp.x ), rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, new Float64Array( rxpyp.y ), rxpyp.strideY, rxpyp.offsetY ); + dspmv( value, rxpyp.uplo, rxpyp.N, rxpyp.alpha, new Float64Array( rxpyp.AP ), rxpyp.offsetAP, new Float64Array( rxpyp.x ), rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, new Float64Array( rxpyp.y ), rxpyp.strideY, rxpyp.offsetY ); }; } }); @@ -128,7 +128,7 @@ tape( 'the function throws an error if provided an invalid second argument', fun function badValue( value ) { return function badValue() { - dspmv( rxpyp.order, value, rxpyp.N, rxpyp.alpha, new Float64Array( rxpyp.AP ), new Float64Array( rxpyp.x ), rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, new Float64Array( rxpyp.y ), rxpyp.strideY, rxpyp.offsetY ); + dspmv( rxpyp.order, value, rxpyp.N, rxpyp.alpha, new Float64Array( rxpyp.AP ), rxpyp.offsetAP, new Float64Array( rxpyp.x ), rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, new Float64Array( rxpyp.y ), rxpyp.strideY, rxpyp.offsetY ); }; } }); @@ -150,12 +150,12 @@ tape( 'the function throws an error if provided an invalid third argument', func function badValue( value ) { return function badValue() { - dspmv( rxpyp.order, rxpyp.uplo, value, rxpyp.alpha, new Float64Array( rxpyp.AP ), new Float64Array( rxpyp.x ), rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, new Float64Array( rxpyp.y ), rxpyp.strideY, rxpyp.offsetY ); + dspmv( rxpyp.order, rxpyp.uplo, value, rxpyp.alpha, new Float64Array( rxpyp.AP ), rxpyp.offsetAP, new Float64Array( rxpyp.x ), rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, new Float64Array( rxpyp.y ), rxpyp.strideY, rxpyp.offsetY ); }; } }); -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 eighth argument', function test( t ) { var values; var i; @@ -170,12 +170,12 @@ tape( 'the function throws an error if provided an invalid seventh argument', fu function badValue( value ) { return function badValue() { - dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, new Float64Array( rxpyp.AP ), new Float64Array( rxpyp.x ), value, rxpyp.offsetX, rxpyp.beta, new Float64Array( rxpyp.y ), rxpyp.strideY, rxpyp.offsetY ); + dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, new Float64Array( rxpyp.AP ), rxpyp.offsetAP, new Float64Array( rxpyp.x ), value, rxpyp.offsetX, rxpyp.beta, new Float64Array( rxpyp.y ), rxpyp.strideY, rxpyp.offsetY ); }; } }); -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 twelfth argument', function test( t ) { var values; var i; @@ -190,7 +190,7 @@ tape( 'the function throws an error if provided an invalid eleventh argument', f function badValue( value ) { return function badValue() { - dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, new Float64Array( rxpyp.AP ), new Float64Array( rxpyp.x ), rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, new Float64Array( rxpyp.y ), value, rxpyp.offsetY ); + dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, new Float64Array( rxpyp.AP ), rxpyp.offsetAP, new Float64Array( rxpyp.x ), rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, new Float64Array( rxpyp.y ), value, rxpyp.offsetY ); }; } }); @@ -208,7 +208,7 @@ tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y expected = new Float64Array( rxpyp.y_out ); - out = dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, ap, x, rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, y, rxpyp.strideY, rxpyp.offsetY ); + out = dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, ap, rxpyp.offsetAP, x, rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, y, rxpyp.strideY, rxpyp.offsetY ); t.strictEqual( out, y, 'returns expected value' ); isApprox( t, y, expected, 2.0 ); @@ -228,7 +228,7 @@ tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y expected = new Float64Array( cxpyp.y_out ); - out = dspmv( cxpyp.order, cxpyp.uplo, cxpyp.N, cxpyp.alpha, ap, x, cxpyp.strideX, cxpyp.offsetX, cxpyp.beta, y, cxpyp.strideY, cxpyp.offsetY ); + out = dspmv( cxpyp.order, cxpyp.uplo, cxpyp.N, cxpyp.alpha, ap, cxpyp.offsetAP, x, cxpyp.strideX, cxpyp.offsetX, cxpyp.beta, y, cxpyp.strideY, cxpyp.offsetY ); t.strictEqual( out, y, 'returns expected value' ); isApprox( t, y, expected, 2.0 ); @@ -248,7 +248,7 @@ tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y expected = new Float64Array( rxoyt.y_out ); - out = dspmv( rxoyt.order, rxoyt.uplo, rxoyt.N, rxoyt.alpha, ap, x, rxoyt.strideX, rxoyt.offsetX, rxoyt.beta, y, rxoyt.strideY, rxoyt.offsetY ); + out = dspmv( rxoyt.order, rxoyt.uplo, rxoyt.N, rxoyt.alpha, ap, rxoyt.offsetAP, x, rxoyt.strideX, rxoyt.offsetX, rxoyt.beta, y, rxoyt.strideY, rxoyt.offsetY ); t.strictEqual( out, y, 'returns expected value' ); isApprox( t, y, expected, 2.0 ); @@ -268,7 +268,7 @@ tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y expected = new Float64Array( cxoyt.y_out ); - out = dspmv( cxoyt.order, cxoyt.uplo, cxoyt.N, cxoyt.alpha, ap, x, cxoyt.strideX, cxoyt.offsetX, cxoyt.beta, y, cxoyt.strideY, cxoyt.offsetY ); + out = dspmv( cxoyt.order, cxoyt.uplo, cxoyt.N, cxoyt.alpha, ap, cxoyt.offsetAP, x, cxoyt.strideX, cxoyt.offsetX, cxoyt.beta, y, cxoyt.strideY, cxoyt.offsetY ); t.strictEqual( out, y, 'returns expected value' ); isApprox( t, y, expected, 2.0 ); @@ -288,7 +288,7 @@ tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y expected = new Float64Array( rxpyn.y_out ); - out = dspmv( rxpyn.order, rxpyn.uplo, rxpyn.N, rxpyn.alpha, ap, x, rxpyn.strideX, rxpyn.offsetX, rxpyn.beta, y, rxpyn.strideY, rxpyn.offsetY ); + out = dspmv( rxpyn.order, rxpyn.uplo, rxpyn.N, rxpyn.alpha, ap, rxpyn.offsetAP, x, rxpyn.strideX, rxpyn.offsetX, rxpyn.beta, y, rxpyn.strideY, rxpyn.offsetY ); t.strictEqual( out, y, 'returns expected value' ); isApprox( t, y, expected, 2.0 ); @@ -308,7 +308,7 @@ tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y expected = new Float64Array( cxpyn.y_out ); - out = dspmv( cxpyn.order, cxpyn.uplo, cxpyn.N, cxpyn.alpha, ap, x, cxpyn.strideX, cxpyn.offsetX, cxpyn.beta, y, cxpyn.strideY, cxpyn.offsetY ); + out = dspmv( cxpyn.order, cxpyn.uplo, cxpyn.N, cxpyn.alpha, ap, cxpyn.offsetAP, x, cxpyn.strideX, cxpyn.offsetX, cxpyn.beta, y, cxpyn.strideY, cxpyn.offsetY ); t.strictEqual( out, y, 'returns expected value' ); isApprox( t, y, expected, 2.0 ); @@ -328,7 +328,7 @@ tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y expected = new Float64Array( rxnyp.y_out ); - out = dspmv( rxnyp.order, rxnyp.uplo, rxnyp.N, rxnyp.alpha, ap, x, rxnyp.strideX, rxnyp.offsetX, rxnyp.beta, y, rxnyp.strideY, rxnyp.offsetY ); + out = dspmv( rxnyp.order, rxnyp.uplo, rxnyp.N, rxnyp.alpha, ap, rxnyp.offsetAP, x, rxnyp.strideX, rxnyp.offsetX, rxnyp.beta, y, rxnyp.strideY, rxnyp.offsetY ); t.strictEqual( out, y, 'returns expected value' ); isApprox( t, y, expected, 2.0 ); @@ -348,7 +348,7 @@ tape( 'the function performs the matrix-vector operation `y = alpha*A*x + beta*y expected = new Float64Array( cxnyp.y_out ); - out = dspmv( cxnyp.order, cxnyp.uplo, cxnyp.N, cxnyp.alpha, ap, x, cxnyp.strideX, cxnyp.offsetX, cxnyp.beta, y, cxnyp.strideY, cxnyp.offsetY ); + out = dspmv( cxnyp.order, cxnyp.uplo, cxnyp.N, cxnyp.alpha, ap, cxnyp.offsetAP, x, cxnyp.strideX, cxnyp.offsetX, cxnyp.beta, y, cxnyp.strideY, cxnyp.offsetY ); t.strictEqual( out, y, 'returns expected value' ); isApprox( t, y, expected, 2.0 ); @@ -365,7 +365,7 @@ tape( 'the function returns a reference to the second input vector', function te x = new Float64Array( rxpyp.x ); y = new Float64Array( rxpyp.y ); - out = dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, ap, x, rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, y, rxpyp.strideY, rxpyp.offsetY ); + out = dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, rxpyp.alpha, ap, rxpyp.offsetAP, x, rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, y, rxpyp.strideY, rxpyp.offsetY ); t.strictEqual( out, y, 'returns expected value' ); t.end(); @@ -384,11 +384,11 @@ tape( 'if `N` is zero or the scalar constants are zero and unity, respectively, expected = new Float64Array( rxpyp.y ); - out = dspmv( rxpyp.order, rxpyp.uplo, 0, rxpyp.alpha, ap, x, rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, y, rxpyp.strideY, rxpyp.offsetY ); + out = dspmv( rxpyp.order, rxpyp.uplo, 0, rxpyp.alpha, ap, rxpyp.offsetAP, x, rxpyp.strideX, rxpyp.offsetX, rxpyp.beta, y, rxpyp.strideY, rxpyp.offsetY ); t.strictEqual( out, y, 'returns expected value' ); t.deepEqual( y, expected, 'returns expected value' ); - out = dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, 0.0, ap, x, rxpyp.strideX, rxpyp.offsetX, 1.0, y, rxpyp.strideY, rxpyp.offsetY ); + out = dspmv( rxpyp.order, rxpyp.uplo, rxpyp.N, 0.0, ap, rxpyp.offsetAP, x, rxpyp.strideX, rxpyp.offsetX, 1.0, y, rxpyp.strideY, rxpyp.offsetY ); t.strictEqual( out, y, 'returns expected value' ); t.deepEqual( y, expected, 'returns expected value' ); @@ -408,11 +408,11 @@ tape( 'if `N` is zero or the scalar constants are zero and unity, respectively, expected = new Float64Array( cxpyp.y ); - out = dspmv( cxpyp.order, cxpyp.uplo, 0, cxpyp.alpha, ap, x, cxpyp.strideX, cxpyp.offsetX, cxpyp.beta, y, cxpyp.strideY, cxpyp.offsetY ); + out = dspmv( cxpyp.order, cxpyp.uplo, 0, cxpyp.alpha, ap, cxpyp.offsetAP, x, cxpyp.strideX, cxpyp.offsetX, cxpyp.beta, y, cxpyp.strideY, cxpyp.offsetY ); t.strictEqual( out, y, 'returns expected value' ); t.deepEqual( y, expected, 'returns expected value' ); - out = dspmv( cxpyp.order, cxpyp.uplo, cxpyp.N, 0.0, ap, x, cxpyp.strideX, cxpyp.offsetX, 1.0, y, cxpyp.strideY, cxpyp.offsetY ); + out = dspmv( cxpyp.order, cxpyp.uplo, cxpyp.N, 0.0, ap, cxpyp.offsetAP, x, cxpyp.strideX, cxpyp.offsetX, 1.0, y, cxpyp.strideY, cxpyp.offsetY ); t.strictEqual( out, y, 'returns expected value' ); t.deepEqual( y, expected, 'returns expected value' ); @@ -432,13 +432,13 @@ tape( 'when `α` is zero, the function scales the second input vector (row-major expected = new Float64Array( [ 5.0, 5.0, 5.0 ] ); - out = dspmv( 'row-major', 'upper', 3, 0.0, ap, x, 1, 0, 5.0, y, 1, 0 ); + out = dspmv( 'row-major', 'upper', 3, 0.0, ap, 0, x, 1, 0, 5.0, y, 1, 0 ); t.strictEqual( out, y, 'returns expected value' ); t.deepEqual( y, expected, 'returns expected value' ); expected = new Float64Array( [ 0.0, 0.0, 0.0 ] ); - out = dspmv( 'row-major', 'upper', 3, 0.0, ap, x, 1, 0, 0.0, y, -1, 2 ); + out = dspmv( 'row-major', 'upper', 3, 0.0, ap, 0, x, 1, 0, 0.0, y, -1, 2 ); t.strictEqual( out, y, 'returns expected value' ); t.deepEqual( y, expected, 'returns expected value' ); @@ -458,13 +458,13 @@ tape( 'when `α` is zero, the function scales the second input vector (row-major expected = new Float64Array( [ 5.0, 5.0, 5.0 ] ); - out = dspmv( 'row-major', 'lower', 3, 0.0, ap, x, 1, 0, 5.0, y, -1, 2 ); + out = dspmv( 'row-major', 'lower', 3, 0.0, ap, 0, x, 1, 0, 5.0, y, -1, 2 ); t.strictEqual( out, y, 'returns expected value' ); t.deepEqual( y, expected, 'returns expected value' ); expected = new Float64Array( [ 0.0, 0.0, 0.0 ] ); - out = dspmv( 'row-major', 'lower', 3, 0.0, ap, x, 1, 0, 0.0, y, 1, 0 ); + out = dspmv( 'row-major', 'lower', 3, 0.0, ap, 0, x, 1, 0, 0.0, y, 1, 0 ); t.strictEqual( out, y, 'returns expected value' ); t.deepEqual( y, expected, 'returns expected value' ); @@ -484,13 +484,13 @@ tape( 'when `α` is zero, the function scales the second input vector (column-ma expected = new Float64Array( [ 5.0, 5.0, 5.0 ] ); - out = dspmv( 'column-major', 'upper', 3, 0.0, ap, x, 1, 0, 5.0, y, 1, 0 ); + out = dspmv( 'column-major', 'upper', 3, 0.0, ap, 0, x, 1, 0, 5.0, y, 1, 0 ); t.strictEqual( out, y, 'returns expected value' ); t.deepEqual( y, expected, 'returns expected value' ); expected = new Float64Array( [ 0.0, 0.0, 0.0 ] ); - out = dspmv( 'column-major', 'upper', 3, 0.0, ap, x, 1, 0, 0.0, y, -1, 2 ); + out = dspmv( 'column-major', 'upper', 3, 0.0, ap, 0, x, 1, 0, 0.0, y, -1, 2 ); t.strictEqual( out, y, 'returns expected value' ); t.deepEqual( y, expected, 'returns expected value' ); @@ -510,13 +510,13 @@ tape( 'when `α` is zero, the function scales the second input vector (column-ma expected = new Float64Array( [ 5.0, 5.0, 5.0 ] ); - out = dspmv( 'column-major', 'lower', 3, 0.0, ap, x, 1, 0, 5.0, y, -1, 2 ); + out = dspmv( 'column-major', 'lower', 3, 0.0, ap, 0, x, 1, 0, 5.0, y, -1, 2 ); t.strictEqual( out, y, 'returns expected value' ); t.deepEqual( y, expected, 'returns expected value' ); expected = new Float64Array( [ 0.0, 0.0, 0.0 ] ); - out = dspmv( 'column-major', 'lower', 3, 0.0, ap, x, 1, 0, 0.0, y, 1, 0 ); + out = dspmv( 'column-major', 'lower', 3, 0.0, ap, 0, x, 1, 0, 0.0, y, 1, 0 ); t.strictEqual( out, y, 'returns expected value' ); t.deepEqual( y, expected, 'returns expected value' ); @@ -536,7 +536,7 @@ tape( 'the function supports complex access patterns (row-major)', function test expected = new Float64Array( rxnyn.y_out ); - out = dspmv( rxnyn.order, rxnyn.uplo, rxnyn.N, rxnyn.alpha, ap, x, rxnyn.strideX, rxnyn.offsetX, rxnyn.beta, y, rxnyn.strideY, rxnyn.offsetY ); + out = dspmv( rxnyn.order, rxnyn.uplo, rxnyn.N, rxnyn.alpha, ap, rxnyn.offsetAP, x, rxnyn.strideX, rxnyn.offsetX, rxnyn.beta, y, rxnyn.strideY, rxnyn.offsetY ); t.strictEqual( out, y, 'returns expected value' ); isApprox( t, y, expected, 2.0 ); @@ -556,7 +556,7 @@ tape( 'the function supports complex access patterns (column-major)', function t expected = new Float64Array( cxnyn.y_out ); - out = dspmv( cxnyn.order, cxnyn.uplo, cxnyn.N, cxnyn.alpha, ap, x, cxnyn.strideX, cxnyn.offsetX, cxnyn.beta, y, cxnyn.strideY, cxnyn.offsetY ); + out = dspmv( cxnyn.order, cxnyn.uplo, cxnyn.N, cxnyn.alpha, ap, cxnyn.offsetAP, x, cxnyn.strideX, cxnyn.offsetX, cxnyn.beta, y, cxnyn.strideY, cxnyn.offsetY ); t.strictEqual( out, y, 'returns expected value' ); isApprox( t, y, expected, 2.0 );