From e15d0190cd453551618a6d3a8756aa398ded63af Mon Sep 17 00:00:00 2001
From: jstark <jstark@mpi-cbg.de>
Date: Thu, 9 Jun 2022 09:07:11 +0200
Subject: [PATCH] Overloading analytical sdf for sphere and disk to enable
center as array.
---
.../tests/analytical_SDF/AnalyticalSDF.hpp | 158 ++++++++++++++++--
1 file changed, 141 insertions(+), 17 deletions(-)
diff --git a/src/level_set/redistancing_Sussman/tests/analytical_SDF/AnalyticalSDF.hpp b/src/level_set/redistancing_Sussman/tests/analytical_SDF/AnalyticalSDF.hpp
index 735e61b5..341b018a 100644
--- a/src/level_set/redistancing_Sussman/tests/analytical_SDF/AnalyticalSDF.hpp
+++ b/src/level_set/redistancing_Sussman/tests/analytical_SDF/AnalyticalSDF.hpp
@@ -43,8 +43,8 @@
* sphere of given radius, where the SDF has positive values inside and negative values outside the sphere.
*/
template <typename point_type, typename space_type>
-space_type get_analytic_sdf_sphere(point_type coords, space_type radius,
- space_type center_x=0, space_type center_y=0, space_type center_z=0)
+space_type get_analytic_sdf_sphere(const point_type coords, const space_type radius,
+ const space_type center_x=0, const space_type center_y=0, const space_type center_z=0)
{
typedef typename std::remove_const_t<std::remove_reference_t<decltype(coords.get(0))>> coord_type;
if(!(std::is_same<space_type, coord_type>::value))
@@ -54,9 +54,35 @@ space_type get_analytic_sdf_sphere(point_type coords, space_type radius,
}
const space_type X = coords.get(0), Y = coords.get(1), Z = coords.get(2);
return (radius -
- sqrt((X - center_x) * (X - center_x)
- + (Y - center_y) * (Y - center_y)
- + (Z - center_z) * (Z - center_z)));
+ sqrt((X - center_x) * (X - center_x)
+ + (Y - center_y) * (Y - center_y)
+ + (Z - center_z) * (Z - center_z)));
+}
+
+/**@brief Based on coordinates, radius and center, returns signed distance function of a sphere.
+ *
+ * @tparam point_type Template type of input coordinates.
+ * @tparam space_type Template type of space.
+ * @param coords Input coordinates of type point_type.
+ * @param radius Radius of type space_type.
+ * @param center Center array of type space_type[3].
+ * @return Value of signed distance function.
+ */
+template <typename point_type, typename space_type>
+space_type get_analytic_sdf_sphere(const point_type coords, const space_type radius, const space_type center[point_type::dims])
+{
+ size_t x = 0, y = 1, z = 2;
+ typedef typename std::remove_const_t<std::remove_reference_t<decltype(coords.get(0))>> coord_type;
+ if(!(std::is_same<space_type, coord_type>::value))
+ {
+ std::cout << "Radius-, Center- and Space-type of grid must be the same! Aborting..." << std::endl;
+ abort();
+ }
+ const space_type X = coords.get(0), Y = coords.get(1), Z = coords.get(2);
+ return (radius -
+ sqrt((X - center[x]) * (X - center[x])
+ + (Y - center[y]) * (Y - center[y])
+ + (Z - center[z]) * (Z - center[z])));
}
/**@brief Initializes the exact solution of the signed distance function of a sphere on an OpenFPM grid.
@@ -71,11 +97,10 @@ space_type get_analytic_sdf_sphere(point_type coords, space_type radius,
* @param center_x Space_type x-coordinate of sphere center.
* @param center_y Space_type y-coordinate of sphere center.
* @param center_z Space_type z-coordinate of sphere center.
-
*/
template <size_t SDF_exact, typename grid_type, typename space_type>
-void init_analytic_sdf_sphere(grid_type & grid, space_type radius, space_type center_x=0, space_type center_y=0,
- space_type center_z=0)
+void init_analytic_sdf_sphere(grid_type & grid, const space_type radius, const space_type center_x=0, const space_type center_y=0,
+ const space_type center_z=0)
{
if(!(std::is_same<typename grid_type::stype, space_type>::value))
{
@@ -94,7 +119,7 @@ void init_analytic_sdf_sphere(grid_type & grid, space_type radius, space_type ce
}
}
-/**@brief Computes the analytical signed distance function of a circle for a given point in space.
+/**@brief Computes the analytical signed distance function of a sphere for a given point in space.
*
* @details At the center of the circle, \a φ_SDF_analytic = \a Radius. Moving outwards from the center on, the
* value for the SDF decreases steadily, eventually becomes 0 at the circle surface and negative beyond the surface. The
@@ -102,19 +127,59 @@ void init_analytic_sdf_sphere(grid_type & grid, space_type radius, space_type ce
*
* @f[ \phi_{SDF}(x, y, z) = R - \sqrt{((x-a)^2 + (y-b)^2 + (z-c)^2)} @f]
*
+ * @tparam SDF_exact Size_t index of property to which analytical SDF of sphere will be saved.
+ * @tparam grid_type Template type of input grid.
+ * @tparam space_type Template type of space.
+ * @param grid Input grid of type grid_type.
+ * @param radius Radius of type space_type.
+ * @param center Array containing 3 elements for center coordinate in x, y and z.
+ */
+template <size_t SDF_exact, typename grid_type, typename space_type>
+void init_analytic_sdf_sphere(grid_type & grid, const space_type radius, const space_type center[grid_type::dims])
+{
+ if(!(std::is_same<typename grid_type::stype, space_type>::value))
+ {
+ std::cout << "Radius-, Center- and Space-type of grid must be the same! Aborting..." << std::endl;
+ abort();
+ }
+
+ auto dom = grid.getDomainIterator();
+ while(dom.isNext())
+ {
+ auto key = dom.get();
+ Point<grid_type::dims, typename grid_type::stype> coords = grid.getPos(key);
+ grid.template getProp<SDF_exact>(key) = get_analytic_sdf_sphere(coords, radius, center);
+ ++dom;
+ }
+}
+
+
+
+
+///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
+
+
+
+/**@brief Computes the analytical signed distance function of a circle for a given point in space.
+ *
+ * @details At the center of the circle, \a φ_SDF_analytic = \a Radius. Moving outwards from the center on, the
+ * value for the SDF decreases steadily, eventually becomes 0 at the circle surface and negative beyond the surface. The
+ * analytic SDF for a circle of radius \a R and center (\a a, \a b) is thus:
+ *
+ * @f[ \phi_{SDF}(x, y, z) = R - \sqrt{((x-a)^2 + (y-b)^2)} @f]
+ *
* @tparam point_type Template type of OpenFPM Point<dimension, type>.
* @tparam space_type Template type of radius.
* @param coords Point_type coordinates of point.
- * @param radius Space_type radius of the sphere.
- * @param center_x Space_type x-coordinate of sphere center.
- * @param center_y Space_type y-coordinate of sphere center.
- * @param center_z Space_type z-coordinate of sphere center.
+ * @param radius Space_type radius of the disk.
+ * @param center_x Space_type x-coordinate of disk center.
+ * @param center_y Space_type y-coordinate of disk center.
* @return Space_type variable that contains the exact solution for the signed distance function of a given point in a
- * sphere of given radius, where the SDF has positive values inside and negative values outside the sphere.
+ * disk of given radius, where the SDF has positive values inside and negative values outside the disk.
*/
template <typename point_type, typename space_type>
-space_type get_analytic_sdf_circle(point_type coords, space_type radius,
- space_type center_x=0, space_type center_y=0)
+space_type get_analytic_sdf_circle(point_type coords, const space_type radius,
+ const space_type center_x=0, const space_type center_y=0)
{
typedef typename std::remove_const_t<std::remove_reference_t<decltype(coords.get(0))>> coord_type;
if(!(std::is_same<space_type, coord_type>::value))
@@ -128,6 +193,31 @@ space_type get_analytic_sdf_circle(point_type coords, space_type radius,
+ (Y - center_y) * (Y - center_y)));
}
+/**@brief Based on coordinates, radius and center, returns signed distance function of a disk.
+ *
+ * @tparam point_type Template type of input coordinates.
+ * @tparam space_type Template type of space.
+ * @param coords Input coordinates of type point_type.
+ * @param radius Radius of type space_type.
+ * @param center Center array of type space_type[2].
+ * @return Value of signed distance function.
+ */
+template <typename point_type, typename space_type>
+space_type get_analytic_sdf_circle(const point_type coords, const space_type radius, const space_type center[point_type::dims])
+{
+ size_t x = 0, y = 1;
+ typedef typename std::remove_const_t<std::remove_reference_t<decltype(coords.get(0))>> coord_type;
+ if(!(std::is_same<space_type, coord_type>::value))
+ {
+ std::cout << "Radius-, Center- and Space-type of grid must be the same! Aborting..." << std::endl;
+ abort();
+ }
+ const space_type X = coords.get(0), Y = coords.get(1);
+ return (radius -
+ sqrt((X - center[x]) * (X - center[x])
+ + (Y - center[y]) * (Y - center[y])));
+}
+
/**@brief Initializes the exact solution of the signed distance function of a circle on an OpenFPM grid.
*
* @details Solves the exact SDF for each grid nodes and writes the solution to a given property.
@@ -141,7 +231,7 @@ space_type get_analytic_sdf_circle(point_type coords, space_type radius,
* @param center_y Y-coordinate of the circle center.
*/
template <size_t SDF_exact, typename grid_type, typename space_type>
-void init_analytic_sdf_circle(grid_type & grid, space_type radius, space_type center_x=0, space_type center_y=0)
+void init_analytic_sdf_circle(grid_type & grid, const space_type radius, const space_type center_x=0, const space_type center_y=0)
{
if(!(std::is_same<typename grid_type::stype, space_type>::value))
{
@@ -159,6 +249,40 @@ void init_analytic_sdf_circle(grid_type & grid, space_type radius, space_type ce
}
}
+/**@brief Computes the analytical signed distance function of a disk for a given point in space.
+ *
+ * @details At the center of the circle, \a φ_SDF_analytic = \a Radius. Moving outwards from the center on, the
+ * value for the SDF decreases steadily, eventually becomes 0 at the circle surface and negative beyond the surface. The
+ * analytic SDF for a circle of radius \a R and center (\a a, \a b) is thus:
+ *
+ * @f[ \phi_{SDF}(x, y, z) = R - \sqrt{((x-a)^2 + (y-b)^2)} @f]
+ *
+ * @tparam SDF_exact Size_t index of property to which analytical SDF of circle will be saved.
+ * @tparam grid_type Template type of input grid.
+ * @tparam space_type Template type of space.
+ * @param grid Input grid of type grid_type.
+ * @param radius Radius of type space_type.
+ * @param center Array containing 2 elements for center coordinate in x and y.
+ */
+template <size_t SDF_exact, typename grid_type, typename space_type>
+void init_analytic_sdf_circle(grid_type & grid, const space_type radius, const space_type center[grid_type::dims])
+{
+ if(!(std::is_same<typename grid_type::stype, space_type>::value))
+ {
+ std::cout << "Radius-, Center- and Space-type of grid must be the same! Aborting..." << std::endl;
+ abort();
+ }
+
+ auto dom = grid.getDomainIterator();
+ while(dom.isNext())
+ {
+ auto key = dom.get();
+ Point<grid_type::dims, typename grid_type::stype> coords = grid.getPos(key);
+ grid.template getProp<SDF_exact>(key) = get_analytic_sdf_circle(coords, radius, center);
+ ++dom;
+ }
+}
+
#endif //ANALYTICAL_SDF_HPP
--
GitLab