Экспоненциальное сглаживание фрактала Ньютона

Я пишу себе генератор фракталов Ньютона. Все изображения выглядели так:

Но на самом деле мне бы хотелось, чтобы это выглядело немного более плавно — наверняка я провел некоторое исследование и переехал http://www.hiddendimension.com/FractalMath/Convergent_Fractals_Main.html и это выглядит довольно правильно, за исключением того, что на краях бассейнов есть некоторые проблемы ..

Это цикл моего поколения:

while (i < 6000 && fabs(z.r) < 10000 && !found){
f = computeFunction(z, params, paramc[0]);
d = computeFunction(z, paramsD, paramc[1]);

iterexp = iterexp + exp(-fabs(z.r) - 0.5 / (fabs(subComplex(zo, z).r)));

zo = z;

z = subComplex(z, divComplex(f, d));

i++;

for (int j = 0; j < paramc[0] - 1; j++){
if (compComplex(z, zeros[j], RESOLUTION)){
resType[x + xRes * y] = j;
result[x + xRes * y] = iterexp;
found = true;
break;
}
}

if (compComplex(z, zo, RESOLUTION/100)){
resType[x + xRes * y] = 12;
break;
}
}

Окраска:

const int xRes = res[0];
const int yRes = res[1];

for (int y = 0; y < fraktal->getHeight(); y++){
for (int x = 0; x < fraktal->getWidth(); x++){
int type, it;
double conDiv;
if (genCL && genCL->err == CL_SUCCESS){
conDiv = genCL->result[x + y * xRes];
type = genCL->typeRes[x + y * xRes];
it = genCL->iterations[x + y * xRes];
} else {
type = 3;
conDiv = runNewton(std::complex<double>((double)((x - (double)(xRes / 2)) / zoom[0]), (double)((y - (double)(yRes / 2)) / zoom[1])), type);
}

if (type < 15){
Color col;
col.setColorHexRGB(colors[type]);
col.setColorHSV(col.getHue(), col.getSaturation(), 1-conDiv);
fraktal->setPixel(x, y, col);
} else {
fraktal->setPixel(x, y, conDiv, conDiv, conDiv, 1);
}
}
}

Я ценю любую помощь, чтобы на самом деле сгладить это 😉
Спасибо,
— Фодинабор