Discrete Hausdorff Distance

Input: data set opts;
Output: Hausdorff distance;
Submodules: -

ClearAll[DiscreteHausdorffDistance];
Options[DiscreteHausdorffDistance] = Options[Nearest];
DiscreteHausdorffDistance[s1_, s2_, opts : OptionsPattern[]] := Max[
	Nearest[s1 -> "Distance", s2, 1, opts],
	Nearest[s2 -> "Distance", s1, 1, opts]
];