因此在使用double作为键时存在一些问题std::map。
首先,NaN比它本身少的是一个问题。如果有NaN插入的机会,请使用以下方法:
struct safe_double_less {
bool operator()(double left, double right) const {
bool leftNaN = std::isnan(left);
bool rightNaN = std::isnan(right);
if (leftNaN != rightNaN)
return leftNaN
return left
}
};
但这可能过于偏执。不要,我再说一次,不要,在您传递给std::set或类似物的比较运算符中包含epsilon阈值:这将违反容器的订购要求,并导致无法预料的不确定行为。
(出于充分的理由,我在排序NaN中将大于doubles,包括在内+inf。小于doubles也可以工作)。
因此,请使用default operator
接下来,我建议您使用std::multimap或std::multiset,因为您应该期望每个查找都有多个值。您也可以使内容管理成为日常工作,而不是简单的日常事务,以增加代码的测试覆盖范围。(我很少会推荐这些容器)加上此块operator[],不建议在使用浮点键时使用。
您要使用epsilon的地方是查询容器时。代替使用直接接口,创建一个如下的帮助函数:
// works on both `const` and non-`const` associative containers:
template
auto my_equal_range( Container&& container, double target, double epsilon = 0.00001 )
-> decltype( container.equal_range(target) )
{
auto lower = container.lower_bound( target-epsilon );
auto upper = container.upper_bound( target+epsilon );
return std::make_pair(lower, upper);
}
适用于std::map和std::set(和multi版本)。
(在更现代的代码库中,我希望range>从equal_range函数返回的对象更好。但是现在,我将使其与兼容equal_range)。
这将查找一系列其关键字“足够接近”您所要求的事物的容器,而容器在内部维护其排序保证,并且不执行未定义的行为。
要测试密钥是否存在,请执行以下操作:
template
bool key_exists( Container const& container, double target, double epsilon = 0.00001 ) {
auto range = my_equal_range(container, target, epsilon);
return range.first != range.second;
}
并且如果您要删除/替换条目,则应处理可能有多个条目命中的可能性。
较短的答案是“不使用浮动点值作为键std::set和std::map”,因为它是一个有点麻烦。
如果确实对std::set或使用浮点键std::map,则几乎可以肯定不要对它们执行.find或[],因为极有可能是错误的来源。只要确切的顺序无关紧要(即,一个特定的1.0在前面或后面,或者与另一个1.0完全在同一位置),就可以将其用于自动排序的东西集合。即使那样,我还是要使用多贴图/多集,因为依靠碰撞或缺乏碰撞并不是我要依靠的东西。
很难对IEEE浮点值的确切值进行推理,并且依赖它的代码的脆弱性很常见。