在我的流体模拟过程中,物理时间正在进行为0,0.001,0.002,…,4.598,4.599,4.6,4.601,4.602 ……现在我想选择时间= 0.1,0.2,4.5,……从这个时间序列再做进一步的分析.所以我编写了以下代码来判断fractpart是否达到零.

但我很惊讶我发现以下两种分割方法得到两种不同的结果,我该怎么办？

double param,fractpart,intpart;
double org = 4.6;
double ddd = 0.1;
 
// This is the correct one I need. I got intpart=46 and fractpart=0
// param = org*(1/ddd);
 
// This is not what I want. I got intpart=45 and fractpart=1
param = org/ddd;
 
fractpart = modf(param,&intpart);
Info<< "\n\nfractpart\t=\t"
    << fractpart
    << "\nAnd intpart\t=\t"
    << intpart
    << endl;

为什么会这样发生？
如果你们容忍我一点点,我可以大声喊道：“C委员会可以对此做点什么吗？因为这很令人困惑.”