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

added return value to TestTriNormal

Browse Source
This commit is contained in:
James E McClure 2018-09-13 23:00:51 -04:00
parent b4886cfb92
commit ad259be66e
1 changed files with 4 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
tests/TestTriNormal.cpp
View File

@ -64,6 +64,7 @@ int main(int argc, char **argv)
double a1,a2,a3; double a1,a2,a3;
double Vx,Vy,Vz,Wx,Wy,Wz,nx,ny,nz,norm; double Vx,Vy,Vz,Wx,Wy,Wz,nx,ny,nz,norm;
double isovalue = 0.f; double isovalue = 0.f;
int count_plus=0; int count_minus=0;
for (int k=1; k<Nz-1; k++){ for (int k=1; k<Nz-1; k++){
for (int j=1; j<Ny-1; j++){ for (int j=1; j<Ny-1; j++){
for (int i=1; i<Nx-1; i++){ for (int i=1; i<Nx-1; i++){
@ -78,14 +79,17 @@ int main(int argc, char **argv)
double dotprod=U.x*nx + U.y*ny + U.z*nz; double dotprod=U.x*nx + U.y*ny + U.z*nz;
if (dotprod < 0){ if (dotprod < 0){
printf("negative %f \n",dotprod); printf("negative %f \n",dotprod);
count_minus++;
} }
else{ else{
printf("positive %f \n",dotprod); printf("positive %f \n",dotprod);
count_plus++;
} }
} }
} }
} }
} }
if (count_minus > 0 && count_plus>0) toReturn=1;
} }
MPI_Barrier(comm); MPI_Barrier(comm);
MPI_Finalize(); MPI_Finalize();