chore: add return statements

---
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: na
  - task: lint_javascript_cli
    status: na
  - task: lint_javascript_examples
    status: na
  - task: lint_javascript_tests
    status: na
  - 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: na
  - 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
---

Author: anandkaranubc

lib/node_modules/@stdlib/math/base/special/ellipj/src/main.c

@@ -101,6 +101,7 @@ void stdlib_base_ellipj( const double u, const double m, double* sn, double* cn,
 		*cn = out[ 1 ] / out[ 2 ];
 		*dn = 1.0 / out[ 2 ];
 		*am = 0.0 / 0.0; // NaN
+		return;
 	} else if ( m > 1.0 ) {
 		// A&S 16.11.1 for reciprocal parameter, mapping m > 1 to 0 < mu < 1:
 		k = stdlib_base_sqrt( m );
@@ -109,19 +110,22 @@ void stdlib_base_ellipj( const double u, const double m, double* sn, double* cn,
 		*cn = out[ 2 ];
 		*dn = out[ 1 ];
 		*am = 0.0 / 0.0; // NaN
+		return;
 	} else if ( m == 0.0 ) {
 		// A&S table 16.6, limiting case m = 0: circular trigonometric functions:
 		stdlib_base_sincos( u, sc, sc+1 );
 		*sn = sc[ 0 ];
 		*cn = sc[ 1 ];
 		*dn = 1.0;
 		*am = u;
+		return;
 	} else if ( m == 1.0 ) {
 		// A&S table 16.6: limiting case m = 1: hyperbolic functions:
 		*sn = stdlib_base_tanh( u );
 		*cn = 1.0 / stdlib_base_cosh( u );
 		*dn = *cn;
 		*am = stdlib_base_atan( stdlib_base_sinh( u ) );
+		return;
 	} else if ( m < STDLIB_CONSTANT_FLOAT64_SQRT_EPS ) {
 		// A&S 16.13.1 for small-u approximations. Additionally, compute K at extra cost in order to ensure returned values are correct outside the range [0, K]:
 		K4 = 4.0 * stdlib_base_ellipk( m );
@@ -134,6 +138,7 @@ void stdlib_base_ellipj( const double u, const double m, double* sn, double* cn,
 		*cn = c + ( f * s );
 		*dn = 1.0 - ( 0.5 * m * s * s );
 		*am = u - ( 0.25 * m * ( u - ( s * c ) ) );
+		return;
 	} else if ( m > 1.0 - STDLIB_CONSTANT_FLOAT64_SQRT_EPS ) {
 		// A&S 16.15.1 - 16.15.4 for near-unity approximations. Additionally, compute K at extra cost so that we may reflect as needed to ensure the returned values are correct...
 
@@ -162,6 +167,7 @@ void stdlib_base_ellipj( const double u, const double m, double* sn, double* cn,
 			*sn = -*sn;
 			*cn = -*cn;
 		}
+		return;
 	} else {
 		// A&S 16.4.1. Compute using the arithmetic-geometric mean...
 		a = 1.0;
@@ -208,5 +214,6 @@ void stdlib_base_ellipj( const double u, const double m, double* sn, double* cn,
 				*dn = *cn / dnDenom;
 			}
 		}
+		return;
 	}
 }