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

SCM Repository

[diderot] Diff of /branches/lamont/test/derivs3.diderot
ViewVC logotype

Diff of /branches/lamont/test/derivs3.diderot

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

trunk/test/derivs3.diderot revision 1119, Thu May 5 05:07:35 2011 UTC branches/lamont/test/derivs3.diderot revision 2081, Mon Nov 5 23:26:06 2012 UTC
# Line 12  Line 12 
12  // pixels due to numerical precision issues)  // pixels due to numerical precision issues)
13    
14  // F: full isotropic resolution  // F: full isotropic resolution
15  image(3)[] Fimg = load ("../data/parab/parab3-150.nrrd");  image(3)[] Fimg = image ("../data/parab/parab3-150.nrrd");
16  field#1(3)[] F0 = Fimg ⊛ ctmr;  field#1(3)[] F0 = Fimg ⊛ ctmr;
17  field#2(3)[] F = Fimg ⊛ bspln3;  field#2(3)[] F = Fimg ⊛ bspln3;
18    
19  // FX: one fifth as many samples along X  // FX: one fifth as many samples along X
20  image(3)[] FXimg = load ("../data/parab/parab3-x30.nrrd");  image(3)[] FXimg = image ("../data/parab/parab3-x30.nrrd");
21  field#1(3)[] F0X = FXimg ⊛ ctmr;  field#1(3)[] F0X = FXimg ⊛ ctmr;
22  field#2(3)[] FX = FXimg ⊛ bspln3;  field#2(3)[] FX = FXimg ⊛ bspln3;
23    
24  // FY: one fifth as many samples along Y  // FY: one fifth as many samples along Y
25  image(3)[] FYimg = load ("../data/parab/parab3-y30.nrrd");  image(3)[] FYimg = image ("../data/parab/parab3-y30.nrrd");
26  field#1(3)[] F0Y = FYimg ⊛ ctmr;  field#1(3)[] F0Y = FYimg ⊛ ctmr;
27  field#2(3)[] FY = FYimg ⊛ bspln3;  field#2(3)[] FY = FYimg ⊛ bspln3;
28    
29  // FZ: one fifth as many samples along Z  // FZ: one fifth as many samples along Z
30  image(3)[] FZimg = load ("../data/parab/parab3-z30.nrrd");  image(3)[] FZimg = image ("../data/parab/parab3-z30.nrrd");
31  field#1(3)[] F0Z = FZimg ⊛ ctmr;  field#1(3)[] F0Z = FZimg ⊛ ctmr;
32  field#2(3)[] FZ = FZimg ⊛ bspln3;  field#2(3)[] FZ = FZimg ⊛ bspln3;
33    
# Line 54  Line 54 
54          // val = [|F0X(p)-f|, |F0Y(p)-f|, |F0Z(p)-f|];          // val = [|F0X(p)-f|, |F0Y(p)-f|, |F0Z(p)-f|];
55    
56          // (2) These are magnitudes of the errors in the gradients          // (2) These are magnitudes of the errors in the gradients
57          val = [|∇FX(p)-g|, |∇FY(p)-g|, |∇FZ(p)-g|];          //val = [|∇FX(p)-g|, |∇FY(p)-g|, |∇FZ(p)-g|];
58    
59          // (3) Magnitudes of errors in Hessians          // (3) Magnitudes of errors in Hessians
60          // val = [|∇(∇FX)(p)-h|, |∇(∇FY)(p)-h|, |∇(∇FZ)(p)-h|]/100.0;          val = [|∇(∇FX)(p)-h|, |∇(∇FY)(p)-h|, |∇(∇FZ)(p)-h|]/100.0;
61    
62          stabilize;          stabilize;
63      }      }

Legend:
Removed from v.1119  
changed lines
  Added in v.2081

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