Skip to content
Snippets Groups Projects
Commit e15d0190 authored by jstark's avatar jstark
Browse files

Overloading analytical sdf for sphere and disk to enable center as array.

parent e3ebd6a3
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
src/level_set/redistancing_Sussman/tests/analytical_SDF/AnalyticalSDF.hpp
+ 141
− 17
View file @ e15d0190
......@@ -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 &phi;_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 &phi;_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 &phi;_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
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment