Last change
on this file was
444,
checked in by satin@…, 6 years ago
|
add newlib,libalmos-mkh, restructure shared_syscalls.h and mini-libc
|
File size:
1.1 KB
|
Line | |
---|
1 | /* Copyright (C) 2002 by Red Hat, Incorporated. All rights reserved. |
---|
2 | * |
---|
3 | * Permission to use, copy, modify, and distribute this software |
---|
4 | * is freely granted, provided that this notice is preserved. |
---|
5 | */ |
---|
6 | |
---|
7 | #include "fdlibm.h" |
---|
8 | |
---|
9 | #ifdef __STDC__ |
---|
10 | float fmaf(float x, float y, float z) |
---|
11 | #else |
---|
12 | float fmaf(x,y,z) |
---|
13 | float x; |
---|
14 | float y; |
---|
15 | float z; |
---|
16 | #endif |
---|
17 | { |
---|
18 | /* NOTE: The floating-point exception behavior of this is not as |
---|
19 | * required. But since the basic function is not really done properly, |
---|
20 | * it is not worth bothering to get the exceptions right, either. */ |
---|
21 | /* Let the implementation handle this. */ /* <= NONSENSE! */ |
---|
22 | /* In floating-point implementations in which double is larger than float, |
---|
23 | * computing as double should provide the desired function. Otherwise, |
---|
24 | * the behavior will not be as specified in the standards. */ |
---|
25 | return (float) (((double) x * (double) y) + (double) z); |
---|
26 | } |
---|
27 | |
---|
28 | #ifdef _DOUBLE_IS_32BITS |
---|
29 | |
---|
30 | #ifdef __STDC__ |
---|
31 | double fma(double x, double y, double z) |
---|
32 | #else |
---|
33 | double fma(x,y,z) |
---|
34 | double x; |
---|
35 | double y; |
---|
36 | double z; |
---|
37 | #endif |
---|
38 | { |
---|
39 | return (double) fmaf((float) x, (float) y, (float) z); |
---|
40 | } |
---|
41 | |
---|
42 | #endif /* defined(_DOUBLE_IS_32BITS) */ |
---|
Note: See
TracBrowser
for help on using the repository browser.