"双倍"不打印小数

Question

我想知道为什么在这个程序中，"pi_estimated"不会打印为带有小数位的数字，尽管变量被声明为"double"。但是，它会打印出一个整数。

double get_pi(double required_accuracy)
{
    double pi_estimation=0.0;
    int x,y;
    double p=0.0,q=0.0,r=0.0;
    int D=0;
    for(int N=1;N<=1e2;N++)
    {
        x = rand()%100;
        p = (x/50.0 - 1.0)/100.0;
        y = rand()%100;
        q = (y/50.0 - 1.0)/100.0;
        r = p*p + q*q;
        if((sqrt(r))<1.0)
        {
            D++;
           pi_estimation = 4.0*(double (D/N));
        }
        if(double (4/(N+1)) < (required_accuracy*pi_estimation/100.0))
        {
            cout<<pi_estimation<<endl;
            return (pi_estimation);
        }
    }
}
int main()
{
    double pi_approx=0.0, a, actual_accuracy=0.0;
    for(a=0.1;a>=1e-14;a/=10)
    {
        pi_approx = get_pi(a);
        actual_accuracy = (fabs((pi_approx - M_PI)/(M_PI)))*100.0;
        cout<<actual_accuracy<<endl;
    }
}

Answer 1

这一行是罪魁祸首：

pi_estimation = 4.0*(double (D/N));

由于D和N都是int s，因此D/N是一个int。将int转换为double并不能神奇地使小数突然出现。

这是固定的行：

pi_estimation = 4.0 * (((double) D) / N));

你也可以先乘法，这样你就不需要那么多参数了：

pi_estimation = 4.0 * D / N;

D乘以4.0，所以它变成了double，因为double * int = double。然后除以N.由于(x * y) / z === x * (y / z)(关联属性(，因此表达式是等效的。

Answer 2

问题就在这里：

pi_estimation = 4.0*(double (D/N));

D 和 N 都是整数，因此D/N是您要转换为双精度数然后乘以 4.0 的整数。

您要执行此操作：

pi_estimation = 4.0 * (static_cast<double>(D) / N));

Answer 3

由于 D

和 N 都是整数类型，因此 D/N 以整数算术执行;转换到双倍发生得太晚，因为在转换之前丢失了精度。

一种解决方法是编写4.0 * D / N。这将确保一切都以浮点计算。(由于*和/具有相同的优先级，因此无需编写(double)。