From 05792f846a9aa1acb80cdc8b454a2d6f05e3f6de Mon Sep 17 00:00:00 2001
From: Julian Stahl <julianstahl.dev@gmail.com>
Date: Tue, 20 Sep 2022 23:38:42 +0200
Subject: [PATCH] Add operator() to grid

---
 Isosurfacing_3/include/CGAL/Cartesian_grid_3.h | 6 +++++-
 1 file changed, 5 insertions(+), 1 deletion(-)

diff --git a/Isosurfacing_3/include/CGAL/Cartesian_grid_3.h b/Isosurfacing_3/include/CGAL/Cartesian_grid_3.h
index 7a110e5fd63..5e7ffd188c0 100644
--- a/Isosurfacing_3/include/CGAL/Cartesian_grid_3.h
+++ b/Isosurfacing_3/include/CGAL/Cartesian_grid_3.h
@@ -12,9 +12,9 @@
 #ifndef CGAL_CARTESIAN_GRID_3_H
 #define CGAL_CARTESIAN_GRID_3_H
 
-#include <CGAL/license/Isosurfacing_3.h>
 #include <CGAL/Bbox_3.h>
 #include <CGAL/Image_3.h>
+#include <CGAL/license/Isosurfacing_3.h>
 
 #include <array>
 #include <type_traits>
@@ -46,6 +46,10 @@ public:
         from_image(image);
     }
 
+    FT operator()(const std::array<std::size_t, 3> &idx) const {
+        return values[linear_index(idx[0], idx[1], idx[2])];
+    }
+
     FT value(const std::size_t x, const std::size_t y, const std::size_t z) const {
         return values[linear_index(x, y, z)];
     }