diff --git a/blas/ext/base/scusumkbn/coverage.ndjson b/blas/ext/base/scusumkbn/coverage.ndjson new file mode 100644 index 000000000..d3777bc72 --- /dev/null +++ b/blas/ext/base/scusumkbn/coverage.ndjson @@ -0,0 +1 @@ +[466,466,100,31,31,100,4,4,100,466,466,100,"9b2ca2ce1281165deebb28310a8e1c73ff409245","2024-03-24 00:46:56 -0700"] diff --git a/blas/ext/base/scusumkbn/index.html b/blas/ext/base/scusumkbn/index.html new file mode 100644 index 000000000..4fea54e5a --- /dev/null +++ b/blas/ext/base/scusumkbn/index.html @@ -0,0 +1,206 @@ + + + + + + Code coverage report for blas/ext/base/scusumkbn/lib + + + + + + + + + +
+
+

All files blas/ext/base/scusumkbn/lib

+
+ +
+ 100% + Statements + 466/466 +
+ + +
+ 100% + Branches + 31/31 +
+ + +
+ 100% + Functions + 4/4 +
+ + +
+ 100% + Lines + 466/466 +
+ + +
+

+ Press n or j to go to the next uncovered block, b, p or k for the previous block. +

+ +
+
+
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +
FileStatementsBranchesFunctionsLines
index.js +
+ 		100%73/73100%3/3100%0/0100%73/73
main.js +
+ 		100%35/35100%1/1100%0/0100%35/35
native.js +
+ 		100%35/35100%1/1100%0/0100%35/35
ndarray.js +
+ 		100%96/96100%8/8100%1/1100%96/96
ndarray.native.js +
+ 		100%69/69100%2/2100%1/1100%69/69
scusumkbn.js +
+ 		100%101/101100%14/14100%1/1100%101/101
scusumkbn.native.js +
+ 		100%57/57100%2/2100%1/1100%57/57
+
+
+
+
+ Code coverage generated by + istanbul + at 2024-03-24T07:48:47.237Z +
+ + + + + + + \ No newline at end of file diff --git a/blas/ext/base/scusumkbn/index.js.html b/blas/ext/base/scusumkbn/index.js.html new file mode 100644 index 000000000..c221f9e5f --- /dev/null +++ b/blas/ext/base/scusumkbn/index.js.html @@ -0,0 +1,304 @@ + + + + + + Code coverage report for blas/ext/base/scusumkbn/lib/index.js + + + + + + + + + +
+
+

All files / blas/ext/base/scusumkbn/lib index.js

+
+ +
+ 100% + Statements + 73/73 +
+ + +
+ 100% + Branches + 3/3 +
+ + +
+ 100% + Functions + 0/0 +
+ + +
+ 100% + Lines + 73/73 +
+ + +
+

+ Press n or j to go to the next uncovered block, b, p or k for the previous block. +

+ +
+
+ 

+
1
+2
+3
+4
+5
+6
+7
+8
+9
+10
+11
+12
+13
+14
+15
+16
+17
+18
+19
+20
+21
+22
+23
+24
+25
+26
+27
+28
+29
+30
+31
+32
+33
+34
+35
+36
+37
+38
+39
+40
+41
+42
+43
+44
+45
+46
+47
+48
+49
+50
+51
+52
+53
+54
+55
+56
+57
+58
+59
+60
+61
+62
+63
+64
+65
+66
+67
+68
+69
+70
+71
+72
+73
+743x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+1x
+3x
+2x
+2x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+ 
/**
+* @license Apache-2.0
+*
+* Copyright (c) 2020 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+ 
+'use strict';
+ 
+/**
+* Compute the cumulative sum of single-precision floating-point strided array elements using an improved Kahan–Babuška algorithm.
+*
+* @module @stdlib/blas/ext/base/scusumkbn
+*
+* @example
+* var Float32Array = require( '@stdlib/array/float32' );
+* var scusumkbn = require( '@stdlib/blas/ext/base/scusumkbn' );
+*
+* var x = new Float32Array( [ 1.0, -2.0, 2.0 ] );
+* var y = new Float32Array( x.length );
+* var N = x.length;
+*
+* scusumkbn( N, 0.0, x, 1, y, 1 );
+* // y => <Float32Array>[ 1.0, -1.0, 1.0 ]
+*
+* @example
+* var Float32Array = require( '@stdlib/array/float32' );
+* var floor = require( '@stdlib/math/base/special/floor' );
+* var scusumkbn = require( '@stdlib/blas/ext/base/scusumkbn' );
+*
+* var x = new Float32Array( [ 2.0, 1.0, 2.0, -2.0, -2.0, 2.0, 3.0, 4.0 ] );
+* var y = new Float32Array( x.length );
+* var N = floor( x.length / 2 );
+*
+* scusumkbn.ndarray( N, 0.0, x, 2, 1, y, 1, 0 );
+* // y => <Float32Array>[ 1.0, -1.0, 1.0, 5.0, 0.0, 0.0, 0.0, 0.0 ]
+*/
+ 
+// MODULES //
+ 
+var join = require( 'path' ).join;
+var tryRequire = require( '@stdlib/utils/try-require' );
+var isError = require( '@stdlib/assert/is-error' );
+var main = require( './main.js' );
+ 
+ 
+// MAIN //
+ 
+var scusumkbn;
+var tmp = tryRequire( join( __dirname, './native.js' ) );
+if ( isError( tmp ) ) {
+	scusumkbn = main;
+} else {
+	scusumkbn = tmp;
+}
+ 
+ 
+// EXPORTS //
+ 
+module.exports = scusumkbn;
+ 
+// exports: { "ndarray": "scusumkbn.ndarray" }
+ 
+ +
+
+
+ Code coverage generated by + istanbul + at 2024-03-24T07:48:47.237Z +
+ + + + + + + \ No newline at end of file diff --git a/blas/ext/base/scusumkbn/main.js.html b/blas/ext/base/scusumkbn/main.js.html new file mode 100644 index 000000000..6ab052914 --- /dev/null +++ b/blas/ext/base/scusumkbn/main.js.html @@ -0,0 +1,190 @@ + + + + + + Code coverage report for blas/ext/base/scusumkbn/lib/main.js + + + + + + + + + +
+
+

All files / blas/ext/base/scusumkbn/lib main.js

+
+ +
+ 100% + Statements + 35/35 +
+ + +
+ 100% + Branches + 1/1 +
+ + +
+ 100% + Functions + 0/0 +
+ + +
+ 100% + Lines + 35/35 +
+ + +
+

+ Press n or j to go to the next uncovered block, b, p or k for the previous block. +

+ +
+
+ 

+
1
+2
+3
+4
+5
+6
+7
+8
+9
+10
+11
+12
+13
+14
+15
+16
+17
+18
+19
+20
+21
+22
+23
+24
+25
+26
+27
+28
+29
+30
+31
+32
+33
+34
+35
+361x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+ 
/**
+* @license Apache-2.0
+*
+* Copyright (c) 2020 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+ 
+'use strict';
+ 
+// MODULES //
+ 
+var setReadOnly = require( '@stdlib/utils/define-nonenumerable-read-only-property' );
+var scusumkbn = require( './scusumkbn.js' );
+var ndarray = require( './ndarray.js' );
+ 
+ 
+// MAIN //
+ 
+setReadOnly( scusumkbn, 'ndarray', ndarray );
+ 
+ 
+// EXPORTS //
+ 
+module.exports = scusumkbn;
+ 
+ +
+
+
+ Code coverage generated by + istanbul + at 2024-03-24T07:48:47.237Z +
+ + + + + + + \ No newline at end of file diff --git a/blas/ext/base/scusumkbn/native.js.html b/blas/ext/base/scusumkbn/native.js.html new file mode 100644 index 000000000..72cb8aae2 --- /dev/null +++ b/blas/ext/base/scusumkbn/native.js.html @@ -0,0 +1,190 @@ + + + + + + Code coverage report for blas/ext/base/scusumkbn/lib/native.js + + + + + + + + + +
+
+

All files / blas/ext/base/scusumkbn/lib native.js

+
+ +
+ 100% + Statements + 35/35 +
+ + +
+ 100% + Branches + 1/1 +
+ + +
+ 100% + Functions + 0/0 +
+ + +
+ 100% + Lines + 35/35 +
+ + +
+

+ Press n or j to go to the next uncovered block, b, p or k for the previous block. +

+ +
+
+ 

+
1
+2
+3
+4
+5
+6
+7
+8
+9
+10
+11
+12
+13
+14
+15
+16
+17
+18
+19
+20
+21
+22
+23
+24
+25
+26
+27
+28
+29
+30
+31
+32
+33
+34
+35
+361x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+1x
+ 
/**
+* @license Apache-2.0
+*
+* Copyright (c) 2020 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+ 
+'use strict';
+ 
+// MODULES //
+ 
+var setReadOnly = require( '@stdlib/utils/define-nonenumerable-read-only-property' );
+var scusumkbn = require( './scusumkbn.native.js' );
+var ndarray = require( './ndarray.native.js' );
+ 
+ 
+// MAIN //
+ 
+setReadOnly( scusumkbn, 'ndarray', ndarray );
+ 
+ 
+// EXPORTS //
+ 
+module.exports = scusumkbn;
+ 
+ +
+
+
+ Code coverage generated by + istanbul + at 2024-03-24T07:48:47.237Z +
+ + + + + + + \ No newline at end of file diff --git a/blas/ext/base/scusumkbn/ndarray.js.html b/blas/ext/base/scusumkbn/ndarray.js.html new file mode 100644 index 000000000..c60a58d14 --- /dev/null +++ b/blas/ext/base/scusumkbn/ndarray.js.html @@ -0,0 +1,373 @@ + + + + + + Code coverage report for blas/ext/base/scusumkbn/lib/ndarray.js + + + + + + + + + +
+
+

All files / blas/ext/base/scusumkbn/lib ndarray.js

+
+ +
+ 100% + Statements + 96/96 +
+ + +
+ 100% + Branches + 8/8 +
+ + +
+ 100% + Functions + 1/1 +
+ + +
+ 100% + Lines + 96/96 +
+ + +
+

+ Press n or j to go to the next uncovered block, b, p or k for the previous block. +

+ +
+
+ 

+
1
+2
+3
+4
+5
+6
+7
+8
+9
+10
+11
+12
+13
+14
+15
+16
+17
+18
+19
+20
+21
+22
+23
+24
+25
+26
+27
+28
+29
+30
+31
+32
+33
+34
+35
+36
+37
+38
+39
+40
+41
+42
+43
+44
+45
+46
+47
+48
+49
+50
+51
+52
+53
+54
+55
+56
+57
+58
+59
+60
+61
+62
+63
+64
+65
+66
+67
+68
+69
+70
+71
+72
+73
+74
+75
+76
+77
+78
+79
+80
+81
+82
+83
+84
+85
+86
+87
+88
+89
+90
+91
+92
+93
+94
+95
+96
+972x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+15x
+15x
+15x
+15x
+15x
+15x
+15x
+15x
+15x
+15x
+2x
+2x
+13x
+13x
+13x
+13x
+15x
+1045x
+1045x
+1045x
+1017x
+1045x
+28x
+28x
+1045x
+1045x
+1045x
+1045x
+1045x
+13x
+15x
+2x
+2x
+2x
+2x
+2x
+ 
/**
+* @license Apache-2.0
+*
+* Copyright (c) 2020 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+ 
+'use strict';
+ 
+// MODULES //
+ 
+var float64ToFloat32 = require( '@stdlib/number/float64/base/to-float32' );
+var abs = require( '@stdlib/math/base/special/abs' );
+ 
+ 
+// MAIN //
+ 
+/**
+* Computes the cumulative sum of single-precision floating-point strided array elements using an improved Kahan–Babuška algorithm.
+*
+* ## Method
+*
+* -   This implementation uses an "improved Kahan–Babuška algorithm", as described by Neumaier (1974).
+*
+* ## References
+*
+* -   Neumaier, Arnold. 1974. "Rounding Error Analysis of Some Methods for Summing Finite Sums." _Zeitschrift Für Angewandte Mathematik Und Mechanik_ 54 (1): 39–51. doi:[10.1002/zamm.19740540106](https://doi.org/10.1002/zamm.19740540106).
+*
+* @param {PositiveInteger} N - number of indexed elements
+* @param {number} sum - initial sum
+* @param {Float32Array} x - input array
+* @param {integer} strideX - `x` stride length
+* @param {NonNegativeInteger} offsetX - starting index for `x`
+* @param {Float32Array} y - output array
+* @param {integer} strideY - `y` stride length
+* @param {NonNegativeInteger} offsetY - starting index for `y`
+* @returns {Float32Array} output array
+*
+* @example
+* var Float32Array = require( '@stdlib/array/float32' );
+* var floor = require( '@stdlib/math/base/special/floor' );
+*
+* var x = new Float32Array( [ 2.0, 1.0, 2.0, -2.0, -2.0, 2.0, 3.0, 4.0 ] );
+* var y = new Float32Array( x.length );
+* var N = floor( x.length / 2 );
+*
+* var v = scusumkbn( N, 0.0, x, 2, 1, y, 1, 0 );
+* // returns <Float32Array>[ 1.0, -1.0, 1.0, 5.0, 0.0, 0.0, 0.0, 0.0 ]
+*/
+function scusumkbn( N, sum, x, strideX, offsetX, y, strideY, offsetY ) {
+	var ix;
+	var iy;
+	var s;
+	var v;
+	var t;
+	var c;
+	var i;
+ 
+	if ( N <= 0 ) {
+		return y;
+	}
+	ix = offsetX;
+	iy = offsetY;
+	s = sum;
+	c = 0.0;
+	for ( i = 0; i < N; i++ ) {
+		v = x[ ix ];
+		t = float64ToFloat32( s + v );
+		if ( abs( s ) >= abs( v ) ) {
+			c = float64ToFloat32( c + float64ToFloat32( float64ToFloat32( s-t ) + v ) ); // eslint-disable-line max-len
+		} else {
+			c = float64ToFloat32( c + float64ToFloat32( float64ToFloat32( v-t ) + s ) ); // eslint-disable-line max-len
+		}
+		s = t;
+		y[ iy ] = float64ToFloat32( s + c );
+		ix += strideX;
+		iy += strideY;
+	}
+	return y;
+}
+ 
+ 
+// EXPORTS //
+ 
+module.exports = scusumkbn;
+ 
+ +
+
+
+ Code coverage generated by + istanbul + at 2024-03-24T07:48:47.237Z +
+ + + + + + + \ No newline at end of file diff --git a/blas/ext/base/scusumkbn/ndarray.native.js.html b/blas/ext/base/scusumkbn/ndarray.native.js.html new file mode 100644 index 000000000..227d2a262 --- /dev/null +++ b/blas/ext/base/scusumkbn/ndarray.native.js.html @@ -0,0 +1,292 @@ + + + + + + Code coverage report for blas/ext/base/scusumkbn/lib/ndarray.native.js + + + + + + + + + +
+
+

All files / blas/ext/base/scusumkbn/lib ndarray.native.js

+
+ +
+ 100% + Statements + 69/69 +
+ + +
+ 100% + Branches + 2/2 +
+ + +
+ 100% + Functions + 1/1 +
+ + +
+ 100% + Lines + 69/69 +
+ + +
+

+ Press n or j to go to the next uncovered block, b, p or k for the previous block. +

+ +
+
+ 

+
1
+2
+3
+4
+5
+6
+7
+8
+9
+10
+11
+12
+13
+14
+15
+16
+17
+18
+19
+20
+21
+22
+23
+24
+25
+26
+27
+28
+29
+30
+31
+32
+33
+34
+35
+36
+37
+38
+39
+40
+41
+42
+43
+44
+45
+46
+47
+48
+49
+50
+51
+52
+53
+54
+55
+56
+57
+58
+59
+60
+61
+62
+63
+64
+65
+66
+67
+68
+69
+702x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+15x
+15x
+15x
+15x
+15x
+15x
+15x
+15x
+15x
+15x
+15x
+2x
+2x
+2x
+2x
+2x
+ 
/**
+* @license Apache-2.0
+*
+* Copyright (c) 2020 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+ 
+'use strict';
+ 
+// MODULES //
+ 
+var minViewBufferIndex = require( '@stdlib/strided/base/min-view-buffer-index' );
+var offsetView = require( '@stdlib/strided/base/offset-view' );
+var addon = require( './scusumkbn.native.js' );
+ 
+ 
+// MAIN //
+ 
+/**
+* Computes the cumulative sum of single-precision floating-point strided array elements using an improved Kahan–Babuška algorithm.
+*
+* @param {PositiveInteger} N - number of indexed elements
+* @param {number} sum - initial sum
+* @param {Float32Array} x - input array
+* @param {integer} strideX - `x` stride length
+* @param {NonNegativeInteger} offsetX - starting index for `x`
+* @param {Float32Array} y - output array
+* @param {integer} strideY - `y` stride length
+* @param {NonNegativeInteger} offsetY - starting index for `y`
+* @returns {Float32Array} output array
+*
+* @example
+* var Float32Array = require( '@stdlib/array/float32' );
+* var floor = require( '@stdlib/math/base/special/floor' );
+*
+* var x = new Float32Array( [ 2.0, 1.0, 2.0, -2.0, -2.0, 2.0, 3.0, 4.0 ] );
+* var y = new Float32Array( x.length );
+* var N = floor( x.length / 2 );
+*
+* var v = scusumkbn( N, 0.0, x, 2, 1, y, 1, 0 );
+* // returns <Float32Array>[ 1.0, -1.0, 1.0, 5.0, 0.0, 0.0, 0.0, 0.0 ]
+*/
+function scusumkbn( N, sum, x, strideX, offsetX, y, strideY, offsetY ) {
+	var viewX;
+	var viewY;
+	offsetX = minViewBufferIndex( N, strideX, offsetX );
+	offsetY = minViewBufferIndex( N, strideY, offsetY );
+ 
+	viewX = offsetView( x, offsetX );
+	viewY = offsetView( y, offsetY );
+	addon( N, sum, viewX, strideX, viewY, strideY );
+	return y;
+}
+ 
+ 
+// EXPORTS //
+ 
+module.exports = scusumkbn;
+ 
+ +
+
+
+ Code coverage generated by + istanbul + at 2024-03-24T07:48:47.237Z +
+ + + + + + + \ No newline at end of file diff --git a/blas/ext/base/scusumkbn/scusumkbn.js.html b/blas/ext/base/scusumkbn/scusumkbn.js.html new file mode 100644 index 000000000..412ddfd7c --- /dev/null +++ b/blas/ext/base/scusumkbn/scusumkbn.js.html @@ -0,0 +1,388 @@ + + + + + + Code coverage report for blas/ext/base/scusumkbn/lib/scusumkbn.js + + + + + + + + + +
+
+

All files / blas/ext/base/scusumkbn/lib scusumkbn.js

+
+ +
+ 100% + Statements + 101/101 +
+ + +
+ 100% + Branches + 14/14 +
+ + +
+ 100% + Functions + 1/1 +
+ + +
+ 100% + Lines + 101/101 +
+ + +
+

+ Press n or j to go to the next uncovered block, b, p or k for the previous block. +

+ +
+
+ 

+
1
+2
+3
+4
+5
+6
+7
+8
+9
+10
+11
+12
+13
+14
+15
+16
+17
+18
+19
+20
+21
+22
+23
+24
+25
+26
+27
+28
+29
+30
+31
+32
+33
+34
+35
+36
+37
+38
+39
+40
+41
+42
+43
+44
+45
+46
+47
+48
+49
+50
+51
+52
+53
+54
+55
+56
+57
+58
+59
+60
+61
+62
+63
+64
+65
+66
+67
+68
+69
+70
+71
+72
+73
+74
+75
+76
+77
+78
+79
+80
+81
+82
+83
+84
+85
+86
+87
+88
+89
+90
+91
+92
+93
+94
+95
+96
+97
+98
+99
+100
+101
+1022x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+2x
+14x
+14x
+14x
+14x
+14x
+14x
+14x
+14x
+14x
+14x
+2x
+2x
+14x
+2x
+14x
+10x
+10x
+14x
+2x
+14x
+10x
+10x
+12x
+12x
+14x
+1040x
+1040x
+1040x
+1016x
+1040x
+24x
+24x
+1040x
+1040x
+1040x
+1040x
+1040x
+12x
+14x
+2x
+2x
+2x
+2x
+2x
+ 
/**
+* @license Apache-2.0
+*
+* Copyright (c) 2020 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+ 
+'use strict';
+ 
+// MODULES //
+ 
+var float64ToFloat32 = require( '@stdlib/number/float64/base/to-float32' );
+var abs = require( '@stdlib/math/base/special/abs' );
+ 
+ 
+// MAIN //
+ 
+/**
+* Computes the cumulative sum of single-precision floating-point strided array elements using an improved Kahan–Babuška algorithm.
+*
+* ## Method
+*
+* -   This implementation uses an "improved Kahan–Babuška algorithm", as described by Neumaier (1974).
+*
+* ## References
+*
+* -   Neumaier, Arnold. 1974. "Rounding Error Analysis of Some Methods for Summing Finite Sums." _Zeitschrift Für Angewandte Mathematik Und Mechanik_ 54 (1): 39–51. doi:[10.1002/zamm.19740540106](https://doi.org/10.1002/zamm.19740540106).
+*
+* @param {PositiveInteger} N - number of indexed elements
+* @param {number} sum - initial sum
+* @param {Float32Array} x - input array
+* @param {integer} strideX - `x` stride length
+* @param {Float32Array} y - output array
+* @param {integer} strideY - `y` stride length
+* @returns {Float32Array} output array
+*
+* @example
+* var Float32Array = require( '@stdlib/array/float32' );
+*
+* var x = new Float32Array( [ 1.0, -2.0, 2.0 ] );
+* var y = new Float32Array( x.length );
+* var N = x.length;
+*
+* var v = scusumkbn( N, 0.0, x, 1, y, 1 );
+* // returns <Float32Array>[ 1.0, -1.0, 1.0 ]
+*/
+function scusumkbn( N, sum, x, strideX, y, strideY ) {
+	var ix;
+	var iy;
+	var s;
+	var v;
+	var t;
+	var c;
+	var i;
+ 
+	if ( N <= 0 ) {
+		return y;
+	}
+	if ( strideX < 0 ) {
+		ix = (1-N) * strideX;
+	} else {
+		ix = 0;
+	}
+	if ( strideY < 0 ) {
+		iy = (1-N) * strideY;
+	} else {
+		iy = 0;
+	}
+	s = sum;
+	c = 0.0;
+	for ( i = 0; i < N; i++ ) {
+		v = x[ ix ];
+		t = float64ToFloat32( s + v );
+		if ( abs( s ) >= abs( v ) ) {
+			c = float64ToFloat32( c + float64ToFloat32( float64ToFloat32( s-t ) + v ) ); // eslint-disable-line max-len
+		} else {
+			c = float64ToFloat32( c + float64ToFloat32( float64ToFloat32( v-t ) + s ) ); // eslint-disable-line max-len
+		}
+		s = t;
+		y[ iy ] = float64ToFloat32( s + c );
+		ix += strideX;
+		iy += strideY;
+	}
+	return y;
+}
+ 
+ 
+// EXPORTS //
+ 
+module.exports = scusumkbn;
+ 
+ +
+
+
+ Code coverage generated by + istanbul + at 2024-03-24T07:48:47.237Z +
+ + + + + + + \ No newline at end of file diff --git a/blas/ext/base/scusumkbn/scusumkbn.native.js.html b/blas/ext/base/scusumkbn/scusumkbn.native.js.html new file mode 100644 index 000000000..eadc72432 --- /dev/null +++ b/blas/ext/base/scusumkbn/scusumkbn.native.js.html @@ -0,0 +1,256 @@ + + + + + + Code coverage report for blas/ext/base/scusumkbn/lib/scusumkbn.native.js + + + + + + + + + +
+
+

All files / blas/ext/base/scusumkbn/lib scusumkbn.native.js

+
+ +
+ 100% + Statements + 57/57 +
+ + +
+ 100% + Branches + 2/2 +
+ + +
+ 100% + Functions + 1/1 +
+ + +
+ 100% + Lines + 57/57 +
+ + +
+

+ Press n or j to go to the next uncovered block, b, p or k for the previous block. +

+ +
+
+ 

+
1
+2
+3
+4
+5
+6
+7
+8
+9
+10
+11
+12
+13
+14
+15
+16
+17
+18
+19
+20
+21
+22
+23
+24
+25
+26
+27
+28
+29
+30
+31
+32
+33
+34
+35
+36
+37
+38
+39
+40
+41
+42
+43
+44
+45
+46
+47
+48
+49
+50
+51
+52
+53
+54
+55
+56
+57
+583x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+3x
+29x
+29x
+29x
+29x
+3x
+3x
+3x
+3x
+3x
+ 
/**
+* @license Apache-2.0
+*
+* Copyright (c) 2020 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+ 
+'use strict';
+ 
+// MODULES //
+ 
+var addon = require( './../src/addon.node' );
+ 
+ 
+// MAIN //
+ 
+/**
+* Computes the cumulative sum of single-precision floating-point strided array elements using an improved Kahan–Babuška algorithm.
+*
+* @param {PositiveInteger} N - number of indexed elements
+* @param {number} sum - initial sum
+* @param {Float32Array} x - input array
+* @param {integer} strideX - `x` stride length
+* @param {Float32Array} y - output array
+* @param {integer} strideY - `y` stride length
+* @returns {Float32Array} output array
+*
+* @example
+* var Float32Array = require( '@stdlib/array/float32' );
+*
+* var x = new Float32Array( [ 1.0, -2.0, 2.0 ] );
+* var y = new Float32Array( x.length );
+* var N = x.length;
+*
+* var v = scusumkbn( N, 0.0, x, 1, y, 1 );
+* // returns <Float32Array>[ 1.0, -1.0, 1.0 ]
+*/
+function scusumkbn( N, sum, x, strideX, y, strideY ) {
+	addon( N, sum, x, strideX, y, strideY );
+	return y;
+}
+ 
+ 
+// EXPORTS //
+ 
+module.exports = scusumkbn;
+ 
+ +
+
+
+ Code coverage generated by + istanbul + at 2024-03-24T07:48:47.237Z +
+ + + + + + + \ No newline at end of file