From 73c870ed3209b68b5c8c350534508cc9d95a6bcb Mon Sep 17 00:00:00 2001 From: Szabolcs Nagy Date: Fri, 11 Apr 2014 17:57:30 +0200 Subject: math: fix aliasing violation in long double wrappers modfl and sincosl were passing long double* instead of double* to the wrapped double precision functions (on archs where long double and double have the same size). This is fixed now by using temporaries (this is not optimized to a single branch so the generated code is a bit bigger). Found by Morten Welinder. --- src/math/sincosl.c | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) (limited to 'src/math/sincosl.c') diff --git a/src/math/sincosl.c b/src/math/sincosl.c index 2c600801..d3ac1c4c 100644 --- a/src/math/sincosl.c +++ b/src/math/sincosl.c @@ -4,7 +4,10 @@ #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 void sincosl(long double x, long double *sin, long double *cos) { - sincos(x, (double *)sin, (double *)cos); + double sind, cosd; + sincos(x, &sind, &cosd); + *sin = sind; + *cos = cosd; } #elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384 void sincosl(long double x, long double *sin, long double *cos) -- cgit v1.2.1