Changeset f8dae2

Timestamp:

Oct 20, 2016, 8:52:32 PM (8 years ago)

Author:

Frederik Heber <heber@…>

Branches:

Fix_IndependentFragmentGrids

Parents:

8d6a44

git-author:

Frederik Heber <heber@…> (09/27/16 08:21:28)

git-committer:

Frederik Heber <heber@…> (10/20/16 20:52:32)

Message:

DROPME: Added lots of debugging to ExportGraph_ToJobs.

File:

• src/Fragmentation/Exporters/ExportGraph_ToJobs.cpp (modified) (9 diffs)

src/Fragmentation/Exporters/ExportGraph_ToJobs.cpp

@@ -99 +99 @@
   /// determine center of fragment
   const Vector center = .5*(_minmax.first + _minmax.second);
-  LOG(4, "DEBUG: center of fragment is at " << center);
+  LOG(2, "DEBUG: center of fragment is at " << center);
 
   /// associate center to its containing grid cell (defined by boundary points)
@@ -111 +111 @@
         fabs(lower_center[i] - higher_center[i]) < std::numeric_limits<double>::epsilon()*1e4;
   }
-  LOG(5, "DEBUG: lower_center is " << lower_center << ", higher_center is " << higher_center);
+  LOG(2, "DEBUG: lower_center is " << lower_center << ", higher_center is " << higher_center);
 
   // fashion min max into cubic box extents of at least extent plus empty
@@ -117 +117 @@
   const Vector extent = _minmax.second - _minmax.first;
   double greatest_extent = extent[extent.GreatestComponent()];
+  LOG(2, "DEBUG: extent of fragment is " << extent << " with greatest extent " << greatest_extent);
   if (greatest_extent > _minimum_empty_boundary)
     greatest_extent *= 3.;
@@ -128 +129 @@
   const double greatest_total = total[total.GreatestComponent()];
   greatest_extent = std::min(greatest_extent, greatest_total);
-  LOG(5, "DEBUG: extent of fragment is " << extent << ", greatest_extent is " << greatest_extent);
+  LOG(2, "DEBUG: Using boundary of " << greatest_extent);
 
   /// increase levels around this center to find the matching window
@@ -136 +137 @@
       _domain.getDeltaPerAxis(2)
   };
-  LOG(6, "DEBUG: delta is " << Vector(delta));
+  LOG(2, "DEBUG: delta is " << Vector(delta));
   const double greatest_delta = std::max(delta[0], std::max(delta[1], delta[2]));
 
@@ -147 +148 @@
   int end_steps[NDIM] = { 0, 0, 0 };
   for (size_t i=0;i<NDIM;++i) {
-    begin_index[i] = round(lower_center[i]/delta[i]) - (half_fragment_gridpoints-1);
-    end_index[i] = round(higher_center[i]/delta[i]) + (half_fragment_gridpoints+(unsigned int)(equal_components[i]));
+    const unsigned int lower_index = round(lower_center[i]/delta[i]);
+    const unsigned int higher_index = round(higher_center[i]/delta[i]);
+    LOG(3, "DEBUG: center is within indices [" << lower_index << ", " << higher_index << "]");
+    begin_index[i] = lower_index - (half_fragment_gridpoints-1);
+    end_index[i] = higher_index + (half_fragment_gridpoints+(unsigned int)(equal_components[i]));
+    LOG(3, "DEBUG: Intermediate begin end indices for axis #" << i << " are ["
+        << begin_index[i] << ", " << end_index[i] << "]");
     if (begin_index[i] < 0)
       begin_steps[i] = -begin_index[i];
@@ -160 +166 @@
       end_steps[i] = 0;
     }
+    LOG(3, "DEBUG: Resulting final indices [" << begin_index[i] << ", "
+        << end_index[i] << "] and step indices for axis #" << i << " are ["
+        << begin_steps[i] << ", " << end_steps[i] << "]");
   }
   for (size_t i=0;i<NDIM;++i) {
@@ -165 +174 @@
     fragment_window.end[i] = (end_index[i]-end_steps[i]+begin_steps[i]) * delta[i];
   }
-    LOG(6, "DEBUG: fragment begin is " << Vector(fragment_window.begin)
+    LOG(2, "DEBUG: fragment begin is " << Vector(fragment_window.begin)
         << ", fragment end is " << Vector(fragment_window.end));
 
@@ -191 +200 @@
   }
 #endif
-  if (DoLog(6))
+  if (DoLog(3))
     for (size_t i=0;i<NDIM;++i)
-      LOG(6, "DEBUG: We have " << (fragment_window.end[i] - fragment_window.begin[i])/delta[i]
+      LOG(3, "DEBUG: We have " << (fragment_window.end[i] - fragment_window.begin[i])/delta[i]
           << " gridpoints on axis " << i);