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 2525 - (download) (annotate)
Tue Jan 21 19:14:22 2014 UTC (5 years, 6 months ago) by cchiw
File size: 1536 byte(s)
eintypes->mid-iltypes

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#0(2)[] F4 = tent⊛ load("../data/ddro-80.nrrd");
//Vector Fields
field#1(2)[2] V = load("../data/turb2d.nrrd") ⊛ ctmr;
field#4(3)[] F3 = img ⊛ c4hexic;

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





//Fields
/*
field#1(2)[] field1=F2*9;
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=9;
*/

//Laplacian and Div has split problems

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

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



//Hessian of a Scalar Field
/*
field#0(3)[3,3] hess1 =∇⊗∇F;
tensor [3,3] M=hess1([5,6,9]);
*/

//Hessian of a Vector Field-ERROR

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

//Curl 2-d
/*field#0(2)[] curl2 =∇×V;
real s=curl2([5,6]);
vec2 gradA=[s,s];
*/

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

strand fieldTest (int row)
{
vec2 gradB=gradA*s;
//tensor [3,3] D=M*9;//Hessian Scalar
//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