Commit 57044b5a by Jakub Jelinek Committed by Jakub Jelinek

rem_pio2q.c (__quadmath_kernel_rem_pio2): Fix up fq to y conversion for prec 3…

rem_pio2q.c (__quadmath_kernel_rem_pio2): Fix up fq to y conversion for prec 3 and __FLT_EVAL_METHOD__ != 0.

	* math/rem_pio2q.c (__quadmath_kernel_rem_pio2): Fix up fq to y
	conversion for prec 3 and __FLT_EVAL_METHOD__ != 0.

From-SVN: r176999
parent 8586635c
2011-08-01 Jakub Jelinek <jakub@redhat.com>
* math/rem_pio2q.c (__quadmath_kernel_rem_pio2): Fix up fq to y
conversion for prec 3 and __FLT_EVAL_METHOD__ != 0.
2011-04-03 Tobias Burnus <burnus@net-b.de> 2011-04-03 Tobias Burnus <burnus@net-b.de>
PR fortran/48311 PR fortran/48311
......
...@@ -282,14 +282,20 @@ recompute: ...@@ -282,14 +282,20 @@ recompute:
break; break;
case 3: /* painful */ case 3: /* painful */
for (i=jz;i>0;i--) { for (i=jz;i>0;i--) {
fw = fq[i-1]+fq[i]; #if __FLT_EVAL_METHOD__ != 0
fq[i] += fq[i-1]-fw; volatile
fq[i-1] = fw; #endif
double fv = (double)(fq[i-1]+fq[i]);
fq[i] += fq[i-1]-fv;
fq[i-1] = fv;
} }
for (i=jz;i>1;i--) { for (i=jz;i>1;i--) {
fw = fq[i-1]+fq[i]; #if __FLT_EVAL_METHOD__ != 0
fq[i] += fq[i-1]-fw; volatile
fq[i-1] = fw; #endif
double fv = (double)(fq[i-1]+fq[i]);
fq[i] += fq[i-1]-fv;
fq[i-1] = fv;
} }
for (fw=0.0,i=jz;i>=2;i--) fw += fq[i]; for (fw=0.0,i=jz;i>=2;i--) fw += fq[i];
if(ih==0) { if(ih==0) {
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment