vla_profile_inversion_target.cpp

Go to the documentation of this file.

 
// C/C++ headers
#include <iomanip>
#include <iostream>
 
// climath
#include <climath/interpolation.h>
#include <climath/linalg.h>
 
// athena
#include <athena/athena.hpp>
#include <athena/mesh/mesh.hpp>
 
// harp
#include <harp/radiation.hpp>
#include <harp/radiation_band.hpp>
 
// application
#include <application/application.hpp>
#include <application/exceptions.hpp>
 
// inversion
#include "profile_inversion.hpp"
 
void VLAProfileInversion::CalculateFitTarget(Radiation const *prad, Real *val,
                                             int nvalue, int k, int j) const {
  std::stringstream msg;
  Application::Logger app("inversion");
 
  app->Log("VLAProfileInversion::CalculateFitTarget");
  app->Log("model = " + std::to_string(j));
 
  if (nvalue != prad->GetNumBands()) {
    throw ValueError("CalculateFitTarget", "nvalue", prad->GetNumBands(),
                     nvalue);
  }
 
  // 11. log likelihood
  std::vector<Real> mus, tbs;
 
  for (int b = 0; b < prad->GetNumBands(); ++b) {
    auto pband = prad->GetBand(b);
 
    // emission angles;
    int ndir = pband->GetNumOutgoingRays();
    mus.resize(ndir);
    tbs.resize(ndir);
 
    for (int n = 0; n < ndir; ++n) mus[n] = pband->GetCosinePolarAngle(n);
 
    // brightness temperature
    val[b] = prad->radiance(b, k, j);
 
    if (fit_differential_) {
      // brightness temperature differential
      val[b] -= pband->btoa(0, k, pmy_block_->js - 1);
    }
  }
 
  // app->Log("foward model results", val, nvalue);
}