initial_guess_value.h

Go to the documentation of this file.

#ifndef INITIAL_GUESS_VALUE_H
#define INITIAL_GUESS_VALUE_H
#include "composite_initial_guess.h"

namespace FullPhysics {
/****************************************************************/
class InitialGuessValue : public InitialGuessBuilder {
public:
  virtual ~InitialGuessValue() {}
  virtual int number_element() const { return apriori_.rows(); }
  virtual void build_initial_value(blitz::Array<double, 1>& v, int index) 
    const;
  virtual void build_apriori(blitz::Array<double, 1>& v, int index) const;
  virtual void build_apriori_covariance(blitz::Array<double, 2>& m, 
                                        int index) const;

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

  const blitz::Array<double, 1>& apriori() const {return apriori_;}

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

  void apriori(const blitz::Array<double, 1>& v) {apriori_.reference(v.copy());}

  void apriori_subset(const blitz::Array<bool, 1>& Flag, 
                      const blitz::Array<double, 1>& V);
  void initial_guess_subset(const blitz::Array<bool, 1>& Flag, 
                            const blitz::Array<double, 1>& V);
  void apriori_covariance_subset(const blitz::Array<bool, 1>& Flag, 
                                 const blitz::Array<double, 2>& V);

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

  const blitz::Array<double, 2>& apriori_covariance() const 
  {return apriori_covariance_;}

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

  void apriori_covariance(const blitz::Array<double, 2>& m) 
  {apriori_covariance_.reference(m.copy());}

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

  const blitz::Array<double, 1>& initial_guess() const 
  {
    if(initial_guess_.rows() > 0)
      return initial_guess_;
    else
      return apriori_;
  }

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

  void initial_guess(const blitz::Array<double, 1>& v) 
  {initial_guess_.reference(v.copy());}

  virtual void print(std::ostream& Os) const;

private:
  blitz::Array<double, 1> apriori_, initial_guess_;
  blitz::Array<double, 2> apriori_covariance_;
};
}
#endif