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. E no repl.itno repl.it. Também coloquei no GitHub para referência futura.
Fonte. Note que a aproximação é arbitrária.