QGpCompatibility/CompatDataPoint.h

Go to the documentation of this file.

/***************************************************************************
**
**  This file is part of QGpCompatibility.
**
**  This file may be distributed and/or modified under the terms of the
**  GNU General Public License version 2 or 3 as published by the Free
**  Software Foundation and appearing in the file LICENSE.GPL included
**  in the packaging of this file.
**
**  This file is distributed in the hope that it will be useful, but WITHOUT
**  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
**  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for
**  more details.
**
**  You should have received a copy of the GNU General Public License
**  along with this program. If not, see <http://www.gnu.org/licenses/>.
**
**  See http://www.geopsy.org for more information.
**
**  Created : 2004-10-21
**  Authors:
**    Marc Wathelet
**    Marc Wathelet (ULg, Liège, Belgium)
**    Marc Wathelet (LGIT, Grenoble, France)
**
***************************************************************************/

#ifndef COMPATDataPoint_H
#define COMPATDataPoint_H

#include <QGpCoreTools.h>
#include "QGpCompatibilityDLLExport.h"

namespace QGpCompatibility {

class QGPCOMPATIBILITY_EXPORT CompatDataPoint : public XMLClass
{
public:
  CompatDataPoint () {}
  inline CompatDataPoint (double mean, double stddev, double _weight);

  inline void operator=(const CompatDataPoint& o);

  double& mean() {return _mean;}
  double& stddev() {return _stddev;}
  double& weight() {return _weight;}

  const double& mean() const {return _mean;}
  const double& stddev() const {return _stddev;}
  const double& weight() const {return _weight;}

  inline double misfitOrder1(double value, double invalid);
protected:
  double _mean;
  double _stddev;
  double _weight;
protected:
  virtual const QString& xml_tagName() const {return xmlDataPointTag;}
  static const QString xmlDataPointTag;
  virtual void xml_writeProperties(XML_WRITEPROPERTIES_ARGS) const;
  virtual XMLMember xml_member(XML_MEMBER_ARGS);
  virtual bool xml_setProperty(XML_SETPROPERTY_ARGS);
};

class QGPCOMPATIBILITY_EXPORT CompatVDataPoint : public CompatDataPoint
{
public:
  CompatVDataPoint () {}
  inline CompatVDataPoint (double mean, double stddev, double _weight);

  inline void operator=(const CompatVDataPoint& o);
  inline void operator=(const CompatDataPoint& o);

  void setValue(double val) {_value=val;}
  const double& value() const {return _value;}

  inline double misfitOrder2(int& nMeans, int& nValues, double invalid);
protected:
  double _value;
};

inline CompatDataPoint::CompatDataPoint(double mean, double stddev, double weight)
{
  TRACE;
  _mean=mean;
  _stddev=stddev;
  _weight=weight;
}

inline void CompatDataPoint::operator=(const CompatDataPoint& o)
{
  TRACE;
  _mean=o._mean;
  _stddev=o._stddev;
  _weight=o._weight;
}

inline double CompatDataPoint::misfitOrder1(double value, double invalid)
{
  TRACE;
  if(value!=invalid) {
    double diff;
    if(_stddev)
      diff=(_mean - value)/_stddev;
    else
      diff=(_mean - value)/_mean;
    return (diff > 0) ? diff : -diff;
  }
  return 1e99;
}

inline CompatVDataPoint::CompatVDataPoint(double mean, double stddev, double weight)
    : CompatDataPoint(mean, stddev, weight)
{
  TRACE;
  _value=mean;
}

inline void CompatVDataPoint::operator=(const CompatVDataPoint& o)
{
  TRACE;
  _value=o._value;
  CompatDataPoint::operator=(o);
}

inline void CompatVDataPoint::operator=(const CompatDataPoint& o)
{
  TRACE;
  CompatDataPoint::operator=(o);
  _value=_mean;
}

inline double CompatVDataPoint::misfitOrder2(int& nMeans, int& nValues, double invalid)
{
  TRACE;
  if(_mean!=invalid) {
    if(_value!=invalid) {
      double diff;
      if(_stddev)
        diff=(_mean - _value)/_stddev;
      else
        diff=(_mean - _value)/_mean;
      return diff*diff;
    } else {
      nValues--;
      return 0.0;
    }
  } else {
    nMeans--;
    nValues--;
    return 0.0;
  }
}

QGPCOMPATIBILITY_EXPORT QDataStream& operator<< (QDataStream& s, const CompatDataPoint& p);
QGPCOMPATIBILITY_EXPORT QDataStream& operator>> (QDataStream& s, CompatDataPoint& p);

class QGPCOMPATIBILITY_EXPORT CompatDataPointVector : public QVector<CompatDataPoint>,
      public XMLClass
{
public:
  CompatDataPointVector() {}
  virtual ~CompatDataPointVector() {}
  QVector<double> * meanVector() const;
  QVector<double> * stddevVector() const;
  QVector<double> * weightVector() const;
protected:
  virtual const QString& xml_tagName() const {return xmlDataPointVectorTag;}
  static const QString xmlDataPointVectorTag;
  virtual void xml_writeChildren(XML_WRITECHILDREN_ARGS) const;
  virtual XMLMember xml_member(XML_MEMBER_ARGS);
};

class QGPCOMPATIBILITY_EXPORT CompatVDataPointVector : public QVector<CompatVDataPoint>,
      public XMLClass
{
public:
  CompatVDataPointVector() {}
  virtual ~CompatVDataPointVector() {}
  QVector<double> * valueVector() const;
  QVector<double> * meanVector() const;
  QVector<double> * stddevVector() const;
  QVector<double> * weightVector() const;
protected:
  virtual const QString& xml_tagName() const {return xmlVDataPointVectorTag;}
  static const QString xmlVDataPointVectorTag;
  virtual void xml_writeChildren(XML_WRITECHILDREN_ARGS) const;
  virtual XMLMember xml_member(XML_MEMBER_ARGS);
};

} // namespace QGpCompatibility

#endif // COMPATDATAPOINT_H