← Back to branch summary

~ubuntu-branches/ubuntu/precise/eglibc/precise

« back to all changes in this revision

Viewing changes to sysdeps/ieee754/ldbl-128ibm/s_fmal.c

  • Committer: Package Import Robot
  • Author(s): Matthias Klose
  • Date: 2011-10-04 17:48:26 UTC
  • mfrom: (216.1.23 oneiric)
  • Revision ID: package-import@ubuntu.com-20111004174826-2cyb9ewn3ucymlsx
Tags: 2.13-20ubuntu5
https://launchpad.net/bugs/853688
libc6-dev: Don't break the current {gnat,gcj}-4.4-base versons. LP: #853688.

Show diffs side-by-side

added added

removed removed

Lines of Context:
sysdeps/ieee754/ldbl-128ibm/s_fmal.c
1
 
/* Compute x * y + z as ternary operation.
2
 
   Copyright (C) 2011 Free Software Foundation, Inc.
3
 
   This file is part of the GNU C Library.
4
 
   Contributed by David Flaherty <flaherty@linux.vnet.ibm.com>.
5
 
 
6
 
   The GNU C Library is free software; you can redistribute it and/or
7
 
   modify it under the terms of the GNU Lesser General Public
8
 
   License as published by the Free Software Foundation; either
9
 
   version 2.1 of the License, or (at your option) any later version.
10
 
 
11
 
   The GNU C Library is distributed in the hope that it will be useful,
12
 
   but WITHOUT ANY WARRANTY; without even the implied warranty of
13
 
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14
 
   Lesser General Public License for more details.
15
 
 
16
 
   You should have received a copy of the GNU Lesser General Public
17
 
   License along with the GNU C Library; if not, write to the Free
18
 
   Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
19
 
   02111-1307 USA.  */
20
 
 
21
 
#include <math.h>
22
 
#include <math_ldbl_opt.h>
23
 
 
24
 
long double
25
 
__fmal (long double x, long double y, long double z)
26
 
{
27
 
        /* An IBM long double 128 is really just 2 IEEE64 doubles, and in
28
 
         * the case of inf/nan only the first double counts. So we use the
29
 
         * (double) cast to avoid any data movement.   */
30
 
       if ((finite ((double)x) && finite ((double)y)) && isinf ((double)z))
31
 
               return (z);
32
 
 
33
 
       return (x * y) + z;
34
 
}
35
 
#ifdef IS_IN_libm
36
 
long_double_symbol (libm, __fmal, fmal);
37
 
#else
38
 
long_double_symbol (libc, __fmal, fmal);
39
 
#endif