Menu
▾
▴
#2009 numerical evaluation of elliptic_ec fails for argument > 1
(%i60) elliptic_ec (1.7);
The assertion (AND (>= X 0) (>= Y 0) (>= Z 0) (PLUSP (+ X Y))
(PLUSP (+ X Z)) (PLUSP (+ Y Z))) failed.
which appears to originate in the function DRF.
Not sure what is the right answer here, but at any rate it shouldn't be a Lisp error.
Agreed. Curiously, elliptic_ec(1.7b0) returns a complex number, which appears to be correct.
The issue is that the algorithm for floats assumes that m <= 1, but the bfloat algorithm is correct for all values of m.
This will be fixed shortly.
Fixed in ellipt.lisp, rev 1.73.