chore: minor clean-up

Signed-off-by: Shabareesh Shetty <[email protected]>

Author: ShabiShett07

lib/node_modules/@stdlib/blas/base/dgemv/src/dgemv_cblas.c

@@ -23,7 +23,7 @@
 #include "stdlib/ndarray/base/min_view_buffer_index.h"
 
 /**
-* Performs one of the matrix-vector operations `Y = α*A*X + β*Y` or `Y = α*A^T*X + β*Y`, where `α` and `β` are scalars, `X` and `Y` are vectors, and `A` is an `M` by `N` matrix.
+* Performs one of the matrix-vector operations `y = α*A*x + β*y` or `y = α*A^T*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
@@ -32,13 +32,58 @@
 * @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 x        first input vector
+* @param strideX  `x` stride length
 * @param beta     scalar constant
-* @param Y        second input vector
-* @param strideY  `Y` stride length
+* @param y        second input vector
+* @param strideY  `y` stride length
 * @return         output value
 */
-double API_SUFFIX(c_dgemv)( const CBLAS_LAYOUT order, const CBLAS_TRANSPOSE trans, const CBLAS_INT M, const CBLAS_INT N, const double alpha, const double *A, const CBLAS_INT LDA, const double *X, const CBLAS_INT strideX, const double beta, const double *Y, const CBLAS_INT strideY ) {
-	return API_SUFFIX(cblas_dgemv)( order, trans, M, N, alpha, A, LDA, X, strideX, beta, Y, strideY );
+double API_SUFFIX(c_dgemv)( const CBLAS_LAYOUT order, const CBLAS_TRANSPOSE trans, const CBLAS_INT M, const CBLAS_INT N, const double alpha, const double *A, const CBLAS_INT LDA, const double *x, const CBLAS_INT strideX, const double beta, const double *y, const CBLAS_INT strideY ) {
+	CBLAS_INT sx = strideX;
+	CBLAS_INT sy = strideY;
+	if ( sx < 0 ) {
+		sx = -sx;
+	}
+	if ( sy < 0 ) {
+		sy = -sy;
+	}
+	return API_SUFFIX(cblas_dgemv)( order, trans, M, N, alpha, A, LDA, x, sx, beta, y, sy );
+}
+
+/**
+* Performs one of the matrix-vector operations `y = α*A*x + β*y` or `y = α*A^T*x + β*y`, using alternative indexing semantics and where `α` and `β` are scalars, `x` and `y` are vectors, and `A` is an `M` by `N` matrix.
+*
+* @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 strideA1  stride of the first dimension of `A`
+* @param strideA1  stride of the second dimension of `A`
+* @param offsetA   starting index for `A`
+* @param x         first input vector
+* @param strideX   `x` stride length
+* @param offsetX   starting index for `x`
+* @param beta      scalar constant
+* @param y         second input vector
+* @param strideY   `y` stride length
+* @param offsetY   starting index for `Y`
+* @return          output value
+*/
+double API_SUFFIX(c_dgemv_ndarray)( const CBLAS_TRANSPOSE trans, const CBLAS_INT M, const CBLAS_INT N, const double alpha, const double *A, const CBLAS_INT strideA1, const CBLAS_INT strideA2, const CBLAS_INT offsetA, const double *x, const CBLAS_INT strideX, const CBLAS_INT offsetX, const double beta, const double *y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) {
+	CBLAS_INT sx = strideX;
+	CBLAS_INT sy = strideY;
+	if ( sx < 0 ) {
+		sx = -sx;
+	}
+	if ( sy < 0 ) {
+		sy = -sy;
+	}
+	x += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); // adjust array pointer
+	y += stdlib_strided_min_view_buffer_index( N, strideY, offsetY ); // adjust array pointer
+	const int64_t shape[] = { M, N };
+	const int64_t strides[] = { strideA1, strideA2 };
+	A += stdlib_ndarray_min_view_buffer_index( 2, shape, strides, offsetA ); // adjust array pointer
+	return API_SUFFIX(cblas_dgemv)( order, trans, M, N, alpha, A, LDA, x, sx, beta, y, sy );
 }