Does sympy give erroneous simplifications?

user_1574054 0
Viewed 47

I have the following calculation for sympy:

import sympy

q, r = sympy.symbols("q r")
equation = (((-q + r) - (q - r)) <= 0).simplify()
print(equation) # q >= r

equation = ((sympy.sqrt(2) * (-q + r) - sympy.sqrt(2) * (q - r)) <= 0).simplify()
print(equation) # q <= r

I don't see why the results should differ. What am I missing?

Edit

I am using version 1.5.1 of sympy and can see this on Python 3.6.6 and 3.7.7.

sympy python
1 Answers

A fix for this is given here. It looks like gcd was assumed to behave like igcd (which gives a nonnegative value). But when dealing with non-integers, gcd currently can give a negative result, thus the error. So SymPy will either modify gcd and the simplify code will work or the simplification routine must account for the sign of the extracted gcd.

user_10891610
Related