LsqSqQuadFit< T > Class Template Reference

#include <LsqQuadraticFit.h>

Classes
struct	Coeff

Public Member Functions
CUDA_SUPPORTED_FUNC	LsqSqQuadFit (uint numPts, const T *xval, const T *yval, const T *weights=0)
CUDA_SUPPORTED_FUNC	LsqSqQuadFit ()
CUDA_SUPPORTED_FUNC void	calculate (uint numPts, const T *X, const T *Y, const T *weights)
CUDA_SUPPORTED_FUNC T	compute (T pos)
CUDA_SUPPORTED_FUNC T	computeDeriv (T pos)
CUDA_SUPPORTED_FUNC T	maxPos ()
CUDA_SUPPORTED_FUNC T	vertexForm ()

Public Attributes
T	a
T	b
T	c
T	d
T	h
T	k
float	xoffset

Private Member Functions
CUDA_SUPPORTED_FUNC Coeff	computeSums (const T *X, const T *Y, const T *weights, uint numPts)

Detailed Description

template<typename T>
class LsqSqQuadFit< T >

Definition at line 9 of file LsqQuadraticFit.h.

Constructor & Destructor Documentation

LsqSqQuadFit() [1/2]

template<typename T>

CUDA_SUPPORTED_FUNC LsqSqQuadFit< T >::LsqSqQuadFit ( uint  numPts, const T *  xval, const T *  yval, const T *  weights = 0  )

inline

Definition at line 27 of file LsqQuadraticFit.h.

    {
        calculate(numPts, xval, yval, weights);
        xoffset =0;
    }

LsqSqQuadFit() [2/2]

template<typename T>

CUDA_SUPPORTED_FUNC LsqSqQuadFit< T >::LsqSqQuadFit ( )

inline

Definition at line 33 of file LsqQuadraticFit.h.

    {
        a=b=c=d=0;
        xoffset =0;
    }

Member Function Documentation

calculate()

template<typename T>

CUDA_SUPPORTED_FUNC void LsqSqQuadFit< T >::calculate ( uint  numPts, const T *  X, const T *  Y, const T *  weights  )

inline

Definition at line 39 of file LsqQuadraticFit.h.

    {
        Coeff co = computeSums(X, Y, weights, numPts);
        co.abc(a,b,c,d);
    }

compute()

template<typename T>

CUDA_SUPPORTED_FUNC T LsqSqQuadFit< T >::compute ( T  pos )

inline

Definition at line 45 of file LsqQuadraticFit.h.

    {
        pos -= xoffset;
        return a*pos*pos + b*pos + c;
    }

computeDeriv()

template<typename T>

CUDA_SUPPORTED_FUNC T LsqSqQuadFit< T >::computeDeriv ( T  pos )

inline

Definition at line 51 of file LsqQuadraticFit.h.

    {
        pos -= xoffset;
        return 2*a*pos + b;
    }

computeSums()

template<typename T>

CUDA_SUPPORTED_FUNC Coeff LsqSqQuadFit< T >::computeSums ( const T *  X, const T *  Y, const T *  weights, uint  numPts  )

inlineprivate

Definition at line 72 of file LsqQuadraticFit.h.

    {
        //notation sjk to mean the sum of x_i^j*y_i^k. 
    /*    s40 = getSx4(); //sum of x^4
        s30 = getSx3(); //sum of x^3
        s20 = getSx2(); //sum of x^2
        s10 = getSx();  //sum of x
        

        s21 = getSx2y(); //sum of x^2*y
        s11 = getSxy();  //sum of x*y
        s01 = getSy();   //sum of y
        */

        if (weights) {
            T Sx = 0, Sy = 0;
            T Sx2 = 0, Sx3 = 0;
            T Sxy = 0, Sx4=0, Sx2y=0;
            T Sw = 0;
            for (uint i=0;i<numPts;i++)
            {
                T x = X[i];
                T y = Y[i];
                Sx += x*weights[i];
                Sy += y*weights[i];
                T sq = x*x;
                Sx2 += x*x*weights[i];
                Sx3 += sq*x*weights[i];
                Sx4 += sq*sq*weights[i];
                Sxy += x*y*weights[i];
                Sx2y += sq*y*weights[i];
                Sw += weights[i];
            }

            Coeff co;
            co.s10 = Sx; co.s20 = Sx2; co.s30 = Sx3; co.s40 = Sx4;
            co.s01 = Sy; co.s11 = Sxy; co.s21 = Sx2y;
            co.s00 = Sw;
            return co;
        } else {
            T Sx = 0, Sy = 0;
            T Sx2 = 0, Sx3 = 0;
            T Sxy = 0, Sx4=0, Sx2y=0;
            for (uint i=0;i<numPts;i++)
            {
                T x = X[i];
                T y = Y[i];
                Sx += x;
                Sy += y;
                T sq = x*x;
                Sx2 += x*x;
                Sx3 += sq*x;
                Sx4 += sq*sq;
                Sxy += x*y;
                Sx2y += sq*y;
            }

            Coeff co;
            co.s10 = Sx; co.s20 = Sx2; co.s30 = Sx3; co.s40 = Sx4;
            co.s01 = Sy; co.s11 = Sxy; co.s21 = Sx2y;
            co.s00 = numPts;
            return co;
        }
    }

maxPos()

template<typename T>

CUDA_SUPPORTED_FUNC T LsqSqQuadFit< T >::maxPos ( )

inline

Definition at line 57 of file LsqQuadraticFit.h.

    {
        return -b/(2*a);
    }

vertexForm()

template<typename T>

CUDA_SUPPORTED_FUNC T LsqSqQuadFit< T >::vertexForm ( )

inline

Definition at line 62 of file LsqQuadraticFit.h.

    {
        // Find equation in form a(x-h)^2+k
        h = -b/(2*a);
        k = c - a*h*h;      
        return k;
    }

Member Data Documentation

a

template<typename T>

T LsqSqQuadFit< T >::a

Definition at line 12 of file LsqQuadraticFit.h.

b

template<typename T>

T LsqSqQuadFit< T >::b

Definition at line 12 of file LsqQuadraticFit.h.

c

template<typename T>

T LsqSqQuadFit< T >::c

Definition at line 12 of file LsqQuadraticFit.h.

d

template<typename T>

T LsqSqQuadFit< T >::d

Definition at line 12 of file LsqQuadraticFit.h.

h

template<typename T>

T LsqSqQuadFit< T >::h

Definition at line 12 of file LsqQuadraticFit.h.

k

template<typename T>

T LsqSqQuadFit< T >::k

Definition at line 12 of file LsqQuadraticFit.h.

xoffset

template<typename T>

float LsqSqQuadFit< T >::xoffset

Definition at line 13 of file LsqQuadraticFit.h.

---

The documentation for this class was generated from the following file:

• cputrack/LsqQuadraticFit.h