chore: initial implementation

---
type: pre_push_report
description: Results of running various checks prior to pushing changes.
report:
  - task: run_javascript_examples
    status: passed
  - task: run_c_examples
    status: na
  - task: run_cpp_examples
    status: na
  - task: run_javascript_readme_examples
    status: na
  - task: run_c_benchmarks
    status: na
  - task: run_cpp_benchmarks
    status: na
  - task: run_fortran_benchmarks
    status: na
  - task: run_javascript_benchmarks
    status: na
  - task: run_julia_benchmarks
    status: na
  - task: run_python_benchmarks
    status: na
  - task: run_r_benchmarks
    status: na
  - task: run_javascript_tests
    status: na
---

Author: ShabiShett07

lib/node_modules/@stdlib/blas/base/sgbmv/examples/index.js

@@ -27,15 +27,14 @@ var opts = {
 
 const A = [0,  9, 10,  // First row (upper band)
     11, 12,  5,  // Second row (main diagonal)
-     6,  7,  8,  // Third row (lower band)
-     0,  1,  2,  // Fourth row (lower band)
-     0,  0,  3]; // Packed storage (row-major)
+    1,  2, 0  // Fourth row (lower band)
+]; // Packed storage (row-major)
 const x = [57, 245, 121];
 const y = [0, 0, 0];
 console.log( x );
 console.log( y );
 
-const order = 'column-major';
+const order = 'row-major';
 const trans = 'no-transpose';
 const LDA = 3;
 const M = 3, N = 3;

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

@@ -112,40 +112,34 @@ function sgbmv(trans, M, N, KL, KU, alpha, A, strideA1, strideA2, offsetA, x, st
         kup1 = KU + 1;
         ix1 = offsetX;
         for (i1 = 0; i1 < xlen; i1++) {
-            tmp = f32(alpha * x[ix1]);
+            tmp = f32(alpha * x[ix1]);  // Compute scaled x[i1]
             oa = offsetA + (sa1 * i1);
             iy1 = offsetY;
+
+            console.log(`Processing x[${i1}] = ${x[ix1]}, scaled tmp = ${tmp}`);
+            console.log(`Offset in A (oa) = ${oa}`);
+
             for (i0 = Math.max(0, i1 - KU); i0 < Math.min(ylen, i1 + KL + 1); i0++) {
-                // Calculate diagonal offset
-                let diag_offset = i0 - i1;
-
-                // Fix the a_idx calculation for banded matrix format
-                // The banded matrix layout appears to have a different structure than expected
-                // The correct index depends on the exact layout of your band matrix
-
-                // Based on the debug output pattern, this appears to be the correct formula:
-                if (diag_offset == -1) {
-                    // Lower diagonal elements
-                    a_idx = 3 + (3 * i1);
-                } else if (diag_offset == 0) {
-                    // Main diagonal elements (9, 12, 8)
-                    if (i1 == 0) a_idx = 1;
-                    else if (i1 == 1) a_idx = 4;
-                    else if (i1 == 2) a_idx = 8; // This was wrong in original - should be 8 not 7
-                } else if (diag_offset == 1) {
-                    // Upper diagonal elements
-                    a_idx = 2 + (3 * i1);
-                }
+                // Fixed indexing for banded matrix
+                a_idx = oa + (kup1 + KU) * sa0;
+                console.log(`i0 = ${i0}, a_idx = ${a_idx}, A[a_idx] = ${A[a_idx]}`);
 
-                if (a_idx >= 0 && a_idx < A.length) {
-                    y[iy1] += f32(A[a_idx] * tmp);
+                if (i0 < ylen && a_idx >= 0) {
+                    y[iy1] += f32(A[a_idx] * tmp);  // Accumulate correctly
+                    console.log(`Updated y[${i0}] = ${y[iy1]}`);
                 }
                 iy1 += strideY;
             }
+
+            // Debugging: Print intermediate y after processing x[i1]
+            console.log(`Intermediate y after processing x[${i1}]: [${y}]`);
             ix1 += strideX;
         }
+        console.log("Expected: [3208, 4236, 2072]");
         return y;
     }
+
+    // Leaving the second part unchanged
     kup1 = KU + 1;
     iy1 = offsetY;
     for (i1 = 0; i1 < ylen; i1++) {

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

@@ -90,7 +90,7 @@ function sgemv( trans, M, N, KL, KU, alpha, A, strideA1, strideA2, offsetA, x, s
 	if ( M === 0 || N === 0 || ( alpha === 0.0 && beta === 1.0 ) ) {
 		return y;
 	}
-	return base( trans, M, N, alpha, A, strideA1, strideA2, offsetA, x, strideX, offsetX, beta, y, strideY, offsetY ); // eslint-disable-line max-len
+	return base( trans, M, N, KL, KU, alpha, A, strideA1, strideA2, offsetA, x, strideX, offsetX, beta, y, strideY, offsetY ); // eslint-disable-line max-len
 }
 
 

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

@@ -36,6 +36,8 @@ var base = require( './base.js' );
 * @param {string} trans - specifies whether `A` should be transposed, conjugate-transposed, or not transposed
 * @param {NonNegativeInteger} M - number of rows in the matrix `A`
 * @param {NonNegativeInteger} N - number of columns in the matrix `A`
+* @param {NonNegativeInteger} KL - number of sub-diagonals of matrix `A`
+* @param {NonNegativeInteger} KU - number of super-diagonals of matrix `A`
 * @param {number} alpha - scalar constant
 * @param {Float32Array} A - input matrix
 * @param {PositiveInteger} LDA - stride of the first dimension of `A` (a.k.a., leading dimension of the matrix `A`)
@@ -118,7 +120,7 @@ function sgbmv( order, trans, M, N, KL, KU, alpha, A, LDA, x, strideX, beta, y,
 		sa1 = LDA;
 		sa2 = 1;
 	}
-	return base( trans, M, N, alpha, A, sa1, sa2, 0, x, strideX, ox, beta, y, strideY, oy ); // eslint-disable-line max-len
+	return base( trans, M, N, KL, KU, alpha, A, sa1, sa2, 0, x, strideX, ox, beta, y, strideY, oy ); // eslint-disable-line max-len
 }