def sample_apoly_points_via_giac_rur(manifold, n):
import giac_rur
I = extended_ptolemy_equations(manifold)
R = I.ring()
p = cyclotomic_polynomial(n, var=R('M'))
I = I + [p]
return giac_rur.rational_univariate_representation(I)
def test_as_cusped(manifold):
import giac_rur
for obs in manifold.ptolemy_generalized_obstruction_classes(2):
I = extended_ptolemy_equations(manifold, obs)
R = I.ring()
M, L = R('M'), R('L')
I = I + [M - 1, L - 1]
if I.dimension() == 0:
print(giac_rur.rational_univariate_representation(I))
def test_rur(manifold):
G = manifold.fundamental_group(True, True, False)
I = character_variety(G)
return giac_rur.rational_univariate_representation(I)
def test_direct(manifold):
import giac_rur
for obs in manifold.ptolemy_generalized_obstruction_classes(2):
I = manifold.ptolemy_variety(2, obs).ideal_with_non_zero_condition
if I.dimension() == 0:
print(giac_rur.rational_univariate_representation(I))
def rur_for_dehn_filling(manifold):
import giac_rur
I = ptolemy_ideal_for_filled(manifold)
return giac_rur.rational_univariate_representation(I)