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

SCM Repository

[diderot] View of /branches/charisee/ertest/test2/traceField3/traceField3.diderot
ViewVC logotype

View of /branches/charisee/ertest/test2/traceField3/traceField3.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: 738 byte(s)
load->image, det(), S*T
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;


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 posv3 =[xx,yy,zz];
    output real out=99;
       update
    {
    if(inside(posv3,G3) && inside(posv3,V3)){

                out=trace(G3⊗V3)(posv3)+1;//lhs
           

        }
  
    stabilize;
  }
  
}


initially [ tensorTest(xi, yi, zi) | zi in 0..15, yi in 0..15, xi in 0..15];

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