[sage-support] Re: Possible new bug of solve()

2017-10-08 Thread Emmanuel Charpentier
This is now Trac#23993 , but deserves further discussion. I'll post something on sage-devel after thinking it out. HTH, -- Emmanuel Charpentier Le dimanche 8 octobre 2017 20:59:26 UTC+2, Emmanuel Charpentier a écrit : > > Yet another datapoint : o

[sage-support] Re: Possible new bug of solve()

2017-10-08 Thread Emmanuel Charpentier
This is now Trac#23992 . HTH, -- Emmanuel Charpentier Le dimanche 8 octobre 2017 20:53:55 UTC+2, Emmanuel Charpentier a écrit : > > I'd vote for a bug : since Maxima, used standalone, *can* solve this > system, not being able to use it correctly *

[sage-support] Re: Possible new bug of solve()

2017-10-08 Thread Emmanuel Charpentier
Yet another datapoint : our interface to mathematica *can* be used to get a solution. The problem I got was that the result can't be translated back to Sage : Mathematica express its results as "rules" (e.g "var -> expression" means "expression is a solution for var"). Manually translating thes

[sage-support] Re: Possible new bug of solve()

2017-10-08 Thread Emmanuel Charpentier
I'd vote for a bug : since Maxima, used standalone, *can* solve this system, not being able to use it correctly *is* a bug. Somewhere in our interface... Le dimanche 8 octobre 2017 14:26:00 UTC+2, Maarten Derickx a écrit : > > It deserves a ticket, the main question is more whether the ticket s

[sage-support] Re: Possible new bug of solve()

2017-10-08 Thread Maarten Derickx
It deserves a ticket, the main question is more whether the ticket should be a bug or an enhancement. But that is something that can still be changed later. -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and sto

[sage-support] Re: Possible new bug of solve()

2017-10-08 Thread Emmanuel Charpentier
Another datapoint : Sympy seems to be able to (slowly) solve this system, via : from sympy.solvers import solve as ssolve ssol=ssolve([e.rhs()-e.lhs() for e in [eq1,eq2,eq3]],[x,y,z]) The output can be converted to Sage via: Ssol=map(lambda S:map(lambda v,s:v==SR(repr(s)).simplify_full(), [x,y,