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

SCM Repository

[diderot] Diff of /branches/vis12-cl/test/iso2d-cbrt.diderot
ViewVC logotype

Diff of /branches/vis12-cl/test/iso2d-cbrt.diderot

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

revision 3225, Thu Jul 23 11:25:15 2015 UTC revision 3226, Thu Jul 23 11:33:34 2015 UTC
# Line 24  Line 24 
24          // non-finite position; these will be filtered later.          // non-finite position; these will be filtered later.
25          // This is really what "die" is for, though (see below).          // This is really what "die" is for, though (see below).
26          if (!inside(pos, F) || steps >= stepsMax) {          if (!inside(pos, F) || steps >= stepsMax) {
27              pos = [nan,nan];              pos = nan[2];
28              // HEY (BUG) we shouldn't even have to use non-finite values;              // HEY (BUG) we shouldn't even have to use non-finite values;
29              // we should be able to just say "die;".  Currently, using              // we should be able to just say "die;".  Currently, using
30              // "die" leads to the program never finishing...              // "die" leads to the program never finishing...
# Line 34  Line 34 
34          real gmsq = grad • grad;          real gmsq = grad • grad;
35          // If |∇F|=0 we can't compute update; we have to bail          // If |∇F|=0 we can't compute update; we have to bail
36          if (gmsq == 0.0) {          if (gmsq == 0.0) {
37              pos = [nan,nan];              pos = nan[2];
38              stabilize;              stabilize;
39          }          }
40          // delta = Newton-Raphson step          // delta = Newton-Raphson step

Legend:
Removed from v.3225  
changed lines
  Added in v.3226

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