AnalyticalSDF.hpp

Go to the documentation of this file.

//
// Created by jstark on 27.05.21.
//
 
#ifndef OPENFPM_PDATA_ANALYTICALSDF_HPP
#define OPENFPM_PDATA_ANALYTICALSDF_HPP
 
//
// Created by jstark on 2020-10-05.
//
#ifndef ANALYTICAL_SDF_HPP
#define ANALYTICAL_SDF_HPP
 
#include <iostream>
#include "Vector/vector_dist.hpp"
#include "Grid/grid_dist_id.hpp"
 
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_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))
    {
        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)));
}
 
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])));
}
 
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_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))
    {
        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_x,
                                                                        center_y, center_z);
        ++dom;
    }
}
 
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;
    }
}
 
 
 
 
 
 
 
template <typename point_type, typename space_type>
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))
    {
        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)));
}
 
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])));
}
 
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_x=0, const space_type center_y=0)
{
    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_x,
                                                                        center_y);
        ++dom;
    }
}
 
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
 
 
#endif //OPENFPM_PDATA_ANALYTICALSDF_HPP