以sprintf科学记数法打印所有有效数字

Question

当R用科学符号将一个大数转换为字符串时，它包括所有有效数字，没有尾随零。有没有可能在C中用Sprintf来完成这个？

> as.character(12345e11)
[1] "1.2345e+15"
> as.character(1234500000e6)
[1] "1.2345e+15"
> as.character(1234500001e6)
[1] "1.234500001e+15"

我尝试了sprintf（buf，“%g”，val），但这似乎最多包含5个十进制数字。我还尝试使用sprintf（buf，“%.18g”，val）设置更高的精度，但这将包括非有效数字和尾随零。

有没有一种方法可以让Sprintf（buf，"%g"， val）的行为增加5位数的限制？

Answer 1

代码可以使用“%.18e”或“%.18g”，但问题是“18”应该有多大？18是最佳值吗？答案在于DBL\u DECIMAL\u DIG。

DBL\u DECIMAL\u DIG是要打印的最小有效位数，以确保所有double的double字符串往返于相同的double。

建议使用格式说明符“%.*e”。

请注意，“%.18e”中的“18”是小数点后的有效位数。因此，“%.18e”打印19个有效数字。

使用printf（"%a"， x）；以十六进制输出输出。
对于十进制输出：

#include <float.h>

//    sign + digit +  dp +       digits          + e + sign + expo + \0
char buf[1 + 1 +      1  + (DBL_DECIMAL_DIG - 1) + 1 + 1    + 5    + 1]; 
sprintf(buf, "%.*e", DBL_DECIMAL_DIG - 1, x);

Ref Printf width specifier以保持浮点值的精度

使用典型的doublebinary64格式的像y=1.0/3.0这样的数字需要看到大约53位十进制数字才能看到其精确值。但成功的往返并不需要许多尾随数字。

现在我们知道了要打印的最多数字，使用下面的数字来摆脱那些讨厌的尾随0数字。

#include <float.h>
#include <math.h>
#include <stdio.h>
#include <string.h>

char *trim0(double x, char *buf) {
  sprintf(buf, "% .*e", DBL_DECIMAL_DIG - 1, x);
  if (isfinite(x)) {
    char *p = &buf[DBL_DECIMAL_DIG + 1];  // address of last significand digit
    char *t;
    for (t=p; *t == '0'; t--);
    memmove(t+1, p+1, strlen(p+1)+1);
  }
  return buf;
}

int main(void) {
  char buf[1 + 1 + 1 + (DBL_DECIMAL_DIG - 1) + 1 + 1 + 5 + 1];
  printf("%s\n", trim0(1.2, buf));
  printf("%s\n", trim0(1.0/7, buf));
  return 0;
}

输出

 1.2e+00
 1.4285714285714285e-01