* math/rem_pio2q.c (__quadmath_kernel_rem_pio2): Fix up fq to y

	conversion for prec 3 and __FLT_EVAL_METHOD__ != 0.

git-svn-id: svn+ssh://gcc.gnu.org/svn/gcc/branches/gcc-4_6-branch@177000 138bc75d-0d04-0410-961f-82ee72b054a4

diff --git a/libquadmath/ChangeLog b/libquadmath/ChangeLog
index a8099a6..fd91570 100644 (file)
--- a/libquadmath/ChangeLog
+++ b/libquadmath/ChangeLog
@@ -1,3 +1,8 @@
+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-06-27  Release Manager
 
        * GCC 4.6.1 released.

diff --git a/libquadmath/math/rem_pio2q.c b/libquadmath/math/rem_pio2q.c
index a1c0b83..47ee8ef 100644 (file)
--- a/libquadmath/math/rem_pio2q.c
+++ b/libquadmath/math/rem_pio2q.c
@@ -282,14 +282,20 @@ recompute:
                break;
            case 3:     /* painful */
                for (i=jz;i>0;i--) {
-                   fw      = fq[i-1]+fq[i];
-                   fq[i]  += fq[i-1]-fw;
-                   fq[i-1] = fw;
+#if __FLT_EVAL_METHOD__ != 0
+                   volatile
+#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--) {
-                   fw      = fq[i-1]+fq[i];
-                   fq[i]  += fq[i-1]-fw;
-                   fq[i-1] = fw;
+#if __FLT_EVAL_METHOD__ != 0
+                   volatile
+#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];
                if(ih==0) {