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

SCM Repository

[diderot] View of /branches/charisee/test/E-fRewrites.diderot
ViewVC logotype

View of /branches/charisee/test/E-fRewrites.diderot

Parent Directory Parent Directory | Revision Log Revision Log

Revision 2517 - (download) (annotate)
Thu Dec 12 20:35:28 2013 UTC (6 years, 10 months ago) by cchiw
File size: 1243 byte(s)
regular check-in

image(3)[] img = load("../data/zimg112.nrrd");
//Scalar Fields
field#2(3)[] F = img ⊛ bspln3;
field#1(2)[] F2 = ctmr ⊛ load("../data/ddro-80.nrrd");
field#1(2)[2] V = load("../data/turb2d.nrrd") ⊛ ctmr;
field#4(2)[2] V2 = load("../data/turb2d.nrrd") ⊛ bspln5;
field#4(3)[] F3 = img ⊛ bspln5;

//real laplac =∇•∇F;

//Hessian of scalar field
field#0(3)[3,3] hess1 =∇⊗∇F;
//Hessian of a vector Field 
field#2(2)[2,2,2] hessian=∇⊗∇⊗V2;
field#2(2)[2,2,2] hessian2=hessian;

field#0(2)[2,2] jacobian =∇⊗V;
// curl
//field#0(2)[3] curl1 = ∇×∇F; //3d curl
//field#0(2)[] curl2 =∇×V;  //2d curl

//Trace of Hessian=>Laplacian

tensor[3,3] t=hess1([0,2,4]);
real r=trace(t);
real test=r;

field#1(3)[3] R1=∇(F+F);
field#1(3)[3] testR1=R1;

field#1(3)[3] R2=∇(9*F);
field#1(3)[3] testR2=R2;

real R3=∇•(V+V);
real testR3=R3;

//Curl Grad=>0
field#0(2)[3] R4=∇×∇F3;
field#0(2)[3] testR4=R4;

// Dot Curl=>0
real R5=∇•∇×∇F3;
real testR5=R5;


int m=6;
real f=m*1;
real s=F3([f,f,f]);
real l=s+1;

strand fieldTest (int row)
real row4=row;
real row1=row4*9;
    vec2 pos=[row,row];

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

ViewVC Help
Powered by ViewVC 1.0.0