00001 /*************************************************************************** 00002 ** 00003 ** This file is part of QGpCoreTools. 00004 ** 00005 ** This library is free software; you can redistribute it and/or 00006 ** modify it under the terms of the GNU Lesser General Public 00007 ** License as published by the Free Software Foundation; either 00008 ** version 2.1 of the License, or (at your option) any later version. 00009 ** 00010 ** This file is distributed in the hope that it will be useful, but WITHOUT 00011 ** ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or 00012 ** FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public 00013 ** License for more details. 00014 ** 00015 ** You should have received a copy of the GNU Lesser General Public 00016 ** License along with this library; if not, write to the Free Software 00017 ** Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA 00018 ** 00019 ** See http://www.geopsy.org for more information. 00020 ** 00021 ** Created : 2004-08-23 00022 ** Authors : 00023 ** Marc Wathelet 00024 ** Marc Wathelet (ULg, Liège, Belgium) 00025 ** Marc Wathelet (LGIT, Grenoble, France) 00026 ** 00027 ***************************************************************************/ 00028 00029 #ifndef GridSearch_H 00030 #define GridSearch_H 00031 00032 #include "Point2D.h" 00033 #include "Function2Search.h" 00034 00035 namespace QGpCoreTools { 00036 00037 class GridFunction2; 00038 00039 class QGPCORETOOLS_EXPORT GridSearch : public Function2Search 00040 { 00041 public: 00042 GridSearch(); 00043 ~GridSearch(); 00044 00045 void setGrid(double minX, double maxX, double dX, 00046 double minY, double maxY, double dY); 00047 00048 int nx() const {return _nx;} 00049 int ny() const {return _ny;} 00050 double x(int i) const {return _x[i];} 00051 double y(int i) const {return _y[i];} 00052 00053 double globalMax(); 00054 Function2MaximaIterator localMax(int nMax, double absThres=0.0, double relThres=0.0); 00055 void setRelativePrecision(double p) {_precisionX=p;_precisionY=p; _relative=true;} 00056 void setAbsolutePrecision(double p) {_precisionX=p;_precisionY=p; _relative=false;} 00057 protected: 00058 int _nx, _ny; 00059 double * _x, * _y; 00060 double _precisionX, _precisionY; 00061 private: 00062 // For 2D grids 00063 double refineMax(double minX, double maxX, double minY, double maxY); 00064 inline void checkForMax(double x, double y, double& maxVal, double& maxX, double& maxY); 00065 // For 1D grids 00066 double refineMax(double minX, double maxX); 00067 inline void checkForMax(double x, double& maxVal, double& maxX); 00068 }; 00069 00070 } // namespace QGpCoreTools 00071 00072 #endif // GRIDSEARCH_H
