DrawDisk.hpp

//
// Created by jstark on 2020-08-18.
//
#ifndef ACCURACY_TESTS_CIRCLE_HPP
#define ACCURACY_TESTS_CIRCLE_HPP

#include <iostream>

template <typename point_type, typename radius_type>
bool inside_disk(point_type coords, radius_type radius, double center_x=0, double center_y=0)
{
        const double EPSILON = std::numeric_limits<double>::epsilon();
        const double X = coords.get(0), Y = coords.get(1);
        return (X - center_x) * (X - center_x)
        + (Y - center_y) * (Y - center_y)
        <= radius * radius + EPSILON;
}

template <size_t Phi_0, typename grid_type, typename radius_type>
void init_grid_with_disk(grid_type & grid, radius_type radius, double center_x=0, double center_y=0)
{
        // assign pixel values to domain. For each pixel get factor_refinement number of grid points with corresponding value
        auto dom = grid.getDomainIterator();
        while(dom.isNext())
        {
                auto key = dom.get();

                Point<grid_type::dims, double> coords = grid.getPos(key);

                if (inside_disk(coords, radius, center_x, center_y))
                {
                        grid.template get<Phi_0> (key) = 1;
                }
                else
                {
                        grid.template get<Phi_0> (key) = -1;
                }
                ++dom;
        }
}


#endif //ACCURACY_TESTS_CIRCLE_HPP