def test_canonical_X11_0(self):
X = self.X11
J = Jacobian(X)
g = X.genus()
P0 = X.base_place
oneforms = X.holomorphic_oneforms()
C = oneforms[0].valuation_divisor()
W = 2 * RiemannConstantVector(P0) + AbelMap(C)
self.assertLess(norm(J(W)), 1e-7)
def is_theta_zero(self, X, W, prec=1e-7):
Omega = X.riemann_matrix()
J = Jacobian(X)
W = J(W)
u = RiemannTheta.oscillatory_part(W, Omega)
self.assertLess(abs(u), prec)