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 : 2007-11-12 00022 ** Authors : 00023 ** Marc Wathelet 00024 ** Marc Wathelet (LGIT, Grenoble, France) 00025 ** 00026 ***************************************************************************/ 00027 00028 #ifndef DIRECTIONALSEARCH_H 00029 #define DIRECTIONALSEARCH_H 00030 00031 #include "QGpCoreToolsDLLExport.h" 00032 #include "Function2Search.h" 00033 00034 namespace QGpCoreTools { 00035 00036 class QGPCORETOOLS_EXPORT DirectionalSearch : public Function2Search 00037 { 00038 public: 00039 DirectionalSearch(); 00040 ~DirectionalSearch(); 00041 00042 void setAzimuth(double azimuth); 00043 void setGrid(double minR, double maxR, double dR); 00044 void setFunction(AbstractFunction2 * f) {_function=f;} 00045 00046 double globalMax(); 00047 Function2MaximaIterator localMax(int nMax, double absThres=0.0, double relThres=0.0); 00048 void setRelativePrecision(double p) {_precision=p; _relative=true;} 00049 void setAbsolutePrecision(double p) {_precision=p; _relative=false;} 00050 00051 double value(double r) const {return _function->value(_cx*r, _sy*r);} 00052 protected: 00053 double _cx, _sy; 00054 int _nr; 00055 double * _r; 00056 double _precision; 00057 private: 00058 double refineMax(double minR, double maxR); 00059 inline void checkForMax(double r, double& maxVal, double& maxR); 00060 }; 00061 00062 } // namespace QGpCoreTools 00063 00064 #endif // DIRECTIONALSEARCH_H
