Line data    Source code

       1             : /**
       2             :    \file measure_iface.cpp
       3             :    Renamed measure_iface.cc 2023
       4             :    
       5             :    Bojan Nikolic <b.nikolic@mrao.cam.ac.uk>, <bojan@bnikolic.co.uk>
       6             : 
       7             : */
       8             : 
       9             : #include "measure_iface.h"
      10             : #include "model_iface.h"
      11             : #include "numalgo.h"
      12             : 
      13             : namespace LibAIR2 {
      14             : 
      15          81 :   ALMAMeasure::ALMAMeasure(WVRAtmoQuantModel &model):
      16          81 :     model(model)
      17             :   {
      18          81 :   }
      19             :   
      20         162 :   void ALMAMeasure::AddParams ( std::vector< Minim::DParamCtr > &pars )
      21             :   {
      22         162 :     model.AddParams(pars);
      23         162 :   }
      24             : 
      25          81 :   NormalNoise::NormalNoise(size_t n):
      26          81 :     thermNoise(n)
      27             :   {
      28          81 :   }
      29             :   
      30             : 
      31          81 :   AbsNormMeasure::AbsNormMeasure(WVRAtmoQuantModel &model):
      32             :     ALMAMeasure(model),
      33             :     NormalNoise(4),
      34          81 :     obs(4)
      35             :   {
      36          81 :   }
      37             : 
      38           0 :   void AbsNormMeasure::modelObs(void)
      39             :   {
      40           0 :     model.eval(obs);
      41           0 :   }
      42             : 
      43    77889705 :   double AbsNormMeasure::lLikely (void) const
      44             :   {
      45    77889705 :     std::vector<double> res;
      46    77889705 :     model.eval(res);
      47             :     
      48   155779410 :     return GaussError( obs.begin(), res, thermNoise);
      49             :     
      50    77889705 :   }
      51             : 
      52             : }
      53             :