colored jones polynomial function does not return

[GGN-2015]

Steps To Reproduce

In sage 10.4 use the following command in an interactive terminal of sage.

from sage.all import *
ans = Knot([[1, 5, 2, 4], [3, 7, 4, 6], [5, 3, 6, 2], [7, 13, 8, 12], [9, 18, 10, 19], [11, 15, 12, 14], [13, 9, 14, 8], [15, 22, 16, 1], [17, 20, 18, 21], [19, 10, 20, 11], [21, 16, 22, 17]]).colored_jones_polynomial(2)
print(ans)

Expected Behavior

according to mathematica, this knot has a 2-colored jones polynomial as 1/q^13 - 2/q^12 - 1/q^11 + 7/q^10 - 6/q^9 - 8/q^8 + 21/q^7 - 10/q^6 - 25/q^5 + 43/q^4 - 10/q^3 - 52/q^2 + 68/q - 82 q + 83 q^2 + 14 q^3 - 101 q^4 + 80 q^5 + 25 q^6 - 98 q^7 + 65 q^8 + 25 q^9 - 74 q^10 + 42 q^11 + 17 q^12 - 40 q^13 + 20 q^14 + 7 q^15 - 14 q^16 + 7 q^17 + q^18 - 3 q^19 + q^20.

Actual Behavior

The function doesn't return in 5 min.

Additional Information

Since the pd_code is not very long and mathematica responses quickly, I think that maybe something is wrong with the implementation of colored jones polynomial in sage.

https://ask.sagemath.org/question/79833/colored-jones-polynomial-function-does-not-return/

Environment

• OS: Ubuntu 24.04 (installed with conda 24.9.0)
• Sage Version: 10.4

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