00001 /*************************************************************************** 00002 ** 00003 ** This file is part of QGpCoreWave. 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 THEORETICALFK_H 00029 #define THEORETICALFK_H 00030 00031 #include <math.h> 00032 00033 #include <QGpCoreTools.h> 00034 #include "QGpCoreWaveDLLExport.h" 00035 00036 namespace QGpCoreWave { 00037 00038 class QGPCOREWAVE_EXPORT TheoreticalFK : public AbstractFunction2 00039 { 00040 public: 00041 TheoreticalFK(const QVector<Point2D>& stations); 00042 00043 inline virtual double value(double kx, double ky) const; 00044 const QVector<Point2D>& stations() const {return _stations;} 00045 private: 00046 QVector<Point2D> _stations; 00047 double _fac; 00048 }; 00049 00050 inline double TheoreticalFK::value(double kx, double ky) const 00051 { 00052 TRACE; 00053 double csum=0.0, ssum=0.0, val; 00054 for(int i=_stations.count()-1; i>=0; i--) { 00055 const Point2D& s=_stations[i]; 00056 val=kx*s.x()+ky*s.y(); 00057 csum+=cos(val); 00058 ssum+=sin(val); 00059 } 00060 return (csum*csum+ssum*ssum)*_fac; 00061 } 00062 00063 } // namespace QGpCoreWave 00064 00065 #endif // THEORETICALFK_H
