heartL <- mesh3dL_to_cgalMeshL(data_heart_obsL)
pairL <- get_mesh_pairs(heartL)
metro <- get_mesh_metro_pair(pairL[[1]], nSamples=500L, silent=TRUE)
# Hausdorff distance
get_HD_max <- function(x) {
max(c(x$ForwardSampling$maxdist,
x$BackwardSampling$maxdist))
}
get_HD_max(metro)
Run the code above in your browser using DataLab