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

SCM Repository

[diderot] View of /branches/charisee/ertest/tests/normalize-vec3/normalize-vec3-x1.diderot
ViewVC logotype

View of /branches/charisee/ertest/tests/normalize-vec3/normalize-vec3-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: 832 byte(s)
load->image, det(), S*T
int imgSizeX = 16;
int imgSizeY = 16;
int imgSizeZ = 16;

field#1(3)[3] V = image("../../../data/const-curl3d.nrrd") ⊛ bspln3;

//normalized helicity 

strand fieldTest (int xi, int yi,int zi) {
  	real xx = lerp(2.00, 17.0, 0.0, real(xi), real(imgSizeX)-1.0);
	real yy = lerp(2.00, 17.0, 0.0, real(yi), real(imgSizeY)-1.0);
	real zz = lerp(2.00, 17.0, 0.0, real(zi), real(imgSizeZ)-1.0);
	vec3 pos = [xx,yy,zz];

        //field#1(3)[] F=  (V/|V|)⋅(∇×V/|∇×V|); 
       // output tensor[] out=  (V(pos)/|V(pos)|)⋅(∇×V(pos)/|∇×V(pos)|);
       
       
        field#1(3)[] F=  normalize(V)⋅normalize(∇×V); 
        output tensor[] out=F(pos);
        

      update
    {
         stabilize;
 
    }
}
initially [ fieldTest(xi, yi, zi) | zi in 0..(imgSizeZ-1), yi in 0..(imgSizeY-1), xi in 0..(imgSizeX-1) ];

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