END OF AN ERA, FRACTALFORUMS.COM IS CONTINUED ON FRACTALFORUMS.ORG

int rand(in out ivec2 prnp)
{
	int	r0 = prnp.y >> 15;
	int r1 = prnp.x ^ r0;
	r0 = r1 << 17;
	prnp.y = r0 ^ r1;
	prnp.x = 69069 * prnp.x;
	r1 = prnp.x ^ prnp.y;
	return r1;
}

ivec2 prnp = ivec2(1257353,135);

int irnd = rand(prnp);
float frnd = float(rand(prnp) & 0xfffffff)/268435455.0;

m1, m2, de1, de2

blendMat(m1, m2, ratio)

// polynomial smooth_ min (k = 0.1); --> fastest
float sminpol(float a, float b, float k) {
    float h = clamp(0.5 + 0.5 * (b - a) / k, 0.0, 1.0);
    return mix(b, a, h) - k * h * (1.0 - h);
}

DE = sminpol(de1, de2, k);

float br = (min(de1, de2) / max(de1, de2)); // this gives 1 at the joints and 0 at non-smoothed areas
float m = de1<de2 ? 1 : 0; // comparing DE to check which one was closer if it was not blended
matBlendRatio = mix(m, 0.5, br); // blend between original object color and perfectly blended color (1/2 ratio) on joints