Diener [60
] describes how the algorithm can be enhanced to also determine (integrate) individual null generators of the
event horizon. This requires interpolating the 4-metric to the generator positions but (still) not taking any derivatives of the
4-metric.