MinimalBoundingSphere.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include <vector>
 
#include "MathLib/Point3d.h"
 
namespace GeoLib
{

class MinimalBoundingSphere final
{
public:
    explicit MinimalBoundingSphere(
        MathLib::Point3d const& p,
        double radius = std::numeric_limits<double>::epsilon());
    MinimalBoundingSphere(MathLib::Point3d const& p, MathLib::Point3d const& q);
    MinimalBoundingSphere(MathLib::Point3d const& p, MathLib::Point3d const& q,
                          MathLib::Point3d const& r);
    MinimalBoundingSphere(MathLib::Point3d const& p,
                          MathLib::Point3d const& q,
                          MathLib::Point3d const& r,
                          MathLib::Point3d const& s);
    explicit MinimalBoundingSphere(
        std::vector<MathLib::Point3d*> const& points);

    MathLib::Point3d getCenter() const { return _center; }

    double getRadius() const { return _radius; }

    double pointDistanceSquared(MathLib::Point3d const& pnt) const;
 
private:
    MinimalBoundingSphere();

    static MinimalBoundingSphere recurseCalculation(
        std::vector<MathLib::Point3d*> sphere_points,
        std::size_t start_idx,
        std::size_t length,
        std::size_t n_boundary_points);
 
    double _radius{-1};
    MathLib::Point3d _center{{std::numeric_limits<double>::max(),
                              std::numeric_limits<double>::max(),
                              std::numeric_limits<double>::max()}};
};
 
}  // namespace GeoLib