wrap up

Author: JoMee

src/FEM/ProjectCurrent.hpp

@@ -18,25 +18,20 @@ inline void assemble_current_whitney1(const Mesh& mesh,
                                       const NedelecSpace& space,
                                       policy_type iteration_policy)
 {
-
   using T = Mesh::value_type;
 
-  const auto org_v = mesh.getOrigin();       // ippl::Vector<T,Dim>
-  const auto h_v   = mesh.getMeshSpacing();  // ippl::Vector<T,Dim>
+  const auto origin = mesh.getOrigin();
+  const auto h   = mesh.getMeshSpacing();
   constexpr unsigned Dim = Mesh::Dimension;
 
-  std::array<T,Dim> origin_a{};
-  std::array<T,Dim> h_a{};
-  for (unsigned d = 0; d < Dim; ++d) {
-    origin_a[d] = org_v[d];
-    h_a[d]      = h_v[d];
-  }
-
   Kokkos::parallel_for("assemble_current_whitney1_make_segments", iteration_policy,
     KOKKOS_LAMBDA(const std::size_t p) {
 
-    std::cout << X0(p) << std::endl;
-    std::cout << X1(p) << std::endl;
+    auto segs = ippl::GridPathSegmenter<Dim, T, ippl::DefaultCellCrossingRule>
+                    ::split(X0(p), X1(p), origin, h);
+    
+     
+
 
   });