chore: clean-up, fix tests, and add FIXMEs

The test tolerances for the C implementation are widely different
from the JavaScript implementation. I attempted to disable various
compiler optimizations. E.g.,

```
CFLAGS="-ffp-contract=on" make install-node-addons NODE_ADDONS_PATTERN="stats/base/dists/triangular/mgf"
```

Marked variables within the implementation as `volatile` in the hopes
of preventing compiler schenigans. Refactored as a single expression.
Replaced use of `exp` and `pow` with `math.h` counterparts. Etc.

Someone will need to dig deeper to figure out what is going on for
some pathological test cases.

---
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: passed
  - task: lint_javascript_benchmarks
    status: na
  - task: lint_python
    status: na
  - task: lint_r
    status: na
  - task: lint_c_src
    status: passed
  - task: lint_c_examples
    status: na
  - task: lint_c_benchmarks
    status: passed
  - 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
---

---
type: pre_push_report
description: Results of running various checks prior to pushing changes.
report:
  - task: run_javascript_examples
    status: na
  - task: run_c_examples
    status: na
  - task: run_cpp_examples
    status: na
  - task: run_javascript_readme_examples
    status: na
  - task: run_c_benchmarks
    status: passed
  - task: run_cpp_benchmarks
    status: na
  - task: run_fortran_benchmarks
    status: na
  - task: run_javascript_benchmarks
    status: na
  - task: run_julia_benchmarks
    status: na
  - task: run_python_benchmarks
    status: na
  - task: run_r_benchmarks
    status: na
  - task: run_javascript_tests
    status: passed
---

lib/node_modules/@stdlib/stats/base/dists/triangular/mgf/benchmark/c/benchmark.c

@@ -106,7 +106,7 @@ static double benchmark( void ) {
 		t[ i ] = random_uniform( 0.0, 5.0 );
 		a[ i ] = random_uniform( 0.0, 10.0 );
 		b[ i ] = random_uniform( a[ i ], 40.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-		c[ i ] = random_uniform( a[ i ], b[ i ]);
+		c[ i ] = random_uniform( a[ i ], b[ i ] );
 	}
 
 	tc = tic();

lib/node_modules/@stdlib/stats/base/dists/triangular/mgf/binding.gyp

@@ -167,4 +167,4 @@
       ], # end actions
     }, # end target copy_addon
   ], # end targets
-}
+}

lib/node_modules/@stdlib/stats/base/dists/triangular/mgf/include.gypi

@@ -50,4 +50,4 @@
       '<!@(node -e "var arr = require(\'@stdlib/utils/library-manifest\')(\'./manifest.json\',{},{\'basedir\':process.cwd(),\'paths\':\'posix\'}).libpath; for ( var i = 0; i < arr.length; i++ ) { console.log( arr[ i ] ); }")',
     ],
   }, # end variables
-}
+}

lib/node_modules/@stdlib/stats/base/dists/triangular/mgf/lib/factory.js

@@ -83,8 +83,8 @@ function factory( a, b, c ) {
 		if ( t === 0.0 ) {
 			return 1.0;
 		}
-		ret = (bmc * exp( a * t )) - (bma * exp( c * t ));
-		ret += cma * exp( b * t );
+		ret = ( bmc * exp( a*t ) ) - ( bma * exp( c*t ) );
+		ret += cma * exp( b*t );
 		ret *= 2.0;
 		ret /= bma * cma * bmc * pow( t, 2.0 );
 		return ret;

lib/node_modules/@stdlib/stats/base/dists/triangular/mgf/lib/main.js

@@ -94,8 +94,8 @@ function mgf( t, a, b, c ) {
 	bmc = b - c;
 	bma = b - a;
 	cma = c - a;
-	ret = (bmc * exp( a * t )) - (bma * exp( c * t ));
-	ret += cma * exp( b * t );
+	ret = ( bmc * exp( a*t ) ) - ( bma * exp( c*t ) );
+	ret += cma * exp( b*t );
 	ret *= 2.0;
 	ret /= bma * cma * bmc * pow( t, 2.0 );
 	return ret;

lib/node_modules/@stdlib/stats/base/dists/triangular/mgf/src/main.c

@@ -31,10 +31,14 @@
 * @return    evaluated MGF
 *
 * @example
-* double y = stdlib_base_dists_geometric_mgf( 0.5, -1.0, 1.0, 0.0 );
+* double y = stdlib_base_dists_triangular_mgf( 0.5, -1.0, 1.0, 0.0 );
 * // returns ~1.021
 */
 double stdlib_base_dists_triangular_mgf( const double t, const double a, const double b, const double c ) {
+	double bmc;
+	double bma;
+	double cma;
+	double ret;
 	if (
 		stdlib_base_is_nan( t ) ||
 		stdlib_base_is_nan( a ) ||
@@ -45,14 +49,14 @@ double stdlib_base_dists_triangular_mgf( const double t, const double a, const d
 	) {
 		return 0.0/0.0; // NaN
 	}
-	if( t == 0 ){
+	if ( t == 0.0 ) {
 		return 1.0;
 	}
-	double bmc = b - c;
-	double bma = b - a;
-	double cma = c - a;
-	double ret = (bmc * stdlib_base_exp( a * t )) - (bma * stdlib_base_exp( c * t ));
-	ret += cma * stdlib_base_exp( b * t );
+	bmc = b - c;
+	bma = b - a;
+	cma = c - a;
+	ret = ( bmc * stdlib_base_exp( a*t ) ) - ( bma * stdlib_base_exp( c*t ) );
+	ret += cma * stdlib_base_exp( b*t );
 	ret *= 2.0;
 	ret /= bma * cma * bmc * stdlib_base_pow( t, 2.0 );
 	return ret;

lib/node_modules/@stdlib/stats/base/dists/triangular/mgf/test/test.native.js

@@ -99,11 +99,13 @@ tape( 'the function evaluates the MGF for `x` given a small range `b - a`', opts
 	c = smallRange.c;
 	for ( i = 0; i < x.length; i++ ) {
 		y = mgf( x[i], a[i], b[i], c[i] );
-		if ( y.toFixed( 6 ) === expected[i].toFixed( 6 ) ) {
-			t.equal( y.toFixed( 6 ), expected[i].toFixed( 6 ), 'x: '+x[i]+', a: '+a[i]+', b: '+b[i]+', c: '+c[i]+', y: '+y+', expected: '+expected[i] );
+		if ( y === expected[i] ) {
+			t.equal( y, expected[i], 'x: '+x[i]+', a: '+a[i]+', b: '+b[i]+', c: '+c[i]+', y: '+y+', expected: '+expected[i] );
 		} else {
 			delta = abs( y - expected[ i ] );
-			tol = 100.0 * EPS * abs( expected[ i ] );
+
+			// FIXME: the tolerance for the C implementation is widely different from the JavaScript implementation, and it is not clear why.
+			tol = 1.0e7 * EPS * abs( expected[ i ] );
 			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. a: '+a[i]+'. b: '+b[i]+'. c: '+c[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
 		}
 	}
@@ -128,11 +130,13 @@ tape( 'the function evaluates the MGF for `x` given a medium range `b - a`', opt
 	c = mediumRange.c;
 	for ( i = 0; i < x.length; i++ ) {
 		y = mgf( x[i], a[i], b[i], c[i] );
-		if ( y.toFixed( 6 ) === expected[i].toFixed( 6 ) ) {
-			t.equal( y.toFixed( 6 ), expected[i].toFixed( 6 ), 'x: '+x[i]+', a: '+a[i]+', b: '+b[i]+', c: '+c[i]+', y: '+y+', expected: '+expected[i] );
+		if ( y === expected[i] ) {
+			t.equal( y, expected[i], 'x: '+x[i]+', a: '+a[i]+', b: '+b[i]+', c: '+c[i]+', y: '+y+', expected: '+expected[i] );
 		} else {
 			delta = abs( y - expected[ i ] );
-			tol = 100.0 * EPS * abs( expected[ i ] );
+
+			// FIXME: the tolerance for the C implementation is widely different from the JavaScript implementation, and it is not clear why.
+			tol = 1.0e7 * EPS * abs( expected[ i ] );
 			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. a: '+a[i]+'. b: '+b[i]+'. c: '+c[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
 		}
 	}
@@ -157,11 +161,13 @@ tape( 'the function evaluates the MGF for `x` given a large range `b - a`', opts
 	c = largeRange.c;
 	for ( i = 0; i < x.length; i++ ) {
 		y = mgf( x[i], a[i], b[i], c[i] );
-		if ( y.toFixed( 6 ) === expected[i].toFixed( 6 ) ) {
-			t.equal( y.toFixed( 6 ), expected[i].toFixed( 6 ), 'x: '+x[i]+', a: '+a[i]+', b: '+b[i]+', c: '+c[i]+', y: '+y+', expected: '+expected[i] );
+		if ( y === expected[i] ) {
+			t.equal( y, expected[i], 'x: '+x[i]+', a: '+a[i]+', b: '+b[i]+', c: '+c[i]+', y: '+y+', expected: '+expected[i] );
 		} else {
 			delta = abs( y - expected[ i ] );
-			tol = 100.0 * EPS * abs( expected[ i ] );
+
+			// FIXME: the tolerance for the C implementation is widely different from the JavaScript implementation, and it is not clear why.
+			tol = 1.0e7 * EPS * abs( expected[ i ] );
 			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. a: '+a[i]+'. b: '+b[i]+'. c: '+c[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
 		}
 	}