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

SCM Repository

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

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

Parent Directory Parent Directory | Revision Log Revision Log

Revision 2558 - (download) (annotate)
Mon Mar 3 22:55:30 2014 UTC (5 years, 5 months ago) by cchiw
File size: 2170 byte(s)
Added to Eqfn

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#3(2)[] F7 =bspln5 ⊛ load("../data/ddro-80.nrrd");

field#0(2)[] F4 = tent⊛ load("../data/ddro-80.nrrd");
//Vector Fields

image(2)[2] img2 =load("../data/turb2d.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 ⊛ c4hexic;

field#1(2)[3] RGB = ctmr  ⊛ load("../data/txf/2d-bow.nrrd");


field#1(2)[] field2=F2+9;

field#1(2)[] field1=field2-field2;
real s= field1([6,7]);
vec2 gradA=[2,3];



field#4(3)[] field1=F3*9;
real s= field1([6,7,9]);
vec2 gradA=[2,3];


field#0(2)[] field1=F4*9;
real s= field1([6,7]);
vec2 gradA=[2,3];

field#0(2)[2] grad =∇F2;
vec2 gradA=grad([6,7]);
real s=9;

//Laplacian and Div has split problems
//laplacian 3-d

field#0(3)[] laplac =∇•∇F;
real s=laplac([5,6,7]);
vec2 gradA=s*[8,9];

//laplacian 2-d
field#0(2)[] laplac =∇•∇F2;
real s=laplac([6,7]);
vec2 gradA=s*[8,9];

field#0(2)[] D =∇•V;
field#0(2)[] div =D*9;
vec2 gradA=[6,7];
real m=div(gradA);
real s=m*9;

//Curl 2-d

field#0(2)[] curl2 =∇×V;
real s=curl2([5,6]);
vec2 gradA=[8,8];

//At tthis point need to change statement in strand method
field#0(2)[2,2] jacobian =∇⊗V2;
tensor [2,2] M=jacobian([6,9]);

//Hessian of a Scalar Field 2d

field#0(2)[2,2] hess1 =∇⊗∇F7;
tensor [2,2] M=hess1([6,9]);

//Hessian of a Scalar Field 3-d

field#0(3)[3,3] hess1 =∇⊗∇F;
tensor [3,3] M=hess1([5,6,9]);

//Hessian of a Vector Field

field#2(2)[2,2,2] A=∇⊗∇⊗V2;
tensor[2,2,2] B=A([0,2]);


//vec4 AA=[3,4,5,6];
//vec4 C=(AA+AA)/9;

//Curl 3-d 
//Can't find a field to test curl 3-d on 
//field#1(2)[3] L=∇×RGB;

strand fieldTest (int row)
//vec4 M=C*9;
vec2 gradB=gradA*s;//default
//tensor [2,2] D=M*9;//Jacobian. HessianScalar  2-d
//tensor [3,3] D=M*9;//Hessian Scalar 3-d
//tensor[2,2,2] C=B;//Hessian Vector



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

ViewVC Help
Powered by ViewVC 1.0.0