Realmente não dá para fazer uma comparação precisa com ponto flutuante. Se necessita desta precisão use outro tipo. Mas se quiser minimizar o problema, faça uma função que realize uma aproximação. Exemplo: bool isEqual(double x, double y) { return std::abs(x - y) <= 1e-5 * std::abs(x); } Veja [funcionando no ideone][1]. [Fonte][2]. Note que a aproximação é arbitrária. [1]: http://ideone.com/NqvONj [2]: https://isocpp.org/wiki/faq/newbie#floating-point-arith