OSDN Git Service

* math/rem_pio2q.c (__quadmath_kernel_rem_pio2): Fix up fq to y
authorjakub <jakub@138bc75d-0d04-0410-961f-82ee72b054a4>
Mon, 1 Aug 2011 08:55:33 +0000 (08:55 +0000)
committerjakub <jakub@138bc75d-0d04-0410-961f-82ee72b054a4>
Mon, 1 Aug 2011 08:55:33 +0000 (08:55 +0000)
conversion for prec 3 and __FLT_EVAL_METHOD__ != 0.

git-svn-id: svn+ssh://gcc.gnu.org/svn/gcc/trunk@176999 138bc75d-0d04-0410-961f-82ee72b054a4

libquadmath/ChangeLog
libquadmath/math/rem_pio2q.c

index 522d683..7b35cc9 100644 (file)
@@ -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-04-03  Tobias Burnus  <burnus@net-b.de>
 
        PR fortran/48311
index a1c0b83..47ee8ef 100644 (file)
@@ -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) {