---

/*{
	"CREDIT": "mojovideotech",
	"CATEGORIES": [
		"generator",
		"2d",
		"fractal",
		"space"
	],
	"DESCRIPTION": "SimplyNebulous",
	"INPUTS": [
		{
			"NAME": "rate",
			"TYPE": "float",
			"DEFAULT": 2.67,
			"MIN": -5,
			"MAX": 5
		},
		{
			"NAME": "rot",
			"TYPE": "float",
			"DEFAULT": 2.82,
			"MIN": 0,
			"MAX": 6.28318
		},
		{
			"NAME": "flash",
			"TYPE": "float",
			"DEFAULT": 0.04,
			"MIN": 0.001,
			"MAX": 1
		},
		{
			"NAME": "zoom",
			"TYPE": "float",
			"DEFAULT": 5.5,
			"MIN": 1,
			"MAX": 10
		},
		{
			"NAME": "morph",
			"TYPE": "float",
			"DEFAULT": 0.64,
			"MIN": -5,
			"MAX": 5
		},
		{
			"NAME": "detail",
			"TYPE": "float",
			"DEFAULT": 19,
			"MIN": 1,
			"MAX": 24
		},
		{
			"NAME": "depth",
			"TYPE": "float",
			"DEFAULT": 12,
			"MIN": 7,
			"MAX": 23
		},
		{
			"NAME": "brightness",
			"TYPE": "float",
			"DEFAULT": 1.9,
			"MIN": 1,
			"MAX": 3
		},
		{
			"NAME": "multiplier",
			"TYPE": "float",
			"DEFAULT": 1,
			"MIN": 1,
			"MAX": 5
		},
		{
			"NAME": "red",
			"TYPE": "float",
			"DEFAULT": 0.2,
			"MIN": 0.1,
			"MAX": 1.5
		},
		{
			"NAME": "green",
			"TYPE": "float",
			"DEFAULT": 0.32,
			"MIN": 0.1,
			"MAX": 1.5
		},
		{
			"NAME": "blue",
			"TYPE": "float",
			"DEFAULT": 0.85,
			"MIN": 0.1,
			"MAX": 1.5
		}
	]
}*/

///////////////////////////////////////////
// SimplyNebulous  by mojovideotech
//
// mod of:
// shadertoy.com\/lslGWr  by Josh P
//
// based on :
// www.fractalforums.com/\new-theories-and-research/\very-simple-formula-for-fractal-patterns
//
// Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License.
///////////////////////////////////////////


mat2 rmat(float t) {
    float c = cos(t), s = sin(t);
    return mat2(c,s,-s,s-c);
}

float field(in vec3 p) {
	float strength = 9.0 + flash * log(1.e-6 + fract(sin(TIME) * 4373.11));
	float accum = 0.0, prev = 0.0, tw = 0.0, b = -5.0;
	for (int i = 0; i < 26; ++i) {
		float mag = dot(p, p);
		p = abs(p) / mag + vec3(-0.5, -0.4, - 1.5);
		float w = exp(-float(i) / depth);
		accum += w * exp(-strength * pow(abs(mag - prev), brightness));
		tw += w;
		prev = mag;
		b += 1.0;
		if (b - detail >= 1.0) break;
	}
	return max(0.0, 5.0 * accum / tw - 0.7);
}

void main() {
	float TT = TIME * rate;
	vec2 uv = 2.0 * isf_FragNormCoord.xy - 1.0;
	vec2 uvs = uv * RENDERSIZE.xy / max(RENDERSIZE.x, RENDERSIZE.y) ;
	vec3 p = vec3(uvs / zoom ,morph) + vec3(1.0, -1.3, -0.5);
	p.xz *= rmat(rot);
	float mu = floor(multiplier);
	p += 0.2 * vec3(sin(TT / 13.0 * mu), sin(TT / 89.0 * mu),  sin(TT / 233.0 * mu));
	float t = field(p);
	float v = (1.0 - exp(abs(uv.x) - 1.0) * 5.0 ) * (1.0 - exp(abs(uv.y) - 1.0));
	vec3 col = mix(2.0,0.1, v) * vec3(1.6 * t * t * t, 1.3 * t * t, 1.1 * t);
	col *= vec3(red,green,blue);
	gl_FragColor = vec4(col, 1.0);
}