k_expl.h (260066) | k_expl.h (275819) |
---|---|
1/* from: FreeBSD: head/lib/msun/ld80/s_expl.c 251343 2013-06-03 19:51:32Z kargl */ 2 3/*- 4 * Copyright (c) 2009-2013 Steven G. Kargl 5 * All rights reserved. 6 * 7 * Redistribution and use in source and binary forms, with or without 8 * modification, are permitted provided that the following conditions --- 15 unchanged lines hidden (view full) --- 24 * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT 25 * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF 26 * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 27 * 28 * Optimized by Bruce D. Evans. 29 */ 30 31#include <sys/cdefs.h> | 1/* from: FreeBSD: head/lib/msun/ld80/s_expl.c 251343 2013-06-03 19:51:32Z kargl */ 2 3/*- 4 * Copyright (c) 2009-2013 Steven G. Kargl 5 * All rights reserved. 6 * 7 * Redistribution and use in source and binary forms, with or without 8 * modification, are permitted provided that the following conditions --- 15 unchanged lines hidden (view full) --- 24 * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT 25 * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF 26 * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 27 * 28 * Optimized by Bruce D. Evans. 29 */ 30 31#include <sys/cdefs.h> |
32__FBSDID("$FreeBSD: head/lib/msun/ld80/k_expl.h 260066 2013-12-30 00:51:25Z kargl $"); | 32__FBSDID("$FreeBSD: head/lib/msun/ld80/k_expl.h 275819 2014-12-16 09:21:56Z ed $"); |
33 34/* 35 * See s_expl.c for more comments about __k_expl(). 36 * 37 * See ../src/e_exp.c and ../src/k_exp.h for precision-independent comments 38 * about the secondary kernels. 39 */ 40 --- 253 unchanged lines hidden (view full) --- 294 expt += k - 16382; 295 296 scale1 = 1; 297 half_expt = expt / 2; 298 SET_LDBL_EXPSIGN(scale1, BIAS + half_expt); 299 scale2 = 1; 300 SET_LDBL_EXPSIGN(scale1, BIAS + expt - half_expt); 301 | 33 34/* 35 * See s_expl.c for more comments about __k_expl(). 36 * 37 * See ../src/e_exp.c and ../src/k_exp.h for precision-independent comments 38 * about the secondary kernels. 39 */ 40 --- 253 unchanged lines hidden (view full) --- 294 expt += k - 16382; 295 296 scale1 = 1; 297 half_expt = expt / 2; 298 SET_LDBL_EXPSIGN(scale1, BIAS + half_expt); 299 scale2 = 1; 300 SET_LDBL_EXPSIGN(scale1, BIAS + expt - half_expt); 301 |
302 return (cpackl(cos(y) * exp_x * scale1 * scale2, | 302 return (CMPLXL(cos(y) * exp_x * scale1 * scale2, |
303 sinl(y) * exp_x * scale1 * scale2)); 304} 305#endif /* _COMPLEX_H */ | 303 sinl(y) * exp_x * scale1 * scale2)); 304} 305#endif /* _COMPLEX_H */ |