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: na
  - 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/gsumors/lib/accessors.js

@@ -58,8 +58,9 @@ function gsumors( N, x, strideX, offsetX ) {
 	if ( strideX === 0 ) {
 		return N * get( xbuf, ix );
 	}
-	s = 0.0;
-	for ( i = 0; i < N; i++ ) {
+	s = get( xbuf, ix );
+	ix += strideX;
+	for ( i = 1; i < N; i++ ) {
 		s += get( xbuf, ix );
 		ix += strideX;
 	}

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

@@ -53,9 +53,8 @@ function gsumors( N, x, strideX, offsetX ) {
 	var o;
 	var i;
 
-	s = 0.0;
 	if ( N <= 0 ) {
-		return s;
+		return 0.0;
 	}
 	o = arraylike2object( x );
 	if ( o.accessorProtocol ) {
@@ -65,10 +64,12 @@ function gsumors( N, x, strideX, offsetX ) {
 	if ( strideX === 0 ) {
 		return N * x[ ix ];
 	}
+	s = x[ ix ];
+	ix += strideX;
 
 	// If the stride is equal to `1`, use unrolled loops...
 	if ( strideX === 1 ) {
-		m = N % M;
+		m = (N-1) % M;
 
 		// If we have a remainder, run a clean-up loop...
 		if ( m > 0 ) {
@@ -80,13 +81,13 @@ function gsumors( N, x, strideX, offsetX ) {
 		if ( N < M ) {
 			return s;
 		}
-		for ( i = m; i < N; i += M ) {
+		for ( i = m; i < N-1; i += M ) {
 			s += x[ix] + x[ix+1] + x[ix+2] + x[ix+3] + x[ix+4] + x[ix+5];
 			ix += M;
 		}
 		return s;
 	}
-	for ( i = 0; i < N; i++ ) {
+	for ( i = 1; i < N; i++ ) {
 		s += x[ ix ];
 		ix += strideX;
 	}