Refactor to wrap specific implementation

Author: kgryte

lib/node_modules/@stdlib/stats/base/sdsmean/lib/ndarray.js

@@ -20,7 +20,7 @@
 
 // MODULES //
 
-var float64ToFloat32 = require( '@stdlib/number/float64/base/to-float32' );
+var sdsmeanors = require( '@stdlib/stats/base/sdsmeanors' ).ndarray;
 
 
 // MAIN //
@@ -45,23 +45,7 @@ var float64ToFloat32 = require( '@stdlib/number/float64/base/to-float32' );
 * // returns 1.25
 */
 function sdsmean( N, x, stride, offset ) {
-	var sum;
-	var ix;
-	var i;
-
-	if ( N <= 0 ) {
-		return NaN;
-	}
-	if ( N === 1 || stride === 0 ) {
-		return x[ offset ];
-	}
-	ix = offset;
-	sum = 0.0;
-	for ( i = 0; i < N; i++ ) {
-		sum += x[ ix ];
-		ix += stride;
-	}
-	return float64ToFloat32( sum / N );
+	return sdsmeanors( N, x, stride, offset );
 }
 
 

lib/node_modules/@stdlib/stats/base/sdsmean/lib/sdsmean.js

@@ -20,7 +20,7 @@
 
 // MODULES //
 
-var float64ToFloat32 = require( '@stdlib/number/float64/base/to-float32' );
+var sdsmeanors = require( '@stdlib/stats/base/sdsmeanors' );
 
 
 // MAIN //
@@ -43,27 +43,7 @@ var float64ToFloat32 = require( '@stdlib/number/float64/base/to-float32' );
 * // returns ~0.3333
 */
 function sdsmean( N, x, stride ) {
-	var sum;
-	var ix;
-	var i;
-
-	if ( N <= 0 ) {
-		return NaN;
-	}
-	if ( N === 1 || stride === 0 ) {
-		return x[ 0 ];
-	}
-	if ( stride < 0 ) {
-		ix = (1-N) * stride;
-	} else {
-		ix = 0;
-	}
-	sum = 0.0;
-	for ( i = 0; i < N; i++ ) {
-		sum += x[ ix ];
-		ix += stride;
-	}
-	return float64ToFloat32( sum / N );
+	return sdsmeanors( N, x, stride );
 }
 
 

lib/node_modules/@stdlib/stats/base/sdsmean/manifest.json

@@ -34,7 +34,9 @@
 				"-lm"
 			],
 			"libpath": [],
-			"dependencies": []
+			"dependencies": [
+				"@stdlib/stats/base/sdsmeanors"
+			]
 		}
 	]
 }

lib/node_modules/@stdlib/stats/base/sdsmean/src/sdsmean.c

@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/sdsmean.h"
+#include "stdlib/stats/base/sdsmeanors.h"
 #include <stdint.h>
 
 /**
@@ -28,25 +29,5 @@
 * @return        output value
 */
 float stdlib_strided_sdsmean( const int64_t N, const float *X, const int64_t stride ) {
-	double sum;
-	int64_t ix;
-	int64_t i;
-
-	if ( N <= 0 ) {
-		return 0.0 / 0.0; // NaN
-	}
-	if ( N == 1 || stride == 0 ) {
-		return X[ 0 ];
-	}
-	if ( stride < 0 ) {
-		ix = (1-N) * stride;
-	} else {
-		ix = 0;
-	}
-	sum = 0.0;
-	for ( i = 0; i < N; i++ ) {
-		sum += (double)X[ ix ];
-		ix += stride;
-	}
-	return sum / N;
+	return stdlib_strided_sdsmeanors( N, X, stride );
 }