Merge branch 'develop' into nd/rot90

Author: headlessNode

docs/migration-guides/numpy/README.md

@@ -0,0 +1,74 @@
+<!--
+
+@license Apache-2.0
+
+Copyright (c) 2026 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.
+
+-->
+
+<!-- lint disable expected-html-sections -->
+
+# NumPy
+
+> Migration guide for NumPy.
+
+## Common Operations
+
+<!-- lint disable table-pipe-alignment -->
+
+| Description | NumPy | stdlib |
+| ----------- | ----- | ------ |
+| Test whether an array contains truthy values | `np.any(x)` | [`any(x)`][@stdlib/ndarray/any] |
+| Broadcast an array to a specified shape | `np.broadcast_to(x, shape)` | [`broadcastArray(x, shape)`][@stdlib/ndarray/broadcast-array] |
+| Broadcast a scalar to a specified shape | `np.broadcast_to(np.array(scalar), shape)` | [`broadcastScalar(scalar, shape)`][@stdlib/ndarray/broadcast-scalar] |
+| Copy an array | `np.copy(x)` | [`copy(x)`][@stdlib/ndarray/copy] |
+| Count the number of falsy values in an array | `x.size-np.count_nonzero(x)` | [`countFalsy(x)`][@stdlib/ndarray/count-falsy] |
+| Count the number of truthy values in an array | `np.count_nonzeros(x)` | [`countTruthy(x)`][@stdlib/ndarray/count-truthy] |
+| Test whether an array includes a specific value | `np.any(np.equal(x,v))` | [`includes(x,v)`][@stdlib/ndarray/includes] |
+| Reverse the elements along a dimension | `np.flip(x, axis=dim)` | [`reverseDimension(x, dim)`][@stdlib/ndarray/reverse-dimension] |
+| Prepend a a specified number of singleton dimensions | `np.reshape(x, (1,)*n + x.shape)` | [`prependSingletonDimensions(x, n)`][@stdlib/ndarray/prepend-singleton-dimensions] |
+| Test whether an array contains at least `n` truthy values | `np.count_nonzero(x) >= n` | [`some(x, n)`][@stdlib/ndarray/some] |
+
+<!-- lint enable table-pipe-alignment -->
+
+## Function Reference
+
+{{table}}
+
+<section class="links">
+
+[@stdlib/ndarray/any]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/ndarray/any
+
+[@stdlib/ndarray/broadcast-array]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/ndarray/broadcast-array
+
+[@stdlib/ndarray/broadcast-scalar]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/ndarray/broadcast-scalar
+
+[@stdlib/ndarray/copy]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/ndarray/copy
+
+[@stdlib/ndarray/count-falsy]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/ndarray/count-falsy
+
+[@stdlib/ndarray/count-truthy]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/ndarray/count-truthy
+
+[@stdlib/ndarray/includes]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/ndarray/includes
+
+[@stdlib/ndarray/prepend-singleton-dimensions]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/ndarray/prepend-singleton-dimensions
+
+[@stdlib/ndarray/reverse-dimension]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/ndarray/reverse-dimension
+
+[@stdlib/ndarray/some]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/ndarray/some
+
+</section>
+
+<!-- /.links -->

lib/node_modules/@stdlib/blas/base/cgemv/lib/base.js

@@ -143,7 +143,7 @@ function cgemv( trans, M, N, alpha, A, strideA1, strideA2, offsetA, x, strideX,
 
 	// y = beta*y
 	if ( rebeta === 0.0 && imbeta === 0.0 ) {
-		cfill( ylen, alpha, y, strideY, offsetY );
+		cfill( ylen, beta, y, strideY, offsetY );
 	} else if ( rebeta !== 1.0 || imbeta !== 0.0 ) {
 		cscal( ylen, beta, y, strideY, offsetY );
 	}
@@ -160,8 +160,8 @@ function cgemv( trans, M, N, alpha, A, strideA1, strideA2, offsetA, x, strideX,
 		sign = 1;
 	}
 	oa = offsetA * 2;
-	sa1 = strideA1 *2;
-	sa2 = strideA2 *2;
+	sa1 = strideA1 * 2;
+	sa2 = strideA2 * 2;
 	ox = offsetX * 2;
 	oy = offsetY * 2;
 	sx = strideX * 2;

lib/node_modules/@stdlib/blas/base/cgemv/lib/ndarray.js

@@ -24,6 +24,9 @@ var isMatrixTranspose = require( '@stdlib/blas/base/assert/is-transpose-operatio
 var format = require( '@stdlib/string/format' );
 var base = require( './base.js' );
 
+
+// MAIN //
+
 /**
 * Performs one of the matrix-vector operations `y = α*A*x + β*y` or `y = α*A^T*x + β*y` or `y = α*A^H*x + β*y`, where `α` and `β` are scalars, `x` and `y` are vectors, and `A` is an `M` by `N` matrix.
 *

lib/node_modules/@stdlib/blas/base/docs/types/index.d.ts

@@ -23,6 +23,7 @@
 import assert = require( '@stdlib/blas/base/assert' );
 import caxpy = require( '@stdlib/blas/base/caxpy' );
 import ccopy = require( '@stdlib/blas/base/ccopy' );
+import cgemv = require( '@stdlib/blas/base/cgemv' );
 import cscal = require( '@stdlib/blas/base/cscal' );
 import csrot = require( '@stdlib/blas/base/csrot' );
 import csscal = require( '@stdlib/blas/base/csscal' );
@@ -200,6 +201,51 @@ interface Namespace {
 	*/
 	ccopy: typeof ccopy;
 
+	/**
+	* Performs one of the matrix-vector operations `y = α*A*x + β*y` or `y = α*A^T*x + β*y` or `y = α*A^H*x + β*y`, where `α` and `β` are scalars, `x` and `y` are vectors, and `A` is an `M` by `N` matrix.
+	*
+	* @param order - storage layout
+	* @param trans - specifies whether `A` should be transposed, conjugate-transposed, or not transposed
+	* @param M - number of rows in the matrix `A`
+	* @param N - number of columns in the matrix `A`
+	* @param alpha - scalar constant
+	* @param A - input matrix
+	* @param LDA - stride of the first dimension of `A` (a.k.a., leading dimension of the matrix `A`)
+	* @param x - first input vector
+	* @param strideX - `x` stride length
+	* @param beta - scalar constant
+	* @param y - second input vector
+	* @param strideY - `y` stride length
+	* @returns `y`
+	*
+	* @example
+	* var Complex64Array = require( '@stdlib/array/complex64' );
+	* var Complex64 = require( '@stdlib/complex/float32/ctor' );
+	*
+	* var A = new Complex64Array( [ 1.0, 1.0, 2.0, 2.0, 3.0, 3.0, 4.0, 4.0, 5.0, 5.0, 6.0, 6.0, 7.0, 7.0, 8.0, 8.0 ] );
+	* var x = new Complex64Array( [ 1.0, 1.0, 2.0, 2.0 ] );
+	* var y = new Complex64Array( [ 1.0, 1.0, 2.0, 2.0, 3.0, 3.0, 4.0, 4.0 ] );
+	* var alpha = new Complex64( 0.5, 0.5 );
+	* var beta = new Complex64( 0.5, -0.5 );
+	*
+	* ns.cgemv( 'column-major', 'no-transpose', 4, 2, alpha, A, 4, x, 1, beta, y, 1 );
+	* // y => <Complex64Array>[ -10.0, 11.0, -12.0, 14.0, -14.0, 17.0, -16.0, 20.0 ]
+	*
+	* @example
+	* var Complex64Array = require( '@stdlib/array/complex64' );
+	* var Complex64 = require( '@stdlib/complex/float32/ctor' );
+	*
+	* var A = new Complex64Array( [ 1.0, 1.0, 2.0, 2.0, 3.0, 3.0, 4.0, 4.0, 5.0, 5.0, 6.0, 6.0, 7.0, 7.0, 8.0, 8.0 ] );
+	* var x = new Complex64Array( [ 1.0, 1.0, 2.0, 2.0 ] );
+	* var y = new Complex64Array( [ 1.0, 1.0, 2.0, 2.0, 3.0, 3.0, 4.0, 4.0 ] );
+	* var alpha = new Complex64( 0.5, 0.5 );
+	* var beta = new Complex64( 0.5, -0.5 );
+	*
+	* ns.cgemv.ndarray( 'no-transpose', 4, 2, alpha, A, 1, 4, 0, x, 1, 0, beta, y, 1, 0 );
+	* // y => <Complex64Array>[ -10.0, 11.0, -12.0, 14.0, -14.0, 17.0, -16.0, 20.0 ]
+	*/
+	cgemv: typeof cgemv;
+
 	/**
 	* Scales a single-precision complex floating-point vector by a single-precision complex floating-point constant.
 	*

lib/node_modules/@stdlib/blas/ext/base/README.md

@@ -62,6 +62,7 @@ var o = ns;
 -   <span class="signature">[`dcusumkbn2( N, sum, x, strideX, y, strideY )`][@stdlib/blas/ext/base/dcusumkbn2]</span><span class="delimiter">: </span><span class="description">calculate the cumulative sum of double-precision floating-point strided array elements using a second-order iterative Kahan–Babuška algorithm.</span>
 -   <span class="signature">[`dcusumors( N, sum, x, strideX, y, strideY )`][@stdlib/blas/ext/base/dcusumors]</span><span class="delimiter">: </span><span class="description">calculate the cumulative sum of double-precision floating-point strided array elements using ordinary recursive summation.</span>
 -   <span class="signature">[`dcusumpw( N, sum, x, strideX, y, strideY )`][@stdlib/blas/ext/base/dcusumpw]</span><span class="delimiter">: </span><span class="description">calculate the cumulative sum of double-precision floating-point strided array elements using pairwise summation.</span>
+-   <span class="signature">[`ddiff( N, k, x, strideX, N1, prepend, strideP, N2, append, strideA, out, strideOut, workspace, strideW )`][@stdlib/blas/ext/base/ddiff]</span><span class="delimiter">: </span><span class="description">calculate the k-th discrete forward difference of a double-precision floating-point strided array.</span>
 -   <span class="signature">[`dfill( N, alpha, x, strideX )`][@stdlib/blas/ext/base/dfill]</span><span class="delimiter">: </span><span class="description">fill a double-precision floating-point strided array with a specified scalar constant.</span>
 -   <span class="signature">[`dindexOfRow( order, M, N, A, LDA, x, strideX, workspace, strideW )`][@stdlib/blas/ext/base/dindex-of-row]</span><span class="delimiter">: </span><span class="description">return the index of the first row in a double-precision floating-point input matrix which has the same elements as a provided search vector.</span>
 -   <span class="signature">[`dindexOf( N, searchElement, x, strideX )`][@stdlib/blas/ext/base/dindex-of]</span><span class="delimiter">: </span><span class="description">return the first index of a specified search element in a double-precision floating-point strided array.</span>
@@ -307,6 +308,8 @@ console.log( objectKeys( ns ) );
 
 [@stdlib/blas/ext/base/dcusumpw]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/dcusumpw
 
+[@stdlib/blas/ext/base/ddiff]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/ddiff
+
 [@stdlib/blas/ext/base/dfill]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/dfill
 
 [@stdlib/blas/ext/base/dindex-of-row]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/dindex-of-row