Realmente não dá para fazer uma comparação precisa com ponto flutuante. Se necessita desta precisa, 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