с плавающей точкой — странная ошибка потолка в переполнении стека

Здравствуйте, у меня есть некоторые трудности с функцией ceil в c ++:

У меня есть регулярная сетка точек, и мне нужно выполнить интерполяцию по ней, чтобы вычислить z-значения набора точек.

Чтобы сделать это, мне нужно для каждой вычисленной точки получить ближайшие точки на сетке. Я делаю это так:

y1 = dy*floor(p.y/dy);
y2 = dy*ceil(p.y/dy);

где dy — пространство между двумя точками сетки. (y1, p.y и dy двойные)
Если я отображу результаты с помощью

cout << static_cast<double>(p.y/dy) << ": " << y1 << ", " << y2 << endl;

Я получил эти странные результаты:

0: 0, 0
1: 0.1, 0.1
2: 0.2, 0.2
3: 0.3, 0.4

Три первых результата в порядке, но последний неверен и делает утверждение ошибочным.

Хотелось бы узнать, откуда появилась эта странная ошибка и как ее избежать.
Благодарю.

Я прошу прощения за мой английский

РЕДАКТИРОВАТЬ

Я вызываю функцию с dy = 0.1, но во время выполнения она принимает следующее значение dy = 0.10000000000000001.
p.y инициализируется так:

 const uint N = round((x2 - x1) / dx2);
const uint M = round((y2 - y1) / dy2);

double p = persistence;
double n = number_of_octaves;

// generation of the points where the perlin noise is generated
std::vector<Vertex3d> ret;
std::vector<Vertex3d> dummy;
for (uint i=0;i<=N;++i)
{
for (uint j=0;j<=M;++j)
{
ret.push_back({.x = i*dx2, .y=j*dy2, .z=0});
dummy.push_back({.x = i*dx2, .y=j*dy2, .z=0});
}
}

где x1 = 0 и x2 = 1 (по данным gdb)