pressure_imp_base.h

Go to the documentation of this file.

#ifndef PRESSURE_IMP_BASE_H
#define PRESSURE_IMP_BASE_H
#include "pressure.h"
#include "sub_state_vector_array.h"

namespace FullPhysics {
/****************************************************************/
class PressureImpBase: public SubStateVectorArray<Pressure> {
public:
  virtual ~PressureImpBase() {}
  virtual ArrayAdWithUnit<double, 1> pressure_grid() const
  { fill_cache(); return pgrid; }
  virtual boost::shared_ptr<Pressure> 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 "PressureImpBase"; }
protected:
//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  mutable bool cache_stale;

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  mutable ArrayAdWithUnit<double, 1> pgrid;

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

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

  PressureImpBase() : cache_stale(true) { }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  PressureImpBase(const blitz::Array<double, 1>& Coeff, 
                  const blitz::Array<bool, 1>& Used_flag)
    : SubStateVectorArray<Pressure>(Coeff, Used_flag),
      cache_stale(true) {}
private:
  void fill_cache() const
  {
    if(cache_stale) {
      pgrid.value.resize_number_variable(coeff.number_variable());
      calc_pressure_grid();
    }
    cache_stale = false;
  }
  
};

}
#endif