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

TestTriNormal works

Browse Source
This commit is contained in:
James E McClure 2018-09-13 22:58:43 -04:00
parent afe64e8076
commit b4886cfb92
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
tests/TestTriNormal.cpp
View File

@ -44,7 +44,7 @@ int main(int argc, char **argv)
DoubleArray SDs(Nx,Ny,Nz);
DoubleArray SDs_x(Nx,Ny,Nz);
DoubleArray SDs_y(Nx,Ny,Nz);
DoubleArray SDs_y(Nx,Ny,Nz);
DoubleArray SDs_z(Nx,Ny,Nz);
printf("Set distance map \n");
for (k=0; k<Nz; k++){
@ -63,13 +63,14 @@ int main(int argc, char **argv)
double s,s1,s2,s3;
double a1,a2,a3;
double Vx,Vy,Vz,Wx,Wy,Wz,nx,ny,nz,norm;
double isovalue = 0.f;
for (int k=1; k<Nz-1; k++){
for (int j=1; j<Ny-1; j++){
for (int i=1; i<Nx-1; i++){
object.LocalIsosurface(SDs,isovalue,i,j,k);
for (unsigned long int idx=0; idx<object.TriangleCount; idx++){
e1 = object.Face(idx);
U = TriNormal(e1);
U = object.TriNormal(e1);
// normal from gradient
nx = SDs_x(i,j,k);
ny = SDs_y(i,j,k);