diff options
| author | Jakub Jelinek <jakub@redhat.com> | 2010-10-13 22:27:03 -0400 |
|---|---|---|
| committer | Ulrich Drepper <drepper@gmail.com> | 2010-10-13 22:27:03 -0400 |
| commit | 5e908464b9603e9cd95047d2edd00de5f78fedaf (patch) | |
| tree | 1f526248f31bee031ce4c92d9b1b2e857d6b671b /sysdeps/s390 | |
| parent | f90681487d90a6eea8f1c60021bc8dd83313c59b (diff) | |
| download | glibc-5e908464b9603e9cd95047d2edd00de5f78fedaf.tar glibc-5e908464b9603e9cd95047d2edd00de5f78fedaf.tar.gz glibc-5e908464b9603e9cd95047d2edd00de5f78fedaf.tar.bz2 glibc-5e908464b9603e9cd95047d2edd00de5f78fedaf.zip | |
Implement accurate fma.
Diffstat (limited to 'sysdeps/s390')
| -rw-r--r-- | sysdeps/s390/fpu/s_fma.c | 37 | ||||
| -rw-r--r-- | sysdeps/s390/fpu/s_fmaf.c | 32 |
2 files changed, 69 insertions, 0 deletions
diff --git a/sysdeps/s390/fpu/s_fma.c b/sysdeps/s390/fpu/s_fma.c new file mode 100644 index 0000000000..f682362313 --- /dev/null +++ b/sysdeps/s390/fpu/s_fma.c @@ -0,0 +1,37 @@ +/* Compute x * y + z as ternary operation. S/390 version. + Copyright (C) 2010 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Jakub Jelinek <jakub@redhat.com>, 2010. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, write to the Free + Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA + 02111-1307 USA. */ + +#include <math.h> + +double +__fma (double x, double y, double z) +{ + double r; + asm ("madbr %0,%1,%2" : "=f" (r) : "%f" (x), "fR" (y), "0" (z)); + return r; +} +#ifndef __fma +weak_alias (__fma, fma) +#endif + +#ifdef NO_LONG_DOUBLE +strong_alias (__fma, __fmal) +weak_alias (__fmal, fmal) +#endif diff --git a/sysdeps/s390/fpu/s_fmaf.c b/sysdeps/s390/fpu/s_fmaf.c new file mode 100644 index 0000000000..030d7b229b --- /dev/null +++ b/sysdeps/s390/fpu/s_fmaf.c @@ -0,0 +1,32 @@ +/* Compute x * y + z as ternary operation. S/390 version. + Copyright (C) 2010 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Jakub Jelinek <jakub@redhat.com>, 2010. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, write to the Free + Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA + 02111-1307 USA. */ + +#include <math.h> + +float +__fmaf (float x, float y, float z) +{ + float r; + asm ("maebr %0,%1,%2" : "=f" (r) : "%f" (x), "fR" (y), "0" (z)); + return r; +} +#ifndef __fmaf +weak_alias (__fmaf, fmaf) +#endif |