simplify_full() aborts on a large expression
I have a large expression and executed simplify_full() on it. It aborts with some error, but I don't understand why. I tried to find out the reason by just shortening the expression until it works. The following fails:
var('C1 L1 L2 k w sgn2 sgn3')
expr = (1/2*sqrt(1/3)*sqrt((3*(3*C1^2*L1 - C1^2*L2)*k^2*w^3 - sqrt(3)*C1*sgn3*w*sqrt(-(48*2^(1/3)*k^2 - 4*(2^(1/3)*C1^4*L1^2*L2^2*k^4*w^7 - 2*2^(1/3)*C1^2*L1*L2*k^2*w^3 + 6*(2^(1/3)*C1*L1 + 2^(1/3)*C1*L2)*k^2*w)*w - 2*2^(2/3)*(9*sqrt(256/3*k^2 - 1/3*(4*(C1^8*L1^5*L2^3 - 2*C1^8*L1^4*L2^4 + C1^8*L1^3*L2^5)*k^6*w^15 - 8*(C1^6*L1^4*L2^2 - 4*C1^6*L1^3*L2^3 + C1^6*L1^2*L2^4)*k^4*w^11 + 4*(9*C1^5*L1^4*L2 - 13*C1^5*L1^3*L2^2 - 13*C1^5*L1^2*L2^3 + 9*C1^5*L1*L2^4)*k^4*w^9 - 4*(C1^3*L1^3 - 41*C1^3*L1^2*L2 - 41*C1^3*L1*L2^2 + C1^3*L2^3)*k^2*w^5 - ((27*C1^4*L1^4 - 36*C1^4*L1^3*L2 + 2*C1^4*L1^2*L2^2 - 36*C1^4*L1*L2^3 + 27*C1^4*L2^4)*k^4 - 4*(C1^4*L1^3*L2 - 10*C1^4*L1^2*L2^2 + C1^4*L1*L2^3)*k^2)*w^7 + 8*(2*C1^2*L1*L2 - 5*(3*C1^2*L1^2 + 14*C1^2*L1*L2 + 3*C1^2*L2^2)*k^2)*w^3 + 16*(24*(C1*L1 + C1*L2)*k^2 - C1*L1 - C1*L2)*w)*w - 16/3)*k^2 + 36*k^2 + (2*C1^6*L1^3*L2^3*k^6*w^11 - 6*C1^4*L1^2*L2^2*k^4*w^7 + 18*(C1^3*L1^2*L2 + C1^3*L1*L2^2)*k^4*w^5 - 18*(C1*L1 + C1*L2)*k^2*w + 3*(2*C1^2*L1*L2*k^2 - 3*(3*C1^2*L1^2 - 2*C1^2*L1*L2 + 3*C1^2*L2^2)*k^4)*w^3)*w - 2)^(2/3) - 4*2^(1/3) - ((3*C1^4*L1^2 - 2*C1^4*L1*L2 + 3*C1^4*L2^2)*k^4*w^7 + 8*C1^2*k^2*w^3)*(9*sqrt(256/3*k^2 - 1/3*(4*(C1^8*L1^5*L2^3 - 2*C1^8*L1^4*L2^4 + C1^8*L1^3*L2^5)*k^6*w^15 - 8*(C1^6*L1^4*L2^2 - 4*C1^6*L1^3*L2^3 + C1^6*L1^2*L2^4)*k^4*w^11 + 4*(9*C1^5*L1^4*L2 - 13*C1^5*L1^3*L2^2 - 13*C1^5*L1^2*L2^3 + 9*C1^5*L1*L2^4)*k^4*w^9 - 4*(C1^3*L1^3 - 41*C1^3*L1^2*L2 - 41*C1^3*L1*L2^2 + C1^3*L2^3)*k^2*w^5 - ((27*C1^4*L1^4 - 36*C1^4*L1^3*L2 + 2*C1^4*L1^2*L2^2 - 36*C1^4*L1*L2^3 + 27*C1^4*L2^4)*k^4 - 4*(C1^4*L1^3*L2 - 10*C1^4*L1^2*L2^2 + C1^4*L1*L2^3)*k^2)*w^7 + 8*(2*C1^2*L1*L2 - 5*(3*C1^2*L1^2 + 14*C1^2*L1*L2 + 3*C1^2*L2^2)*k^2)*w^3 + 16*(24*(C1*L1 + C1*L2)*k^2 - C1*L1 - C1*L2)*w)*w - 16/3)*k^2 + 36*k^2 + (2*C1^6*L1^3*L2^3*k^6*w^11 - 6*C1^4*L1^2*L2^2*k^4*w^7 + 18*(C1^3*L1^2*L2 + C1^3*L1*L2^2)*k^4*w^5 - 18*(C1*L1 + C1*L2)*k^2*w + 3*(2*C1^2*L1*L2*k^2 - 3*(3*C1^2*L1^2 - 2*C1^2*L1*L2 + 3*C1^2*L2^2)*k^4)*w^3)*w - 2)^(1/3)/(C1^2*k^2*w^3))/(9*sqrt(256/3*k^2 - 1/3*(4*(C1^8*L1^5*L2^3 - 2*C1^8*L1^4*L2^4 + C1^8*L1^3*L2^5)*k^6*w^15 - 8*(C1^6*L1^4*L2^2 - 4*C1^6*L1^3*L2^3 + C1^6*L1^2*L2^4)*k^4*w^11 + 4*(9*C1^5*L1^4*L2 - 13*C1^5*L1^3*L2^2 - 13*C1^5*L1^2*L2^3 + 9*C1^5*L1*L2^4)*k^4*w^9 - 4*(C1^3*L1^3 - 41*C1^3*L1^2*L2 - 41*C1^3*L1*L2^2 + C1^3*L2^3)*k^2*w^5 - ((27*C1^4*L1^4 - 36*C1^4*L1^3*L2 + 2*C1^4*L1^2*L2^2 - 36*C1^4*L1*L2^3 + 27*C1^4*L2^4)*k^4 - 4*(C1^4*L1^3*L2 - 10*C1^4*L1^2*L2^2 + C1^4*L1*L2^3)*k^2)*w^7 + 8*(2*C1^2*L1*L2 - 5*(3*C1^2*L1^2 + 14*C1^2*L1*L2 + 3*C1^2*L2^2)*k^2)*w^3 + 16*(24*(C1*L1 + C1*L2)*k^2 - C1*L1 - C1*L2)*w)*w - 16/3)*k^2 + 36*k^2 + (2*C1^6*L1^3*L2^3*k^6*w^11 - 6*C1^4*L1^2*L2^2*k^4*w^7 + 18*(C1^3*L1^2*L2 + C1^3*L1*L2^2)*k^4*w^5 - 18*(C1*L1 + C1*L2)*k^2*w + 3*(2*C1^2*L1*L2*k^2 - 3*(3*C1^2*L1^2 - 2*C1^2*L1*L2 + 3*C1^2*L2^2)*k^4)*w^3)*w - 2)^(1/3))*abs(k) + C1^2*k^2*sgn2*sqrt(sqrt(3)*2^(2/3)*(sqrt(3)*(24*2^(2/3)*k^2 - 2*(2^(2/3)*C1^4*L1^2*L2^2*k^4*w^7 - 2*2^(2/3)*C1^2*L1*L2*k^2*w^3 + 6*(2^(2/3)*C1*L1 + 2^(2/3)*C1*L2)*k^2*w)*w - 2*2^(2/3) - 2*(9*sqrt(256/3*k^2 - 1/3*(4*(C1^8*L1^5*L2^3 - 2*C1^8*L1^4*L2^4 + C1^8*L1^3*L2^5)*k^6*w^15 - 8*(C1^6*L1^4*L2^2 - 4*C1^6*L1^3*L2^3 + C1^6*L1^2*L2^4)*k^4*w^11 + 4*(9*C1^5*L1^4*L2 - 13*C1^5*L1^3*L2^2 - 13*C1^5*L1^2*L2^3 + 9*C1^5*L1*L2^4)*k^4*w^9 - 4*(C1^3*L1^3 - 41*C1^3*L1^2*L2 - 41*C1^3*L1*L2^2 + C1^3*L2^3)*k^2*w^5 - ((27*C1^4*L1^4 - 36*C1^4*L1^3*L2 + 2*C1^4*L1^2*L2^2 - 36*C1^4*L1*L2^3 + 27*C1^4*L2^4)*k^4 - 4*(C1^4*L1^3*L2 - 10*C1^4*L1^2*L2^2 + C1^4*L1*L2^3)*k^2)*w^7 + 8*(2*C1^2*L1*L2 - 5*(3*C1^2*L1^2 + 14*C1^2*L1*L2 + 3*C1^2*L2^2)*k^2)*w^3 + 16*(24*(C1*L1 + C1*L2)*k^2 - C1*L1 - C1*L2)*w)*w - 16/3)*k^2 + 36*k^2 + (2*C1^6*L1^3*L2^3*k^6*w^11 - 6*C1^4*L1^2*L2^2*k^4*w^7 + 18*(C1^3*L1^2*L2 + C1^3*L1*L2^2)*k^4*w^5 - 18*(C1*L1 + C1*L2)*k^2*w + 3*(2*C1^2*L1*L2*k^2 - 3*(3*C1^2*L1^2 - 2*C1^2*L1*L2 + 3*C1^2*L2^2)*k^4)*w^3)*w - 2)^(2/3) + ((3*2^(1/3)*C1^4*L1^2 - 2*2^(1/3)*C1^4*L1*L2 + 3*2^(1/3)*C1^4*L2^2)*k^4*w^7 + 8*2^(1/3)*C1^2*k^2*w^3)*(9*sqrt(256/3*k^2 - 1/3*(4*(C1^8*L1^5*L2^3 - 2*C1^8*L1^4*L2^4 + C1^8*L1^3*L2^5)*k^6*w^15 - 8*(C1^6*L1^4*L2^2 - 4*C1^6*L1^3*L2^3 + C1^6*L1^2*L2^4)*k^4*w^11 + 4*(9*C1^5*L1^4*L2 - 13*C1^5*L1^3*L2^2 - 13*C1^5*L1^2*L2^3 + 9*C1^5*L1*L2^4)*k^4*w^9 - 4*(C1^3*L1^3 - 41*C1^3*L1^2*L2 - 41*C1^3*L1*L2^2 + C1^3*L2^3)*k^2*w^5 - ((27*C1^4*L1^4 - 36*C1^4*L1^3*L2 + 2*C1^4*L1^2*L2^2 - 36*C1^4*L1*L2^3 + 27*C1^4*L2^4)*k^4 - 4*(C1^4*L1^3*L2 - 10*C1^4*L1^2*L2^2 + C1^4*L1*L2^3)*k^2)*w^7 + 8*(2*C1^2*L1*L2 - 5*(3*C1^2*L1^2 + 14*C1^2*L1*L2 + 3*C1^2*L2^2)*k^2)*w^3 + 16*(24*(C1*L1 + C1*L2)*k^2 - C1*L1 - C1*L2)*w)*w - 16/3)*k^2 + 36*k^2 + (2*C1^6*L1^3*L2^3*k^6*w^11 - 6*C1^4*L1^2*L2^2*k^4*w^7 + 18*(C1^3*L1^2*L2 + C1^3*L1*L2^2)*k^4*w^5 - 18*(C1*L1 + C1*L2)*k^2*w + 3*(2*C1^2*L1*L2*k^2 - 3*(3*C1^2*L1^2 - 2*C1^2*L1*L2 + 3*C1^2*L2^2)*k^4)*w^3)*w - 2)^(1/3)/(C1^2*k^2*w^3))*C1*w^2*sqrt(-(48*2^(1/3)*k^2 - 4*(2^(1/3)*C1^4*L1^2*L2^2*k^4*w^7 - 2*2^(1/3)*C1^2*L1*L2*k^2*w^3 + 6*(2^(1/3)*C1*L1 + 2^(1/3)*C1*L2)*k^2*w)*w - 2*2^(2/3)*(9*sqrt(256/3*k^2 - 1/3*(4*(C1^8*L1^5*L2^3 - 2*C1^8*L1^4*L2^4 + C1^8*L1^3*L2^5)*k^6*w^15 - 8*(C1^6*L1^4*L2^2 - 4*C1^6*L1^3*L2^3 + C1^6*L1^2*L2^4)*k^4*w^11 + 4*(9*C1^5*L1^4*L2 - 13*C1^5*L1^3*L2^2 - 13*C1^5*L1^2*L2^3 + 9*C1^5*L1*L2^4)*k^4*w^9 - 4*(C1^3*L1^3 - 41*C1^3*L1^2*L2 - 41*C1^3*L1*L2^2 + C1^3*L2^3)*k^2*w^5 - ((27*C1^4*L1^4 - 36*C1^4*L1^3*L2 + 2*C1^4*L1^2*L2^2 - 36*C1^4*L1*L2^3 + 27*C1^4*L2^4)*k^4 - 4*(C1^4*L1^3*L2 - 10*C1^4*L1^2*L2^2 + C1^4*L1*L2^3)*k^2)*w^7 + 8*(2*C1^2*L1*L2 - 5*(3*C1^2*L1^2 + 14*C1^2*L1*L2 + 3*C1^2*L2^2)*k^2)*w^3 + 16*(24*(C1*L1 + C1*L2)*k^2 - C1*L1 - C1*L2)*w)*w - 16/3)*k^2 + 36*k^2 + (2*C1^6*L1^3*L2^3*k^6*w^11 - 6*C1^4*L1^2*L2^2*k^4*w^7 + 18*(C1^3*L1^2*L2 + C1^3*L1*L2^2)*k^4*w^5 - 18*(C1*L1 + C1*L2)*k^2*w + 3*(2*C1^2*L1*L2*k^2 - 3*(3*C1^2*L1^2 - 2*C1^2*L1*L2 + 3*C1^2*L2^2)*k^4)*w^3)*w - 2)^(2/3) - 4*2^(1/3) - ((3*C1^4*L1^2 - 2*C1^4*L1*L2 + 3*C1^4*L2^2)*k^4*w^7 + 8*C1^2*k^2*w^3)*(9*sqrt(256/3*k^2 - 1/3*(4*(C1^8*L1^5*L2^3 - 2*C1^8*L1^4*L2^4 + C1^8*L1^3*L2^5)*k^6*w^15 - 8*(C1^6*L1^4*L2^2 - 4*C1^6*L1^3*L2^3 + C1^6*L1^2*L2^4)*k^4*w^11 + 4*(9*C1^5*L1^4*L2 - 13*C1^5*L1^3*L2^2 - 13*C1^5*L1^2*L2^3 + 9*C1^5*L1*L2^4)*k^4*w^9 - 4*(C1^3*L1^3 - 41*C1^3*L1^2*L2 - 41*C1^3*L1*L2^2 + C1^3*L2^3)*k^2*w^5 - ((27*C1^4*L1^4 - 36*C1^4*L1^3*L2 + 2*C1^4*L1^2*L2^2 - 36*C1^4*L1*L2^3 + 27*C1^4*L2^4)*k^4 - 4*(C1^4*L1^3*L2 - 10*C1^4*L1^2*L2^2 + C1^4*L1*L2^3)*k^2)*w^7 + 8*(2*C1^2*L1*L2 - 5*(3*C1^2*L1^2 + 14*C1^2*L1*L2 + 3*C1^2*L2^2)*k^2)*w^3 + 16*(24*(C1*L1 + C1*L2)*k^2 - C1*L1 - C1*L2)*w)*w - 16/3)*k^2 + 36*k^2 + (2*C1^6*L1^3*L2^3*k^6*w^11 - 6*C1^4*L1^2*L2^2*k^4*w^7 + 18*(C1^3*L1^2*L2 + C1^3*L1*L2^2)*k^4*w^5 - 18*(C1*L1 + C1*L2)*k^2*w + 3*(2*C1^2*L1*L2*k^2 - 3*(3*C1^2*L1^2 - 2*C1^2*L1*L2 + 3*C1^2*L2^2)*k^4)*w^3)*w - 2)^(1/3)/(C1^2*k^2*w^3))/(9*sqrt(256/3*k^2 - 1/3*(4*(C1^8*L1^5*L2^3 - 2*C1^8*L1^4*L2^4 + C1^8*L1^3*L2^5)*k^6*w^15 - 8*(C1^6*L1^4*L2^2 - 4*C1^6*L1^3*L2^3 + C1^6*L1^2*L2^4)*k^4*w^11 + 4*(9*C1^5*L1^4*L2 - 13*C1^5*L1^3*L2^2 - 13*C1^5*L1^2*L2^3 + 9*C1^5*L1*L2^4)*k^4*w^9 - 4*(C1^3*L1^3 - 41*C1^3*L1^2*L2 - 41*C1^3*L1*L2^2 + C1^3*L2^3)*k^2*w^5 - ((27*C1^4*L1^4 - 36*C1^4*L1^3*L2 + 2*C1^4*L1^2*L2^2 - 36*C1^4*L1*L2^3 + 27*C1^4*L2^4)*k^4 - 4*(C1^4*L1^3*L2 - 10*C1^4*L1^2*L2^2 + C1^4*L1*L2^3)*k^2)*w^7 + 8*(2*C1^2*L1*L2 - 5*(3*C1^2*L1^2 + 14*C1^2*L1*L2 + 3*C1^2*L2^2)*k^2)*w^3 + 16*(24*(C1*L1 + C1*L2)*k^2 - C1*L1 - C1*L2)*w)*w - 16/3)*k^2 + 36*k^2 + (2*C1^6*L1^3*L2^3*k^6*w^11 - 6*C1^4*L1^2*L2^2*k^4*w^7 + 18*(C1^3*L1^2*L2 + C1^3*L1*L2^2)*k^4*w^5 - 18*(C1*L1 + C1*L2)*k^2*w + 3*(2*C1^2*L1*L2*k^2 - 3*(3*C1^2*L1^2 - 2*C1^2*L1*L2 + 3*C1^2*L2^2)*k^4)*w^3)*w - 2)^(1/3))*abs(k) + 9*2^(1/3))))))
expr.simplify_full()
The output is:
RuntimeError Traceback (most recent call last)
/opt/sagemath-9.3/local/lib/python3.7/site-packages/sage/interfaces/interface.py in __init__(self, parent, value, is_name, name)
717 try:
--> 718 self._name = parent._create(value, name=name)
719 except (TypeError, RuntimeError, ValueError) as x:
/opt/sagemath-9.3/local/lib/python3.7/site-packages/sage/interfaces/maxima_lib.py in _create(self, value, name)
603 else:
--> 604 self.set(name, value)
605 except RuntimeError as error:
/opt/sagemath-9.3/local/lib/python3.7/site-packages/sage/interfaces/maxima_lib.py in set(self, var, value)
511 cmd = '%s : %s$'%(var, value.rstrip(';'))
--> 512 self.eval(cmd)
513
/opt/sagemath-9.3/local/lib/python3.7/site-packages/sage/interfaces/maxima_lib.py in _eval_line(self, line, locals, reformat, **kwds)
457 if statement:
--> 458 maxima_eval("#$%s$" % statement)
459 if not reformat:
/opt/sagemath-9.3/local/lib/python3.7/site-packages/sage/libs/ecl.pyx in sage.libs.ecl.EclObject.__call__ (build/cythonized/sage/libs/ecl.c:8623)()
853 lispargs = EclObject(list(args))
--> 854 return ecl_wrap(ecl_safe_apply(self.obj,(<EclObject>lispargs).obj))
855
/opt/sagemath-9.3/local/lib/python3.7/site-packages/sage/libs/ecl.pyx in sage.libs.ecl.ecl_safe_apply (build/cythonized/sage/libs/ecl.c:5827)()
359
--> 360 ecl_sig_on()
361 ret = safe_cl_apply(&error,func,args)
RuntimeError: Aborted
During handling of the above exception, another exception occurred:
TypeError Traceback (most recent call last)
<ipython-input-1-e47b1b19c48f> in <module>
3 expr = (Integer(1)/Integer(2)*sqrt(Integer(1)/Integer(3))*sqrt((Integer(3)*(Integer(3)*C1**Integer(2)*L1 - C1**Integer(2)*L2)*k**Integer(2)*w**Integer(3) - sqrt(Integer(3))*C1*sgn3*w*sqrt(-(Integer(48)*Integer(2)**(Integer(1)/Integer(3))*k**Integer(2) - Integer(4)*(Integer(2)**(Integer(1)/Integer(3))*C1**Integer(4)*L1**Integer(2)*L2**Integer(2)*k**Integer(4)*w**Integer(7) - Integer(2)*Integer(2)**(Integer(1)/Integer(3))*C1**Integer(2)*L1*L2*k**Integer(2)*w**Integer(3) + Integer(6)*(Integer(2)**(Integer(1)/Integer(3))*C1*L1 + Integer(2)**(Integer(1)/Integer(3))*C1*L2)*k**Integer(2)*w)*w - Integer(2)*Integer(2)**(Integer(2)/Integer(3))*(Integer(9)*sqrt(Integer(256)/Integer(3)*k**Integer(2) - Integer(1)/Integer(3)*(Integer(4)*(C1**Integer(8)*L1**Integer(5)*L2**Integer(3) - Integer(2)*C1**Integer(8)*L1**Integer(4)*L2**Integer(4) + C1**Integer(8)*L1**Integer(3)*L2**Integer(5))*k**Integer(6)*w**Integer(15) - Integer(8)*(C1**Integer(6)*L1**Integer(4)*L2**Integer(2) - Integer(4)*C1**Integer(6)*L1**Integer(3)*L2**Integer(3) + C1**Integer(6)*L1**Integer(2)*L2**Integer(4))*k**Integer(4)*w**Integer(11) + Integer(4)*(Integer(9)*C1**Integer(5)*L1**Integer(4)*L2 - Integer(13)*C1**Integer(5)*L1**Integer(3)*L2**Integer(2) - Integer(13)*C1**Integer(5)*L1**Integer(2)*L2**Integer(3) + Integer(9)*C1**Integer(5)*L1*L2**Integer(4))*k**Integer(4)*w**Integer(9) - Integer(4)*(C1**Integer(3)*L1**Integer(3) - Integer(41)*C1**Integer(3)*L1**Integer(2)*L2 - Integer(41)*C1**Integer(3)*L1*L2**Integer(2) + C1**Integer(3)*L2**Integer(3))*k**Integer(2)*w**Integer(5) - ((Integer(27)*C1**Integer(4)*L1**Integer(4) - Integer(36)*C1**Integer(4)*L1**Integer(3)*L2 + Integer(2)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2) - Integer(36)*C1**Integer(4)*L1*L2**Integer(3) + Integer(27)*C1**Integer(4)*L2**Integer(4))*k**Integer(4) - Integer(4)*(C1**Integer(4)*L1**Integer(3)*L2 - Integer(10)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2) + C1**Integer(4)*L1*L2**Integer(3))*k**Integer(2))*w**Integer(7) + Integer(8)*(Integer(2)*C1**Integer(2)*L1*L2 - Integer(5)*(Integer(3)*C1**Integer(2)*L1**Integer(2) + Integer(14)*C1**Integer(2)*L1*L2 + Integer(3)*C1**Integer(2)*L2**Integer(2))*k**Integer(2))*w**Integer(3) + Integer(16)*(Integer(24)*(C1*L1 + C1*L2)*k**Integer(2) - C1*L1 - C1*L2)*w)*w - Integer(16)/Integer(3))*k**Integer(2) + Integer(36)*k**Integer(2) + (Integer(2)*C1**Integer(6)*L1**Integer(3)*L2**Integer(3)*k**Integer(6)*w**Integer(11) - Integer(6)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2)*k**Integer(4)*w**Integer(7) + Integer(18)*(C1**Integer(3)*L1**Integer(2)*L2 + C1**Integer(3)*L1*L2**Integer(2))*k**Integer(4)*w**Integer(5) - Integer(18)*(C1*L1 + C1*L2)*k**Integer(2)*w + Integer(3)*(Integer(2)*C1**Integer(2)*L1*L2*k**Integer(2) - Integer(3)*(Integer(3)*C1**Integer(2)*L1**Integer(2) - Integer(2)*C1**Integer(2)*L1*L2 + Integer(3)*C1**Integer(2)*L2**Integer(2))*k**Integer(4))*w**Integer(3))*w - Integer(2))**(Integer(2)/Integer(3)) - Integer(4)*Integer(2)**(Integer(1)/Integer(3)) - ((Integer(3)*C1**Integer(4)*L1**Integer(2) - Integer(2)*C1**Integer(4)*L1*L2 + Integer(3)*C1**Integer(4)*L2**Integer(2))*k**Integer(4)*w**Integer(7) + Integer(8)*C1**Integer(2)*k**Integer(2)*w**Integer(3))*(Integer(9)*sqrt(Integer(256)/Integer(3)*k**Integer(2) - Integer(1)/Integer(3)*(Integer(4)*(C1**Integer(8)*L1**Integer(5)*L2**Integer(3) - Integer(2)*C1**Integer(8)*L1**Integer(4)*L2**Integer(4) + C1**Integer(8)*L1**Integer(3)*L2**Integer(5))*k**Integer(6)*w**Integer(15) - Integer(8)*(C1**Integer(6)*L1**Integer(4)*L2**Integer(2) - Integer(4)*C1**Integer(6)*L1**Integer(3)*L2**Integer(3) + C1**Integer(6)*L1**Integer(2)*L2**Integer(4))*k**Integer(4)*w**Integer(11) + Integer(4)*(Integer(9)*C1**Integer(5)*L1**Integer(4)*L2 - Integer(13)*C1**Integer(5)*L1**Integer(3)*L2**Integer(2) - Integer(13)*C1**Integer(5)*L1**Integer(2)*L2**Integer(3) + Integer(9)*C1**Integer(5)*L1*L2**Integer(4))*k**Integer(4)*w**Integer(9) - Integer(4)*(C1**Integer(3)*L1**Integer(3) - Integer(41)*C1**Integer(3)*L1**Integer(2)*L2 - Integer(41)*C1**Integer(3)*L1*L2**Integer(2) + C1**Integer(3)*L2**Integer(3))*k**Integer(2)*w**Integer(5) - ((Integer(27)*C1**Integer(4)*L1**Integer(4) - Integer(36)*C1**Integer(4)*L1**Integer(3)*L2 + Integer(2)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2) - Integer(36)*C1**Integer(4)*L1*L2**Integer(3) + Integer(27)*C1**Integer(4)*L2**Integer(4))*k**Integer(4) - Integer(4)*(C1**Integer(4)*L1**Integer(3)*L2 - Integer(10)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2) + C1**Integer(4)*L1*L2**Integer(3))*k**Integer(2))*w**Integer(7) + Integer(8)*(Integer(2)*C1**Integer(2)*L1*L2 - Integer(5)*(Integer(3)*C1**Integer(2)*L1**Integer(2) + Integer(14)*C1**Integer(2)*L1*L2 + Integer(3)*C1**Integer(2)*L2**Integer(2))*k**Integer(2))*w**Integer(3) + Integer(16)*(Integer(24)*(C1*L1 + C1*L2)*k**Integer(2) - C1*L1 - C1*L2)*w)*w - Integer(16)/Integer(3))*k**Integer(2) + Integer(36)*k**Integer(2) + (Integer(2)*C1**Integer(6)*L1**Integer(3)*L2**Integer(3)*k**Integer(6)*w**Integer(11) - Integer(6)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2)*k**Integer(4)*w**Integer(7) + Integer(18)*(C1**Integer(3)*L1**Integer(2)*L2 + C1**Integer(3)*L1*L2**Integer(2))*k**Integer(4)*w**Integer(5) - Integer(18)*(C1*L1 + C1*L2)*k**Integer(2)*w + Integer(3)*(Integer(2)*C1**Integer(2)*L1*L2*k**Integer(2) - Integer(3)*(Integer(3)*C1**Integer(2)*L1**Integer(2) - Integer(2)*C1**Integer(2)*L1*L2 + Integer(3)*C1**Integer(2)*L2**Integer(2))*k**Integer(4))*w**Integer(3))*w - Integer(2))**(Integer(1)/Integer(3))/(C1**Integer(2)*k**Integer(2)*w**Integer(3)))/(Integer(9)*sqrt(Integer(256)/Integer(3)*k**Integer(2) - Integer(1)/Integer(3)*(Integer(4)*(C1**Integer(8)*L1**Integer(5)*L2**Integer(3) - Integer(2)*C1**Integer(8)*L1**Integer(4)*L2**Integer(4) + C1**Integer(8)*L1**Integer(3)*L2**Integer(5))*k**Integer(6)*w**Integer(15) - Integer(8)*(C1**Integer(6)*L1**Integer(4)*L2**Integer(2) - Integer(4)*C1**Integer(6)*L1**Integer(3)*L2**Integer(3) + C1**Integer(6)*L1**Integer(2)*L2**Integer(4))*k**Integer(4)*w**Integer(11) + Integer(4)*(Integer(9)*C1**Integer(5)*L1**Integer(4)*L2 - Integer(13)*C1**Integer(5)*L1**Integer(3)*L2**Integer(2) - Integer(13)*C1**Integer(5)*L1**Integer(2)*L2**Integer(3) + Integer(9)*C1**Integer(5)*L1*L2**Integer(4))*k**Integer(4)*w**Integer(9) - Integer(4)*(C1**Integer(3)*L1**Integer(3) - Integer(41)*C1**Integer(3)*L1**Integer(2)*L2 - Integer(41)*C1**Integer(3)*L1*L2**Integer(2) + C1**Integer(3)*L2**Integer(3))*k**Integer(2)*w**Integer(5) - ((Integer(27)*C1**Integer(4)*L1**Integer(4) - Integer(36)*C1**Integer(4)*L1**Integer(3)*L2 + Integer(2)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2) - Integer(36)*C1**Integer(4)*L1*L2**Integer(3) + Integer(27)*C1**Integer(4)*L2**Integer(4))*k**Integer(4) - Integer(4)*(C1**Integer(4)*L1**Integer(3)*L2 - Integer(10)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2) + C1**Integer(4)*L1*L2**Integer(3))*k**Integer(2))*w**Integer(7) + Integer(8)*(Integer(2)*C1**Integer(2)*L1*L2 - Integer(5)*(Integer(3)*C1**Integer(2)*L1**Integer(2) + Integer(14)*C1**Integer(2)*L1*L2 + Integer(3)*C1**Integer(2)*L2**Integer(2))*k**Integer(2))*w**Integer(3) + Integer(16)*(Integer(24)*(C1*L1 + C1*L2)*k**Integer(2) - C1*L1 - C1*L2)*w)*w - Integer(16)/Integer(3))*k**Integer(2) + Integer(36)*k**Integer(2) + (Integer(2)*C1**Integer(6)*L1**Integer(3)*L2**Integer(3)*k**Integer(6)*w**Integer(11) - Integer(6)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2)*k**Integer(4)*w**Integer(7) + Integer(18)*(C1**Integer(3)*L1**Integer(2)*L2 + C1**Integer(3)*L1*L2**Integer(2))*k**Integer(4)*w**Integer(5) - Integer(18)*(C1*L1 + C1*L2)*k**Integer(2)*w + Integer(3)*(Integer(2)*C1**Integer(2)*L1*L2*k**Integer(2) - Integer(3)*(Integer(3)*C1**Integer(2)*L1**Integer(2) - Integer(2)*C1**Integer(2)*L1*L2 + Integer(3)*C1**Integer(2)*L2**Integer(2))*k**Integer(4))*w**Integer(3))*w - Integer(2))**(Integer(1)/Integer(3)))*abs(k) + C1**Integer(2)*k**Integer(2)*sgn2*sqrt(sqrt(Integer(3))*Integer(2)**(Integer(2)/Integer(3))*(sqrt(Integer(3))*(Integer(24)*Integer(2)**(Integer(2)/Integer(3))*k**Integer(2) - Integer(2)*(Integer(2)**(Integer(2)/Integer(3))*C1**Integer(4)*L1**Integer(2)*L2**Integer(2)*k**Integer(4)*w**Integer(7) - Integer(2)*Integer(2)**(Integer(2)/Integer(3))*C1**Integer(2)*L1*L2*k**Integer(2)*w**Integer(3) + Integer(6)*(Integer(2)**(Integer(2)/Integer(3))*C1*L1 + Integer(2)**(Integer(2)/Integer(3))*C1*L2)*k**Integer(2)*w)*w - Integer(2)*Integer(2)**(Integer(2)/Integer(3)) - Integer(2)*(Integer(9)*sqrt(Integer(256)/Integer(3)*k**Integer(2) - Integer(1)/Integer(3)*(Integer(4)*(C1**Integer(8)*L1**Integer(5)*L2**Integer(3) - Integer(2)*C1**Integer(8)*L1**Integer(4)*L2**Integer(4) + C1**Integer(8)*L1**Integer(3)*L2**Integer(5))*k**Integer(6)*w**Integer(15) - Integer(8)*(C1**Integer(6)*L1**Integer(4)*L2**Integer(2) - Integer(4)*C1**Integer(6)*L1**Integer(3)*L2**Integer(3) + C1**Integer(6)*L1**Integer(2)*L2**Integer(4))*k**Integer(4)*w**Integer(11) + Integer(4)*(Integer(9)*C1**Integer(5)*L1**Integer(4)*L2 - Integer(13)*C1**Integer(5)*L1**Integer(3)*L2**Integer(2) - Integer(13)*C1**Integer(5)*L1**Integer(2)*L2**Integer(3) + Integer(9)*C1**Integer(5)*L1*L2**Integer(4))*k**Integer(4)*w**Integer(9) - Integer(4)*(C1**Integer(3)*L1**Integer(3) - Integer(41)*C1**Integer(3)*L1**Integer(2)*L2 - Integer(41)*C1**Integer(3)*L1*L2**Integer(2) + C1**Integer(3)*L2**Integer(3))*k**Integer(2)*w**Integer(5) - ((Integer(27)*C1**Integer(4)*L1**Integer(4) - Integer(36)*C1**Integer(4)*L1**Integer(3)*L2 + Integer(2)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2) - Integer(36)*C1**Integer(4)*L1*L2**Integer(3) + Integer(27)*C1**Integer(4)*L2**Integer(4))*k**Integer(4) - Integer(4)*(C1**Integer(4)*L1**Integer(3)*L2 - Integer(10)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2) + C1**Integer(4)*L1*L2**Integer(3))*k**Integer(2))*w**Integer(7) + Integer(8)*(Integer(2)*C1**Integer(2)*L1*L2 - Integer(5)*(Integer(3)*C1**Integer(2)*L1**Integer(2) + Integer(14)*C1**Integer(2)*L1*L2 + Integer(3)*C1**Integer(2)*L2**Integer(2))*k**Integer(2))*w**Integer(3) + Integer(16)*(Integer(24)*(C1*L1 + C1*L2)*k**Integer(2) - C1*L1 - C1*L2)*w)*w - Integer(16)/Integer(3))*k**Integer(2) + Integer(36)*k**Integer(2) + (Integer(2)*C1**Integer(6)*L1**Integer(3)*L2**Integer(3)*k**Integer(6)*w**Integer(11) - Integer(6)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2)*k**Integer(4)*w**Integer(7) + Integer(18)*(C1**Integer(3)*L1**Integer(2)*L2 + C1**Integer(3)*L1*L2**Integer(2))*k**Integer(4)*w**Integer(5) - Integer(18)*(C1*L1 + C1*L2)*k**Integer(2)*w + Integer(3)*(Integer(2)*C1**Integer(2)*L1*L2*k**Integer(2) - Integer(3)*(Integer(3)*C1**Integer(2)*L1**Integer(2) - Integer(2)*C1**Integer(2)*L1*L2 + Integer(3)*C1**Integer(2)*L2**Integer(2))*k**Integer(4))*w**Integer(3))*w - Integer(2))**(Integer(2)/Integer(3)) + ((Integer(3)*Integer(2)**(Integer(1)/Integer(3))*C1**Integer(4)*L1**Integer(2) - Integer(2)*Integer(2)**(Integer(1)/Integer(3))*C1**Integer(4)*L1*L2 + Integer(3)*Integer(2)**(Integer(1)/Integer(3))*C1**Integer(4)*L2**Integer(2))*k**Integer(4)*w**Integer(7) + Integer(8)*Integer(2)**(Integer(1)/Integer(3))*C1**Integer(2)*k**Integer(2)*w**Integer(3))*(Integer(9)*sqrt(Integer(256)/Integer(3)*k**Integer(2) - Integer(1)/Integer(3)*(Integer(4)*(C1**Integer(8)*L1**Integer(5)*L2**Integer(3) - Integer(2)*C1**Integer(8)*L1**Integer(4)*L2**Integer(4) + C1**Integer(8)*L1**Integer(3)*L2**Integer(5))*k**Integer(6)*w**Integer(15) - Integer(8)*(C1**Integer(6)*L1**Integer(4)*L2**Integer(2) - Integer(4)*C1**Integer(6)*L1**Integer(3)*L2**Integer(3) + C1**Integer(6)*L1**Integer(2)*L2**Integer(4))*k**Integer(4)*w**Integer(11) + Integer(4)*(Integer(9)*C1**Integer(5)*L1**Integer(4)*L2 - Integer(13)*C1**Integer(5)*L1**Integer(3)*L2**Integer(2) - Integer(13)*C1**Integer(5)*L1**Integer(2)*L2**Integer(3) + Integer(9)*C1**Integer(5)*L1*L2**Integer(4))*k**Integer(4)*w**Integer(9) - Integer(4)*(C1**Integer(3)*L1**Integer(3) - Integer(41)*C1**Integer(3)*L1**Integer(2)*L2 - Integer(41)*C1**Integer(3)*L1*L2**Integer(2) + C1**Integer(3)*L2**Integer(3))*k**Integer(2)*w**Integer(5) - ((Integer(27)*C1**Integer(4)*L1**Integer(4) - Integer(36)*C1**Integer(4)*L1**Integer(3)*L2 + Integer(2)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2) - Integer(36)*C1**Integer(4)*L1*L2**Integer(3) + Integer(27)*C1**Integer(4)*L2**Integer(4))*k**Integer(4) - Integer(4)*(C1**Integer(4)*L1**Integer(3)*L2 - Integer(10)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2) + C1**Integer(4)*L1*L2**Integer(3))*k**Integer(2))*w**Integer(7) + Integer(8)*(Integer(2)*C1**Integer(2)*L1*L2 - Integer(5)*(Integer(3)*C1**Integer(2)*L1**Integer(2) + Integer(14)*C1**Integer(2)*L1*L2 + Integer(3)*C1**Integer(2)*L2**Integer(2))*k**Integer(2))*w**Integer(3) + Integer(16)*(Integer(24)*(C1*L1 + C1*L2)*k**Integer(2) - C1*L1 - C1*L2)*w)*w - Integer(16)/Integer(3))*k**Integer(2) + Integer(36)*k**Integer(2) + (Integer(2)*C1**Integer(6)*L1**Integer(3)*L2**Integer(3)*k**Integer(6)*w**Integer(11) - Integer(6)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2)*k**Integer(4)*w**Integer(7) + Integer(18)*(C1**Integer(3)*L1**Integer(2)*L2 + C1**Integer(3)*L1*L2**Integer(2))*k**Integer(4)*w**Integer(5) - Integer(18)*(C1*L1 + C1*L2)*k**Integer(2)*w + Integer(3)*(Integer(2)*C1**Integer(2)*L1*L2*k**Integer(2) - Integer(3)*(Integer(3)*C1**Integer(2)*L1**Integer(2) - Integer(2)*C1**Integer(2)*L1*L2 + Integer(3)*C1**Integer(2)*L2**Integer(2))*k**Integer(4))*w**Integer(3))*w - Integer(2))**(Integer(1)/Integer(3))/(C1**Integer(2)*k**Integer(2)*w**Integer(3)))*C1*w**Integer(2)*sqrt(-(Integer(48)*Integer(2)**(Integer(1)/Integer(3))*k**Integer(2) - Integer(4)*(Integer(2)**(Integer(1)/Integer(3))*C1**Integer(4)*L1**Integer(2)*L2**Integer(2)*k**Integer(4)*w**Integer(7) - Integer(2)*Integer(2)**(Integer(1)/Integer(3))*C1**Integer(2)*L1*L2*k**Integer(2)*w**Integer(3) + Integer(6)*(Integer(2)**(Integer(1)/Integer(3))*C1*L1 + Integer(2)**(Integer(1)/Integer(3))*C1*L2)*k**Integer(2)*w)*w - Integer(2)*Integer(2)**(Integer(2)/Integer(3))*(Integer(9)*sqrt(Integer(256)/Integer(3)*k**Integer(2) - Integer(1)/Integer(3)*(Integer(4)*(C1**Integer(8)*L1**Integer(5)*L2**Integer(3) - Integer(2)*C1**Integer(8)*L1**Integer(4)*L2**Integer(4) + C1**Integer(8)*L1**Integer(3)*L2**Integer(5))*k**Integer(6)*w**Integer(15) - Integer(8)*(C1**Integer(6)*L1**Integer(4)*L2**Integer(2) - Integer(4)*C1**Integer(6)*L1**Integer(3)*L2**Integer(3) + C1**Integer(6)*L1**Integer(2)*L2**Integer(4))*k**Integer(4)*w**Integer(11) + Integer(4)*(Integer(9)*C1**Integer(5)*L1**Integer(4)*L2 - Integer(13)*C1**Integer(5)*L1**Integer(3)*L2**Integer(2) - Integer(13)*C1**Integer(5)*L1**Integer(2)*L2**Integer(3) + Integer(9)*C1**Integer(5)*L1*L2**Integer(4))*k**Integer(4)*w**Integer(9) - Integer(4)*(C1**Integer(3)*L1**Integer(3) - Integer(41)*C1**Integer(3)*L1**Integer(2)*L2 - Integer(41)*C1**Integer(3)*L1*L2**Integer(2) + C1**Integer(3)*L2**Integer(3))*k**Integer(2)*w**Integer(5) - ((Integer(27)*C1**Integer(4)*L1**Integer(4) - Integer(36)*C1**Integer(4)*L1**Integer(3)*L2 + Integer(2)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2) - Integer(36)*C1**Integer(4)*L1*L2**Integer(3) + Integer(27)*C1**Integer(4)*L2**Integer(4))*k**Integer(4) - Integer(4)*(C1**Integer(4)*L1**Integer(3)*L2 - Integer(10)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2) + C1**Integer(4)*L1*L2**Integer(3))*k**Integer(2))*w**Integer(7) + Integer(8)*(Integer(2)*C1**Integer(2)*L1*L2 - Integer(5)*(Integer(3)*C1**Integer(2)*L1**Integer(2) + Integer(14)*C1**Integer(2)*L1*L2 + Integer(3)*C1**Integer(2)*L2**Integer(2))*k**Integer(2))*w**Integer(3) + Integer(16)*(Integer(24)*(C1*L1 + C1*L2)*k**Integer(2) - C1*L1 - C1*L2)*w)*w - Integer(16)/Integer(3))*k**Integer(2) + Integer(36)*k**Integer(2) + (Integer(2)*C1**Integer(6)*L1**Integer(3)*L2**Integer(3)*k**Integer(6)*w**Integer(11) - Integer(6)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2)*k**Integer(4)*w**Integer(7) + Integer(18)*(C1**Integer(3)*L1**Integer(2)*L2 + C1**Integer(3)*L1*L2**Integer(2))*k**Integer(4)*w**Integer(5) - Integer(18)*(C1*L1 + C1*L2)*k**Integer(2)*w + Integer(3)*(Integer(2)*C1**Integer(2)*L1*L2*k**Integer(2) - Integer(3)*(Integer(3)*C1**Integer(2)*L1**Integer(2) - Integer(2)*C1**Integer(2)*L1*L2 + Integer(3)*C1**Integer(2)*L2**Integer(2))*k**Integer(4))*w**Integer(3))*w - Integer(2))**(Integer(2)/Integer(3)) - Integer(4)*Integer(2)**(Integer(1)/Integer(3)) - ((Integer(3)*C1**Integer(4)*L1**Integer(2) - Integer(2)*C1**Integer(4)*L1*L2 + Integer(3)*C1**Integer(4)*L2**Integer(2))*k**Integer(4)*w**Integer(7) + Integer(8)*C1**Integer(2)*k**Integer(2)*w**Integer(3))*(Integer(9)*sqrt(Integer(256)/Integer(3)*k**Integer(2) - Integer(1)/Integer(3)*(Integer(4)*(C1**Integer(8)*L1**Integer(5)*L2**Integer(3) - Integer(2)*C1**Integer(8)*L1**Integer(4)*L2**Integer(4) + C1**Integer(8)*L1**Integer(3)*L2**Integer(5))*k**Integer(6)*w**Integer(15) - Integer(8)*(C1**Integer(6)*L1**Integer(4)*L2**Integer(2) - Integer(4)*C1**Integer(6)*L1**Integer(3)*L2**Integer(3) + C1**Integer(6)*L1**Integer(2)*L2**Integer(4))*k**Integer(4)*w**Integer(11) + Integer(4)*(Integer(9)*C1**Integer(5)*L1**Integer(4)*L2 - Integer(13)*C1**Integer(5)*L1**Integer(3)*L2**Integer(2) - Integer(13)*C1**Integer(5)*L1**Integer(2)*L2**Integer(3) + Integer(9)*C1**Integer(5)*L1*L2**Integer(4))*k**Integer(4)*w**Integer(9) - Integer(4)*(C1**Integer(3)*L1**Integer(3) - Integer(41)*C1**Integer(3)*L1**Integer(2)*L2 - Integer(41)*C1**Integer(3)*L1*L2**Integer(2) + C1**Integer(3)*L2**Integer(3))*k**Integer(2)*w**Integer(5) - ((Integer(27)*C1**Integer(4)*L1**Integer(4) - Integer(36)*C1**Integer(4)*L1**Integer(3)*L2 + Integer(2)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2) - Integer(36)*C1**Integer(4)*L1*L2**Integer(3) + Integer(27)*C1**Integer(4)*L2**Integer(4))*k**Integer(4) - Integer(4)*(C1**Integer(4)*L1**Integer(3)*L2 - Integer(10)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2) + C1**Integer(4)*L1*L2**Integer(3))*k**Integer(2))*w**Integer(7) + Integer(8)*(Integer(2)*C1**Integer(2)*L1*L2 - Integer(5)*(Integer(3)*C1**Integer(2)*L1**Integer(2) + Integer(14)*C1**Integer(2)*L1*L2 + Integer(3)*C1**Integer(2)*L2**Integer(2))*k**Integer(2))*w**Integer(3) + Integer(16)*(Integer(24)*(C1*L1 + C1*L2)*k**Integer(2) - C1*L1 - C1*L2)*w)*w - Integer(16)/Integer(3))*k**Integer(2) + Integer(36)*k**Integer(2) + (Integer(2)*C1**Integer(6)*L1**Integer(3)*L2**Integer(3)*k**Integer(6)*w**Integer(11) - Integer(6)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2)*k**Integer(4)*w**Integer(7) + Integer(18)*(C1**Integer(3)*L1**Integer(2)*L2 + C1**Integer(3)*L1*L2**Integer(2))*k**Integer(4)*w**Integer(5) - Integer(18)*(C1*L1 + C1*L2)*k**Integer(2)*w + Integer(3)*(Integer(2)*C1**Integer(2)*L1*L2*k**Integer(2) - Integer(3)*(Integer(3)*C1**Integer(2)*L1**Integer(2) - Integer(2)*C1**Integer(2)*L1*L2 + Integer(3)*C1**Integer(2)*L2**Integer(2))*k**Integer(4))*w**Integer(3))*w - Integer(2))**(Integer(1)/Integer(3))/(C1**Integer(2)*k**Integer(2)*w**Integer(3)))/(Integer(9)*sqrt(Integer(256)/Integer(3)*k**Integer(2) - Integer(1)/Integer(3)*(Integer(4)*(C1**Integer(8)*L1**Integer(5)*L2**Integer(3) - Integer(2)*C1**Integer(8)*L1**Integer(4)*L2**Integer(4) + C1**Integer(8)*L1**Integer(3)*L2**Integer(5))*k**Integer(6)*w**Integer(15) - Integer(8)*(C1**Integer(6)*L1**Integer(4)*L2**Integer(2) - Integer(4)*C1**Integer(6)*L1**Integer(3)*L2**Integer(3) + C1**Integer(6)*L1**Integer(2)*L2**Integer(4))*k**Integer(4)*w**Integer(11) + Integer(4)*(Integer(9)*C1**Integer(5)*L1**Integer(4)*L2 - Integer(13)*C1**Integer(5)*L1**Integer(3)*L2**Integer(2) - Integer(13)*C1**Integer(5)*L1**Integer(2)*L2**Integer(3) + Integer(9)*C1**Integer(5)*L1*L2**Integer(4))*k**Integer(4)*w**Integer(9) - Integer(4)*(C1**Integer(3)*L1**Integer(3) - Integer(41)*C1**Integer(3)*L1**Integer(2)*L2 - Integer(41)*C1**Integer(3)*L1*L2**Integer(2) + C1**Integer(3)*L2**Integer(3))*k**Integer(2)*w**Integer(5) - ((Integer(27)*C1**Integer(4)*L1**Integer(4) - Integer(36)*C1**Integer(4)*L1**Integer(3)*L2 + Integer(2)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2) - Integer(36)*C1**Integer(4)*L1*L2**Integer(3) + Integer(27)*C1**Integer(4)*L2**Integer(4))*k**Integer(4) - Integer(4)*(C1**Integer(4)*L1**Integer(3)*L2 - Integer(10)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2) + C1**Integer(4)*L1*L2**Integer(3))*k**Integer(2))*w**Integer(7) + Integer(8)*(Integer(2)*C1**Integer(2)*L1*L2 - Integer(5)*(Integer(3)*C1**Integer(2)*L1**Integer(2) + Integer(14)*C1**Integer(2)*L1*L2 + Integer(3)*C1**Integer(2)*L2**Integer(2))*k**Integer(2))*w**Integer(3) + Integer(16)*(Integer(24)*(C1*L1 + C1*L2)*k**Integer(2) - C1*L1 - C1*L2)*w)*w - Integer(16)/Integer(3))*k**Integer(2) + Integer(36)*k**Integer(2) + (Integer(2)*C1**Integer(6)*L1**Integer(3)*L2**Integer(3)*k**Integer(6)*w**Integer(11) - Integer(6)*C1**Integer(4)*L1**Integer(2)*L2**Integer(2)*k**Integer(4)*w**Integer(7) + Integer(18)*(C1**Integer(3)*L1**Integer(2)*L2 + C1**Integer(3)*L1*L2**Integer(2))*k**Integer(4)*w**Integer(5) - Integer(18)*(C1*L1 + C1*L2)*k**Integer(2)*w + Integer(3)*(Integer(2)*C1**Integer(2)*L1*L2*k**Integer(2) - Integer(3)*(Integer(3)*C1**Integer(2)*L1**Integer(2) - Integer(2)*C1**Integer(2)*L1*L2 + Integer(3)*C1**Integer(2)*L2**Integer(2))*k**Integer(4))*w**Integer(3))*w - Integer(2))**(Integer(1)/Integer(3)))*abs(k) + Integer(9)*Integer(2)**(Integer(1)/Integer(3)))))))
4
----> 5 expr.simplify_full()
/opt/sagemath-9.3/local/lib/python3.7/site-packages/sage/symbolic/expression.pyx in sage.symbolic.expression.Expression.simplify_full (build/cythonized/sage/symbolic/expression.cpp:53294)()
10196 x = self
10197 x = x.simplify_factorial()
> 10198 x = x.simplify_rectform()
10199 x = x.simplify_trig()
10200 x = x.simplify_rational()
/opt/sagemath-9.3/local/lib/python3.7/site-packages/sage/symbolic/expression.pyx in sage.symbolic.expression.Expression.simplify_rectform (build/cythonized/sage/symbolic/expression.cpp:54760)()
10347
10348 """
> 10349 simplified_expr = self.rectform()
10350
10351 if complexity_measure is None:
/opt/sagemath-9.3/local/lib/python3.7/site-packages/sage/symbolic/expression.pyx in sage.symbolic.expression.Expression.rectform (build/cythonized/sage/symbolic/expression.cpp:52806)()
10028
10029 """
> 10030 return self.maxima_methods().rectform()
10031
10032 def unhold(self, exclude=None):
/opt/sagemath-9.3/local/lib/python3.7/site-packages/sage/symbolic/maxima_wrapper.py in __call__(self, *args, **kwds)
30 """
31 return super(MaximaFunctionElementWrapper, self).__call__(*args,
---> 32 **kwds).sage()
33
34 class MaximaWrapper(SageObject):
/opt/sagemath-9.3/local/lib/python3.7/site-packages/sage/interfaces/interface.py in __call__(self, *args, **kwds)
677
678 def __call__(self, *args, **kwds):
--> 679 return self._obj.parent().function_call(self._name, [self._obj] + list(args), kwds)
680
681 def help(self):
/opt/sagemath-9.3/local/lib/python3.7/site-packages/sage/interfaces/interface.py in function_call(self, function, args, kwds)
598 [s.name() for s in args],
599 ['%s=%s'%(key,value.name()) for key, value in kwds.items()])
--> 600 return self.new(s)
601
602 def _function_call_string(self, function, args, kwds):
/opt/sagemath-9.3/local/lib/python3.7/site-packages/sage/interfaces/interface.py in new(self, code)
367
368 def new(self, code):
--> 369 return self(code)
370
371 ###################################################################
/opt/sagemath-9.3/local/lib/python3.7/site-packages/sage/interfaces/interface.py in __call__(self, x, name)
293
294 if isinstance(x, str):
--> 295 return cls(self, x, name=name)
296 try:
297 # Special methods do not and should not have an option to
/opt/sagemath-9.3/local/lib/python3.7/site-packages/sage/interfaces/interface.py in __init__(self, parent, value, is_name, name)
718 self._name = parent._create(value, name=name)
719 except (TypeError, RuntimeError, ValueError) as x:
--> 720 raise TypeError(x)
721
722 def _latex_(self):
TypeError: Aborted
When I leave just the last simple term + 9*2^(1/3) away, it works:
var('C1 L1 L2 k w sgn2 sgn3')
expr = (1/2*sqrt(1/3)*sqrt((3*(3*C1^2*L1 - C1^2*L2)*k^2*w^3 - sqrt(3)*C1*sgn3*w*sqrt(-(48*2^(1/3)*k^2 - 4*(2^(1/3)*C1^4*L1^2*L2^2*k^4*w^7 - 2*2^(1/3)*C1^2*L1*L2*k^2*w^3 + 6*(2^(1/3)*C1*L1 + 2^(1/3)*C1*L2)*k^2*w)*w - 2*2^(2/3)*(9*sqrt(256/3*k^2 - 1/3*(4*(C1^8*L1^5*L2^3 - 2*C1^8*L1^4*L2^4 + C1^8*L1^3*L2^5)*k^6*w^15 - 8*(C1^6*L1^4*L2^2 - 4*C1^6*L1^3*L2^3 + C1^6*L1^2*L2^4)*k^4*w^11 + 4*(9*C1^5*L1^4*L2 - 13*C1^5*L1^3*L2^2 - 13*C1^5*L1^2*L2^3 + 9*C1^5*L1*L2^4)*k^4*w^9 - 4*(C1^3*L1^3 - 41*C1^3*L1^2*L2 - 41*C1^3*L1*L2^2 + C1^3*L2^3)*k^2*w^5 - ((27*C1^4*L1^4 - 36*C1^4*L1^3*L2 + 2*C1^4*L1^2*L2^2 - 36*C1^4*L1*L2^3 + 27*C1^4*L2^4)*k^4 - 4*(C1^4*L1^3*L2 - 10*C1^4*L1^2*L2^2 + C1^4*L1*L2^3)*k^2)*w^7 + 8*(2*C1^2*L1*L2 - 5*(3*C1^2*L1^2 + 14*C1^2*L1*L2 + 3*C1^2*L2^2)*k^2)*w^3 + 16*(24*(C1*L1 + C1*L2)*k^2 - C1*L1 - C1*L2)*w)*w - 16/3)*k^2 + 36*k^2 + (2*C1^6*L1^3*L2^3*k^6*w^11 - 6*C1^4*L1^2*L2^2*k^4*w^7 + 18*(C1^3*L1^2*L2 + C1^3*L1*L2^2)*k^4*w^5 - 18*(C1*L1 + C1*L2)*k^2*w + 3*(2*C1^2*L1*L2*k^2 - 3*(3*C1^2*L1^2 - 2*C1^2*L1*L2 + 3*C1^2*L2^2)*k^4)*w^3)*w - 2)^(2/3) - 4*2^(1/3) - ((3*C1^4*L1^2 - 2*C1^4*L1*L2 + 3*C1^4*L2^2)*k^4*w^7 + 8*C1^2*k^2*w^3)*(9*sqrt(256/3*k^2 - 1/3*(4*(C1^8*L1^5*L2^3 - 2*C1^8*L1^4*L2^4 + C1^8*L1^3*L2^5)*k^6*w^15 - 8*(C1^6*L1^4*L2^2 - 4*C1^6*L1^3*L2^3 + C1^6*L1^2*L2^4)*k^4*w^11 + 4*(9*C1^5*L1^4*L2 - 13*C1^5*L1^3*L2^2 - 13*C1^5*L1^2*L2^3 + 9*C1^5*L1*L2^4)*k^4*w^9 - 4*(C1^3*L1^3 - 41*C1^3*L1^2*L2 - 41*C1^3*L1*L2^2 + C1^3*L2^3)*k^2*w^5 - ((27*C1^4*L1^4 - 36*C1^4*L1^3*L2 + 2*C1^4*L1^2*L2^2 - 36*C1^4*L1*L2^3 + 27*C1^4*L2^4)*k^4 - 4*(C1^4*L1^3*L2 - 10*C1^4*L1^2*L2^2 + C1^4*L1*L2^3)*k^2)*w^7 + 8*(2*C1^2*L1*L2 - 5*(3*C1^2*L1^2 + 14*C1^2*L1*L2 + 3*C1^2*L2^2)*k^2)*w^3 + 16*(24*(C1*L1 + C1*L2)*k^2 - C1*L1 - C1*L2)*w)*w - 16/3)*k^2 + 36*k^2 + (2*C1^6*L1^3*L2^3*k^6*w^11 - 6*C1^4*L1^2*L2^2*k^4*w^7 + 18*(C1^3*L1^2*L2 + C1^3*L1*L2^2)*k^4*w^5 - 18*(C1*L1 + C1*L2)*k^2*w + 3*(2*C1^2*L1*L2*k^2 - 3*(3*C1^2*L1^2 - 2*C1^2*L1*L2 + 3*C1^2*L2^2)*k^4)*w^3)*w - 2)^(1/3)/(C1^2*k^2*w^3))/(9*sqrt(256/3*k^2 - 1/3*(4*(C1^8*L1^5*L2^3 - 2*C1^8*L1^4*L2^4 + C1^8*L1^3*L2^5)*k^6*w^15 - 8*(C1^6*L1^4*L2^2 - 4*C1^6*L1^3*L2^3 + C1^6*L1^2*L2^4)*k^4*w^11 + 4*(9*C1^5*L1^4*L2 - 13*C1^5*L1^3*L2^2 - 13*C1^5*L1^2*L2^3 + 9*C1^5*L1*L2^4)*k^4*w^9 - 4*(C1^3*L1^3 - 41*C1^3*L1^2*L2 - 41*C1^3*L1*L2^2 + C1^3*L2^3)*k^2*w^5 - ((27*C1^4*L1^4 - 36*C1^4*L1^3*L2 + 2*C1^4*L1^2*L2^2 - 36*C1^4*L1*L2^3 + 27*C1^4*L2^4)*k^4 - 4*(C1^4*L1^3*L2 - 10*C1^4*L1^2*L2^2 + C1^4*L1*L2^3)*k^2)*w^7 + 8*(2*C1^2*L1*L2 - 5*(3*C1^2*L1^2 + 14*C1^2*L1*L2 + 3*C1^2*L2^2)*k^2)*w^3 + 16*(24*(C1*L1 + C1*L2)*k^2 - C1*L1 - C1*L2)*w)*w - 16/3)*k^2 + 36*k^2 + (2*C1^6*L1^3*L2^3*k^6*w^11 - 6*C1^4*L1^2*L2^2*k^4*w^7 + 18*(C1^3*L1^2*L2 + C1^3*L1*L2^2)*k^4*w^5 - 18*(C1*L1 + C1*L2)*k^2*w + 3*(2*C1^2*L1*L2*k^2 - 3*(3*C1^2*L1^2 - 2*C1^2*L1*L2 + 3*C1^2*L2^2)*k^4)*w^3)*w - 2)^(1/3))*abs(k) + C1^2*k^2*sgn2*sqrt(sqrt(3)*2^(2/3)*(sqrt(3)*(24*2^(2/3)*k^2 - 2*(2^(2/3)*C1^4*L1^2*L2^2*k^4*w^7 - 2*2^(2/3)*C1^2*L1*L2*k^2*w^3 + 6*(2^(2/3)*C1*L1 + 2^(2/3)*C1*L2)*k^2*w)*w - 2*2^(2/3) - 2*(9*sqrt(256/3*k^2 - 1/3*(4*(C1^8*L1^5*L2^3 - 2*C1^8*L1^4*L2^4 + C1^8*L1^3*L2^5)*k^6*w^15 - 8*(C1^6*L1^4*L2^2 - 4*C1^6*L1^3*L2^3 + C1^6*L1^2*L2^4)*k^4*w^11 + 4*(9*C1^5*L1^4*L2 - 13*C1^5*L1^3*L2^2 - 13*C1^5*L1^2*L2^3 + 9*C1^5*L1*L2^4)*k^4*w^9 - 4*(C1^3*L1^3 - 41*C1^3*L1^2*L2 - 41*C1^3*L1*L2^2 + C1^3*L2^3)*k^2*w^5 - ((27*C1^4*L1^4 - 36*C1^4*L1^3*L2 + 2*C1^4*L1^2*L2^2 - 36*C1^4*L1*L2^3 + 27*C1^4*L2^4)*k^4 - 4*(C1^4*L1^3*L2 - 10*C1^4*L1^2*L2^2 + C1^4*L1*L2^3)*k^2)*w^7 + 8*(2*C1^2*L1*L2 - 5*(3*C1^2*L1^2 + 14*C1^2*L1*L2 + 3*C1^2*L2^2)*k^2)*w^3 + 16*(24*(C1*L1 + C1*L2)*k^2 - C1*L1 - C1*L2)*w)*w - 16/3)*k^2 + 36*k^2 + (2*C1^6*L1^3*L2^3*k^6*w^11 - 6*C1^4*L1^2*L2^2*k^4*w^7 + 18*(C1^3*L1^2*L2 + C1^3*L1*L2^2)*k^4*w^5 - 18*(C1*L1 + C1*L2)*k^2*w + 3*(2*C1^2*L1*L2*k^2 - 3*(3*C1^2*L1^2 - 2*C1^2*L1*L2 + 3*C1^2*L2^2)*k^4)*w^3)*w - 2)^(2/3) + ((3*2^(1/3)*C1^4*L1^2 - 2*2^(1/3)*C1^4*L1*L2 + 3*2^(1/3)*C1^4*L2^2)*k^4*w^7 + 8*2^(1/3)*C1^2*k^2*w^3)*(9*sqrt(256/3*k^2 - 1/3*(4*(C1^8*L1^5*L2^3 - 2*C1^8*L1^4*L2^4 + C1^8*L1^3*L2^5)*k^6*w^15 - 8*(C1^6*L1^4*L2^2 - 4*C1^6*L1^3*L2^3 + C1^6*L1^2*L2^4)*k^4*w^11 + 4*(9*C1^5*L1^4*L2 - 13*C1^5*L1^3*L2^2 - 13*C1^5*L1^2*L2^3 + 9*C1^5*L1*L2^4)*k^4*w^9 - 4*(C1^3*L1^3 - 41*C1^3*L1^2*L2 - 41*C1^3*L1*L2^2 + C1^3*L2^3)*k^2*w^5 - ((27*C1^4*L1^4 - 36*C1^4*L1^3*L2 + 2*C1^4*L1^2*L2^2 - 36*C1^4*L1*L2^3 + 27*C1^4*L2^4)*k^4 - 4*(C1^4*L1^3*L2 - 10*C1^4*L1^2*L2^2 + C1^4*L1*L2^3)*k^2)*w^7 + 8*(2*C1^2*L1*L2 - 5*(3*C1^2*L1^2 + 14*C1^2*L1*L2 + 3*C1^2*L2^2)*k^2)*w^3 + 16*(24*(C1*L1 + C1*L2)*k^2 - C1*L1 - C1*L2)*w)*w - 16/3)*k^2 + 36*k^2 + (2*C1^6*L1^3*L2^3*k^6*w^11 - 6*C1^4*L1^2*L2^2*k^4*w^7 + 18*(C1^3*L1^2*L2 + C1^3*L1*L2^2)*k^4*w^5 - 18*(C1*L1 + C1*L2)*k^2*w + 3*(2*C1^2*L1*L2*k^2 - 3*(3*C1^2*L1^2 - 2*C1^2*L1*L2 + 3*C1^2*L2^2)*k^4)*w^3)*w - 2)^(1/3)/(C1^2*k^2*w^3))*C1*w^2*sqrt(-(48*2^(1/3)*k^2 - 4*(2^(1/3)*C1^4*L1^2*L2^2*k^4*w^7 - 2*2^(1/3)*C1^2*L1*L2*k^2*w^3 + 6*(2^(1/3)*C1*L1 + 2^(1/3)*C1*L2)*k^2*w)*w - 2*2^(2/3)*(9*sqrt(256/3*k^2 - 1/3*(4*(C1^8*L1^5*L2^3 - 2*C1^8*L1^4*L2^4 + C1^8*L1^3*L2^5)*k^6*w^15 - 8*(C1^6*L1^4*L2^2 - 4*C1^6*L1^3*L2^3 + C1^6*L1^2*L2^4)*k^4*w^11 + 4*(9*C1^5*L1^4*L2 - 13*C1^5*L1^3*L2^2 - 13*C1^5*L1^2*L2^3 + 9*C1^5*L1*L2^4)*k^4*w^9 - 4*(C1^3*L1^3 - 41*C1^3*L1^2*L2 - 41*C1^3*L1*L2^2 + C1^3*L2^3)*k^2*w^5 - ((27*C1^4*L1^4 - 36*C1^4*L1^3*L2 + 2*C1^4*L1^2*L2^2 - 36*C1^4*L1*L2^3 + 27*C1^4*L2^4)*k^4 - 4*(C1^4*L1^3*L2 - 10*C1^4*L1^2*L2^2 + C1^4*L1*L2^3)*k^2)*w^7 + 8*(2*C1^2*L1*L2 - 5*(3*C1^2*L1^2 + 14*C1^2*L1*L2 + 3*C1^2*L2^2)*k^2)*w^3 + 16*(24*(C1*L1 + C1*L2)*k^2 - C1*L1 - C1*L2)*w)*w - 16/3)*k^2 + 36*k^2 + (2*C1^6*L1^3*L2^3*k^6*w^11 - 6*C1^4*L1^2*L2^2*k^4*w^7 + 18*(C1^3*L1^2*L2 + C1^3*L1*L2^2)*k^4*w^5 - 18*(C1*L1 + C1*L2)*k^2*w + 3*(2*C1^2*L1*L2*k^2 - 3*(3*C1^2*L1^2 - 2*C1^2*L1*L2 + 3*C1^2*L2^2)*k^4)*w^3)*w - 2)^(2/3) - 4*2^(1/3) - ((3*C1^4*L1^2 - 2*C1^4*L1*L2 + 3*C1^4*L2^2)*k^4*w^7 + 8*C1^2*k^2*w^3)*(9*sqrt(256/3*k^2 - 1/3*(4*(C1^8*L1^5*L2^3 - 2*C1^8*L1^4*L2^4 + C1^8*L1^3*L2^5)*k^6*w^15 - 8*(C1^6*L1^4*L2^2 - 4*C1^6*L1^3*L2^3 + C1^6*L1^2*L2^4)*k^4*w^11 + 4*(9*C1^5*L1^4*L2 - 13*C1^5*L1^3*L2^2 - 13*C1^5*L1^2*L2^3 + 9*C1^5*L1*L2^4)*k^4*w^9 - 4*(C1^3*L1^3 - 41*C1^3*L1^2*L2 - 41*C1^3*L1*L2^2 + C1^3*L2^3)*k^2*w^5 - ((27*C1^4*L1^4 - 36*C1^4*L1^3*L2 + 2*C1^4*L1^2*L2^2 - 36*C1^4*L1*L2^3 + 27*C1^4*L2^4)*k^4 - 4*(C1^4*L1^3*L2 - 10*C1^4*L1^2*L2^2 + C1^4*L1*L2^3)*k^2)*w^7 + 8*(2*C1^2*L1*L2 - 5*(3*C1^2*L1^2 + 14*C1^2*L1*L2 + 3*C1^2*L2^2)*k^2)*w^3 + 16*(24*(C1*L1 + C1*L2)*k^2 - C1*L1 - C1*L2)*w)*w - 16/3)*k^2 + 36*k^2 + (2*C1^6*L1^3*L2^3*k^6*w^11 - 6*C1^4*L1^2*L2^2*k^4*w^7 + 18*(C1^3*L1^2*L2 + C1^3*L1*L2^2)*k^4*w^5 - 18*(C1*L1 + C1*L2)*k^2*w + 3*(2*C1^2*L1*L2*k^2 - 3*(3*C1^2*L1^2 - 2*C1^2*L1*L2 + 3*C1^2*L2^2)*k^4)*w^3)*w - 2)^(1/3)/(C1^2*k^2*w^3))/(9*sqrt(256/3*k^2 - 1/3*(4*(C1^8*L1^5*L2^3 - 2*C1^8*L1^4*L2^4 + C1^8*L1^3*L2^5)*k^6*w^15 - 8*(C1^6*L1^4*L2^2 - 4*C1^6*L1^3*L2^3 + C1^6*L1^2*L2^4)*k^4*w^11 + 4*(9*C1^5*L1^4*L2 - 13*C1^5*L1^3*L2^2 - 13*C1^5*L1^2*L2^3 + 9*C1^5*L1*L2^4)*k^4*w^9 - 4*(C1^3*L1^3 - 41*C1^3*L1^2*L2 - 41*C1^3*L1*L2^2 + C1^3*L2^3)*k^2*w^5 - ((27*C1^4*L1^4 - 36*C1^4*L1^3*L2 + 2*C1^4*L1^2*L2^2 - 36*C1^4*L1*L2^3 + 27*C1^4*L2^4)*k^4 - 4*(C1^4*L1^3*L2 - 10*C1^4*L1^2*L2^2 + C1^4*L1*L2^3)*k^2)*w^7 + 8*(2*C1^2*L1*L2 - 5*(3*C1^2*L1^2 + 14*C1^2*L1*L2 + 3*C1^2*L2^2)*k^2)*w^3 + 16*(24*(C1*L1 + C1*L2)*k^2 - C1*L1 - C1*L2)*w)*w - 16/3)*k^2 + 36*k^2 + (2*C1^6*L1^3*L2^3*k^6*w^11 - 6*C1^4*L1^2*L2^2*k^4*w^7 + 18*(C1^3*L1^2*L2 + C1^3*L1*L2^2)*k^4*w^5 - 18*(C1*L1 + C1*L2)*k^2*w + 3*(2*C1^2*L1*L2*k^2 - 3*(3*C1^2*L1^2 - 2*C1^2*L1*L2 + 3*C1^2*L2^2)*k^4)*w^3)*w - 2)^(1/3))*abs(k))))))
expr.simplify_full()
Here the output comes without an error:
1/6*sqrt(3)*sqrt(sqrt(3)*sqrt(2)*sqrt(-(2*2^(1/3)*C1^5*L1^2*L2^2*k^4*w^10*abs(k) - 4*2^(1/3)*C1^3 ...
Please try to ask a question in a minimal form. A lot of work is needed to trace back the error, and so far i see no reason for med to do the work. Try to isolate the inner most pieces like square roots. If there is no simplification inside them, and there is none for some square root i found somewhere, which is the expectation in the simplification? More or less, the machine has to do the same. It would make sense especially in the case of such a question to mention why this question is interesting or important, how the expression was generated, and what did you try to do on the road for special pieces inside it.
I have updated the expression. I shortened it until it works (see above). One expression works (the second one I posted), but when adding the simple term + 9*2^(1/3) at the end, it doesn't work any more. Maybe some system limitations lead to the error (memory, stack (number of nested brackets, ...). It would also be interesting if on other systems it's like on my system (maybe others have different system settings and both expressions work ?)