diff --git a/IMG_5493.JPG b/IMG_5493.JPG new file mode 100644 index 0000000..f47aa3c Binary files /dev/null and b/IMG_5493.JPG differ diff --git a/resolucionlarga.ipynb b/resolucionlarga.ipynb new file mode 100644 index 0000000..c2a5380 --- /dev/null +++ b/resolucionlarga.ipynb @@ -0,0 +1,279 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": 3, + "metadata": { + "collapsed": false + }, + "outputs": [ + ], + "source": [ + "S. = PolynomialRing(QQ, order='neglex')\n", + "f=((y^2+x^3)^2 + x^5*y)*((x^2+y^3)^2 + y^5*x)\n", + "g=(y^6+x^7)*(x^6+y^7)\n", + "h = f^3 - g^2\n", + "F= z^3 - 3 * f * z + 2 * g" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "data": { + "text/plain": [ + "0" + ] + }, + "execution_count": 4, + "metadata": { + }, + "output_type": "execute_result" + } + ], + "source": [ + "h = f^3 - g^2\n", + "F.discriminant(z) - 108 * h" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "data": { + "text/plain": [ + "False" + ] + }, + "execution_count": 5, + "metadata": { + }, + "output_type": "execute_result" + } + ], + "source": [ + "D=F.discriminant(z)\n", + "load('ZetaFunctionsNewtonND.py')\n", + "ZD = ZetaFunctions(D)\n", + "ZD.is_newton_degenerated(local=True)" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Newton's polyhedron of -108*y^26 + 108*y^30 + 324*x*y^29 + 648*x^2*y^27 + 324*x^2*y^28 + 1296*x^3*y^26 + 108*x^3*y^27 + 648*x^3*y^28 + 1620*x^4*y^24 + 648*x^4*y^25 + 1944*x^4*y^27 + 1944*x^5*y^23 + 3888*x^5*y^25 + 1944*x^5*y^26 + 324*x^5*y^27 - 216*x^6*y^19 + 2160*x^6*y^21 + 324*x^6*y^22 + 7776*x^6*y^24 + 648*x^6*y^25 + 2592*x^6*y^26 + 1080*x^7*y^20 + 9720*x^7*y^22 + 3888*x^7*y^23 + 1944*x^7*y^24 + 5832*x^7*y^25 + 1620*x^8*y^18 + 11664*x^8*y^21 + 13608*x^8*y^23 + 5184*x^8*y^24 + 1296*x^8*y^25 + 324*x^9*y^17 + 12960*x^9*y^19 + 1944*x^9*y^20 + 4860*x^9*y^21 + 21384*x^9*y^22 + 1620*x^9*y^23 + 6048*x^9*y^24 + 648*x^10*y^15 + 7776*x^10*y^18 + 30132*x^10*y^20 + 9720*x^10*y^21 + 7776*x^10*y^22 + 10368*x^10*y^23 + 324*x^10*y^24 + 9720*x^11*y^16 + 6480*x^11*y^18 + 30132*x^11*y^19 + 28512*x^11*y^21 + 7776*x^11*y^22 + 2916*x^11*y^23 + 1944*x^12*y^15 + 36288*x^12*y^17 + 4860*x^12*y^18 + 19440*x^12*y^19 + 33696*x^12*y^20 + 4104*x^12*y^21 + 8424*x^12*y^22 + 3456*x^13*y^13 + 4860*x^13*y^15 + 19440*x^13*y^16 + 55728*x^13*y^18 + 12960*x^13*y^19 + 15552*x^13*y^20 + 11016*x^13*y^21 + 648*x^13*y^22 + 25164*x^14*y^14 + 25920*x^14*y^16 + 42768*x^14*y^17 + 4860*x^14*y^18 + 34992*x^14*y^19 + 6804*x^14*y^20 + 3240*x^14*y^21 + 648*x^15*y^10 + 1944*x^15*y^12 + 4860*x^15*y^13 + 58752*x^15*y^15 + 6480*x^15*y^16 + 34992*x^15*y^17 + 31104*x^15*y^18 + 5508*x^15*y^19 + 6480*x^15*y^20 + 108*x^15*y^21 + 9720*x^16*y^11 + 19440*x^16*y^13 + 25920*x^16*y^14 + 6480*x^16*y^15 + 60264*x^16*y^16 + 9720*x^16*y^17 + 15552*x^16*y^18 + 6480*x^16*y^19 + 648*x^16*y^20 + 324*x^17*y^9 + 36288*x^17*y^12 + 42768*x^17*y^14 + 34992*x^17*y^15 + 9720*x^17*y^16 + 23328*x^17*y^17 + 3888*x^17*y^18 + 1620*x^17*y^19 + 1620*x^18*y^8 + 7776*x^18*y^10 + 6480*x^18*y^11 + 4860*x^18*y^12 + 55728*x^18*y^13 + 4860*x^18*y^14 + 31104*x^18*y^15 + 15552*x^18*y^16 + 3888*x^18*y^17 + 2160*x^18*y^18 - 216*x^19*y^6 + 12960*x^19*y^9 + 30132*x^19*y^11 + 19440*x^19*y^12 + 12960*x^19*y^13 + 34992*x^19*y^14 + 5508*x^19*y^15 + 6480*x^19*y^16 + 1620*x^19*y^17 + 1080*x^20*y^7 + 1944*x^20*y^9 + 30132*x^20*y^10 + 33696*x^20*y^12 + 15552*x^20*y^13 + 6804*x^20*y^14 + 6480*x^20*y^15 + 648*x^20*y^16 + 2160*x^21*y^6 + 11664*x^21*y^8 + 4860*x^21*y^9 + 9720*x^21*y^10 + 28512*x^21*y^11 + 4104*x^21*y^12 + 11016*x^21*y^13 + 3240*x^21*y^14 + 108*x^21*y^15 + 324*x^22*y^6 + 9720*x^22*y^7 + 21384*x^22*y^9 + 7776*x^22*y^10 + 7776*x^22*y^11 + 8424*x^22*y^12 + 648*x^22*y^13 + 1944*x^23*y^5 + 3888*x^23*y^7 + 13608*x^23*y^8 + 1620*x^23*y^9 + 10368*x^23*y^10 + 2916*x^23*y^11 + 1620*x^24*y^4 + 7776*x^24*y^6 + 1944*x^24*y^7 + 5184*x^24*y^8 + 6048*x^24*y^9 + 324*x^24*y^10 + 648*x^25*y^4 + 3888*x^25*y^5 + 648*x^25*y^6 + 5832*x^25*y^7 + 1296*x^25*y^8 - 108*x^26 + 1296*x^26*y^3 + 1944*x^26*y^5 + 2592*x^26*y^6 + 648*x^27*y^2 + 108*x^27*y^3 + 1944*x^27*y^4 + 324*x^27*y^5 + 324*x^28*y^2 + 648*x^28*y^3 + 324*x^29*y + 108*x^30:\n", + " support points = [(0, 26, 0), (0, 30, 0), (1, 29, 0), (2, 27, 0), (2, 28, 0), (3, 26, 0), (3, 27, 0), (3, 28, 0), (4, 24, 0), (4, 25, 0), (4, 27, 0), (5, 23, 0), (5, 25, 0), (5, 26, 0), (5, 27, 0), (6, 19, 0), (6, 21, 0), (6, 22, 0), (6, 24, 0), (6, 25, 0), (6, 26, 0), (7, 20, 0), (7, 22, 0), (7, 23, 0), (7, 24, 0), (7, 25, 0), (8, 18, 0), (8, 21, 0), (8, 23, 0), (8, 24, 0), (8, 25, 0), (9, 17, 0), (9, 19, 0), (9, 20, 0), (9, 21, 0), (9, 22, 0), (9, 23, 0), (9, 24, 0), (10, 15, 0), (10, 18, 0), (10, 20, 0), (10, 21, 0), (10, 22, 0), (10, 23, 0), (10, 24, 0), (11, 16, 0), (11, 18, 0), (11, 19, 0), (11, 21, 0), (11, 22, 0), (11, 23, 0), (12, 15, 0), (12, 17, 0), (12, 18, 0), (12, 19, 0), (12, 20, 0), (12, 21, 0), (12, 22, 0), (13, 13, 0), (13, 15, 0), (13, 16, 0), (13, 18, 0), (13, 19, 0), (13, 20, 0), (13, 21, 0), (13, 22, 0), (14, 14, 0), (14, 16, 0), (14, 17, 0), (14, 18, 0), (14, 19, 0), (14, 20, 0), (14, 21, 0), (15, 10, 0), (15, 12, 0), (15, 13, 0), (15, 15, 0), (15, 16, 0), (15, 17, 0), (15, 18, 0), (15, 19, 0), (15, 20, 0), (15, 21, 0), (16, 11, 0), (16, 13, 0), (16, 14, 0), (16, 15, 0), (16, 16, 0), (16, 17, 0), (16, 18, 0), (16, 19, 0), (16, 20, 0), (17, 9, 0), (17, 12, 0), (17, 14, 0), (17, 15, 0), (17, 16, 0), (17, 17, 0), (17, 18, 0), (17, 19, 0), (18, 8, 0), (18, 10, 0), (18, 11, 0), (18, 12, 0), (18, 13, 0), (18, 14, 0), (18, 15, 0), (18, 16, 0), (18, 17, 0), (18, 18, 0), (19, 6, 0), (19, 9, 0), (19, 11, 0), (19, 12, 0), (19, 13, 0), (19, 14, 0), (19, 15, 0), (19, 16, 0), (19, 17, 0), (20, 7, 0), (20, 9, 0), (20, 10, 0), (20, 12, 0), (20, 13, 0), (20, 14, 0), (20, 15, 0), (20, 16, 0), (21, 6, 0), (21, 8, 0), (21, 9, 0), (21, 10, 0), (21, 11, 0), (21, 12, 0), (21, 13, 0), (21, 14, 0), (21, 15, 0), (22, 6, 0), (22, 7, 0), (22, 9, 0), (22, 10, 0), (22, 11, 0), (22, 12, 0), (22, 13, 0), (23, 5, 0), (23, 7, 0), (23, 8, 0), (23, 9, 0), (23, 10, 0), (23, 11, 0), (24, 4, 0), (24, 6, 0), (24, 7, 0), (24, 8, 0), (24, 9, 0), (24, 10, 0), (25, 4, 0), (25, 5, 0), (25, 6, 0), (25, 7, 0), (25, 8, 0), (26, 0, 0), (26, 3, 0), (26, 5, 0), (26, 6, 0), (27, 2, 0), (27, 3, 0), (27, 4, 0), (27, 5, 0), (28, 2, 0), (28, 3, 0), (29, 1, 0), (30, 0, 0)]\n", + " vertices = [(0, 26, 0), (26, 0, 0), (19, 6, 0), (6, 19, 0)]\n", + " number of proper faces = 19\n", + " Facet 0: -156 + 7*y + 6*x >= 0\n", + " Facet 1: -25 + y + x >= 0\n", + " Facet 2: -156 + 6*y + 7*x >= 0\n", + " Facet 3: x >= 0\n", + " Facet 4: y >= 0\n", + " Facet 5: z >= 0\n" + ] + } + ], + "source": [ + "ZD.give_info_newton(faces=False, cones=False, compact = False)" + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": { + "collapsed": false + }, + "outputs": [ + ], + "source": [ + "def estrato(ff, chi, pesos):\n", + " Z1 = ZetaFunctions(ff)\n", + " zetamon1 = Z1.monodromy_zeta().value()^chi\n", + " zetatop1 = chi * Z1.topological_zeta(local=True, weights=pesos).value()\n", + " return zetamon1, zetatop1" + ] + }, + { + "cell_type": "markdown", + "metadata": { + "collapsed": false + }, + "source": [ + "Veamos que la resolución larga no aporta\n", + "\n" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "data": { + "text/plain": [ + "(0, 0)" + ] + }, + "execution_count": 8, + "metadata": { + }, + "output_type": "execute_result" + } + ], + "source": [ + "M3, Z3 = estrato(z^3+x^26,-1, [5,1,1])\n", + "M9, Z9 = estrato(z^3+x^52,-1, [11,1,1])\n", + "M8, Z8 = estrato(z^3+x^26,1, [6,1,1])\n", + "M39, Z39 = estrato(z^3+x^26*y^52,1, [5,11,1])\n", + "M98, Z98 = estrato(z^3+x^52*y^26, 1, [11,6,1])\n", + "M9f, Z9f = estrato(z^3+x^50*y*z + x^52, 1, [11,1,1])\n", + "Z3 + Z39 + Z98 + Z8, Z9 + Z9f\n" + ] + }, + { + "cell_type": "markdown", + "metadata": { + "collapsed": false + }, + "source": [ + "Cálculo completo\n", + "\n" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "ename": "SyntaxError", + "evalue": "invalid syntax (3357930205.py, line 3)", + "output_type": "error", + "traceback": [ + "\u001b[0;36m Cell \u001b[0;32mIn [9], line 3\u001b[0;36m\u001b[0m\n\u001b[0;31m M17, Z17 = estrato(z**Integer(3)+, Integer(2), [Integer(2),Integer(13),Integer(1)])\u001b[0m\n\u001b[0m ^\u001b[0m\n\u001b[0;31mSyntaxError\u001b[0m\u001b[0;31m:\u001b[0m invalid syntax\n" + ] + } + ], + "source": [ + "M1, Zt1 = estrato(z^3+x^4*z^2 + x^13, -13, [2,1,1])\n", + "M1d, Zt1d = estrato(z^3+x^4*z^2 + x^13*y, 13, [2,1,1])\n", + "M17, Z17 = estrato(z^3+x^4*y^28*z^2 + x^13*y^78 , 2, [2,13,1])\n", + "M7, Z7 = estrato(z^3+x^4*z + x^13*y, -4, [13,1,1])\n", + "M7d, Z7d = estrato(z^3+x^4*z + x^13*y, -4, [13,1,1])\n", + "M7g, Z7g\n" + ] + }, + { + "cell_type": "code", + "execution_count": 55, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Newton's polyhedron of z^3 + x^4*z + x^13:\n", + " support points = [(0, 0, 3), (4, 0, 1), (13, 0, 0)]\n", + " vertices = [(0, 0, 3), (4, 0, 1), (13, 0, 0)]\n", + " number of proper faces = 15\n", + " Facet 0: z >= 0\n", + " Facet 1: -13 + 9*z + x >= 0\n", + " Facet 2: -6 + 2*z + x >= 0\n", + " Facet 3: x >= 0\n", + " Facet 4: y >= 0\n" + ] + }, + { + "data": { + "text/plain": [ + "(-1/6) * (s + 2/3)^-1 * (s + 1)^-1 * (s - 2)" + ] + }, + "execution_count": 55, + "metadata": { + }, + "output_type": "execute_result" + } + ], + "source": [ + "G=z^3+x^4*z + x^13\n", + "ZG=ZetaFunctions(G)\n", + "ZG.give_info_newton(faces=False, cones=False, compact = False)\n", + "ZG.topological_zeta(local= True, weights = [2,1,1])" + ] + }, + { + "cell_type": "code", + "execution_count": 0, + "metadata": { + "collapsed": false + }, + "outputs": [ + ], + "source": [ + ] + }, + { + "cell_type": "code", + "execution_count": 0, + "metadata": { + "collapsed": false + }, + "outputs": [ + ], + "source": [ + ] + } + ], + "metadata": { + "kernelspec": { + "argv": [ + "sage-10.0", + "--python", + "-m", + "sage.repl.ipython_kernel", + "--matplotlib=inline", + "-f", + "{connection_file}" + ], + "display_name": "SageMath 10.0", + "env": { + }, + "language": "sagemath", + "metadata": { + "cocalc": { + "description": "Open-source mathematical software system", + "priority": 10, + "url": "https://www.sagemath.org/" + } + }, + "name": "sage-10.0", + "resource_dir": "/ext/jupyter/kernels/sage-10.0" + } + }, + "nbformat": 4, + "nbformat_minor": 4 +} \ No newline at end of file