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

SCM Repository

[diderot] Diff of /branches/pure-cfg/test/vr-curv-quad.diderot
ViewVC logotype

Diff of /branches/pure-cfg/test/vr-curv-quad.diderot

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

revision 1071, Tue May 3 21:29:54 2011 UTC revision 1078, Wed May 4 05:05:07 2011 UTC
# Line 58  Line 58 
58                vec3 grad = -∇F(pos);                vec3 grad = -∇F(pos);
59                vec3 norm = normalize(grad);                vec3 norm = normalize(grad);
60                // begin curvature computation                // begin curvature computation
61                tensor[3,3] H = -∇(∇F)(pos);                tensor[3,3] H = ∇(∇F)(pos);
62                tensor[3,3] P = identity[3] - norm⊗norm;                tensor[3,3] P = identity[3] - norm⊗norm;
63                tensor[3,3] G = (P•H•P)/|grad|;                tensor[3,3] G = -(P•H•P)/|grad|;
64                real disc = max(0.0, sqrt(2.0*|G|^2 - trace(G)^2));                real disc = max(0.0, sqrt(2.0*|G|^2 - trace(G)^2));
65                real k1 = (trace(G) + disc)/2.0;                real k1 = (trace(G) + disc)/2.0;
66                real k2 = (trace(G) - disc)/2.0;                real k2 = (trace(G) - disc)/2.0;

Legend:
Removed from v.1071  
changed lines
  Added in v.1078

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