PARP Research Group University of Murcia, Spain


src/qvmath/qvmultidimensionalminimizer.h

Go to the documentation of this file.
00001 /*
00002  *      Copyright (C) 2007, 2008. PARP Research Group.
00003  *      <http://perception.inf.um.es>
00004  *      University of Murcia, Spain.
00005  *
00006  *      This file is part of the QVision library.
00007  *
00008  *      QVision is free software: you can redistribute it and/or modify
00009  *      it under the terms of the GNU Lesser General Public License as
00010  *      published by the Free Software Foundation, version 3 of the License.
00011  *
00012  *      QVision is distributed in the hope that it will be useful,
00013  *      but WITHOUT ANY WARRANTY; without even the implied warranty of
00014  *      MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
00015  *      GNU Lesser General Public License for more details.
00016  *
00017  *      You should have received a copy of the GNU Lesser General Public
00018  *      License along with QVision. If not, see <http://www.gnu.org/licenses/>.
00019  */
00020 
00024 
00025 #ifndef QVMULTIDIMENSIONALFUNCTIONMINIMIZER_H
00026 #define QVMULTIDIMENSIONALFUNCTIONMINIMIZER_H
00027 
00028 #include <qvmath.h>
00029 #include <QVVector>
00030 
00031 #include <gsl/gsl_errno.h>
00032 #include <gsl/gsl_math.h>
00033 #include <gsl/gsl_multimin.h>
00034 
00035 #ifndef DOXYGEN_IGNORE_THIS
00036 
00044 class QVMultidimensionalMinimizer
00045         {
00046         private:
00047                 int iterations, vectorSize;
00048                 QVVector actual, step;
00049                 double minimumValue;
00050 
00051                 static double static_function (const gsl_vector *v, void *multidimensionalOptimizerPtr);
00052 
00053         public: 
00054                 virtual double function (const QVVector &v) = 0;
00055 
00056                 QVMultidimensionalMinimizer(const int vectorSize, const double stepSize):
00057                         iterations(0), vectorSize(vectorSize), actual(vectorSize,0), step(vectorSize, stepSize) { }
00058 
00059                 void setStartingVector(const QVVector &vector)  { actual = vector; }
00060                 void setStepVector(const QVVector &vector)      { step = vector; }
00061                 void setMinimum(const QVVector &vector)         { actual = vector; }
00062                 void setMinimumValue(const double value)        { minimumValue = value; }
00063 
00064                 const QVVector & getMinimum()           const   { return actual; }
00065                 const double & getMinimumValue()        const   { return minimumValue; }
00066                 const int & getIterations()             const   { return iterations; }
00067 
00068                 const QList<QVVector> gridIterate(const QVVector &firstVertex, const QVVector &secondVertex, const int divisions = 5, const int maxIterations = 100, const double maxError = 1e-5);
00069                 bool iterate(const int maxIterations = 100, const double maxError = 1e-5);
00070         };
00071 
00072 #endif
00073 #endif
00074 



QVision framework. PARP research group, copyright 2007, 2008.