Jacobian Hess element comparison is not symmetric

[vincentmacri]

Steps To Reproduce

K = GF(7)
Kx.<x> = FunctionField(K)
_.<v> = K[]
F.<y> = Kx.extension(v^2 - (x^6 + 3))
O = F.maximal_order()

D3 = O.ideal(x + 3, y + 5).place() + O.ideal (x + 4, y + 5).place() + O.ideal (x + 5, y + 5).place() + O.ideal(x + 6, y + 5).place()
D3_reduced = O.ideal(x+1, y+2).place() + O.ideal(x+2, y+2).place()

J = F.jacobian(model='hess')
D3_jacobian = J(D3)
D3_reduced_jacobian = J(D3_reduced)

print(D3_jacobian == D3_reduced_jacobian)  # True
print(D3_reduced_jacobian == D3_jacobian)  # False

Expected Behavior

The == operator should be symmetric. I believe in the example above the comparison result should be True.

Actual Behavior

The == operator is not symmetric.

Additional Information

I believe the relevant code here was added by @kwankyu in #35467.

Environment

- **OS**: Fedora 40
- **Sage Version**: 10.4, compiled with conda

Checklist

• I have searched the existing issues for a bug report that matches the one I want to file, without success.
• I have read the documentation and troubleshoot guide