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

make sure halo is filled when measuring object

Browse Source
This commit is contained in:
James McClure
2022-03-04 18:59:11 -05:00
parent eccebcd95a
commit b39dd1a474
1 changed files with 7 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
analysis/Minkowski.cpp
View File

@@ -67,6 +67,7 @@ Minkowski::~Minkowski() {
void Minkowski::ComputeScalar(const DoubleArray &Field, const double isovalue) { void Minkowski::ComputeScalar(const DoubleArray &Field, const double isovalue) {
PROFILE_START("ComputeScalar"); PROFILE_START("ComputeScalar");
Xi = Ji = Ai = 0.0; Xi = Ji = Ai = 0.0;
DCEL object; DCEL object;
int e1, e2, e3; int e1, e2, e3;
@@ -160,6 +161,10 @@ void Minkowski::MeasureObject() {
* 1 - labels the rest of the * 1 - labels the rest of the
*/ */
//DoubleArray smooth_distance(Nx,Ny,Nz); //DoubleArray smooth_distance(Nx,Ny,Nz);
fillHalo<double> fillData(Dm->Comm, Dm->rank_info, {Nx-2,Ny-2,Nz-2}, {1, 1, 1}, 0, 1);
for (int k = 0; k < Nz; k++) { for (int k = 0; k < Nz; k++) {
for (int j = 0; j < Ny; j++) { for (int j = 0; j < Ny; j++) {
for (int i = 0; i < Nx; i++) { for (int i = 0; i < Nx; i++) {
@@ -168,6 +173,8 @@ void Minkowski::MeasureObject() {
} }
} }
CalcDist(distance, id, *Dm); CalcDist(distance, id, *Dm);
fillData.fill(distance);
//Mean3D(distance,smooth_distance); //Mean3D(distance,smooth_distance);
//Eikonal(distance, id, *Dm, 20, {true, true, true}); //Eikonal(distance, id, *Dm, 20, {true, true, true});
ComputeScalar(distance, 0.0); ComputeScalar(distance, 0.0);