fix: update implementation to preserve signed zeros

---
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: na
  - 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/ssumpw/lib/ndarray.js

@@ -20,7 +20,7 @@
 
 // MODULES //
 
-var float64ToFloat32 = require( '@stdlib/number/float64/base/to-float32' );
+var f32 = require( '@stdlib/number/float64/base/to-float32' );
 var floor = require( '@stdlib/math/base/special/floor' );
 
 
@@ -81,9 +81,10 @@ function ssumpw( N, x, strideX, offsetX ) {
 	}
 	if ( N < 8 ) {
 		// Use simple summation...
-		s = 0.0;
-		for ( i = 0; i < N; i++ ) {
-			s = float64ToFloat32( s + x[ ix ] );
+		s = x[ ix ];
+		ix += strideX;
+		for ( i = 1; i < N; i++ ) {
+			s = f32( s + x[ ix ] );
 			ix += strideX;
 		}
 		return s;
@@ -102,30 +103,30 @@ function ssumpw( N, x, strideX, offsetX ) {
 
 		M = N % 8;
 		for ( i = 8; i < N-M; i += 8 ) {
-			s0 = float64ToFloat32( s0 + x[ ix ] );
-			s1 = float64ToFloat32( s1 + x[ ix+strideX ] );
-			s2 = float64ToFloat32( s2 + x[ ix+(2*strideX) ] );
-			s3 = float64ToFloat32( s3 + x[ ix+(3*strideX) ] );
-			s4 = float64ToFloat32( s4 + x[ ix+(4*strideX) ] );
-			s5 = float64ToFloat32( s5 + x[ ix+(5*strideX) ] );
-			s6 = float64ToFloat32( s6 + x[ ix+(6*strideX) ] );
-			s7 = float64ToFloat32( s7 + x[ ix+(7*strideX) ] );
+			s0 = f32( s0 + x[ ix ] );
+			s1 = f32( s1 + x[ ix+strideX ] );
+			s2 = f32( s2 + x[ ix+(2*strideX) ] );
+			s3 = f32( s3 + x[ ix+(3*strideX) ] );
+			s4 = f32( s4 + x[ ix+(4*strideX) ] );
+			s5 = f32( s5 + x[ ix+(5*strideX) ] );
+			s6 = f32( s6 + x[ ix+(6*strideX) ] );
+			s7 = f32( s7 + x[ ix+(7*strideX) ] );
 			ix += 8 * strideX;
 		}
 		// Pairwise sum the accumulators:
-		s = float64ToFloat32( float64ToFloat32( float64ToFloat32(s0+s1) + float64ToFloat32(s2+s3) ) + float64ToFloat32( float64ToFloat32(s4+s5) + float64ToFloat32(s6+s7) ) ); // eslint-disable-line max-len
+		s = f32( f32( f32(s0+s1) + f32(s2+s3) ) + f32( f32(s4+s5) + f32(s6+s7) ) ); // eslint-disable-line max-len
 
 		// Clean-up loop...
 		for ( i; i < N; i++ ) {
-			s = float64ToFloat32( s + x[ ix ] );
+			s = f32( s + x[ ix ] );
 			ix += strideX;
 		}
 		return s;
 	}
 	// Recurse by dividing by two, but avoiding non-multiples of unroll factor...
 	n = floor( N/2 );
 	n -= n % 8;
-	return float64ToFloat32( ssumpw( n, x, strideX, ix ) + ssumpw( N-n, x, strideX, ix+(n*strideX) ) ); // eslint-disable-line max-len
+	return f32( ssumpw( n, x, strideX, ix ) + ssumpw( N-n, x, strideX, ix+(n*strideX) ) ); // eslint-disable-line max-len
 }
 
 

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

@@ -22,6 +22,7 @@
 
 var tape = require( 'tape' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
+var isNegativeZero = require( '@stdlib/math/base/assert/is-negative-zero' );
 var Float32Array = require( '@stdlib/array/float32' );
 var ssumpw = require( './../lib/ndarray.js' );
 
@@ -74,6 +75,17 @@ tape( 'the function calculates the sum of all strided array elements', function
 	t.end();
 });
 
+tape( 'the function preserves the sign of zero', function test( t ) {
+	var x;
+	var v;
+
+	x = new Float32Array( [ -0.0, -0.0, -0.0, -0.0, -0.0 ] );
+	v = ssumpw( x.length, x, 1, 0 );
+	t.strictEqual( isNegativeZero( v ), true, 'returns expected value' );
+
+	t.end();
+});
+
 tape( 'if provided an `N` parameter less than or equal to `0`, the function returns `0.0`', function test( t ) {
 	var x;
 	var v;

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

@@ -23,6 +23,7 @@
 var resolve = require( 'path' ).resolve;
 var tape = require( 'tape' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
+var isNegativeZero = require( '@stdlib/math/base/assert/is-negative-zero' );
 var Float32Array = require( '@stdlib/array/float32' );
 var tryRequire = require( '@stdlib/utils/try-require' );
 
@@ -83,6 +84,17 @@ tape( 'the function calculates the sum of all strided array elements', opts, fun
 	t.end();
 });
 
+tape( 'the function preserves the sign of zero', opts, function test( t ) {
+	var x;
+	var v;
+
+	x = new Float32Array( [ -0.0, -0.0, -0.0, -0.0, -0.0 ] );
+	v = ssumpw( x.length, x, 1, 0 );
+	t.strictEqual( isNegativeZero( v ), true, 'returns expected value' );
+
+	t.end();
+});
+
 tape( 'if provided an `N` parameter less than or equal to `0`, the function returns `0.0`', opts, function test( t ) {
 	var x;
 	var v;