8000 feat: add C implementation for `math/base/special/ellipj` by anandkaranubc · Pull Request #6965 · stdlib-js/stdlib · GitHub
[go: up one dir, main page]

Skip to content

feat: add C implementation for math/base/special/ellipj #6965

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 10 commits into from
Jun 12, 2025
Prev Previous commit
Next Next commit
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
---
  • Loading branch information
anandkaranubc committed May 11, 2025
commit a0af388a3ffb9ca8e8f8296c7243f9d6c2e929db
7 changes: 7 additions & 0 deletions lib/node_modules/@stdlib/math/base/special/ellipj/src/main.c
Original file line number Diff line number Diff line change
Expand Up @@ -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 );
Expand All @@ -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 );
Expand All @@ -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...

Expand Down Expand Up @@ -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;
Expand Down Expand Up @@ -208,5 +214,6 @@ void stdlib_base_ellipj( const double u, const double m, double* sn, double* cn,
*dn = *cn / dnDenom;
}
}
return;
}
}
Loading
0