This is the mail archive of the libc-hacker@sourceware.org mailing list for the glibc project.

Note that libc-hacker is a closed list. You may look at the archives of this list, but subscription and posting are not open.

Index Nav:	[Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav:	[Date Prev] [Date Next]	[Thread Prev] [Thread Next]
Other format:	[Raw text]

[PATCH] bits/mathdef.h changes for s390 and ia64, ldbl-128 fmal fix

From: Jakub Jelinek <jakub at redhat dot com>
To: Ulrich Drepper <drepper at gmail dot com>
Cc: Glibc hackers <libc-hacker at sources dot redhat dot com>
Date: Wed, 20 Oct 2010 13:29:57 +0200
Subject: [PATCH] bits/mathdef.h changes for s390 and ia64, ldbl-128 fmal fix
Reply-to: Jakub Jelinek <jakub at redhat dot com>

Hi!

This provides FP_FAST* definitions for ia64 and s390{,x} and fixes a buglet
in ldbl-128 fmal.

2010-10-20  Jakub Jelinek  <jakub@redhat.com>

	* sysdeps/ia64/bits/mathdef.h (FP_FAST_FMA, FP_FAST_FMAF,
	FP_FAST_FMAL): Define.
	* sysdeps/s390/bits/mathdef.h (FP_FAST_FMA, FP_FAST_FMAF): Likewise.

	* sysdeps/ieee754/ldbl-128/s_fmal.c (__fmal): Fix up inline asm.

--- libc/sysdeps/ia64/bits/mathdef.h
+++ libc/sysdeps/ia64/bits/mathdef.h
@@ -1,4 +1,4 @@
-/* Copyright (C) 2000, 2001, 2004 Free Software Foundation, Inc.
+/* Copyright (C) 2000, 2001, 2004, 2010 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
 
    The GNU C Library is free software; you can redistribute it and/or
@@ -34,4 +34,9 @@ typedef double double_t;	/* `double' expressions are evaluated as
 # define FP_ILOGB0	(-2147483647 - 1)
 # define FP_ILOGBNAN	2147483647
 
+/* The IA-64 architecture has fused multiply/add instructions.  */
+# define FP_FAST_FMA 1
+# define FP_FAST_FMAF 1
+# define FP_FAST_FMAL 1
+
 #endif	/* ISO C99 */
--- libc/sysdeps/ieee754/ldbl-128/s_fmal.c
+++ libc/sysdeps/ieee754/ldbl-128/s_fmal.c
@@ -176,7 +176,7 @@ __fmal (long double x, long double y, long double z)
 	u.ieee.mantissa3 |= fetestexcept (FE_INEXACT) != 0;
       v.d = a1 + u.d;
       /* Ensure the addition is not scheduled after fetestexcept call.  */
-      asm volatile ("" : "m" (v));
+      asm volatile ("" : : "m" (v));
       int j = fetestexcept (FE_INEXACT) != 0;
       feupdateenv (&env);
       /* Ensure the following computations are performed in default rounding
--- libc/sysdeps/s390/bits/mathdef.h
+++ libc/sysdeps/s390/bits/mathdef.h
@@ -1,4 +1,4 @@
-/* Copyright (C) 1997, 1998, 1999, 2000, 2004, 2006
+/* Copyright (C) 1997, 1998, 1999, 2000, 2004, 2006, 2010
    Free Software Foundation, Inc.
    This file is part of the GNU C Library.
 
@@ -35,4 +35,8 @@ typedef double double_t;	/* `double' expressions are evaluated as
 # define FP_ILOGB0	(-2147483647)
 # define FP_ILOGBNAN	2147483647
 
+/* The S/390 architecture has fused multiply/add instructions.  */
+# define FP_FAST_FMA 1
+# define FP_FAST_FMAF 1
+
 #endif	/* ISO C99 */

	Jakub

Index Nav:	[Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav:	[Date Prev] [Date Next]	[Thread Prev] [Thread Next]