JACOUSTICS::JKatoomba< JEstimator > Struct Template Reference

Template specialisation of fit function of acoustic model based on linear approximation. More...

#include <JKatoomba.hh>

Inheritance diagram for JACOUSTICS::JKatoomba< JEstimator >:

Public Member Functions
	JKatoomba (const JDetector &detector, const JSoundVelocity &velocity)
	Constructor. More...
template<class T >
const JModel &	operator() (T __begin, T __end) const
	Get start values of string parameters. More...
template<class JPDF_t >
double	getToA (const JModel &model, const JHit< JPDF_t > &hit) const
	Get estimated time-of-arrival for given hit. More...

Public Attributes
const JDetector &	detector
const JSoundVelocity &	velocity
JLANG::JSharedPointer < JMEstimator >	estimator
	M-Estimator function. More...

Private Attributes
JModel	value

Detailed Description

template<>
struct JACOUSTICS::JKatoomba< JEstimator >

Template specialisation of fit function of acoustic model based on linear approximation.

Definition at line 153 of file JKatoomba.hh.

Constructor & Destructor Documentation

JACOUSTICS::JKatoomba< JEstimator >::JKatoomba ( const JDetector &  detector, const JSoundVelocity &  velocity  )

inline

Constructor.

Parameters

detector	detector
velocity	sound velocity

Definition at line 162 of file JKatoomba.hh.

                                              :
      JKatoomba_t(detector, velocity)
    {};

Member Function Documentation

template<class T >

const JModel& JACOUSTICS::JKatoomba< JEstimator >::operator() ( T  __begin, T  __end  ) const

inline

Get start values of string parameters.

Note that this method may throw an exception.

Parameters

__begin	begin of hits
__end	end of hits

Returns

model

Definition at line 177 of file JKatoomba.hh.

    {
      using namespace std;
      using namespace JPP;

      value = JModel(__begin, __end);                 // set up global model according data

      struct H_t :                                    // H-equation as per hit
            public JModel::JString,
            public JModel::JEmitter
      {} H;

      struct { int t1; int tx; int ty; } i;           // indices of H-equation in global model 

      JMatrixNS      V(value.getN());                 // V = H^T x H
      vector<double> Y(value.getN(), 0.0);            // Y = H^T x y
      

      for (T hit = __begin; hit != __end; ++hit) {

        const JString&    string   = detector[hit->getString()];
        const JPosition3D position = string.getPosition(hit->getFloor());
        const double      Vi       = velocity.getInverseVelocity(hit->getDistance(position), hit->getZ(), position.getZ());

        const double      h1 = string[hit->getFloor()].getHeight();
        const JPosition3D p1 = string.getPosition()  -  hit->getPosition();
        const double      ds = sqrt(p1.getLengthSquared()  +  h1*h1  +  2.0*p1.getZ()*h1);
        const double      y  = hit->getValue()  -  Vi*ds;

        H.t1  =  1.0;
        H.tx  =  Vi * p1.getX() * h1 / ds;
        H.ty  =  Vi * p1.getY() * h1 / ds;
        
        i.t1  =  value.getIndex(hit->getEKey(),   &H_t::t1);
        i.tx  =  value.getIndex(hit->getString(), &H_t::tx);
        i.ty  =  value.getIndex(hit->getString(), &H_t::ty);

        V(i.t1, i.t1)  +=  H.t1 * H.t1;
        V(i.t1, i.tx)  +=  H.t1 * H.tx;
        V(i.t1, i.ty)  +=  H.t1 * H.ty;

        V(i.tx, i.t1)  +=  H.tx * H.t1;
        V(i.ty, i.t1)  +=  H.ty * H.t1;

        V(i.tx, i.tx)  +=  H.tx * H.tx;
        V(i.tx, i.ty)  +=  H.tx * H.ty;

        V(i.ty, i.tx)  +=  H.ty * H.tx;

        V(i.ty, i.ty)  +=  H.ty * H.ty;

        Y[i.t1]  +=  H.t1 * y;
        Y[i.tx]  +=  H.tx * y;
        Y[i.ty]  +=  H.ty * y;
      }

      // evaluate (H^T x H)^-1 x H^T x y 

      V.invert();

      for (size_t row = 0; row != value.getN(); ++row) {
        for (size_t col = 0; col != value.getN(); ++col) {
          value[row] += V(row,col) * Y[col];
        }
      }
      
      return value;
    }

template<class JPDF_t >

double JACOUSTICS::JKatoomba_t::getToA ( const JModel &  model, const JHit< JPDF_t > &  hit  ) const

inlineinherited

Get estimated time-of-arrival for given hit.

Parameters

model	model
hit	hit

Returns

time-of-arrival

Definition at line 62 of file JKatoomba.hh.

    {
      const JString&         string     = detector[hit.getString()];
      const JModel::JString& parameters = model.string[hit.getString()];
      const JPosition3D      position   = string.getPosition(parameters, hit.getFloor());

      const double D     = hit.getDistance(position);
      const double Vi    = velocity.getInverseVelocity(D, hit.getZ(), position.getZ());

      return model.emitter[hit.getEKey()].t1  +  D * Vi;
    }

Member Data Documentation

JModel JACOUSTICS::JKatoomba< JEstimator >::value

mutableprivate

Definition at line 247 of file JKatoomba.hh.

const JDetector& JACOUSTICS::JKatoomba_t::detector

inherited

Definition at line 75 of file JKatoomba.hh.

const JSoundVelocity& JACOUSTICS::JKatoomba_t::velocity

inherited

Definition at line 76 of file JKatoomba.hh.

JLANG::JSharedPointer<JMEstimator> JACOUSTICS::JKatoomba_t::estimator

inherited

M-Estimator function.

Definition at line 78 of file JKatoomba.hh.

---

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

• JKatoomba.hh