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

SCM Repository

[diderot] View of /branches/charisee/ertest/vispaper/volume2-x1/volume2-x1.diderot
ViewVC logotype

View of /branches/charisee/ertest/vispaper/volume2-x1/volume2-x1.diderot

Parent Directory Parent Directory | Revision Log Revision Log


Revision 3066 - (download) (annotate)
Sat Mar 14 16:34:25 2015 UTC (4 years, 6 months ago) by cchiw
File size: 888 byte(s)
load->image, det(), S*T
field#3(2)[] f= image ("../../../data/ddro.nrrd") ⊛ bspln5;
int imgSize = 300;

field#2(2)[2] g = ∇f;
field#1(2)[2,2] H = ∇⊗∇f;
field#2(2)[] m=|g|;
field#2(2)[2] n=g/m;

//Various ways to write ∇M.
field#1(2)[2] M1=∇(g•g/2);
field#1(2)[2] M2=∇(m*m/2);
field#1(2)[2] M3= ((H•g)+(H•g))/2;

//Various ways to write C.
field#1(2)[] C1= g•M3/m;
field#1(2)[] C2 = g•((g•H+g•H)/(2*m)); //use to make vis.nrrd
field#1(2)[] C3 = g•H•g/|g|;

//field#1(2)[] D5= n•M3;//Use b.nrrd

strand fieldTest (int xi, int yi) {
    vec2 pos = [lerp(0.0, 1.0, -0.5, real(xi), real(imgSize)-0.5),
            lerp(0.0, 1.0, -0.5, real(yi), real(imgSize)-0.5)];

    //real a = alpha(2, m(pos));
    output tensor[] t2=C3(pos);

     update
        {
                  

        stabilize;
    }
}
initially [ fieldTest(xi, yi) |  yi in 0..(imgSize-1), xi in 0..(imgSize-1)];


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