Home My Page Projects Code Snippets Project Openings diderot

# SCM Repository

[diderot] Diff of /branches/charisee/ertest/vispaper/tiso3d/faiso3d.diderot
 [diderot] / branches / charisee / ertest / vispaper / tiso3d / faiso3d.diderot

# Diff of /branches/charisee/ertest/vispaper/tiso3d/faiso3d.diderot

revision 3109, Wed Mar 18 16:27:22 2015 UTC revision 3110, Wed Mar 18 16:32:05 2015 UTC
# Line 20  Line 20
20  */  */
21
22  field#2(3)[3,3] V = bspln3 ⊛ image("tball.nrrd");  field#2(3)[3,3] V = bspln3 ⊛ image("tball.nrrd");
23  field#2(3)[3,3] E = V - trace(V)*identity[3,3]/3;  field#2(3)[3,3] E = V - trace(V)*identity[3]/3;
24  field#2(3)[] F = sqrt(3.0/2.0)*|E|/|V| - 0.4;  field#2(3)[] F = sqrt(3.0/2.0)*|E|/|V| - 0.4;
25
26  int grid = 30;  int grid = 30;
# Line 32  Line 32
32     int steps = 0;     int steps = 0;
33     update {     update {
34        // Stop if we're no longer inside or taken too many steps.        // Stop if we're no longer inside or taken too many steps.
35        if (!inside(x, F) || steps >= stepsMax)        if (!inside(x, V) || steps >= stepsMax)
36           stabilize;           stabilize;
37        // subsequent expressions are undefined if |∇F| is zero        // subsequent expressions are undefined if |∇F| is zero
38        if (|∇F(x)| == 0.0)        if (|∇F(x)| == 0.0)

Legend:
 Removed from v.3109 changed lines Added in v.3110