Weierstrass special cases

paulmasson · Mar 14, 2024 · 00957bd · 00957bd
1 parent f0a4e40
commit 00957bd
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 0 deletions.
diff --git a/build/math.js b/build/math.js
@@ -1401,6 +1401,8 @@ function weierstrassP( x, g2, g3 ) {
 
   if ( !isComplex(x) ) x = complex(x);
 
+  if ( isZero(g2) && isZero(g3) ) return pow( x, -2 );
+
   var [ e1, e2, e3 ] = weierstrassRoots( g2, g3 );
 
   // Whittaker & Watson, Section 22.351
@@ -1415,6 +1417,8 @@ function weierstrassPPrime( x, g2, g3 ) {
 
   if ( !isComplex(x) ) x = complex(x);
 
+  if ( isZero(g2) && isZero(g3) ) return mul( -2, pow( x, -3 ) );
+
   var [ e1, e2, e3 ] = weierstrassRoots( g2, g3 );
 
   // Whittaker & Watson, Section 22.351

diff --git a/src/functions/elliptic-functions.js b/src/functions/elliptic-functions.js
@@ -445,6 +445,8 @@ function weierstrassP( x, g2, g3 ) {
 
   if ( !isComplex(x) ) x = complex(x);
 
+  if ( isZero(g2) && isZero(g3) ) return pow( x, -2 );
+
   var [ e1, e2, e3 ] = weierstrassRoots( g2, g3 );
 
   // Whittaker & Watson, Section 22.351
@@ -459,6 +461,8 @@ function weierstrassPPrime( x, g2, g3 ) {
 
   if ( !isComplex(x) ) x = complex(x);
 
+  if ( isZero(g2) && isZero(g3) ) return mul( -2, pow( x, -3 ) );
+
   var [ e1, e2, e3 ] = weierstrassRoots( g2, g3 );
 
   // Whittaker & Watson, Section 22.351