Moving least squares interpolation

CMakeLists.txt

@@ -69,11 +69,15 @@ target_include_directories(ArborX INTERFACE
   $<BUILD_INTERFACE:${CMAKE_CURRENT_SOURCE_DIR}/src/details>
   $<BUILD_INTERFACE:${CMAKE_CURRENT_SOURCE_DIR}/src/geometry>
   $<BUILD_INTERFACE:${CMAKE_CURRENT_SOURCE_DIR}/src/kokkos_ext>
+  $<BUILD_INTERFACE:${CMAKE_CURRENT_SOURCE_DIR}/src/interpolation>
+  $<BUILD_INTERFACE:${CMAKE_CURRENT_SOURCE_DIR}/src/interpolation/details>
   $<BUILD_INTERFACE:${CMAKE_CURRENT_BINARY_DIR}/include>
   $<INSTALL_INTERFACE:include/ArborX>
   $<INSTALL_INTERFACE:include/ArborX/details>
   $<INSTALL_INTERFACE:include/ArborX/geometry>
   $<INSTALL_INTERFACE:include/ArborX/kokkos_ext>
+  $<INSTALL_INTERFACE:include/ArborX/interpolation>
+  $<INSTALL_INTERFACE:include/ArborX/interpolation/details>
 )
 
 install(TARGETS ArborX

examples/CMakeLists.txt

@@ -18,4 +18,5 @@ if(Boost_FOUND)
   add_subdirectory(viz)
   add_subdirectory(raytracing)
   add_subdirectory(brute_force)
+  add_subdirectory(moving_least_squares)
 endif()

examples/moving_least_squares/CMakeLists.txt

@@ -0,0 +1,3 @@
+add_executable(ArborX_Example_MovingLeastSquares.exe moving_least_squares.cpp)
+target_link_libraries(ArborX_Example_MovingLeastSquares.exe ArborX::ArborX)
+add_test(NAME ArborX_Example_MovingLeastSquares COMMAND ArborX_Example_MovingLeastSquares.exe)

examples/moving_least_squares/moving_least_squares.cpp

@@ -0,0 +1,104 @@
+/****************************************************************************
+ * Copyright (c) 2023 by the ArborX authors                                 *
+ * All rights reserved.                                                     *
+ *                                                                          *
+ * This file is part of the ArborX library. ArborX is                       *
+ * distributed under a BSD 3-clause license. For the licensing terms see    *
+ * the LICENSE file in the top-level directory.                             *
+ *                                                                          *
+ * SPDX-License-Identifier: BSD-3-Clause                                    *
+ ****************************************************************************/
+
+#include <ArborX.hpp>
+#include <ArborX_InterpMovingLeastSquares.hpp>
+
+#include <Kokkos_Core.hpp>
+
+#include <iostream>
+
+using Point = ArborX::ExperimentalHyperGeometry::Point<2, double>;
+
+KOKKOS_FUNCTION double functionToApproximate(Point const &p)
+{
+  return Kokkos::cos(p[0] + p[1] / 4);
+}
+
+int main(int argc, char *argv[])
+{
+  Kokkos::ScopeGuard guard(argc, argv);
+
+  using ExecutionSpace = Kokkos::DefaultExecutionSpace;
+  using MemorySpace = typename ExecutionSpace::memory_space;
+  ExecutionSpace space{};
+
+  // Source space is a 3x3 grid
+  // Target space is a single point
+  //
+  //  ^
+  //  |
+  //  S     S     S
+  //  |
+  //  |   T
+  //  S     S  T  S
+  //  |
+  //  | T
+  // -S-----S-----S->
+  //  |
+
+  // Set up points
+  Kokkos::View<Point *, MemorySpace> src_points("Example::src_points", 9);
+  Kokkos::View<Point *, MemorySpace> tgt_points("Example::tgt_points", 3);
+  Kokkos::parallel_for(
+      "Example::make_points", Kokkos::RangePolicy<ExecutionSpace>(space, 0, 1),
+      KOKKOS_LAMBDA(int const) {
+        src_points(0) = {0., 0.};
+        src_points(1) = {1., 0.};
+        src_points(2) = {2., 0.};
+        src_points(3) = {0., 1.};
+        src_points(4) = {1., 1.};
+        src_points(5) = {2., 1.};
+        src_points(6) = {0., 2.};
+        src_points(7) = {1., 2.};
+        src_points(8) = {2., 2.};
+        tgt_points(0) = {4. / 6., 4. / 3.};
+        tgt_points(1) = {9. / 6., 3. / 3.};
+        tgt_points(2) = {2. / 6., 1. / 3.};
+      });
+
+  // Set up values
+  Kokkos::View<double *, MemorySpace> src_values("Example::src_values", 9);
+  Kokkos::View<double *, MemorySpace> app_values("Example::app_values", 3);
+  Kokkos::View<double *, MemorySpace> tgt_values("Example::tgt_values", 3);
+  Kokkos::parallel_for(
+      "Example::make_values", Kokkos::RangePolicy<ExecutionSpace>(space, 0, 9),
+      KOKKOS_LAMBDA(int const i) {
+        src_values(i) = functionToApproximate(src_points(i));
+        if (i < 3)
+          tgt_values(i) = functionToApproximate(tgt_points(i));
+      });
+
+  // Build the moving least squares coefficients
+  ArborX::Interpolation::MovingLeastSquares<MemorySpace> mls(space, src_points,
+                                                             tgt_points);
+
+  // Interpolate
+  mls.interpolate(space, src_values, app_values);
+
+  // Show results
+  auto app_values_host =
+      Kokkos::create_mirror_view_and_copy(Kokkos::HostSpace{}, app_values);
+  auto tgt_values_host =
+      Kokkos::create_mirror_view_and_copy(Kokkos::HostSpace{}, tgt_values);
+  auto diff = [=](int const i) {
+    return Kokkos::abs(app_values_host(i) - tgt_values_host(i));
+  };
+
+  std::cout << "Approximated values: " << app_values_host(0) << ' '
+            << app_values_host(1) << ' ' << app_values_host(2) << '\n';
+  std::cout << "Real values        : " << tgt_values_host(0) << ' '
+            << tgt_values_host(1) << ' ' << tgt_values_host(2) << '\n';
+  std::cout << "Differences        : " << diff(0) << ' ' << diff(1) << ' '
+            << diff(2) << '\n';
+
+  return 0;
+}

src/geometry/ArborX_GeometryTraits.hpp

@@ -90,7 +90,7 @@ struct is_triangle : std::is_same<typename tag<Geometry>::type, TriangleTag>
 {};
 
 template <typename Geometry>
-void check_valid_geometry_traits(Geometry const &)
+KOKKOS_FUNCTION constexpr void check_valid_geometry_traits(Geometry const &)
 {
   static_assert(
       !Kokkos::is_detected<DimensionNotSpecializedArchetypeAlias, Geometry>{},