nlls_max_a_posteriori.h

Go to the documentation of this file.

#ifndef NLLS_MAX_A_POSTERIORI_H
#define NLLS_MAX_A_POSTERIORI_H
#include <max_a_posteriori.h>
#include <nlls_problem.h>
#include <nlls_problem_state.h>
#include <boost/shared_ptr.hpp>

namespace FullPhysics {
/******************************************************************
  This class casts a maximum a posteriori problem into the
  Nonlinear Least Squares problem.
*******************************************************************/
class NLLSMaxAPosteriori: 
    public NLLSProblem, public NLLSProblemState {

public:

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  NLLSMaxAPosteriori(const boost::shared_ptr<MaxAPosteriori>& map)
    : NLLSProblem(), NLLSProblemState(), MAP(map)
  {}

  virtual ~NLLSMaxAPosteriori() {}


//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  virtual blitz::Array<double, 1> residual();


//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  virtual blitz::Array<double, 2> jacobian();


//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  virtual int residual_size() const
  { return MAP->measurement_size()+MAP->expected_parameter_size(); }


//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  virtual int expected_parameter_size() const
  { return MAP->expected_parameter_size(); }


//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  virtual void parameters(const blitz::Array<double, 1>& x);


//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  virtual blitz::Array<double, 1> parameters() const
  { return NLLSProblemState::parameters(); }



  boost::shared_ptr<MaxAPosteriori> max_a_posteriori()
  { return MAP; }


//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  virtual void print(std::ostream& Os) const 
  { Os << "NLLSMaxAPosteriori"; }


protected:

  const boost::shared_ptr<MaxAPosteriori> MAP;

};
}
#endif