// lic-turb.diderot
// 
// demo of line integral convolution (LIC) on 2D turbulent flow
//
// process output with:
// unu reshape -i lic-turb2d.txt -s 3 1020 561 | unu quantize -b 8  -o lic-turb2d.ppm

int imgSizeX = 500;
int imgSizeY = 500;
real h = 0.005;    // step size of integration
int stepNum = 10; // take this many steps both upstream and downstream

// field#2(2)[2] V = load("imp-u.nrrd") ⊛ bspln3;
field#1(2)[2] V = load("imp-u.nrrd") ⊛ ctmr;
field#0(2)[] R = load("../../data/R.nrrd") ⊛ tent;

strand LIC (int xi, int yi) {
	real xx = lerp(0.0, 8.0, 0.0, real(xi), real(imgSizeX)-1.0);
	real yy = lerp(0.0, 8.0, 0.0, real(yi), real(imgSizeY)-1.0);
	vec2 pos0 = [xx,yy];
	output vec3 RGB = [0.0, 0.0, 0.0];

	update {
                tensor[2,2] J = ∇⊗V(pos0);
                RGB = [J[0,0], J[0,1], J[0,0]];
                stabilize;
	}
}

initially [ LIC(xi, yi) | yi in 0..(imgSizeY-1), xi in 0..(imgSizeX-1) ];

<script type="text/javascript">//<![CDATA[
		function toggle_ffErrors() {
			var errorsblock = document.getElementById("ffErrorsBlock");
			var errorsgroup = document.getElementById("ffErrors");
			if (errorsblock.style.display == "none") {
				errorsblock.style.display = "block";
				errorsgroup.style.right = "10px";
			} else {
				errorsblock.style.display = "none";
				errorsgroup.style.right = "300px";
			}
		}
//]]></script>
<div id="ffErrors">
<a href="javascript:toggle_ffErrors();">Click to toggle</a>
<div id="ffErrorsBlock">
	<div class="error">does not end with &lt;/html&gt; tag</div>
	<div class="error">does not end with &lt;/body&gt; tag</div>
	<div class="info">The output has ended thus:           stabilize;
	}
}

initially [ LIC(xi, yi) | yi in 0..(imgSizeY-1), xi in 0..(imgSizeX-1) ];</div>
</div></div>
</body></html>