stokes_coefficient_imp_base.h

Go to the documentation of this file.

#ifndef STOKES_COEFFICIENT_IMP_BASE_H
#define STOKES_COEFFICIENT_IMP_BASE_H
#include "stokes_coefficient.h"
#include "sub_state_vector_array.h"

namespace FullPhysics {
/****************************************************************/
class StokesCoefficientImpBase: public SubStateVectorArray<StokesCoefficient> {
public:
  virtual ~StokesCoefficientImpBase() {}
  virtual ArrayAd<double, 2> stokes_coefficient() const
  { fill_cache(); return stokes_coeff; }
  virtual boost::shared_ptr<StokesCoefficient> clone() const = 0;
  virtual void update_sub_state_hook() 
  { cache_stale = true; }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  virtual void print(std::ostream& Os) const { Os << desc(); }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  virtual std::string desc() const { return "StokesCoefficientImpBase"; }
protected:
//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  mutable bool cache_stale;

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  mutable ArrayAd<double, 2> stokes_coeff;

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  virtual void calc_stokes_coeff() const = 0;

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

  StokesCoefficientImpBase() : cache_stale(true) { }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  StokesCoefficientImpBase(const blitz::Array<double, 1>& Coeff, 
                           const blitz::Array<bool, 1>& Used_flag)
    : SubStateVectorArray<StokesCoefficient>(Coeff, Used_flag),
      cache_stale(true) {}
private:
  void fill_cache() const
  {
    if(cache_stale) {
      stokes_coeff.resize_number_variable(coeff.number_variable());
      calc_stokes_coeff();
    }
    cache_stale = false;
  }
};

}
#endif