粒子测试用例容差问题：EXPECT_NEAR精度解析与通用测试方法

1. 对EXPECT_NEAR的理解是否正确？ #

不正确。EXPECT_NEAR(val1, val2, tolerance)的核心逻辑是判断两个值的绝对差值是否小于等于容差，即|val1 - val2| ≤ tolerance，并非仅看小数点后五位。容差是一个绝对值阈值，和数值的量级直接相关——比如对29.9625来说，0.00001的容差允许的范围是29.96249到29.96251；但对0.00044463402来说，同样的容差允许的范围是-0.00045463402到-0.00043463402。

2. 为什么TestParticleKineticEnergy用0.00001会失败？ #

结合动能公式和浮点数特性，原因主要有两点：

• 浮点数累积误差：动能计算公式为0.5 * 质量 * 速度平方，如果你的real类型是float（仅6-7位十进制精度），多步运算后可能产生累积误差，导致实际计算出的kinetic与预期值29.9625的绝对差值超过0.00001。
• 量级差异导致容错空间不同：29.9625的数值量级远大于TestParticleForce中的力值（~0.0004），相同的绝对容差对大数值来说苛刻得多。比如力值的0.00001容差占自身的2.2%左右，而动能的0.00001容差仅占自身的0.00003%，容错空间极小，稍有计算偏差就会触发失败。

3. 如何使用通用容差对测试值进行验证？ #

可以通过以下几种方式实现通用容差验证：

• 定义统一的容差常量：在测试文件开头定义全局或局部常量，方便统一修改和复用：

  const real KINETIC_TOL = 1e-5; // 0.00001
  const real FORCE_TOL = 1e-5;
  
  TEST(TestParticle, TestParticleKineticEnergy) {
      // ...
      EXPECT_NEAR(29.9625, kinetic, KINETIC_TOL);
  }
  

• 使用相对容差（适配跨量级数值）：如果测试的数值量级差异大，绝对容差容易失效，可自己实现相对误差判断：

  // 判断相对误差是否小于等于1e-5（万分之一）
  EXPECT_TRUE(std::abs((kinetic - 29.9625) / 29.9625) <= 1e-5);
  

• 封装自定义断言宏：如果需要频繁使用相对容差，可封装成宏简化调用：

  #define EXPECT_REL_NEAR(val1, val2, rel_tol) \
      EXPECT_TRUE(std::abs((val1 - val2) / val2) <= rel_tol)
  
  // 使用示例
  EXPECT_REL_NEAR(29.9625, kinetic, 1e-5);
  

• 测试套件级统一容差：如果多个测试属于同一套件，可使用TEST_F在SetUp中初始化容差成员变量，实现全套件统一使用。

内容的提问来源于stack exchange，提问作者user366312