Merge branch 'develop' into feat/lanczos-sum

Author: anandkaranubc

lib/node_modules/@stdlib/blas/ext/base/README.md

@@ -44,6 +44,8 @@ var o = ns;
 <div class="namespace-toc">
 
 -   <span class="signature">[`cfill( N, alpha, x, strideX )`][@stdlib/blas/ext/base/cfill]</span><span class="delimiter">: </span><span class="description">fill a single-precision complex floating-point strided array with a specified scalar constant.</span>
+-   <span class="signature">[`csum( N, x, strideX )`][@stdlib/blas/ext/base/csum]</span><span class="delimiter">: </span><span class="description">calculate the sum of single-precision complex floating-point strided array elements.</span>
+-   <span class="signature">[`csumkbn( N, x, strideX )`][@stdlib/blas/ext/base/csumkbn]</span><span class="delimiter">: </span><span class="description">calculate the sum of single-precision complex floating-point strided array elements using an improved Kahan–Babuška algorithm.</span>
 -   <span class="signature">[`dapx( N, alpha, x, strideX )`][@stdlib/blas/ext/base/dapx]</span><span class="delimiter">: </span><span class="description">add a scalar constant to each element in a double-precision floating-point strided array.</span>
 -   <span class="signature">[`dapxsum( N, alpha, x, strideX )`][@stdlib/blas/ext/base/dapxsum]</span><span class="delimiter">: </span><span class="description">add a scalar constant to each double-precision floating-point strided array element and compute the sum.</span>
 -   <span class="signature">[`dapxsumkbn( N, alpha, x, strideX )`][@stdlib/blas/ext/base/dapxsumkbn]</span><span class="delimiter">: </span><span class="description">add a scalar constant to each double-precision floating-point strided array element and compute the sum using an improved Kahan–Babuška algorithm.</span>
@@ -105,6 +107,7 @@ var o = ns;
 -   <span class="signature">[`gfillBy( N, x, strideX, clbk[, thisArg] )`][@stdlib/blas/ext/base/gfill-by]</span><span class="delimiter">: </span><span class="description">fill a strided array according to a provided callback function.</span>
 -   <span class="signature">[`gfill( N, alpha, x, strideX )`][@stdlib/blas/ext/base/gfill]</span><span class="delimiter">: </span><span class="description">fill a strided array with a specified scalar constant.</span>
 -   <span class="signature">[`gnannsumkbn( N, x, strideX, out, strideOut )`][@stdlib/blas/ext/base/gnannsumkbn]</span><span class="delimiter">: </span><span class="description">calculate the sum of strided array elements, ignoring `NaN` values and using an improved Kahan–Babuška algorithm.</span>
+-   <span class="signature">[`gnannsumpw( N, x, strideX, out, strideOut )`][@stdlib/blas/ext/base/gnannsumpw]</span><span class="delimiter">: </span><span class="description">calculate the sum of strided array elements, ignoring `NaN` values and using pairwise summation.</span>
 -   <span class="signature">[`gnansum( N, x, strideX )`][@stdlib/blas/ext/base/gnansum]</span><span class="delimiter">: </span><span class="description">calculate the sum of strided array elements, ignoring `NaN` values.</span>
 -   <span class="signature">[`gnansumkbn( N, x, strideX )`][@stdlib/blas/ext/base/gnansumkbn]</span><span class="delimiter">: </span><span class="description">calculate the sum of strided array elements, ignoring `NaN` values and using an improved Kahan–Babuška algorithm.</span>
 -   <span class="signature">[`gnansumkbn2( N, x, strideX )`][@stdlib/blas/ext/base/gnansumkbn2]</span><span class="delimiter">: </span><span class="description">calculate the sum of strided array elements, ignoring `NaN` values and using a second-order iterative Kahan–Babuška algorithm.</span>
@@ -161,6 +164,8 @@ var o = ns;
 -   <span class="signature">[`ssumpw( N, x, strideX )`][@stdlib/blas/ext/base/ssumpw]</span><span class="delimiter">: </span><span class="description">calculate the sum of single-precision floating-point strided array elements using pairwise summation.</span>
 -   <span class="signature">[`wasm`][@stdlib/blas/ext/base/wasm]</span><span class="delimiter">: </span><span class="description">extensions to basic linear algebra subprograms (BLAS) compiled to WebAssembly.</span>
 -   <span class="signature">[`zfill( N, alpha, x, strideX )`][@stdlib/blas/ext/base/zfill]</span><span class="delimiter">: </span><span class="description">fill a double-precision complex floating-point strided array with a specified scalar constant.</span>
+-   <span class="signature">[`zsum( N, x, strideX )`][@stdlib/blas/ext/base/zsum]</span><span class="delimiter">: </span><span class="description">calculate the sum of double-precision complex floating-point strided array elements.</span>
+-   <span class="signature">[`zsumkbn( N, x, strideX )`][@stdlib/blas/ext/base/zsumkbn]</span><span class="delimiter">: </span><span class="description">calculate the sum of double-precision complex floating-point strided array elements using an improved Kahan–Babuška algorithm.</span>
 
 </div>
 
@@ -213,6 +218,10 @@ console.log( objectKeys( ns ) );
 
 [@stdlib/blas/ext/base/cfill]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/cfill
 
+[@stdlib/blas/ext/base/csum]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/csum
+
+[@stdlib/blas/ext/base/csumkbn]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/csumkbn
+
 [@stdlib/blas/ext/base/dapx]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/dapx
 
 [@stdlib/blas/ext/base/dapxsum]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/dapxsum
@@ -335,6 +344,8 @@ console.log( objectKeys( ns ) );
 
 [@stdlib/blas/ext/base/gnannsumkbn]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/gnannsumkbn
 
+[@stdlib/blas/ext/base/gnannsumpw]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/gnannsumpw
+
 [@stdlib/blas/ext/base/gnansum]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/gnansum
 
 [@stdlib/blas/ext/base/gnansumkbn]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/gnansumkbn
@@ -447,6 +458,10 @@ console.log( objectKeys( ns ) );
 
 [@stdlib/blas/ext/base/zfill]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/zfill
 
+[@stdlib/blas/ext/base/zsum]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/zsum
+
+[@stdlib/blas/ext/base/zsumkbn]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/zsumkbn
+
 <!-- </toc-links> -->
 
 </section>

lib/node_modules/@stdlib/blas/ext/base/ndarray/README.md

@@ -45,12 +45,14 @@ The namespace exposes the following APIs:
 
 <div class="namespace-toc">
 
+-   <span class="signature">[`csum( arrays )`][@stdlib/blas/ext/base/ndarray/csum]</span><span class="delimiter">: </span><span class="description">compute the sum of all elements in a one-dimensional single-precision complex floating-point ndarray.</span>
 -   <span class="signature">[`dcusum( arrays )`][@stdlib/blas/ext/base/ndarray/dcusum]</span><span class="delimiter">: </span><span class="description">compute the cumulative sum of a one-dimensional double-precision floating-point ndarray.</span>
 -   <span class="signature">[`dsum( arrays )`][@stdlib/blas/ext/base/ndarray/dsum]</span><span class="delimiter">: </span><span class="description">compute the sum of all elements in a one-dimensional double-precision floating-point ndarray.</span>
 -   <span class="signature">[`gcusum( arrays )`][@stdlib/blas/ext/base/ndarray/gcusum]</span><span class="delimiter">: </span><span class="description">compute the cumulative sum of a one-dimensional ndarray.</span>
 -   <span class="signature">[`gsum( arrays )`][@stdlib/blas/ext/base/ndarray/gsum]</span><span class="delimiter">: </span><span class="description">compute the sum of all elements in a one-dimensional ndarray.</span>
 -   <span class="signature">[`scusum( arrays )`][@stdlib/blas/ext/base/ndarray/scusum]</span><span class="delimiter">: </span><span class="description">compute the cumulative sum of a one-dimensional single-precision floating-point ndarray.</span>
 -   <span class="signature">[`ssum( arrays )`][@stdlib/blas/ext/base/ndarray/ssum]</span><span class="delimiter">: </span><span class="description">compute the sum of all elements in a one-dimensional single-precision floating-point ndarray.</span>
+-   <span class="signature">[`zsum( arrays )`][@stdlib/blas/ext/base/ndarray/zsum]</span><span class="delimiter">: </span><span class="description">compute the sum of all elements in a one-dimensional double-precision complex floating-point ndarray.</span>
 
 </div>
 
@@ -93,6 +95,8 @@ console.log( objectKeys( ns ) );
 
 <!-- <toc-links> -->
 
+[@stdlib/blas/ext/base/ndarray/csum]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/ndarray/csum
+
 [@stdlib/blas/ext/base/ndarray/dcusum]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/ndarray/dcusum
 
 [@stdlib/blas/ext/base/ndarray/dsum]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/ndarray/dsum
@@ -105,6 +109,8 @@ console.log( objectKeys( ns ) );
 
 [@stdlib/blas/ext/base/ndarray/ssum]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/ndarray/ssum
 
+[@stdlib/blas/ext/base/ndarray/zsum]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/ndarray/zsum
+
 <!-- </toc-links> -->
 
 </section>

lib/node_modules/@stdlib/constants/float32/gamma-lanczos-g/README.md

@@ -41,7 +41,7 @@ var FLOAT32_GAMMA_LANCZOS_G = require( '@stdlib/constants/float32/gamma-lanczos-
 Arbitrary constant `g` to be used in [Lanczos approximation][lanczos-approximation] functions.
 
 ```javascript
-var bool = ( FLOAT32_GAMMA_LANCZOS_G === 10.900510787963867 );
+var bool = ( FLOAT32_GAMMA_LANCZOS_G === 1.42845618724823 );
 // returns true
 ```
 
@@ -59,7 +59,7 @@ var bool = ( FLOAT32_GAMMA_LANCZOS_G === 10.900510787963867 );
 var FLOAT32_GAMMA_LANCZOS_G = require( '@stdlib/constants/float32/gamma-lanczos-g' );
 
 console.log( FLOAT32_GAMMA_LANCZOS_G );
-// => 10.900510787963867
+// => 1.42845618724823
 ```
 
 </section>
@@ -126,7 +126,7 @@ Macro for the arbitrary constant `g` to be used in [Lanczos approximation][lancz
 
 ## References
 
--   Pugh, Glendon R. 2004. "An analysis of the Lanczos gamma approximation." PhD thesis, University of British Columbia. [&lt;https://web.viu.ca/pughg/phdThesis/phdThesis.pdf>][@pugh:2004a].
+-   Boost C++ Libraries. "The Lanczos Approximation." [https://beta.boost.org/doc/libs/1_82_0/libs/math/doc/html/math_toolkit/lanczos.html][@boost:lanczos].
 
 </section>
 
@@ -146,7 +146,7 @@ Macro for the arbitrary constant `g` to be used in [Lanczos approximation][lancz
 
 [lanczos-approximation]: https://en.wikipedia.org/wiki/Lanczos_approximation
 
-[@pugh:2004a]: https://web.viu.ca/pughg/phdThesis/phdThesis.pdf
+[@boost:lanczos]: https://beta.boost.org/doc/libs/1_82_0/libs/math/doc/html/math_toolkit/lanczos.html
 
 </section>
 

lib/node_modules/@stdlib/constants/float32/gamma-lanczos-g/docs/repl.txt

@@ -5,7 +5,7 @@
     Examples
     --------
     > {{alias}}
-    10.900510787963867
+    1.42845618724823
 
     See Also
     --------

lib/node_modules/@stdlib/constants/float32/gamma-lanczos-g/docs/types/index.d.ts

@@ -23,7 +23,7 @@
 *
 * @example
 * var g = FLOAT32_GAMMA_LANCZOS_G;
-* // returns 10.900510787963867
+* // returns 1.42845618724823
 */
 declare const FLOAT32_GAMMA_LANCZOS_G: number;
 

lib/node_modules/@stdlib/constants/float32/gamma-lanczos-g/examples/index.js

@@ -21,4 +21,4 @@
 var FLOAT32_GAMMA_LANCZOS_G = require( './../lib' );
 
 console.log( FLOAT32_GAMMA_LANCZOS_G );
-// => 10.900510787963867
+// => 1.42845618724823

lib/node_modules/@stdlib/constants/float32/gamma-lanczos-g/include/stdlib/constants/float32/gamma_lanczos_g.h

@@ -22,6 +22,6 @@
 /**
 * Macro for the arbitrary constant `g` to be used in Lanczos approximation functions.
 */
-#define STDLIB_CONSTANT_FLOAT32_GAMMA_LANCZOS_G 10.900511f
+#define STDLIB_CONSTANT_FLOAT32_GAMMA_LANCZOS_G 1.428456135094165802001953125f
 
 #endif // !STDLIB_CONSTANTS_FLOAT32_GAMMA_LANCZOS_G_H

lib/node_modules/@stdlib/constants/float32/gamma-lanczos-g/lib/index.js

@@ -26,7 +26,7 @@
 *
 * @example
 * var FLOAT32_GAMMA_LANCZOS_G = require( '@stdlib/constants/float32/gamma-lanczos-g' );
-* // returns 10.900510787963867
+* // returns 1.42845618724823
 */
 
 // MODULES //
@@ -41,10 +41,10 @@ var float64ToFloat32 = require( '@stdlib/number/float64/base/to-float32' );
 *
 * @constant
 * @type {number}
-* @default 10.900510787963867
+* @default 1.42845618724823
 * @see [Lanczos Approximation]{@link https://en.wikipedia.org/wiki/Lanczos_approximation}
 */
-var FLOAT32_GAMMA_LANCZOS_G = float64ToFloat32( 10.900511 );
+var FLOAT32_GAMMA_LANCZOS_G = float64ToFloat32( 1.428456135094165802001953125 );
 
 
 // EXPORTS //

lib/node_modules/@stdlib/constants/float32/gamma-lanczos-g/test/test.js

@@ -32,8 +32,8 @@ tape( 'main export is a number', function test( t ) {
 	t.end();
 });
 
-tape( 'the exported value equals `10.900510787963867`', function test( t ) {
-	var expected = 10.900510787963867;
+tape( 'the exported value equals `1.42845618724823`', function test( t ) {
+	var expected = 1.42845618724823;
 	t.strictEqual( FLOAT32_GAMMA_LANCZOS_G, expected, 'returns expected value' );
 	t.end();
 });

lib/node_modules/@stdlib/repl/code-blocks/data/data.csv

@@ -1850,7 +1850,7 @@ base.strided.dnannsumkbn2.ndarray,"var x = new Float64Array( [ 1.0, -2.0, NaN, 2
 base.strided.dnannsumors,"var x = new Float64Array( [ 1.0, -2.0, NaN, 2.0 ] );\nvar out = new Float64Array( 2 );\nbase.strided.dnannsumors( x.length, x, 1, out, 1 )\nx = new Float64Array( [ -2.0, 1.0, 1.0, -5.0, 2.0, -1.0, NaN, NaN ] );\nout = new Float64Array( 2 );\nbase.strided.dnannsumors( 4, x, 2, out, 1 )\nvar x0 = new Float64Array( [ 1.0, -2.0, 3.0, 2.0, 5.0, 1.0, NaN, NaN ] );\nvar x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 );\nout = new Float64Array( 2 );\nbase.strided.dnannsumors( 4, x1, 2, out, 1 )\n"
 base.strided.dnannsumors.ndarray,"var x = new Float64Array( [ 1.0, -2.0, NaN, 2.0 ] );\nvar out = new Float64Array( 2 );\nbase.strided.dnannsumors.ndarray( x.length, x, 1, 0, out, 1, 0 )\nvar x = new Float64Array( [ 1.0, -2.0, 3.0, 2.0, 5.0, 1.0, NaN, NaN ] );\nout = new Float64Array( 2 );\nbase.strided.dnannsumors.ndarray( 4, x, 2, 1, out, 1, 0 )\n"
 base.strided.dnannsumpw,"var x = new Float64Array( [ 1.0, -2.0, NaN, 2.0 ] );\nvar out = new Float64Array( 2 );\nbase.strided.dnannsumpw( x.length, x, 1, out, 1 )\nx = new Float64Array( [ -2.0, 1.0, 1.0, -5.0, 2.0, -1.0, NaN, NaN ] );\nout = new Float64Array( 2 );\nbase.strided.dnannsumpw( 4, x, 2, out, 1 )\nvar x0 = new Float64Array( [ 1.0, -2.0, 3.0, 2.0, 5.0, 1.0, NaN, NaN ] );\nvar x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 );\nout = new Float64Array( 2 );\nbase.strided.dnannsumpw( 4, x1, 2, out, 1 )\n"
-base.strided.dnannsumpw.ndarray,"var x = new Float64Array( [ 1.0, -2.0, NaN, 2.0 ] );\nvar out = new Float64Array( 2 );\nbase.strided.dnannsumpw.ndarray( x.length, x, 1, 0, out, 1, 0 )\nvar x = new Float64Array( [ 1.0, -2.0, 3.0, 2.0, 5.0, 1.0, NaN, NaN ] );\nout = new Float64Array( 2 );\nbase.strided.dnannsumpw.ndarray( 4, x, 2, 1, out, 1, 0 )\n"
+base.strided.dnannsumpw.ndarray,"var x = new Float64Array( [ 1.0, -2.0, NaN, 2.0 ] );\nvar out = new Float64Array( 2 );\nbase.strided.dnannsumpw.ndarray( x.length, x, 1, 0, out, 1, 0 )\nx = new Float64Array( [ 1.0, -2.0, 3.0, 2.0, 5.0, 1.0, NaN, NaN ] );\nout = new Float64Array( 2 );\nbase.strided.dnannsumpw.ndarray( 4, x, 2, 1, out, 1, 0 )\n"
 base.strided.dnanrange,"var x = new Float64Array( [ 1.0, -2.0, NaN, 2.0 ] );\nbase.strided.dnanrange( x.length, x, 1 )\nx = new Float64Array( [ -2.0, 1.0, 1.0, -5.0, 2.0, -1.0, NaN ] );\nbase.strided.dnanrange( 3, x, 2 )\nvar x0 = new Float64Array( [ 1.0, -2.0, 3.0, 2.0, 5.0, -1.0, NaN ] );\nvar x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 );\nbase.strided.dnanrange( 3, x1, 2 )\n"
 base.strided.dnanrange.ndarray,"var x = new Float64Array( [ 1.0, -2.0, 2.0, NaN ] );\nbase.strided.dnanrange.ndarray( x.length, x, 1, 0 )\nvar x = new Float64Array( [ 1.0, -2.0, 3.0, 2.0, 5.0, -1.0, NaN ] );\nbase.strided.dnanrange.ndarray( 3, x, 2, 1 )\n"
 base.strided.dnanstdev,"var x = new Float64Array( [ 1.0, -2.0, NaN, 2.0 ] );\nbase.strided.dnanstdev( x.length, 1, x, 1 )\nx = new Float64Array( [ -2.0, 1.0, 1.0, -5.0, 2.0, -1.0 ] );\nbase.strided.dnanstdev( 3, 1, x, 2 )\nvar x0 = new Float64Array( [ 1.0, -2.0, 3.0, 2.0, 5.0, -1.0 ] );\nvar x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 );\nbase.strided.dnanstdev( 3, 1, x1, 2 )\n"