Use mpfr function to get meaningful decimal digits in "format"

d-torrance · Nov 7, 2024 · a74365c · a74365c
1 parent 6a1ee7c
commit a74365c
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/M2/Macaulay2/d/gmp1.d b/M2/Macaulay2/d/gmp1.d
@@ -73,7 +73,8 @@ export format(
      ng := sign(x);
      if isinf(x) then return array(string)(if ng then "-" else "","infinity");
      if isnan(x) then return array(string)(if ng then "-" else "","NotANumber");
-     meaningful := int(floor(precision(x) / log2ten)) + 1;
+     meaningful := Ccode(int,
+	 "mpfr_get_str_ndigits(", base, ", ", precision(x), ")");
      if s == 0 || s > meaningful then s = meaningful; -- print at most the "meaningful" digits
      sgn := "";
      if ng then (