Vinzent Steinberg

bc.. Currently the interface is not stateless, see mp.dps for instance (or is this "high level"?). Are you planning to change this? Why do you prefer a stateless interface? The...

bc.. Something like "inexact=True" could be stored inside the result. This would be a better place IMHO. The high-level functions could check the result and raise exceptions if wanted, while...

bc.. The solver exists and works, your starting point just has to be close enough to the root. Depending on the system, "close" can mean +-0.001. It converges only locally,...

bc.. The result is not exactly the same, I suspect there is a typo or bug somewhere. Anyway, have fun! For the impatient: $ python sofa.py [...] step 5: x0:...

bc.. Got this when trying to start from (0., 0., 0., 0.): Traceback (most recent call last): File "sofa.py", line 17, in print msolve((A, B, phi, theta), system, (0., 0.,...

bc.. In fact, it's based on a numerical function. The only symbolic math is done to calculate the Jacobian, and this can be done numerically (it's even faster for really...

bc.. This example should work now in mpmath, using findroot(). **Status:** Started p. Original comment: "http://code.google.com/p/mpmath/issues/detail?id=56#c7":http://code.google.com/p/mpmath/issues/detail?id=56#c7 p. Original author: "https://code.google.com/u/[email protected]/":https://code.google.com/u/[email protected]/

bc.. http://www.dd.chalmers.se/~frejohl/code/libarith/ Why not add the rationals from there? **Status:** Started p. Original comment: "http://code.google.com/p/mpmath/issues/detail?id=53#c2":http://code.google.com/p/mpmath/issues/detail?id=53#c2 p. Original author: "https://code.google.com/u/[email protected]/":https://code.google.com/u/[email protected]/

bc.. Which interface issues are to be sorted out? sqrt(mpq) etc.? sqrt(4)? p. Original comment: "http://code.google.com/p/mpmath/issues/detail?id=53#c4":http://code.google.com/p/mpmath/issues/detail?id=53#c4 p. Original author: "https://code.google.com/u/[email protected]/":https://code.google.com/u/[email protected]/

bc.. I have no problems with sqrt(4) being an integer. I'm not sure about mpq simplifying to integer. Would it affect speed dramatically? Please note that matrix([[1]]) does not simplify...