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

SCM Repository

[diderot] View of /branches/charisee/ertest/test2/curl3-x5/curl3-x5.diderot
ViewVC logotype

View of /branches/charisee/ertest/test2/curl3-x5/curl3-x5.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: 939 byte(s)
load->image, det(), S*T
field#2(3)[] s3 = image("../../../data/vfrhand-nohip.nhdr")⊛ bspln3;
field#2(3)[] t3 = bspln3 ⊛ image("../../../data/txs.nrrd");
field#2(3)[3] V3= image("../../../data/const-curl3d.nrrd") ⊛ bspln3;
field#2(3)[3] G3= 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(0.3, 6.75, 0.0, real(xi), real(imgSizeX)-1.0);
    real yy = lerp(0.6, 3.69, 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 tensor[3] out=[99,99,99];
    
       update
    {
        if (inside(posv3,V3))
        { if (inside(posv3,G3))
        {
            out=∇×(G3×V3)(posv3);//lhs
            }
        }
  
    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