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

SCM Repository

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

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

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

revision 3037, Wed Mar 11 13:42:30 2015 UTC revision 3038, Wed Mar 11 14:36:06 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.3037  
changed lines
  Added in v.3038

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