[Bug math/15734] New: IBM long double erfl and erfcl very inaccurate

[Bug math/15734] New: IBM long double erfl and erfcl very inaccurate

[amodra at gmail dot com]

http://sourceware.org/bugzilla/show_bug.cgi?id=15734

            Bug ID: 15734
           Summary: IBM long double erfl and erfcl very inaccurate
           Product: glibc
           Version: 2.18
            Status: NEW
          Severity: normal
          Priority: P2
         Component: math
          Assignee: unassigned at sourceware dot org
          Reporter: amodra at gmail dot com

The IBM long double implementation of erfl() and erfcl() both give inaccurate
results when given a negative argument such that the component doubles have
differing sign.  Also, erfl() with a large negative argument returns +1.0
rather than -1.0.

The following compiled with -O0 shows one such bad result.  (Compare with
results from -O2 where gcc evalates erfl and erfl.)

#include <stdio.h>
#include <math.h>

int
main (void)
{
  long double x, z;
  union { long double ld; double d[2]; } u;

  x = -0x1.fffffffffffff8p-2L;
  z = erfl (x);

  u.ld = x;
  printf ("%24.21Lf %36La %21a %21a\n", u.ld, u.ld, u.d[0], u.d[1]);
  u.ld = z;
  printf ("%24.21Lf %36La %21a %21a\n", u.ld, u.ld, u.d[0], u.d[1]);

  z = erfcl (x);
  u.ld = z;
  printf ("%24.21Lf %36La %21a %21a\n", u.ld, u.ld, u.d[0], u.d[1]);

  return 0;
}

-- 
You are receiving this mail because:
You are on the CC list for the bug.

[Bug math/15734] IBM long double erfl and erfcl very inaccurate

[amodra at gmail dot com]

http://sourceware.org/bugzilla/show_bug.cgi?id=15734

Alan Modra <amodra at gmail dot com> changed:

           What    |Removed                     |Added
----------------------------------------------------------------------------
             Status|NEW                         |ASSIGNED
           Assignee|unassigned at sourceware dot org   |amodra at gmail dot com

-- 
You are receiving this mail because:
You are on the CC list for the bug.

[Bug math/15734] IBM long double erfl and erfcl very inaccurate

[amodra at gmail dot com]

https://sourceware.org/bugzilla/show_bug.cgi?id=15734

Alan Modra <amodra at gmail dot com> changed:

           What    |Removed                     |Added
----------------------------------------------------------------------------
             Status|ASSIGNED                    |RESOLVED
                URL|                            |http://sourceware.org/ml/li
                   |                            |bc-alpha/2013-08/msg00084.h
                   |                            |tml
         Resolution|---                         |FIXED
   Target Milestone|---                         |2.19

--- Comment #1 from Alan Modra <amodra at gmail dot com> ---
Fixed

-- 
You are receiving this mail because:
You are on the CC list for the bug.

[Bug math/15734] IBM long double erfl and erfcl very inaccurate

[fweimer at redhat dot com]

https://sourceware.org/bugzilla/show_bug.cgi?id=15734

Florian Weimer <fweimer at redhat dot com> changed:

           What    |Removed                     |Added
----------------------------------------------------------------------------
              Flags|                            |security-

-- 
You are receiving this mail because:
You are on the CC list for the bug.