Ticket #9958: trac_9958-fix_transcendental.patch

File trac_9958-fix_transcendental.patch, 1.0 KB (added by fbissey, 10 years ago)

fixing numerical noise in transcendental.py- split from numerical noise pt1 - updated 2011092

  • sage/functions/transcendental.py

    # HG changeset patch
    # User Francois Bissey <francois.bissey@canterbury.ac.nz>
    # Date 1314647533 -43200
    # Node ID 846f00ceda1ad84481fae9813219a48e06ae54eb
    # Parent  a48af224839e4eb2555175d15d703ecd81d0b4c1
    trac 9958: patch numerical noise in sage/functions/transcendental.py
    
    diff --git a/sage/functions/transcendental.py b/sage/functions/transcendental.py
    a b  
    7979    EXAMPLES::
    8080   
    8181        sage: exponential_integral_1(2)
    82         0.048900510708061118
     82        0.04890051070806112
    8383        sage: w = exponential_integral_1(2,4); w
    84         [0.048900510708061118, 0.0037793524098489067, 0.00036008245216265873, 3.7665622843924751e-05] # 32-bit
    85         [0.048900510708061118, 0.0037793524098489063, 0.00036008245216265873, 3.7665622843924534e-05] # 64-bit
     84        [0.04890051070806112, 0.003779352409848906..., 0.00036008245216265873, 3.7665622843924...e-05]
    8685   
    8786    IMPLEMENTATION: We use the PARI C-library functions eint1 and
    8887    veceint1.