00001 /*************************************************************************** 00002 ** 00003 ** This file is part of warand. 00004 ** 00005 ** This file may be distributed and/or modified under the terms of the 00006 ** GNU General Public License version 2 or 3 as published by the Free 00007 ** Software Foundation and appearing in the file LICENSE.GPL included 00008 ** in the packaging of this file. 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 General Public License for 00013 ** more details. 00014 ** 00015 ** You should have received a copy of the GNU General Public License 00016 ** along with this program. If not, see <http://www.gnu.org/licenses/>. 00017 ** 00018 ** See http://www.geopsy.org for more information. 00019 ** 00020 ** Created : 2008-07-13 00021 ** Authors: 00022 ** Marc Wathelet 00023 ** Marc Wathelet (LGIT, Grenoble, France) 00024 ** 00025 ***************************************************************************/ 00026 00027 #ifndef PHASESHIFTER_H 00028 #define PHASESHIFTER_H 00029 00030 #include <QGpCoreTools.h> 00031 00032 class PhaseShifter : public AbstractFunction2 00033 { 00034 public: 00035 PhaseShifter(ArrayStations * array); 00036 ~PhaseShifter(); 00037 00038 void setMaximumSlowness(double slow) {_squaredSlowMax=slow*slow;} 00039 double maximumSlowness() const {return sqrt(_squaredSlowMax);} 00040 00041 void setMaximumWavenumber(double k) {_squaredKmaxGrid=k*k;} 00042 00043 void setOmega(double omega); 00044 void setCrossSpectrum(ComplexMatrix m) {_crossSpectrum=m;} 00045 00046 virtual void initGrid(int n); 00047 virtual void initGrid(double kx, double ky, int index); 00048 virtual double value(double kx, double ky) const; 00049 virtual double value(double kx, double ky, int index) const; 00050 private: 00051 // Maximum K^2 from grid limits 00052 double _squaredKmaxGrid; 00053 // Maximum slowness^2 00054 double _squaredSlowMax; 00055 // Maximum K^2 from grid limits and maximum slowness 00056 double _squaredKmaxTotal; 00057 00058 QVector<Point2D> _relativeCoordinates; 00059 QVector< ComplexMatrix > _gridShifts; 00060 QVector< ComplexMatrix > _gridShiftConjugates; 00061 ComplexMatrix _crossSpectrum; 00062 }; 00063 00064 #endif // PHASESHIFTER_H
