Merge branch 'develop' into feat/heavisidef

Author: anandkaranubc

lib/node_modules/@stdlib/array/typed/test/test.js

@@ -809,6 +809,8 @@ tape( 'the function returns a typed array (default, complex typed array)', funct
 	var out;
 
 	arr = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0 ] );
+
+	// Note: this is effectively a reinterpretation of the complex number array as an array of interleaved real and imaginary components
 	out = typedarray( arr );
 	t.strictEqual( instanceOf( out, Float64Array ), true, 'returns expected value' );
 	t.strictEqual( out.length, 4, 'returns expected value' );
@@ -828,6 +830,8 @@ tape( 'the function returns a typed array (dtype=float64, complex typed array)',
 	var out;
 
 	arr = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0 ] );
+
+	// Note: this is effectively a reinterpretation of the complex number array as an array of interleaved real and imaginary components
 	out = typedarray( arr, 'float64' );
 	t.strictEqual( instanceOf( out, Float64Array ), true, 'returns expected value' );
 	t.strictEqual( out.length, 4, 'returns expected value' );
@@ -847,6 +851,8 @@ tape( 'the function returns a typed array (dtype=float32, complex typed array)',
 	var out;
 
 	arr = new Complex64Array( [ 1.0, 2.0, 3.0, 4.0 ] );
+
+	// Note: this is effectively a reinterpretation of the complex number array as an array of interleaved real and imaginary components
 	out = typedarray( arr, 'float32' );
 	t.strictEqual( instanceOf( out, Float32Array ), true, 'returns expected value' );
 	t.strictEqual( out.length, 4, 'returns expected value' );
@@ -908,6 +914,8 @@ tape( 'the function returns a typed array (dtype=int32, complex typed array)', f
 	var out;
 
 	arr = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0 ] );
+
+	// Note: this is effectively a reinterpretation of the complex number array as an array of interleaved real and imaginary components
 	out = typedarray( arr, 'int32' );
 	t.strictEqual( instanceOf( out, Int32Array ), true, 'returns expected value' );
 	t.strictEqual( out.length, 4, 'returns expected value' );
@@ -927,6 +935,8 @@ tape( 'the function returns a typed array (dtype=int16, complex typed array)', f
 	var out;
 
 	arr = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0 ] );
+
+	// Note: this is effectively a reinterpretation of the complex number array as an array of interleaved real and imaginary components
 	out = typedarray( arr, 'int16' );
 	t.strictEqual( instanceOf( out, Int16Array ), true, 'returns expected value' );
 	t.strictEqual( out.length, 4, 'returns expected value' );
@@ -946,6 +956,8 @@ tape( 'the function returns a typed array (dtype=int8, complex typed array)', fu
 	var out;
 
 	arr = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0 ] );
+
+	// Note: this is effectively a reinterpretation of the complex number array as an array of interleaved real and imaginary components
 	out = typedarray( arr, 'int8' );
 	t.strictEqual( instanceOf( out, Int8Array ), true, 'returns expected value' );
 	t.strictEqual( out.length, 4, 'returns expected value' );
@@ -965,6 +977,8 @@ tape( 'the function returns a typed array (dtype=uint32, complex typed array)',
 	var out;
 
 	arr = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0 ] );
+
+	// Note: this is effectively a reinterpretation of the complex number array as an array of interleaved real and imaginary components
 	out = typedarray( arr, 'uint32' );
 	t.strictEqual( instanceOf( out, Uint32Array ), true, 'returns expected value' );
 	t.strictEqual( out.length, 4, 'returns expected value' );
@@ -984,6 +998,8 @@ tape( 'the function returns a typed array (dtype=uint16, complex typed array)',
 	var out;
 
 	arr = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0 ] );
+
+	// Note: this is effectively a reinterpretation of the complex number array as an array of interleaved real and imaginary components
 	out = typedarray( arr, 'uint16' );
 	t.strictEqual( instanceOf( out, Uint16Array ), true, 'returns expected value' );
 	t.strictEqual( out.length, 4, 'returns expected value' );
@@ -1003,6 +1019,8 @@ tape( 'the function returns a typed array (dtype=uint8, complex typed array)', f
 	var out;
 
 	arr = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0 ] );
+
+	// Note: this is effectively a reinterpretation of the complex number array as an array of interleaved real and imaginary components
 	out = typedarray( arr, 'uint8' );
 	t.strictEqual( instanceOf( out, Uint8Array ), true, 'returns expected value' );
 	t.strictEqual( out.length, 4, 'returns expected value' );
@@ -1022,6 +1040,8 @@ tape( 'the function returns a typed array (dtype=uint8c, complex typed array)',
 	var out;
 
 	arr = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0 ] );
+
+	// Note: this is effectively a reinterpretation of the complex number array as an array of interleaved real and imaginary components
 	out = typedarray( arr, 'uint8c' );
 	t.strictEqual( instanceOf( out, Uint8ClampedArray ), true, 'returns expected value' );
 	t.strictEqual( out.length, 4, 'returns expected value' );

lib/node_modules/@stdlib/math/base/special/gammasgn/README.md

@@ -26,17 +26,12 @@ limitations under the License.
 
 The sign of the [gamma-function][@stdlib/math/base/special/gamma] is defined as
 
-<!-- <equation class="equation" label="eq:gamma_sign_function" align="center" raw="\operatorname{gammasgn} ( x ) = \begin{cases} 1 & \textrm{if}\ \Gamma > 1 \\ -1 & \textrm{if}\ \Gamma < 1 \\ 0 & \textrm{otherwise}\ \end{cases}" alt="Sign of the gamma function"> -->
+<!-- <equation class="equation" label="eq:gamma_sign_function" align="center" raw="\operatorname{gammasgn} ( x ) = \begin{cases} 1 & \textrm{if}\ \Gamma > 0 \\ -1 & \textrm{if}\ \Gamma < 0 \\ 0 & \textrm{otherwise}\ \end{cases}" alt="Sign of the gamma function"> -->
 
 ```math
-\mathop{\mathrm{gammasgn}} ( x ) = \begin{cases} 1 & \textrm{if}\ \Gamma > 1 \\ -1 & \textrm{if}\ \Gamma < 1 \\ 0 & \textrm{otherwise}\ \end{cases}
+\mathop{\mathrm{gammasgn}} ( x ) = \begin{cases} 1 & \textrm{if}\ \Gamma > 0 \\ -1 & \textrm{if}\ \Gamma < 0 \\ 0 & \textrm{otherwise}\ \end{cases}
 ```
 
-<!-- <div class="equation" align="center" data-raw-text="\operatorname{gammasgn} ( x ) = \begin{cases} 1 &amp; \textrm{if}\ \Gamma &gt; 1 \\ -1 &amp; \textrm{if}\ \Gamma &lt; 1 \\ 0 &amp; \textrm{otherwise}\ \end{cases}" data-equation="eq:gamma_sign_function">
-    <img src="https://cdn.jsdelivr.net/gh/stdlib-js/stdlib@50b141156b147529227e2eb4247eda81c781dec9/lib/node_modules/@stdlib/math/base/special/gammasgn/docs/img/equation_gamma_sign_function.svg" alt="Sign of the gamma function">
-    <br>
-</div> -->
-
 <!-- </equation> -->
 
 The [gamma function][@stdlib/math/base/special/gamma] can be computed as the product of `gammasgn(x)` and `exp(gammaln(x))`.
@@ -185,7 +180,7 @@ int main( void ) {
     for ( i = 0; i < 100; i++ ) {
         x = ( (double)rand() / (double)RAND_MAX ) * 100.0;
         v = stdlib_base_gammasgn( x );
-        printf( "gammasgn%lf = %lf\n", x, v );
+        printf( "gammasgn(%lf) = %lf\n", x, v );
     }
 }
 ```

lib/node_modules/@stdlib/math/base/special/gammasgn/benchmark/benchmark.js

@@ -34,7 +34,7 @@ bench( pkg, function benchmark( b ) {
 	var y;
 	var i;
 
-	x = uniform( 100, 0.0, 171.0 );
+	x = uniform( 100, -10.0, 10.0 );
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {

lib/node_modules/@stdlib/math/base/special/gammasgn/benchmark/benchmark.native.js

@@ -43,7 +43,7 @@ bench( pkg+'::native', opts, function benchmark( b ) {
 	var y;
 	var i;
 
-	x = uniform( 100, 0.0, 171.0 );
+	x = uniform( 100, -10.0, 10.0 );
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {

lib/node_modules/@stdlib/math/base/special/gammasgn/benchmark/c/native/benchmark.c

@@ -97,7 +97,7 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		x[ i ] = ( 171.0 * rand_double() ) - 0.0;
+		x[ i ] = ( 100.0 * rand_double() ) - 50.0;
 	}
 
 	t = tic();

lib/node_modules/@stdlib/math/base/special/gammasgn/benchmark/python/scipy/benchmark.py

@@ -73,7 +73,7 @@ def print_results(elapsed):
 def benchmark():
     """Run the benchmark and print benchmark results."""
     setup = "from scipy.special import gammasgn; from random import random;"
-    stmt = "y = gammasgn(171.0*random() - 0.0)"
+    stmt = "y = gammasgn(100.0*random() - 50.0)"
 
     t = timeit.Timer(stmt, setup=setup)