Алгоритм Форчуна

Псевдокод алгоритма[4].

Пусть ${\displaystyle *(z)}$ будет преобразованием ${\displaystyle *(z)=(z_{x},z_{y}+d(z))}$,
  где ${\displaystyle d(z)}$ — евклидово расстояние между  z и ближайшей точкой
Пусть T будет «береговой линией»
Пусть ${\displaystyle R_{p}}$ будет областью, покрывающей точку p.
Пусть ${\displaystyle C_{pq}}$ будет граничным лучом между точками p и q.
Пусть ${\displaystyle p_{1},p_{2},...,p_{m}}$ будут точками с минимальной y-координатой, упорядоченными по x-координате
${\displaystyle Q\gets S-{p_{1},p_{2},...,p_{m}}}$
создаёт начальные вертикальные граничные лучи ${\displaystyle C_{p_{1},p_{2}}^{0},C_{p_{2},p_{3}}^{0},...C_{p_{m-1},p_{m}}^{0}}$
${\displaystyle T\gets *(R_{p_{1}}),C_{p_{1},p_{2}}^{0},*(R_{p_{2}}),C_{p_{2},p_{3}}^{0},...,*(R_{p_{m-1}}),C_{p_{m-1},p_{m}}^{0},*(R_{p_{m}})}$
пока не IsEmpty(Q) выполнить
    ${\displaystyle p\gets \mathrm {DeleteMin} (Q)}$
    в случае если 
    p является точкой в ${\displaystyle *(V)}$:
        Находим область ${\displaystyle *(R_{q})}$ в T, содержащую p,
          ограниченную кривой ${\displaystyle C_{rq}}$ слева и кривой ${\displaystyle C_{qs}}$ справа
        создаём новые граничные лучи ${\displaystyle C_{pq}^{-}}$ и ${\displaystyle C_{pq}^{+}}$ с основанием p
        заменяем ${\displaystyle *(R_{q})}$ на ${\displaystyle *(R_{q}),C_{pq}^{-},*(R_{p}),C_{pq}^{+},*(R_{q})}$ в T
        удаляем из Q любое пересечение между ${\displaystyle C_{rq}}$ и ${\displaystyle C_{qs}}$
        вставляем в Q любое пересечение ${\displaystyle C_{rq}}$ и ${\displaystyle C_{pq}^{-}}$
        вставляем в Q любое пересечение ${\displaystyle C_{pq}^{+}}$ и ${\displaystyle C_{qs}}$
    p является вершиной Вороного в ${\displaystyle *(V)}$:
        пусть p будет пересечением ${\displaystyle C_{qr}}$ слева и ${\displaystyle C_{rs}}$ справа
        пусть ${\displaystyle C_{uq}}$ будет левым соседом ${\displaystyle C_{qr}}$ и
          пусть ${\displaystyle C_{sv}}$ будет правым соседом ${\displaystyle C_{rs}}$ в T
        создаём новый граничный луч ${\displaystyle C_{qs}^{0}}$, если ${\displaystyle q_{y}=s_{y}}$,
          или создаём ${\displaystyle C_{qs}^{+}}$, если p правая часть более высокой  из q и s,
          в противном случае создаём ${\displaystyle C_{qs}^{-}}$
        заменяем ${\displaystyle C_{qr},*(R_{r}),C_{rs}}$ на вновь созданную ${\displaystyle C_{qs}}$ в T
        удаляем из Q любое пересечение ${\displaystyle C_{uq}}$ и ${\displaystyle C_{qr}}$
        удаляем из Q любое пересечение ${\displaystyle C_{rs}}$ и ${\displaystyle C_{sv}}$
        вставляем в Q любое пересечение ${\displaystyle C_{uq}}$ и ${\displaystyle C_{qs}}$
        вставляем в Q любое пересечение ${\displaystyle C_{qs}}$ и ${\displaystyle C_{sv}}$
        записываем p как вершину ${\displaystyle C_{qr}}$ и ${\displaystyle C_{rs}}$ и основание ${\displaystyle C_{qs}}$
        выводим сегменты границ ${\displaystyle C_{qr}}$ и ${\displaystyle C_{rs}}$
    конец в случае
конец пока
выводим оставшиеся граничные лучи из T