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

pn

Browse Source
This commit is contained in:
James E McClure
2018-03-23 09:25:52 -04:00
parent 3c809bfb29
commit 4a60e4d560
1 changed files with 3 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
tests/lbpm_porenetwork_pp.cpp
View File

@@ -209,19 +209,19 @@ int main(int argc, char **argv)
// value of s along center line {x=alpha*s, y = beta*s, z=gamma*s} that is closest to i,j,k
double s = (alpha*xi + beta*yj + gamma*zk)/(alpha*alpha + beta*beta + gamma*gamma);
double distance=0.f;
/*if (s > length){
if (s > length){
distance = radius - sqrt((xi-X)*(xi-X) + (yj-Y)*(yj-Y) + (zk-Z)*(zk-Z));
}
else if (s<0){
distance = radius - sqrt((xi-x)*(xi-x) + (yj-y)*(yj-y) + (zk-z)*(zk-z));
}
else{*/
else{
double xs = alpha*s;
double ys = beta*s;
double zs = gamma*s;
distance = radius - sqrt((xi-xs)*(xi-xs) + (yj-ys)*(yj-ys) + (zk-zs)*(zk-zs));
//if (distance>0)printf("s=%f,alpha=%f,beta=%f,gamma=%f,distance=%f\n",s,alpha,beta,gamma,distance);
//}
}
if (distance > Averages.SDs(i,j,k)) Averages.SDs(i,j,k) = distance;
}