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(point_type coords, space_type radius,
                                   space_type center_x=0, space_type center_y=0, 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 <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)
{
    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 <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)
{
    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, space_type radius, space_type center_x=0, 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;
    }
}

#endif //ANALYTICAL_SDF_HPP


#endif //OPENFPM_PDATA_ANALYTICALSDF_HPP