W8WD1 "controller design for the brockett-integrator" does not finish calculation

[JuliusFiedler]

steps to reproduce:
ackrep -c W8WD1 in the appropriate environment
I canceled the calculation after ~7 min

Traceback (most recent call last):
  File "execscript.py", line 45, in <module>
    solution_data = solution.solve(ps.ProblemSpecification)
  File "/code/ackrep_data/problem_solutions/control_brockett_integrator/solution.py", line 56, in solve
    r1_opt_interp = mgic.shortest_curve(problem_spec)
  File "/code/ackrep_data/method_packages/method_geometry_inspired_control/_build/method_geometry_inspired_control.py", line 80, in shortest_curve
    sol = sp.solve(Le.diff(r1), r1)
  File "/opt/conda/lib/python3.8/site-packages/sympy/solvers/solvers.py", line 1106, in solve
    solution = _solve(f[0], *symbols, **flags)
  File "/opt/conda/lib/python3.8/site-packages/sympy/solvers/solvers.py", line 1732, in _solve
    result = [s for s in result if
  File "/opt/conda/lib/python3.8/site-packages/sympy/solvers/solvers.py", line 1733, in <listcomp>
    not any(checksol(d, {symbol: s}, **flags)
  File "/opt/conda/lib/python3.8/site-packages/sympy/solvers/solvers.py", line 1733, in <genexpr>
    not any(checksol(d, {symbol: s}, **flags)
  File "/opt/conda/lib/python3.8/site-packages/sympy/solvers/solvers.py", line 301, in checksol
    if not val.is_constant(*list(sol.keys()), simplify=not minimal):
  File "/opt/conda/lib/python3.8/site-packages/sympy/core/power.py", line 1937, in is_constant
    bz = b.equals(0)
  File "/opt/conda/lib/python3.8/site-packages/sympy/core/expr.py", line 736, in equals
    fac_zero = [fac.equals(0) for fac in factors]
  File "/opt/conda/lib/python3.8/site-packages/sympy/core/expr.py", line 736, in <listcomp>
    fac_zero = [fac.equals(0) for fac in factors]
  File "/opt/conda/lib/python3.8/site-packages/sympy/core/expr.py", line 740, in equals
    constant = diff.is_constant(simplify=False, failing_number=True)
  File "/opt/conda/lib/python3.8/site-packages/sympy/core/power.py", line 1943, in is_constant
    bcon = b.is_constant(*wrt)
  File "/opt/conda/lib/python3.8/site-packages/sympy/core/expr.py", line 660, in is_constant
    if b is not None and b is not S.NaN and b.equals(a) is False:
  File "/opt/conda/lib/python3.8/site-packages/sympy/core/expr.py", line 736, in equals
    fac_zero = [fac.equals(0) for fac in factors]
  File "/opt/conda/lib/python3.8/site-packages/sympy/core/expr.py", line 736, in <listcomp>
    fac_zero = [fac.equals(0) for fac in factors]
  File "/opt/conda/lib/python3.8/site-packages/sympy/core/expr.py", line 817, in equals
    mp = minimal_polynomial(diff)
  File "/opt/conda/lib/python3.8/site-packages/sympy/polys/numberfields/minpoly.py", line 701, in minimal_polynomial
    result = _minpoly_compose(ex, x, domain)
  File "/opt/conda/lib/python3.8/site-packages/sympy/polys/numberfields/minpoly.py", line 581, in _minpoly_compose
    res = _minpoly_add(x, dom, *ex.args)
  File "/opt/conda/lib/python3.8/site-packages/sympy/polys/numberfields/minpoly.py", line 373, in _minpoly_add
    mp = _minpoly_op_algebraic_element(Add, p, px, x, dom, mp1=mp)
  File "/opt/conda/lib/python3.8/site-packages/sympy/polys/numberfields/minpoly.py", line 286, in _minpoly_op_algebraic_element
    _, factors = r.factor_list()
  File "/opt/conda/lib/python3.8/site-packages/sympy/polys/polytools.py", line 3353, in factor_list
    coeff, factors = f.rep.factor_list()
  File "/opt/conda/lib/python3.8/site-packages/sympy/polys/polyclasses.py", line 823, in factor_list
    coeff, factors = dmp_factor_list(f.rep, f.lev, f.dom)
  File "/opt/conda/lib/python3.8/site-packages/sympy/polys/factortools.py", line 1381, in dmp_factor_list
    return dup_factor_list(f, K0)
  File "/opt/conda/lib/python3.8/site-packages/sympy/polys/factortools.py", line 1330, in dup_factor_list
    coeff, factors = dup_zz_factor(f, K)
  File "/opt/conda/lib/python3.8/site-packages/sympy/polys/factortools.py", line 676, in dup_zz_factor
    H = dup_zz_zassenhaus(g, K)
  File "/opt/conda/lib/python3.8/site-packages/sympy/polys/factortools.py", line 394, in dup_zz_zassenhaus
    H = dup_mul(H, g[i], K)
  File "/opt/conda/lib/python3.8/site-packages/sympy/polys/densearith.py", line 783, in dup_mul
    lo, hi = dup_mul(fl, gl, K), dup_mul(fh, gh, K)
  File "/opt/conda/lib/python3.8/site-packages/sympy/polys/densearith.py", line 783, in dup_mul
    lo, hi = dup_mul(fl, gl, K), dup_mul(fh, gh, K)
  File "/opt/conda/lib/python3.8/site-packages/sympy/polys/densearith.py", line 766, in dup_mul
    for j in range(max(0, i - dg), min(df, i) + 1):`
KeyboardInterrupt