Home My Page Projects Code Snippets Project Openings diderot
Summary Activity Tracker Tasks SCM

SCM Repository

[diderot] Diff of /branches/pure-cfg/test/iso2d.diderot
ViewVC logotype

Diff of /branches/pure-cfg/test/iso2d.diderot

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 975, Mon Apr 25 23:18:20 2011 UTC revision 976, Tue Apr 26 00:26:08 2011 UTC
# Line 9  Line 9 
9    
10  int gridSize = 100;  int gridSize = 100;
11  real isoval = 30.0;  real isoval = 30.0;
12  image(2)[] img = load("../data/ddro-64-unit.nrrd");  //image(2)[] img = load("../data/ddro-64-unit.nrrd");
13    image(2)[] img = load("../data/tmp.nrrd");
14  field#1(2)[] F = img ⊛ ctmr;  field#1(2)[] F = img ⊛ ctmr;
15  int stepsMax = 10;  int stepsMax = 10;
16  real epsilon = 0.001;  real epsilon = 0.001;
# Line 33  Line 34 
34          vec2 grad = ∇F(pos);          vec2 grad = ∇F(pos);
35          real gmsq = grad • grad;          real gmsq = grad • grad;
36          // If |∇F|=0 we can't compute update; we have to bail          // If |∇F|=0 we can't compute update; we have to bail
37          if (gmsq < 0.00001) {          if (gmsq == 0.0) {
38              pos = [nan,nan];              pos = [nan,nan];
39              stabilize;              stabilize;
40          }          }

Legend:
Removed from v.975  
changed lines
  Added in v.976

root@smlnj-gforge.cs.uchicago.edu
ViewVC Help
Powered by ViewVC 1.0.0