move assertion that SD don't move by more than 1 domain from adjust_x to unpack

Author: pdziekan

src/impl/particles_impl_pack.ipp

@@ -99,16 +99,6 @@ namespace libcloudphxx
         thrust::make_permutation_iterator(x.begin(), lft_id.begin()), // in place
         detail::remote<real_t>(x0, x1)
       );
-#if !defined(NDEBUG)
-      auto min_it = thrust::min_element(
-        thrust::make_permutation_iterator(x.begin(), lft_id.begin()),
-        thrust::make_permutation_iterator(x.begin(), lft_id.begin()) + lft_count);
-      if(*min_it < x0)
-      {
-        std::cerr << "x (" << *min_it << ")  < x0 (" << x0 << ") after adjustment for distmem copy, potentially SD moved by more than one process/GPU domain size" << std::endl;
-        assert(0);
-      }
-#endif
     }
 
     template <typename real_t, backend_t device>
@@ -120,16 +110,6 @@ namespace libcloudphxx
         thrust::make_permutation_iterator(x.begin(), rgt_id.begin()), // in place
         detail::remote<real_t>(x1, x0)
       );
-#if !defined(NDEBUG)
-      auto max_it = thrust::max_element(
-        thrust::make_permutation_iterator(x.begin(), lft_id.begin()),
-        thrust::make_permutation_iterator(x.begin(), lft_id.begin()) + lft_count);
-      if(*max_it >= x1)
-      {
-        std::cerr << "x (" << *max_it << ")  >= x1 (" << x1 << ") after adjustment for distmem copy, potentially SD moved by more than one process/GPU domain size" << std::endl;
-        assert(0);
-      }
-#endif
     }
   };
 };

src/impl/particles_impl_unpack.ipp

@@ -31,6 +31,20 @@ namespace libcloudphxx
         thrust::copy( in_real_bfr.begin() + distance * n_copied, in_real_bfr.begin() + (distance+1) * n_copied, (*it)->begin() + n_part_old);
         it++;
       }
+#if !defined(NDEBUG)
+      auto min_it = thrust::min_element(x.begin() + n_part_old, x.end());
+      if(*min_it < opts_init.x0)
+      {
+        std::cerr << "x (" << *min_it << ")  < opts_init.x0 (" << opts_init.x0 << ") after unpacking, potentially SD moved by more than one process/GPU domain size" << std::endl;
+        assert(0);
+      }
+      auto max_it = thrust::max_element(x.begin() + n_part_old, x.end());
+      if(*max_it >= opts_init.x1)
+      {
+        std::cerr << "x (" << *max_it << ")  >= opts_init.x1 (" << opts_init.x1 << ") after unpacking, potentially SD moved by more than one process/GPU domain size" << std::endl;
+        assert(0);
+      }
+#endif
     }
 
     template <typename real_t, backend_t device>