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

SCM Repository

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

Diff of /branches/pure-cfg/test/zslice-k1k2.diderot

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

revision 1100, Wed May 4 20:03:36 2011 UTC revision 1101, Wed May 4 20:52:56 2011 UTC
# Line 28  Line 28 
28             if (g[2] > 0.10) {             if (g[2] > 0.10) {
29                vec3 n = normalize(g);                vec3 n = normalize(g);
30                tensor[3,3] H = -∇(∇F)(pos);                tensor[3,3] H = -∇(∇F)(pos);
31                // HEY (BUG) replace with: identity[3] - n⊗n                tensor[3,3] P = identity[3] - n⊗n;
               // when 3x3 matrix subtraction is working  
               tensor[3,3] P = identity[3] + (-1.0*n⊗n);  
32                tensor[3,3] G = (P•H•P)/|g|;                tensor[3,3] G = (P•H•P)/|g|;
33                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));
34                real k1 = (trace(G) + disc)/2.0;                real k1 = (trace(G) + disc)/2.0;

Legend:
Removed from v.1100  
changed lines
  Added in v.1101

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