Description: Support ARM64 ("aarch64") platform
Author: Jakub Bogusz
--- a/scl.c	2015/01/03 20:50:15	1.142
+++ b/scl.c	2016/08/10 02:33:05	1.143
@@ -166,33 +166,34 @@
     /* try first with starved precision */
     {
       num = scm_ash(mant, MAKINUM(e2 - point));
-      bigrecy(mant);
+      if (num != mant) bigrecy(mant);
       quo = scm_round_quotient(num, VELTS(pows5)[(long) point]);
       if (pmantexp2dbl(quo, point) != f) {
-	bigrecy(quo); quo = num;
+	if (quo != num) { bigrecy(quo); quo = num; }
 	num = scm_ash(quo, MAKINUM(1L));
-	bigrecy(quo);
+	if (num != quo) bigrecy(quo);
 	quo = scm_round_quotient(num, VELTS(pows5)[(long) --point]);
       }
+      if (num != quo) bigrecy(num);
     }
   } else {   /* e2 <= 0 */
     /* try first with starved precision */
     {
       SCM den = scm_ash(MAKINUM(1L), MAKINUM(point - e2));
       num = product(mant, VELTS(pows5)[- (long) point]);
-      bigrecy(mant);
+      if (num != mant) { bigrecy(mant); if (mant != MAKINUM(1)) mant = INUM0; }
       quo = scm_round_quotient(num, den);
       if (pmantexp2dbl(quo, point) != f) {
-	bigrecy(quo); quo = num;
+	if (quo != num) { bigrecy(quo); quo = num; }
 	point--;
 	num = product(quo, MAKINUM(10));
 	if (mant != MAKINUM(1)) bigrecy(quo);
 	quo = scm_round_quotient(num, den);
-      }
+	if (quo != num) bigrecy(num);
+      } else if ((mant != MAKINUM(1)) && (quo != num)) bigrecy(num);
       bigrecy(den);
     }
   }
-  bigrecy(num);
   a[ch++] = '.';
   /* if (sizeof(UBIGLONG)>=sizeof(double)) /\* Is ulong larger than mantissa? *\/ */
   /*   ch += iulong2str(num2ulong(quo, (char *)ARG1, s_number2string), 10, &a[ch]); */
@@ -680,7 +681,7 @@
       ans = ldexp(ans, bex + point);
     }
     if (num != quo) bigrecy(quo);
-    if (bmant != MAKINUM(1L)) bigrecy(num);
+    if ((num != bmant) && (bmant != MAKINUM(1L))) bigrecy(num);
     return ans;
   } else {
     int maxpow = LENGTH(pows5) - 1;
@@ -696,7 +697,7 @@
       quo = scm_round_quotient(num, scm_ash(scl, MAKINUM(1L)));
     }
     if (-point > maxpow) bigrecy(scl);
-    bigrecy(num);
+    if (num != quo) bigrecy(num);
     ans = ldexp(int2dbl(quo), bex + point);
     bigrecy(quo);
     return ans;
@@ -2478,7 +2479,7 @@
 	tmp = acc;
 	acc = product(tmp, tz1);
 #ifdef BIGDIG
-	bigrecy(tmp);
+	if (acc != tmp) bigrecy(tmp);
 #endif
 	break;
       }
@@ -2486,7 +2487,7 @@
 	tmp = acc;
 	acc = product(tmp, tz1);
 #ifdef BIGDIG
-	bigrecy(tmp);
+	if (acc != tmp) bigrecy(tmp);
 #endif
       }
       tmp = tz1;
@@ -2841,8 +2842,8 @@
       bex += j * BITSPERDIG;
       if (bex > 0) ans = ldexp(ans, bex);
     }
+    if (quo != num) bigrecy(quo);
     if (num != b) bigrecy(num);
-    if (quo != b) bigrecy(quo);
     if (tc16_bigneg==TYP16(b)) return -ans;
     return ans;
   }