Ticket #2303: 2303-scipy-bessel.patch

File 2303-scipy-bessel.patch, 1.4 kB (added by robertwb, 5 months ago)
  • a/sage/functions/special.py

    old new  
    394394    INPUT: 
    395395        nu -- a real (or complex, for pari) number 
    396396        z  -- a real (positive) 
    397         algorithm - "pari" or "maxima" or "scipy" 
     397        algorithm - "pari" or "maxima" 
    398398        prec - real precision (for Pari only) 
    399399     
    400400    DEFINITION: 
     
    437437        1.32616018371265... 
    438438        sage: bessel_I(0,1,"maxima")    
    439439        1.26606587775200... 
    440         sage: bessel_I(1,1,"scipy") 
    441         0.565159103992... 
    442440 
    443441    """ 
    444442    if algorithm=="pari": 
     
    447445        b = R(pari(nu).besseli(z)) 
    448446        pari.set_real_precision(a) 
    449447        return b 
    450     elif algorithm=="scipy": 
     448    elif algorithm=="scipy" and False: 
     449        # Broken x86-Linux box (pentium4-fc6) using gcc-4.3.1 
    451450        if prec != 53: 
    452451            raise ValueError, "for the scipy algorithm the precision must be 53" 
    453452        import scipy.special