diff --git a/src/Base/Precision.h b/src/Base/Precision.h new file mode 100644 index 0000000000..6e38a9a039 --- /dev/null +++ b/src/Base/Precision.h @@ -0,0 +1,175 @@ +/*************************************************************************** + * Copyright (c) 2022 Werner Mayer * + * * + * This file is part of the FreeCAD CAx development system. * + * * + * This library is free software; you can redistribute it and/or * + * modify it under the terms of the GNU Library General Public * + * License as published by the Free Software Foundation; either * + * version 2 of the License, or (at your option) any later version. * + * * + * This library 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 * + * GNU Library General Public License for more details. * + * * + * You should have received a copy of the GNU Library General Public * + * License along with this library; see the file COPYING.LIB. If not, * + * write to the Free Software Foundation, Inc., 59 Temple Place, * + * Suite 330, Boston, MA 02111-1307, USA * + * * + ***************************************************************************/ + + +#ifndef BASE_PRECISION_H +#define BASE_PRECISION_H + +#include + +namespace Base { + +// The methods are copied from OCC's Precision class +class Precision { +public: + /*! + * \brief Angular + * Returns the recommended precision value when checking the equality of two angles (given in radians). + * \return + */ + static double Angular() { return 1.e-12; } + + /*! + * \brief Confusion + * Returns the recommended precision value when checking coincidence of two points in real space. + * \return + */ + static double Confusion() { return 1.e-7; } + + /*! + * \brief SquareConfusion + * Returns square of \ref Confusion. + * \return + */ + static double SquareConfusion() { return Confusion() * Confusion(); } + + /*! + * \brief Intersection + * Returns the precision value in real space, frequently + * used by intersection algorithms to decide that a solution is reached. + * \return + */ + static double Intersection() { return Confusion() * 0.01; } + + /*! + * \brief Approximation + * Returns the precision value in real space, frequently used + * by approximation algorithms. + * \return + */ + static double Approximation() { return Confusion() * 10.0; } + + /*! + * \brief Parametric + * Convert a real space precision to a parametric space precision. + * \param P + * \param T + * \return + */ + static double Parametric (const double P, const double T) { return P / T; } + + /*! + * \brief PConfusion + * Returns a precision value in parametric space. + * \param T + * \return + */ + static double PConfusion (const double T) { return Parametric (Confusion(), T); } + + /*! + * \brief PConfusion + * Used to test distances in parametric space on a default curve. + * \return + */ + static double PConfusion() { return Parametric (Confusion()); } + + /*! + * \brief SquarePConfusion + * Returns square of \ref PConfusion. + * \return + */ + static double SquarePConfusion() { return PConfusion() * PConfusion(); } + + /*! + * \brief PIntersection + * Returns a precision value in parametric space, which may be used by intersection algorithms, + * to decide that a solution is reached. + * \param T + * \return + */ + static double PIntersection (const double T) { return Parametric(Intersection(),T); } + + /*! + * \brief PApproximation + * Returns a precision value in parametric space, which may be used by approximation algorithms. + * \param T + * \return + */ + static double PApproximation (const double T) { return Parametric(Approximation(),T); } + + /*! + * \brief Parametric + * Convert a real space precision to a parametric space precision on a default curve. + * \param P + * \return + */ + static double Parametric (const double P) { return Parametric (P, 100.0); } + + /*! + * \brief PIntersection + * Used for Intersections in parametric space on a default curve. + * \return + */ + static double PIntersection() { return Parametric (Intersection()); } + + /*! + * \brief PApproximation + * Used for Approximations in parametric space on a default curve. + * \return + */ + static double PApproximation() { return Parametric (Approximation()); } + + /*! + * \brief IsInfinite + * Returns True if R may be considered as an infinite number. Currently Abs(R) > 1e100 + * \param R + * \return + */ + static double IsInfinite (const double R) { return std::fabs (R) >= (0.5 * Precision::Infinite()); } + + /*! + * \brief IsPositiveInfinite + * Returns True if R may be considered as a positive infinite number. Currently R > 1e100 + * \param R + * \return + */ + static double IsPositiveInfinite (const double R) { return R >= (0.5 * Precision::Infinite()); } + + /*! + * \brief IsNegativeInfinite + * Returns True if R may be considered as a negative infinite number. Currently R < -1e100 + * \param R + * \return + */ + static bool IsNegativeInfinite (const double R) { return R <= -(0.5 * Precision::Infinite()); } + + /*! + * \brief Infinite + * Returns a big number that can be considered as infinite. Use -Infinite() for a negative big number. + * \return + */ + static double Infinite() { return 2.e+100; } +}; + +} + +#endif // BASE_PRECISION_H