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

add geometric test for 3D topo

Browse Source
This commit is contained in:
James McClure 2019-05-22 17:33:10 -06:00
parent 06ecb0db4e
commit f360eb1ecb
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
tests/TestTopo3D.cpp
View File

@ -127,7 +127,7 @@ int main(int argc, char **argv)
//.............................................................................. //..............................................................................
if (x <= 0 || y<=0) { if (x <= 0 || y<=0) {
// Single torus // Single torus
Object.distance(i,j,k) = sqrt((sqrt(x*x+y*y) - R1)*(sqrt(x*x+y*y) - R1) + z*z) - R2; Object.distance(i,j,k) = R2 - sqrt((sqrt(x*x+y*y) - R1)*(sqrt(x*x+y*y) - R1) + z*z);
} }
else { else {
double d1 = R2-sqrt(x*x +(y-R1)*(y-R1) + z*z); double d1 = R2-sqrt(x*x +(y-R1)*(y-R1) + z*z);
@ -201,7 +201,7 @@ int main(int argc, char **argv)
//.............................................................................. //..............................................................................
// Bowl // Bowl
if (z > 0 ){ if (z > 0 ){
Object.distance(i,j,k) = sqrt((sqrt(x*x+y*y) - R1)*(sqrt(x*x+y*y) - R1) + z*z) - R2; Object.distance(i,j,k) = R2-sqrt((sqrt(x*x+y*y) - R1)*(sqrt(x*x+y*y) - R1) + z*z);
} }
else else
{ {