chore: minor clean-up

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

Author: ShabiShett07

lib/node_modules/@stdlib/blas/base/dgemv/include/stdlib/blas/base/dgemv_cblas.h

@@ -16,28 +16,74 @@
 * limitations under the License.
 */
 
-/**
-* Header file containing function declarations for the C interface to the CBLAS Level 2 routine `cblas_dgemv`.
-*/
-#ifndef DGEMV_CBLAS_H
-#define DGEMV_CBLAS_H
-
+#include "stdlib/blas/base/dgemv.h"
+#include "stdlib/blas/base/dgemv_cblas.h"
 #include "stdlib/blas/base/shared.h"
+#include "stdlib/strided/base/min_view_buffer_index.h"
+#include "stdlib/ndarray/base/min_view_buffer_index.h"
 
-/*
-* If C++, prevent name mangling so that the compiler emits a binary file having undecorated names, thus mirroring the behavior of a C compiler.
+/**
+* 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
+* @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
+* @return         output value
 */
-#ifdef __cplusplus
-extern "C" {
-#endif
+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`, 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`, 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
 */
-void API_SUFFIX(cblas_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, double *Y, const CBLAS_INT strideY );
-
-#ifdef __cplusplus
+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 );
 }
-#endif
-
-#endif // !DGEMV_CBLAS_H