PISM, A Parallel Ice Sheet Model
stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800
|
◆ 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 sub-domain was updated and update ghosts in those cases only. FIXME: replace this with a better algorithm. Definition at line 828 of file PicoGeometry.cc. References pism::array::Scalar::as_int(), pism::GlobalMax(), pism::array::Array::grid(), pism::array::Array2D< T >::star(), pism::array::Array::stencil_width(), and pism::array::Array::update_ghosts(). Referenced by pism::ocean::PicoGeometry::compute_distances_cf(), and pism::ocean::PicoGeometry::compute_distances_gl(). |