Dear Sage-devel,
I recently stumbled upon the following odd behavior in the new release of
sage. The issue does not occur in sage-7.4
Am I missing something ?
sage: var('a00, a10, a01, a11, u00, u10, u01, u11')
(a00, a10, a01, a11, u00, u10, u01, u11)
sage: f = (((a01*a10 - a00*a11)*u01*u10 - (a01*a10 - a00*a11)*u00*u11)*u10*(
1/u00 - u01*u10/(u00^2*(u01*u10/u00 - u11)))/(a10*u00*u01 - a00*u
....: 01*u10 + (a11*u00 - a01*u10)*u11) - a10 - ((a01*a10 -
a00*a11)*u01*u10 - (a01*a10 - a00*a11)*u00*u11)*u10*u11/((a10*u00*u01 +
a11*u01*u
....: 10 - (a00*u00 + a01*u10)*u11)*u00*(u01*u10/u00 - u11)))*(a01 -
((a01*a10 - a00*a11)*u01*u10 - (a01*a10 - a00*a11)*u00*u11)*u01/((a10*u0
....: 0*u01 - a00*u01*u10 + (a11*u00 - a01*u10)*u11)*(u01*u10/u00 - u11)) -
((a01*a10 - a00*a11)*u01*u10 - (a01*a10 - a00*a11)*u00*u11)*u01/(
....: (a10*u00*u01 + a11*u01*u10 - (a00*u00 + a01*u10)*u11)*(u01*u10/u00 -
u11))) - (((a01*a10 - a00*a11)*u01*u10 - (a01*a10 - a00*a11)*u00*u
....: 11)*u00*(1/u00 - u01*u10/(u00^2*(u01*u10/u00 - u11)))/(a10*u00*u01 -
a00*u01*u10 + (a11*u00 - a01*u10)*u11) - a00 - ((a01*a10 - a00*a11
....: )*u01*u10 - (a01*a10 - a00*a11)*u00*u11)*u01*u10/((a10*u00*u01 +
a11*u01*u10 - (a00*u00 + a01*u10)*u11)*u00*(u01*u10/u00 - u11)))*(a11
....: - ((a01*a10 - a00*a11)*u01*u10 - (a01*a10 -
a00*a11)*u00*u11)*u01*u10/((a10*u00*u01 - a00*u01*u10 + (a11*u00 -
a01*u10)*u11)*u00*(u01*u
....: 10/u00 - u11)) - ((a01*a10 - a00*a11)*u01*u10 - (a01*a10 -
a00*a11)*u00*u11)*u11/((a10*u00*u01 + a11*u01*u10 - (a00*u00 + a01*u10)*u11)
....: *(u01*u10/u00 - u11)))
sage: f.numerator()
n = 1
is_pyobject = 1
---------------------------------------------------------------------------
RuntimeError Traceback (most recent call last)
<ipython-input-3-d636c277ed2c> in <module>()
----> 1 f.numerator()
/Applications/sage-7.5/src/sage/symbolic/expression.pyx in
sage.symbolic.expression.Expression.numerator
(/Applications/sage-7.5/src/build/cythonized/sage/symbolic/expression.cpp:45612)
()
8581 cdef GEx oper, power
8582 if normalize:
-> 8583 return new_Expression_from_GEx(self._parent, self._gobj.
numer())
8584 elif is_a_mul(self._gobj):
8585 for i from 0 <= i < self._gobj.nops():
RuntimeError: num2canonical: can't happen
--
You received this message because you are subscribed to the Google Groups
"sage-devel" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at https://groups.google.com/group/sage-devel.
For more options, visit https://groups.google.com/d/optout.
