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 2611 - (download) (annotate)
Mon May 5 21:21:12 2014 UTC (5 years, 2 months ago) by cchiw
File size: 2278 byte(s)
InnerProduct, DoubleDot:

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");

vec2 v2=[4,5];


/*
//Fields

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

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




//Fields

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


//Fields-tent

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


*/
//grad
field#0(2)[2] grad =∇F2;
vec2 gradA=grad([6,7]);
real s=(v2•v2);
/*
//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 gradA2=s*[8,9];
vec2 gradA=gradA2+[9,4];


//Div
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
//Jacobian
/*
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;
vec2 vv2=[8,9];
tensor [2,2] vv3=[vv2,vv2];
strand fieldTest (int row)
{
//vec4 M=C*9;
vec2 gradB=gradA*s;//default
//tensor [2,2] D=M+vv3;//Jacobian. HessianScalar  2-d
//tensor [3,3] D=M+M*9;//Hessian Scalar 3-d
//tensor[2,2,2] C=B;//Hessian Vector

    update
    {
     

 
    }
}

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

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