refactor: simplify implementation

---
type: pre_commit_static_analysis_report
description: Results of running static analysis checks when committing changes.
report:
  - task: lint_filenames
    status: passed
  - task: lint_editorconfig
    status: passed
  - task: lint_markdown
    status: na
  - task: lint_package_json
    status: na
  - task: lint_repl_help
    status: na
  - task: lint_javascript_src
    status: passed
  - task: lint_javascript_cli
    status: na
  - task: lint_javascript_examples
    status: na
  - task: lint_javascript_tests
    status: na
  - task: lint_javascript_benchmarks
    status: na
  - task: lint_python
    status: na
  - task: lint_r
    status: na
  - task: lint_c_src
    status: passed
  - task: lint_c_examples
    status: na
  - task: lint_c_benchmarks
    status: na
  - task: lint_c_tests_fixtures
    status: na
  - task: lint_shell
    status: na
  - task: lint_typescript_declarations
    status: na
  - task: lint_typescript_tests
    status: na
  - task: lint_license_headers
    status: passed
---

Author: kgryte

lib/node_modules/@stdlib/blas/ext/base/sapxsumkbn2/lib/ndarray.js

@@ -20,8 +20,8 @@
 
 // MODULES //
 
-var float64ToFloat32 = require( '@stdlib/number/float64/base/to-float32' );
-var abs = require( '@stdlib/math/base/special/abs' );
+var f32 = require( '@stdlib/number/float64/base/to-float32' );
+var ssumkbn2 = require( '@stdlib/blas/ext/base/ssumkbn2' ).ndarray;
 
 
 // MAIN //
@@ -53,46 +53,7 @@ var abs = require( '@stdlib/math/base/special/abs' );
 * // returns 25.0
 */
 function sapxsumkbn2( N, alpha, x, strideX, offsetX ) {
-	var sum;
-	var ccs;
-	var ix;
-	var cs;
-	var cc;
-	var v;
-	var t;
-	var c;
-	var i;
-
-	if ( N <= 0 ) {
-		return 0.0;
-	}
-	ix = offsetX;
-	if ( strideX === 0 ) {
-		return float64ToFloat32( N * float64ToFloat32( alpha + x[ ix ] ) );
-	}
-	sum = 0.0;
-	ccs = 0.0; // second order correction term for lost low order bits
-	cs = 0.0; // first order correction term for lost low order bits
-	for ( i = 0; i < N; i++ ) {
-		v = float64ToFloat32( alpha + x[ ix ] );
-		t = float64ToFloat32( sum+v );
-		if ( abs( sum ) >= abs( v ) ) {
-			c = float64ToFloat32( float64ToFloat32( sum-t ) + v );
-		} else {
-			c = float64ToFloat32( float64ToFloat32( v-t ) + sum );
-		}
-		sum = t;
-		t = float64ToFloat32( cs+c );
-		if ( abs( cs ) >= abs( c ) ) {
-			cc = float64ToFloat32( float64ToFloat32( cs-t ) + c );
-		} else {
-			cc = float64ToFloat32( float64ToFloat32( c-t ) + cs );
-		}
-		cs = t;
-		ccs = float64ToFloat32( ccs+cc );
-		ix += strideX;
-	}
-	return float64ToFloat32( sum + float64ToFloat32( cs+ccs ) );
+	return f32( f32( N * alpha ) + ssumkbn2( N, x, strideX, offsetX ) );
 }
 
 

lib/node_modules/@stdlib/blas/ext/base/sapxsumkbn2/manifest.json

@@ -44,7 +44,7 @@
         "@stdlib/napi/create-double",
         "@stdlib/strided/base/stride2offset",
         "@stdlib/blas/base/shared",
-        "@stdlib/math/base/special/absf"
+        "@stdlib/blas/ext/base/ssumkbn2"
       ]
     },
     {
@@ -60,7 +60,7 @@
       "dependencies": [
         "@stdlib/strided/base/stride2offset",
         "@stdlib/blas/base/shared",
-        "@stdlib/math/base/special/absf"
+        "@stdlib/blas/ext/base/ssumkbn2"
       ]
     },
     {
@@ -76,7 +76,7 @@
       "dependencies": [
         "@stdlib/strided/base/stride2offset",
         "@stdlib/blas/base/shared",
-        "@stdlib/math/base/special/absf"
+        "@stdlib/blas/ext/base/ssumkbn2"
       ]
     }
   ]

lib/node_modules/@stdlib/blas/ext/base/sapxsumkbn2/src/main.c

@@ -19,7 +19,7 @@
 #include "stdlib/blas/ext/base/sapxsumkbn2.h"
 #include "stdlib/strided/base/stride2offset.h"
 #include "stdlib/blas/base/shared.h"
-#include "stdlib/math/base/special/absf.h"
+#include "stdlib/blas/ext/base/ssumkbn2.h"
 
 /**
 * Adds a scalar constant to each single-precision floating-point strided array element and computes the sum using a second-order iterative Kahan–Babuška algorithm.
@@ -62,44 +62,5 @@ float API_SUFFIX(stdlib_strided_sapxsumkbn2)( const CBLAS_INT N, const float alp
 * @return         output value
 */
 float API_SUFFIX(stdlib_strided_sapxsumkbn2_ndarray)( const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) {
-	CBLAS_INT ix;
-	CBLAS_INT i;
-	float sum;
-	float ccs;
-	float cs;
-	float cc;
-	float v;
-	float t;
-	float c;
-
-	if ( N <= 0 ) {
-		return 0.0f;
-	}
-	ix = offsetX;
-	if ( strideX == 0 ) {
-		return N * ( alpha + X[ ix ] );
-	}
-	sum = 0.0f;
-	ccs = 0.0f; // second order correction term for lost lower order bits
-	cs = 0.0f; // first order correction term for lost low order bits
-	for ( i = 0; i < N; i++ ) {
-		v = alpha + X[ ix ];
-		t = sum + v;
-		if ( stdlib_base_absf( sum ) >= stdlib_base_absf( v ) ) {
-			c = (sum-t) + v;
-		} else {
-			c = (v-t) + sum;
-		}
-		sum = t;
-		t = cs + c;
-		if ( stdlib_base_absf( cs ) >= stdlib_base_absf( c ) ) {
-			cc = (cs-t) + c;
-		} else {
-			cc = (c-t) + cs;
-		}
-		cs = t;
-		ccs += cc;
-		ix += strideX;
-	}
-	return sum + cs + ccs;
+	return ( N * alpha ) + API_SUFFIX(stdlib_strided_ssumkbn2_ndarray)( N, X, strideX, offsetX );
 }