From mboxrd@z Thu Jan  1 00:00:00 1970
Return-Path: <glibc-bugs-return-18833-listarch-glibc-bugs=sources.redhat.com@sourceware.org>
Received: (qmail 11896 invoked by alias); 3 Jun 2013 17:01:18 -0000
Mailing-List: contact glibc-bugs-help@sourceware.org; run by ezmlm
Precedence: bulk
List-Id: <glibc-bugs.sourceware.org>
List-Subscribe: <mailto:glibc-bugs-subscribe@sourceware.org>
List-Post: <mailto:glibc-bugs@sourceware.org>
List-Help: <mailto:glibc-bugs-help@sourceware.org>, <http://sourceware.org/lists.html#faqs>
Sender: glibc-bugs-owner@sourceware.org
Received: (qmail 11843 invoked by uid 55); 3 Jun 2013 17:01:12 -0000
From: "joseph at codesourcery dot com" <sourceware-bugzilla@sourceware.org>
To: glibc-bugs@sourceware.org
Subject: [Bug math/15563] sincos() is incorrect for long double and large inputs on x86_64
Date: Mon, 03 Jun 2013 17:01:00 -0000
X-Bugzilla-Reason: CC
X-Bugzilla-Type: changed
X-Bugzilla-Watch-Reason: None
X-Bugzilla-Product: glibc
X-Bugzilla-Component: math
X-Bugzilla-Version: 2.18
X-Bugzilla-Keywords:
X-Bugzilla-Severity: normal
X-Bugzilla-Who: joseph at codesourcery dot com
X-Bugzilla-Status: NEW
X-Bugzilla-Priority: P2
X-Bugzilla-Assigned-To: unassigned at sourceware dot org
X-Bugzilla-Target-Milestone: ---
X-Bugzilla-Flags:
X-Bugzilla-Changed-Fields:
Message-ID: <bug-15563-131-Z7b75mNCSZ@http.sourceware.org/bugzilla/>
In-Reply-To: <bug-15563-131@http.sourceware.org/bugzilla/>
References: <bug-15563-131@http.sourceware.org/bugzilla/>
Content-Type: text/plain; charset="UTF-8"
Content-Transfer-Encoding: 7bit
X-Bugzilla-URL: http://sourceware.org/bugzilla/
Auto-Submitted: auto-generated
MIME-Version: 1.0
X-SW-Source: 2013-06/txt/msg00012.txt.bz2

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

--- Comment #1 from joseph at codesourcery dot com <joseph at codesourcery dot com> ---
I can't tell what you think the bug is here.  You give a testcase for 
double but talk about long double, and with the obvious substitutions

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

int main (void)
{
  volatile long double x = 1.0e22;
  long double s1, s2, c1;

  sincosl (x, &s1, &c1);
  s2 = sinl (x);
  printf ("s1 = %.17Lg\n", s1);
  printf ("s2 = %.17Lg\n", s2);
  return 0;
}

I get on x86_64 (with current glibc)

s1 = -0.8522008497671888
s2 = -0.8522008497671888

which appears to be what you'd expect.  So what program are you building, 
with what options, and running with current glibc, and getting what 
results, and what would you like instead?

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