Skip to content

Instantly share code, notes, and snippets.

@edgarcosta
Created December 4, 2025 01:18
Show Gist options
  • Select an option

  • Save edgarcosta/40e19c0ee5b5694213cc222d2e79b0dc to your computer and use it in GitHub Desktop.

Select an option

Save edgarcosta/40e19c0ee5b5694213cc222d2e79b0dc to your computer and use it in GitHub Desktop.
128
2025-11-18 09:55:05.955863-05:00
189
110
1 ----------------------------------------------------------------------------------------------------
1 2025-11-18 09:55:29.444685-05:00
1 hs 0:03:46.636426
1 ----------------------------------------------------------------------------------------------------
Caught: NotSmoothError()
Traceback (most recent call last):
File "sage/misc/cachefunc.pyx", line 1962, in sage.misc.cachefunc.CachedMethodCaller.__call__
return cache[k]
KeyError: ((1,), ())
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
File "/home/edgarcosta/projects/five-nomial-quartics/code/debug.py", line 27, in compute_period
Y.period_matrix
File "/home/sage/sage-10.3/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/lefschetz_family/hypersurface.py", line 175, in period_matrix
self._period_matrix = block_matrix([[self.period_matrix_modification*matrix(self.homology).transpose()], [matrix(self.intersection_product*self.lift_modification(self.fibre_class))]])
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/sage/sage-10.3/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/lefschetz_family/hypersurface.py", line 153, in period_matrix_modification
integrated_thimbles = self.integrated_thimbles
^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/sage/sage-10.3/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/lefschetz_family/hypersurface.py", line 659, in integrated_thimbles
transition_matrices = self.transition_matrices
^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/sage/sage-10.3/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/lefschetz_family/hypersurface.py", line 533, in transition_matrices
rat_coefs = self.family.coordinates([self._restrict_form(w) for w in self.cohomology_internal])
^^^^^^^^^^^
File "/home/sage/sage-10.3/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/lefschetz_family/hypersurface.py", line 245, in family
self._family = Family(RtoS(self.P))
^^^^^^^^^^^^^^^^^^^^
File "/home/sage/sage-10.3/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/lefschetz_family/numperiods/family.py", line 74, in __init__
self.coho1 = self.cohomologyAt(basepoint)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "sage/misc/cachefunc.pyx", line 1967, in sage.misc.cachefunc.CachedMethodCaller.__call__
w = self._instance_call(*args, **kwds)
File "sage/misc/cachefunc.pyx", line 1842, in sage.misc.cachefunc.CachedMethodCaller._instance_call
return self.f(self._instance, *args, **kwds)
File "/home/sage/sage-10.3/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/lefschetz_family/numperiods/family.py", line 87, in cohomologyAt
return cohomology.Cohomology(self.pol(self.pol.base_ring()(t))/self.denom(self.pol.base_ring()(t)), shift=self.shift)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/sage/sage-10.3/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/lefschetz_family/numperiods/cohomology.py", line 48, in __init__
raise NotSmoothError()
lefschetz_family.exceptions.NotSmoothError
Y.ctx = <lefschetz_family.context.Context object at 0x7fafe8b58b10>
Y._P = y^4 + x^3*z + z^4 + x^2*w^2 + y*w^3
Y._dim = 2
Y._fibration = [(-7, -5, -10, 8), (-7, -7, -3, -5), (-2, -9, -5, -7)]
Y._critical_values = [-1.774316362609199? - 0.5412934630328315?*I, -1.345235496239728? - 0.6491629577976294?*I, -0.4763396119472564? - 0.9432433608841687?*I, -0.3121467129517910? - 0.4064313368762971?*I, -0.2881668639585622? - 0.3728829111687313?*I, -0.005433465395361827? - 1.018161881193612?*I, -0.001647546977624030? - 1.831305204880877?*I, 0.2783588215676724? - 1.445835998166678?*I, 1.049429649178422? - 1.282966302206489?*I, 1.337293237935781? - 1.244878165228499?*I, 3.371812801154517? - 0.08280196151425278?*I, 0.729736281072912? - 0.9778100201290040?*I, 0.7336216350561964? - 0.7034732213316668?*I, 0.10793173671158141? - 0.3055165567742057?*I, 0.989236709935345? - 0.8487005326886322?*I, 1.653884067559782? - 0.4758221693673788?*I, 2.374618950430569? - 0.1793761399560670?*I, 2.898754513324695?, 2.374618950430569? + 0.1793761399560670?*I, 1.653884067559782? + 0.4758221693673788?*I, 0.989236709935345? + 0.8487005326886322?*I, 0.7336216350561964? + 0.7034732213316668?*I, 0.10793173671158141? + 0.3055165567742057?*I, 0.729736281072912? + 0.9778100201290040?*I, 3.371812801154517? + 0.08280196151425278?*I, 1.337293237935781? + 1.244878165228499?*I, 1.049429649178422? + 1.282966302206489?*I, 0.2783588215676724? + 1.445835998166678?*I, -0.001647546977624030? + 1.831305204880877?*I, -0.005433465395361827? + 1.018161881193612?*I, -0.10137295341624369?, -0.2881668639585622? + 0.3728829111687313?*I, -0.3121467129517910? + 0.4064313368762971?*I, -0.4763396119472564? + 0.9432433608841687?*I, -1.345235496239728? + 0.6491629577976294?*I, -1.774316362609199? + 0.5412934630328315?*I]
Y._basepoint = -7
Y._fundamental_group = <lefschetz_family.voronoi.FundamentalGroupVoronoi object at 0x7fafe8b94990>
Y._Hypersurface__restriction_variable = 0
Y._fibre = Hypersurface of dimension 1 and degree 4
1 ----------------------------------------------------------------------------------------------------
2 ----------------------------------------------------------------------------------------------------
2 2025-11-19 12:38:47.779909-05:00
2 hs 0:02:50.974552
2 4 days, 0:47:04.620555 2025-11-23 13:25:52.400464-05:00
2 ----------------------------------------------------------------------------------------------------
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment