I'm sorry that the function in question is rather messy, but the problem is rather specific. I'm setting a constant 'c' and evaluating the local maximum of a function like so:
c = 1E6 find_local_maximum(lambda f: abs(1000*(2*pi*f*i)/(0.1*(2*pi*f*i)^2+1000*(2*pi*f*i)+c)), 1000, 3000)
In this case, I get the expected output. If I set
c = 1E7, or a nearby value, like
c = 1.123E7 or
c = 3E7, I get the following error:
RuntimeError: ECL says: THROW: The catch MACSYMA-QUIT is undefined.
Now, if I set
c = 5E7or
c = 1E8, I get no error.
It's also worth noting that
find_local_minimum with the exact same set of values and function does not have this same issue. Finally, the type of c is
sage.rings.real_mpfr.RealLiteral in each case.
What could be causing this problem? Is there a workaround or a potential bug fix?