OilfieldToolsNet/LBPM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Added TriLinearPoly to common/pmmc.h

Browse Source
This commit is contained in:
James E McClure 2014-10-16 10:24:39 -04:00
parent a45e23d9d5
commit 39f0ac6394
1 changed files with 14 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
common/pmmc.h
View File

@ -306,12 +306,19 @@ class TriLinearPoly{
* Values are provided at the corners in CubeValues * Values are provided at the corners in CubeValues
* x,y,z must be defined on [0,1] where the length of the cube edge is one * x,y,z must be defined on [0,1] where the length of the cube edge is one
*/ */
int ic,jc,kc,
double a,b,c,d,e,f,g,h; double a,b,c,d,e,f,g,h;
double x,y,z; double x,y,z;
DoubleArray CubeValues; DoubleArray CubeValues;
public: public:
TriLinearPoly(DoubleArray &A, int i, int j, int k){ TriLinearPoly(){
CubeValues.New(2,2,2); CubeValues.New(2,2,2);
}
void assign(DoubleArray &A, int i, int j, int k){
ic=i; jc=j; kc=k
CubeValues(0,0,0) = A(i,j,k); CubeValues(0,0,0) = A(i,j,k);
CubeValues(1,0,0) = A(i+1,j,k); CubeValues(1,0,0) = A(i+1,j,k);
CubeValues(0,1,0) = A(i,j+1,k); CubeValues(0,1,0) = A(i,j+1,k);
@ -320,7 +327,7 @@ public:
CubeValues(1,0,1) = A(i+1,j,k+1); CubeValues(1,0,1) = A(i+1,j,k+1);
CubeValues(0,1,1) = A(i,j+1,k+1); CubeValues(0,1,1) = A(i,j+1,k+1);
CubeValues(1,1,1) = A(i+1,j+1,k+1); CubeValues(1,1,1) = A(i+1,j+1,k+1);
a = CubeValues(0,0,0); a = CubeValues(0,0,0);
b = CubeValues(1,0,0)-a; b = CubeValues(1,0,0)-a;
c = CubeValues(0,1,0)-a; c = CubeValues(0,1,0)-a;
@ -330,9 +337,12 @@ public:
g = CubeValues(0,1,1)-a-c-d; g = CubeValues(0,1,1)-a-c-d;
h = CubeValues(1,1,1)-a-b-c-d-e-f-g; h = CubeValues(1,1,1)-a-b-c-d-e-f-g;
} }
double eval(double x, double y, double z){ double eval(Point P){
double returnValue; double returnValue;
x = P.x - 1.0*ic;
y = P.y - 1.0*jc;
z = P.z - 1.0*kc;
returnValue = a + b*x + c*y+d*z + e*x*y + f*x*z + g*y*z + h*x*y*z; returnValue = a + b*x + c*y+d*z + e*x*y + f*x*z + g*y*z + h*x*y*z;
return returnValue; return returnValue;
} }