Ticket #5739: trac_5739.patch

File trac_5739.patch, 1.1 KB (added by mhansen, 10 years ago)
  • sage/rings/complex_double.pyx

    # HG changeset patch
    # User Mike Hansen <mhansen@gmail.com>
    # Date 1263711833 28800
    # Node ID b5bf25c4ee6e6446ec7a903e2e12d59d4c2c6cc9
    # Parent  5c04ecef65b3f8c31d7a9200fa28284b8844244b
    Trac #5739: zeta(CDF(1)) go boom!
    
    diff --git a/sage/rings/complex_double.pyx b/sage/rings/complex_double.pyx
    a b  
    20172017       
    20182018            sage: z = CDF(1, 1)
    20192019            sage: z.zeta()
    2020             0.582158059752 - 0.926848564331*I
     2020            0.582158059752004 - 0.926848564330807*I
    20212021            sage: zeta(z)
    2022             0.582158059752 - 0.926848564331*I
     2022            0.582158059752004 - 0.926848564330807*I
     2023            sage: zeta(CDF(1))
     2024            Traceback (most recent call last):
     2025            ...
     2026            ValueError: zeta(1) pole
    20232027        """
    2024         cdef pari_sp sp
    2025         sp = avma
    2026         return self._new_from_gen_c(  gzeta(self._gen(), PREC),   sp)
     2028        import mpmath
     2029        from sage.libs.mpmath.utils import call
     2030        return CDF(call(mpmath.zeta, self))
    20272031
    20282032    def algdep(self, long n):
    20292033        """