MinimalBoundingSphere.h

Go to the documentation of this file.

 
#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