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 : 2004-10-22 00022 ** Authors : 00023 ** Marc Wathelet 00024 ** Marc Wathelet (ULg, Liège, Belgium) 00025 ** Marc Wathelet (LGIT, Grenoble, France) 00026 ** 00027 ***************************************************************************/ 00028 00029 #ifndef MODALELLIPTICITY_H 00030 #define MODALELLIPTICITY_H 00031 00032 #include <QGpCoreTools.h> 00033 00034 #include "QGpCoreWaveDLLExport.h" 00035 #include "ModalStorage.h" 00036 #include "Rayleigh.h" 00037 00038 namespace QGpCoreWave { 00039 00040 class Dispersion; 00041 00042 class QGPCOREWAVE_EXPORT Ellipticity : public ModalStorage 00043 { 00044 TRANSLATIONS("EllipticityFactory") 00045 public: 00046 Ellipticity(int nModes, const QVector<double>* x) : ModalStorage(nModes, x) {} 00047 00048 QList<double> peaks(int modeIndex, Dispersion& disp, Rayleigh * waveModel); 00049 double peakMisfit(int modeIndex, const RealStatisticalValue& val, Dispersion& disp, Rayleigh * waveModel); 00050 private: 00051 inline double refineMax(int modeIndex, int iOmega, Dispersion& disp, 00052 RootSolver<Rayleigh>& , bool& badSampling); 00053 }; 00054 00055 } // namespace QGpCoreWave 00056 00057 #endif // ELLIPTICITY_H