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

SCM Repository

[diderot] View of /branches/charisee/ertest/test2/divergence3-x2/divergence3-x2.diderot
ViewVC logotype

View of /branches/charisee/ertest/test2/divergence3-x2/divergence3-x2.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: 913 byte(s)
load->image, det(), S*T
field#2(3)[] S3= image("../../../data/zimg112.nrrd") ⊛ bspln3;
field#2(3)[3] G3= image("../../../data/const-curl3d.nrrd") ⊛ bspln3;
field#2(3)[3] V3= image("../../../data/const-curl3d.nrrd") ⊛ bspln5;
int imgSizeX = 16;
int imgSizeY = 16;
int imgSizeZ = 16;

field#2(3)[3] F3=S3*G3;

strand tensorTest (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 pos3 = [xx,yy,zz];
        output real out=99;
    
       update
    {
        if (inside(pos3,S3))
        {
        out=∇•F3(pos3);//lhs
       //  out=((∇S3(pos3))•G3(pos3) + S3(pos3)*(∇•G3(pos3)));//rhs

        }
  
    stabilize;
  }
}

initially [ tensorTest(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