PISM, A Parallel Ice Sheet Model
stable v1.2 committed by Constantine Khrulev on 20200211 20:24:05 0900

◆ eikonal_equation()
Find an approximate solution of the Eikonal equation on a given domain. To specify the problem, the input field (mask) should be filled with
For example, to compute distances from the grounding line within ice shelves, fill generic ice shelf locations with zeros, set neighbors of the grounding line to 1, and the rest of the grid with 1 or some other negative number. Note that this implementation updates ghosts every iteration. We could speed this up by checking if a point at a boundary of the processor subdomain was updated and update ghosts in those cases only. Definition at line 579 of file PicoGeometry.cc. References pism::IceModelVec2Int::as_int(), pism::GlobalMax(), pism::IceModelVec::grid(), pism::IceModelVec2Int::int_star(), pism::PointsWithGhosts::next(), pism::IceModelVec::stencil_width(), and pism::IceModelVec::update_ghosts(). Referenced by pism::ocean::PicoGeometry::compute_distances_cf(), and pism::ocean::PicoGeometry::compute_distances_gl(). 