fix: preserve the sign of zero

---
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: passed
  - 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/scusumkbn/lib/ndarray.js

@@ -71,8 +71,23 @@ function scusumkbn( N, sum, x, strideX, offsetX, y, strideY, offsetY ) {
 	ix = offsetX;
 	iy = offsetY;
 	s = sum;
+
+	// In order to preserve the sign of zero which can be lost during compensated summation below, find the first non-zero element...
+	if ( s === 0.0 ) {
+		for ( i = 0; i < N; i++ ) {
+			v = x[ ix ];
+			if ( v !== 0.0 ) {
+				break;
+			}
+			y[ iy ] = float64ToFloat32( s + v );
+			ix += strideX;
+			iy += strideY;
+		}
+	} else {
+		i = 0;
+	}
 	c = 0.0;
-	for ( i = 0; i < N; i++ ) {
+	for ( ; i < N; i++ ) {
 		v = x[ ix ];
 		t = float64ToFloat32( s + v );
 		if ( abs( s ) >= abs( v ) ) {

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

@@ -80,8 +80,23 @@ void API_SUFFIX(stdlib_strided_scusumkbn_ndarray)( const CBLAS_INT N, const floa
 	ix = offsetX;
 	iy = offsetY;
 	s = sum;
+
+	// In order to preserve the sign of zero which can be lost during compensated summation below, find the first non-zero element...
+	if ( s == 0.0 ) {
+		for ( i = 0; i < N; i++ ) {
+			v = X[ ix ];
+			if ( v != 0.0f ) {
+				break;
+			}
+			Y[ iy ] = s + v;
+			ix += strideX;
+			iy += strideY;
+		}
+	} else {
+		i = 0;
+	}
 	c = 0.0f;
-	for ( i = 0; i < N; i++ ) {
+	for ( ; i < N; i++ ) {
 		v = X[ ix ];
 		t = s + v;
 		if ( stdlib_base_absf( s ) >= stdlib_base_absf( v ) ) {

lib/node_modules/@stdlib/blas/ext/base/scusumkbn/test/test.ndarray.js

@@ -22,6 +22,7 @@
 
 var tape = require( 'tape' );
 var isnanf = require( '@stdlib/math/base/assert/is-nanf' );
+var isSameFloat32Array = require( '@stdlib/assert/is-same-float32array' );
 var Float32Array = require( '@stdlib/array/float32' );
 var scusumkbn = require( './../lib/ndarray.js' );
 
@@ -126,6 +127,26 @@ tape( 'the function calculates the cumulative sum', function test( t ) {
 	t.end();
 });
 
+tape( 'the function preserves the sign of zero', function test( t ) {
+	var expected;
+	var x;
+	var y;
+
+	x = new Float32Array( [ -0.0, -0.0, -0.0, 0.0, 1.0 ] );
+	y = new Float32Array( x.length );
+
+	scusumkbn( x.length, -0.0, x, 1, 0, y, 1, 0 );
+	expected = new Float32Array([
+		-0.0,
+		-0.0,
+		-0.0,
+		0.0,
+		1.0
+	]);
+	t.strictEqual( isSameFloat32Array( y, expected ), true, 'returns expected value' );
+	t.end();
+});
+
 tape( 'the function returns a reference to the output array', function test( t ) {
 	var out;
 	var x;

lib/node_modules/@stdlib/blas/ext/base/scusumkbn/test/test.ndarray.native.js

@@ -23,6 +23,7 @@
 var resolve = require( 'path' ).resolve;
 var tape = require( 'tape' );
 var isnanf = require( '@stdlib/math/base/assert/is-nanf' );
+var isSameFloat32Array = require( '@stdlib/assert/is-same-float32array' );
 var Float32Array = require( '@stdlib/array/float32' );
 var tryRequire = require( '@stdlib/utils/try-require' );
 
@@ -135,6 +136,26 @@ tape( 'the function calculates the cumulative sum', opts, function test( t ) {
 	t.end();
 });
 
+tape( 'the function preserves the sign of zero', opts, function test( t ) {
+	var expected;
+	var x;
+	var y;
+
+	x = new Float32Array( [ -0.0, -0.0, -0.0, 0.0, 1.0 ] );
+	y = new Float32Array( x.length );
+
+	scusumkbn( x.length, -0.0, x, 1, 0, y, 1, 0 );
+	expected = new Float32Array([
+		-0.0,
+		-0.0,
+		-0.0,
+		0.0,
+		1.0
+	]);
+	t.strictEqual( isSameFloat32Array( y, expected ), true, 'returns expected value' );
+	t.end();
+});
+
 tape( 'the function returns a reference to the output array', opts, function test( t ) {
 	var out;
 	var x;