Last change
on this file since 559 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
|
Rev | Line | |
---|
[444] | 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.