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 : 2002-10-11 00022 ** Authors : 00023 ** Marc Wathelet 00024 ** Marc Wathelet (ULg, Liège, Belgium) 00025 ** Marc Wathelet (LGIT, Grenoble, France) 00026 ** 00027 ***************************************************************************/ 00028 00029 #ifndef GaussDistribution_H 00030 #define GaussDistribution_H 00031 00032 #include "QGpCoreToolsDLLExport.h" 00033 00034 namespace QGpCoreTools { 00035 00036 class QGPCORETOOLS_EXPORT GaussDistribution 00037 { 00038 public: 00043 GaussDistribution(double precision=1e-4); 00050 GaussDistribution(double mean, double stddev,double precision=1e-4); 00052 ~GaussDistribution(); 00057 double prob(double a, double b) { 00058 return (_calcErrorFunction((b-_mean)*_inv_stddev) 00059 -_calcErrorFunction((a-_mean)*_inv_stddev)); 00060 } 00062 double part(double a) { 00063 return _calcErrorFunction((a-_mean)*_inv_stddev); 00064 } 00066 double partSigma(double a) {return _calcErrorFunction(a);} 00068 double rho(double val); 00069 void setMean(double m) {_mean=m;} 00070 double mean() {return _mean;} 00071 void setStddev(double s) {_inv_stddev=1.0/s;} 00072 double stddev() {return 1.0/_inv_stddev;} 00073 bool isValid() {return _inv_stddev!=0;} 00074 protected: 00076 double _calcErrorFunction(double z); 00077 double _mean; 00078 double _inv_stddev; 00080 double _precision; 00081 }; 00082 00083 } // namespace QGpCoreTools 00084 00085 #endif // GAUSSDISTRIBUTION_H
