PARP Research Group

Universidad de Murcia

---

src/qvmath/qvmath.h Go to the documentation of this file. /* * Copyright (C) 2007, 2008, 2009, 2010, 2011, 2012. PARP Research Group. * <http://perception.inf.um.es> * University of Murcia, Spain. * * This file is part of the QVision library. * * QVision is free software: you can redistribute it and/or modify * it under the terms of the GNU Lesser General Public License as * published by the Free Software Foundation, version 3 of the License. * * QVision 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 Lesser General Public License for more details. * * You should have received a copy of the GNU Lesser General Public * License along with QVision. If not, see <http://www.gnu.org/licenses/>. */ #ifndef QVMATH_H #define QVMATH_H #ifndef M_SQRT1_2 #define M_SQRT1_2 7.0710678118654752440E-1 #endif #ifndef M_PI #define M_PI 3.14159265358979323846264338327950288 /* pi */ #endif #ifndef M_PI_2 #define M_PI_2 1.5707963267948966192E0 #endif #ifndef M_PI_4 #define M_PI_4 0.785398163397448309615660845819875721 /* pi/4 */ #endif #ifdef max #undef max #endif #ifdef min #undef min #endif #ifdef _MSC_VER #include <float.h> // for _isnan() on VC++ #define isnan(x) _isnan(x) // VC++ uses _isnan() instead of isnan() #define isinf(x) _isinf(x) #else #include <cmath> // for isnan() and isinf() everywhere else extern "C" int isnan(double); extern "C" int isinf(double); #endif #include <iostream> #include <math.h> #define _USE_MATH_DEFINES #include <cmath> #include <qvdefines.h> // Get BLAS version #ifdef MKL_AVAILABLE #include <mkl.h> #include <mkl_lapack.h> #include <mkl_dss.h> #include <mkl_rci.h> #include <mkl_spblas.h> #elif GSL_AVAILABLE // If the MKL is not available, use the BLAS implementation from the GSL. #include <gsl/gsl_blas.h> // #include <gsl/gsl_blas.h> // #include <gsl/gsl_math.h> // #include <gsl/gsl_eigen.h> // #include <gsl/gsl_linalg.h> #endif #ifdef GSL_AVAILABLE #include <gsl/gsl_math.h> #include <gsl/gsl_eigen.h> #include <gsl/gsl_linalg.h> #include <gsl/gsl_errno.h> #include <gsl/gsl_math.h> #include <gsl/gsl_min.h> #include <gsl/gsl_multimin.h> #include <gsl/gsl_multifit_nlin.h> #endif #include <QPair> #include <QPointF> #include <QVComplex> #include <QVVector> #include <qvmath/qvbitcount.h> #ifndef POW2 #define POW2(X) ((X)*(X)) #endif #ifndef SIGN #define SIGN(X) (((X)>=0)?1:-1) #endif #ifndef ABS #define ABS(X) (((X)>0)?(X):(-(X))) #endif #ifndef MIN #define MIN(X,Y) (((X)<(Y))?(X):(Y)) #endif #ifndef MAX #define MAX(X,Y) (((X)>(Y))?(X):(Y)) #endif #ifndef RANDOM #define RANDOM(MIN, MAX) (rand()%((MAX)-(MIN)+1) + (MIN)) #endif #ifndef PI #define PI (3.14159265358979323846) #endif #ifndef EPSILON #define EPSILON (10E-6) #endif int qvFactorial(const int n); double qvCombination(const int setRange, const int subsetRange); double qvAngle(const QPointF &point1, const QPointF &point2); #ifndef DOXYGEN_IGNORE_THIS double qvAngle(const QPointF &); double qvClockWiseAngle(const QPointF &, const QPointF &); int qvRandom(const int minValue, const int maxValue); double norm2(const QPointF &p); double random(const double min, const double max); #endif //#ifndef DOXYGEN_IGNORE_THIS //void qvSolveCubicPolynomialAux(const double a0, const double a1, const double a2, const double a3, double cg[6]); //#endif QVector<QVComplex> qvSolveCubicPolynomial(const double a0, const double a1, const double a2, const double a3); double qvPointLineDistance(const QVVector &, const QPointF &); double qvSymmetricFloor( const double value ); #ifndef DOXYGEN_IGNORE_THIS double relativeEuclideanDistance(const double &v1, const double &v2); double relativeEuclideanDistance(const QVVector &v1, const QVVector &v2); #endif // DOXYGEN_IGNORE_THIS #endif // QVMATH_H

---

QVision framework. PARP research group. Copyright © 2007, 2008, 2009, 2010, 2011.