C中的双截断和数学问题

double ballsFit(double r,double l,double w,double h,double Boxthick)
{
double ballsInL,ballsInW,ballsInH;
int ballsinBox;

ballsInL= (int)((l-(2*Boxthick))/(r*2));
ballsInW= (int)((w-(2*Boxthick))/(r*2));
ballsInH= (int)((h-(2*Boxthick))/(r*2));
ballsinBox=(ballsInL*ballsInW*ballsInH);
printf("LENGTH=%f\nWidth=%f\nHight=%f\nBALLS=%d\n",ballsInL,ballsInH,ballsinBox);
return ballsinBox;
}

基本问题是浮点数学是不精确的. @H_696_5@ @H_696_5@例如,数字0.1 – 您在有问题的示例中提到的厚度值 – 不能完全表示为double.将0.1分配给变量时,存储的内容是0.1的近似值.

@H_696_5@我建议你阅读What Every Computer Scientist Should Know About Floating-Point Arithmetic.

@H_696_5@ @H_696_5@几乎可以肯定的是,被乘数(至少其中一些)不是它们的样子.如果它们恰好是4.0,4.0和2.0,那么它们的产品就是32.0.如果你打印出双打能够代表的所有数字,我很确定你会看到很多9,如3.99999999999 ……等等.因此,产品有点小于32. double-to-int转换只是简单地删除小数部分,所以最终得到31.

@H_696_5@当然,如果计算精确的话,你并不总是得到比它们更小的数字;您还可以获得比您预期的更大的数字.