functions.h

/*****************************************************************************
 * VCLib                                                                     *
 * Visual Computing Library                                                  *
 *                                                                           *
 * Copyright(C) 2021-2025                                                    *
 * Visual Computing Lab                                                      *
 * ISTI - Italian National Research Council                                  *
 *                                                                           *
 * All rights reserved.                                                      *
 *                                                                           *
 * This program is free software; you can redistribute it and/or modify      *
 * it under the terms of the Mozilla Public License Version 2.0 as published *
 * by the Mozilla Foundation; either version 2 of the License, or            *
 * (at your option) any later version.                                       *
 *                                                                           *
 * This program is distributed in the hope that it will be useful,           *
 * but WITHOUT ANY WARRANTY; without even the implied warranty of            *
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the              *
 * Mozilla Public License Version 2.0                                        *
 * (https://www.mozilla.org/en-US/MPL/2.0/) for more details.                *
 ****************************************************************************/
 
#ifndef VCL_ALGORITHMS_CORE_DISTANCE_FUNCTIONS_H
#define VCL_ALGORITHMS_CORE_DISTANCE_FUNCTIONS_H
 
#include "element.h"
#include "misc.h"
 
#include <vclib/concepts/mesh.h>
 
namespace vcl {
 
namespace detail {
 
template<typename Obj1, typename Obj2, typename ST>
concept BoundedDistFunctionExists =
    requires (Obj1 o1, Obj2 o2, ST v) { boundedDistance(o1, o2, v); };
 
} // namespace detail
 
template<typename Obj1, typename Obj2>
auto distFunction()
{
    auto f = [](const Obj1& o1, const Obj2& o2) {
        return distance(o1, o2);
    };
 
    return f;
}
 
template<typename Obj1, typename Obj2, typename ScalarType = double>
auto boundedDistFunction()
{
    if constexpr (detail::BoundedDistFunctionExists<Obj1, Obj2, ScalarType>) {
        auto f = [](const Obj1& o1, const Obj2& o2, ScalarType s) {
            return boundedDistance(o1, o2, s);
        };
 
        return f;
    }
    else {
        auto f = [](const Obj1& o1, const Obj2& o2, ScalarType) {
            return distance(o1, o2);
        };
 
        return f;
    }
}
 
} // namespace vcl
 
#endif // VCL_ALGORITHMS_CORE_DISTANCE_FUNCTIONS_H