linear.h

Go to the documentation of this file.

#ifndef linear_h
#define linear_h

/*+
________________________________________________________________________

 (C) dGB Beheer B.V.; (LICENSE) http://opendtect.org/OpendTect_license.txt
 Author:        A.H.Bril
 Date:          Jan 2005
 RCS:           $Id$
________________________________________________________________________

-*/


#include "algomod.h"
#include "mathfunc.h"
namespace Geom { template <class T> class Point2D; }


template <class T>
mClass(Algo) LineParameters : public MathFunction<T,T>
{
public:
                LineParameters( T i0=0, T i1=0 )
                : a0(i0), ax(i1)                {}
 
    inline T    getValue( T x ) const
                        { return a0 + ax * x; }
    inline T    getXValue( T y ) const
                        { return ax ? (y - a0) / ax : 0; }
    inline T    getProjectedX( T x, T y ) const
                        { return (x + ax * (y - a0)) / (1 + ax * ax); }
    inline T    getValue( const T* x ) const
                        { return getValue(*x); }
 
    T           a0, ax;
};

typedef LineParameters<float> LinePars;


template <class T>
mClass(Algo) PlaneParameters : public MathXYFunction<T,T>
{
public:
                PlaneParameters( T i0=0, T i1=0, T i2=0 )
                : a0(i0), ax(i1), ay(i2)        {}

    inline T    getValue( T x, T y ) const
                { return a0 + ax * x + ay * y; }
    inline T    getValue( const T* x ) const
                        { return getValue(x[0],x[1]); }

    T           a0, ax, ay;

};

typedef PlaneParameters<float> PlanePars;


mExpClass(Algo) LinStats2D
{
public:
                LinStats2D() : corrcoeff(0)     {}

    LinePars    lp;             
    LinePars    sd;             
    float       corrcoeff;      

    void        use(const float*,const float*,int nrpts);
    void        use(const Geom::Point2D<float>*,int nrpts);
};


mExpClass(Algo) LinStats3D
{
public:
                LinStats3D() : corrcoeff(0)     {}

    PlanePars   pp;             
    PlanePars   sd;             
    float       corrcoeff;      

};


#endif