From 430ad68dfb183bbadc9610cca6d3d4f50379df07 Mon Sep 17 00:00:00 2001
From: abhiyanpaudel <abhiyanp93@gmail.com>
Date: Tue, 7 Apr 2026 14:00:13 -0400
Subject: [PATCH] fix conversion from flat array to rank 2 array

---
 src/pcms/transfer/mesh_intersection.cpp | 9 +++++++--
 1 file changed, 7 insertions(+), 2 deletions(-)

diff --git a/src/pcms/transfer/mesh_intersection.cpp b/src/pcms/transfer/mesh_intersection.cpp
index 85e41967..6d638814 100644
--- a/src/pcms/transfer/mesh_intersection.cpp
+++ b/src/pcms/transfer/mesh_intersection.cpp
@@ -24,8 +24,13 @@ void FindIntersections::adjBasedIntersectSearch(
 
   const auto flat_centroids =
     pcms::get_entity_centroids(target_mesh_, Omega_h::FACE);
-  auto centroids = Kokkos::View<const Omega_h::Real* [2]>(
-    flat_centroids.data(), target_mesh_.nfaces());
+
+  Kokkos::View<Omega_h::Real* [2]> centroids("centroids", target_mesh_.nfaces());
+
+  Omega_h::parallel_for(target_mesh_.nfaces(), OMEGA_H_LAMBDA(const Omega_h::LO i) {
+		  centroids(i,0) = flat_centroids[2*i+0];
+		  centroids(i,1) = flat_centroids[2*i+1];
+  });
 
   pcms::GridPointSearch2D search_cell(source_mesh_, 20, 20);
   auto results = search_cell(centroids);