/* ../../../bin/diderotc iso3d.diderot \ && ./iso3d \ && unu crop -i iso3d.txt -min 0 0 -max 1 M | \ unu jhisto -b 600 600 -min -1.2 -1.2 -max 1.2 1.2 -t float | \ unu 1op log1p | \ unu quantize -b 8 -o iso3d.png \ && open iso3d.png */ field#2(3)[] V= bspln3 ⊛ load("parab0.nrrd"); field#2(3)[] U= bspln3 ⊛ load("parab1.nrrd"); // isosurface of V itself: works //field#1(3)[] F = V - 0.5; // isosurface of |∇V|: works //field#1(3)[] F = |∇V| - 1; // lifting vector op to constant vector and field: does not work? field#1(3)[] F = ∇V • [1,0,0]; //field#1(3)[] F = |∇V - [1,0,0]|; // lifting vector op to two fields // field#1(3)[] F = ∇V • ∇U - 1; // does work //field#1(3)[] F = |∇V - ∇U| - 0.01; // never finishes compiling? int grid = 30; int stepsMax = 10; real epsilon = 0.000001; strand RootFind(real xi,real xj, real xk) { output vec3 x = [xi,xj,xk]; int steps = 0; update { // Stop if we're no longer inside or taken too many steps. if (!inside(x, V) || steps > stepsMax) die; // subsequent expressions are undefined if |∇F| is zero if (|∇F(x)| == 0.0) die; // the Newton-Raphson step vec3 delta = normalize(∇F(x)) * F(x)/|∇F(x)|; // we've converged if the change is small enough if (|delta| < epsilon) stabilize; x -= delta; steps += 1; } } initially { RootFind(lerp(-2, 2, -0.5, ui, grid-0.5), lerp(-2, 2, -0.5, vi, grid-0.5), lerp(-2, 2, -0.5, wi, grid-0.5)) | wi in 0..(grid-1), vi in 0..(grid-1), ui in 0..(grid-1) };
Click to toggle
does not end with </html> tag
does not end with </body> tag
The output has ended thus: -2, 2, -0.5, wi, grid-0.5)) | wi in 0..(grid-1), vi in 0..(grid-1), ui in 0..(grid-1) };