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

SCM Repository

[diderot] View of /branches/charisee/test/t2.diderot
ViewVC logotype

View of /branches/charisee/test/t2.diderot

Parent Directory Parent Directory | Revision Log Revision Log


Revision 2499 - (download) (annotate)
Mon Nov 4 16:10:36 2013 UTC (6 years ago) by cchiw
File size: 1167 byte(s)
remove generic rep, change Summation index, and add tests
// t3-test translate Fields

vec3 x = [0.0, 0.0, 5.0];
vec3 y = [0.0, 1.0, 0.0];
vec3 z = [1.0, 0.0, 0.0];

image(3)[] img = load("../data/zimg112.nrrd");
field#1(2)[2] V = load("../data/turb2d.nrrd") ⊛ ctmr;
field#1(2)[] F2 = ctmr ⊛ load("../data/ddro-80.nrrd");
field#2(3)[] F = img ⊛ bspln3;
field#0(2)[3] B = tent ⊛ load("../data/txf/2d-bow.nrrd");



//Simple Field Operations

/*
field#2(3)[] G=F+1;

field#2(3)[] G2=8+F;
field#2(3)[] J=1-F;
field#2(3)[] tree=F-1;
field#2(3)[] K=F*3;
field#2(3)[] K2=4*F;
field#1(2)[2] V2=V+V;
field#0(2)[3] V3=B-B;
field#2(3)[] R=-F;



//grad
field#1(3)[3] L =∇F;

//divergence
field#0(2)[2,2] N =∇⊗V;
field#0(2)[2,2] N1=N+N;
field#0(3)[3,3] N2 =∇⊗L;

// curl
field#0(2)[3] P2 =∇×L;

//don't know divide fields
//field#2(3)[] L=F/5;

//probe
vec2 pos2=[8,9];
vec3 pos3=[8,9,10];
//K(3)[]
real p2=K(pos3);
//P (2)[]
//real p3=P(pos2);
//B(2)[3]
vec3 p1=B(pos2);
//(2)[2] V
vec2 bird=V(pos2);
tensor [2,2] peanut=N(pos2);
*/

strand fieldTest (int row)
{//vec3 u=x+y ;
real row1=row*9;
    update
    {
     
    vec2 pos=[row,row1];
   // real del=P(pos);
    }
}

initially [ fieldTest( r) | r in 0..1];

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