k_expl.h (260066) | k_expl.h (275819) |
---|---|
1/* from: FreeBSD: head/lib/msun/ld128/s_expl.c 251345 2013-06-03 20:09:22Z 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/ld128/s_expl.c 251345 2013-06-03 20:09:22Z 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/ld128/k_expl.h 260066 2013-12-30 00:51:25Z kargl $"); | 32__FBSDID("$FreeBSD: head/lib/msun/ld128/k_expl.h 275819 2014-12-16 09:21:56Z ed $"); |
33 34/* 35 * ld128 version of k_expl.h. See ../ld80/s_expl.c for most comments. 36 * 37 * See ../src/e_exp.c and ../src/k_exp.h for precision-independent comments 38 * about the secondary kernels. 39 */ 40 --- 276 unchanged lines hidden (view full) --- 317 expt += k - 16382; 318 319 scale1 = 1; 320 half_expt = expt / 2; 321 SET_LDBL_EXPSIGN(scale1, BIAS + half_expt); 322 scale2 = 1; 323 SET_LDBL_EXPSIGN(scale1, BIAS + expt - half_expt); 324 | 33 34/* 35 * ld128 version of k_expl.h. See ../ld80/s_expl.c for most comments. 36 * 37 * See ../src/e_exp.c and ../src/k_exp.h for precision-independent comments 38 * about the secondary kernels. 39 */ 40 --- 276 unchanged lines hidden (view full) --- 317 expt += k - 16382; 318 319 scale1 = 1; 320 half_expt = expt / 2; 321 SET_LDBL_EXPSIGN(scale1, BIAS + half_expt); 322 scale2 = 1; 323 SET_LDBL_EXPSIGN(scale1, BIAS + expt - half_expt); 324 |
325 return (cpackl(cos(y) * exp_x * scale1 * scale2, | 325 return (CMPLXL(cos(y) * exp_x * scale1 * scale2, |
326 sinl(y) * exp_x * scale1 * scale2)); 327} 328#endif /* _COMPLEX_H */ | 326 sinl(y) * exp_x * scale1 * scale2)); 327} 328#endif /* _COMPLEX_H */ |