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

debugging...

Browse Source
This commit is contained in:
James McClure
2014-04-21 19:51:18 -04:00
parent f60243ea93
commit 07fe4f65e9
1 changed files with 10 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
common/pmmc.h
View File

@@ -3498,6 +3498,7 @@ inline void pmmc_CubeTrimSurfaceInterpValues(DoubleArray &CubeValues, DoubleArra
Point A,B,C;
int p;
double vA,vB,vC;
double dA,dB,dC;
double x,y,z;
double s,s1,s2,s3,temp;
double a,b,c,d,e,f,g,h;
@@ -3564,26 +3565,32 @@ inline void pmmc_CubeTrimSurfaceInterpValues(DoubleArray &CubeValues, DoubleArra
A = Points(Triangles(0,r));
B = Points(Triangles(1,r));
C = Points(Triangles(2,r));
vA = SurfaceValues(Triangles(0,r));
vB = SurfaceValues(Triangles(1,r));
vC = SurfaceValues(Triangles(2,r));
dA = DistanceValues(Triangles(0,r));
dB = DistanceValues(Triangles(1,r));
dC = DistanceValues(Triangles(2,r));
// Compute length of sides (assume dx=dy=dz)
s1 = sqrt((A.x-B.x)*(A.x-B.x)+(A.y-B.y)*(A.y-B.y)+(A.z-B.z)*(A.z-B.z));
s2 = sqrt((A.x-C.x)*(A.x-C.x)+(A.y-C.y)*(A.y-C.y)+(A.z-C.z)*(A.z-C.z));
s3 = sqrt((B.x-C.x)*(B.x-C.x)+(B.y-C.y)*(B.y-C.y)+(B.z-C.z)*(B.z-C.z));
s = 0.5*(s1+s2+s3);
temp = s*(s-s1)*(s-s2)*(s-s3);
if (temp > 0.0){
if (DistanceValues(Triangles(0,r)) > mindist){
if (dA > mindist){
integral += sqrt(temp)*0.33333333333333333*(vA);
area += sqrt(temp)*0.33333333333333333;
}
if (DistanceValues(Triangles(1,r)) > mindist){
if (dB > mindist){
integral += sqrt(temp)*0.33333333333333333*(vB);
area += sqrt(temp)*0.33333333333333333;
}
if (DistanceValues(Triangles(2,r)) > mindist){
if (dC > mindist){
integral += sqrt(temp)*0.33333333333333333*(vC);
area += sqrt(temp)*0.33333333333333333;
}