Thanks Marc for your answer! It would be very useful to have explicitely
mentioned in the documentation of these functions that the result is not
constant over time (ie explain what consistency the user is supposed to
expect). Looking at the name "is_ordinary", I do not suspect it to be
depende
The function ps_modsym_from_simple_modsym_space calls internally a
randomized method from modular symbols (dual_eigenvector). Here is a
snippet of code that also returns different answers depending on the seed
(on a fresh session, since dual_eigenvector gets cached).
f = Newforms(32, 8, names='
