prey.h

Go to the documentation of this file.

#ifndef prey_h
#define prey_h
 
#include "popinfovector.h"
#include "conversionindex.h"
#include "commentstream.h"
#include "charptrvector.h"
#include "hasname.h"
#include "livesonareas.h"
#include "doublematrix.h"
#include "areatime.h"
#include "agebandmatrix.h"
#include "keeper.h"
#include "gadget.h"
#include "modelvariable.h"
 
enum PreyType { STOCKPREY = 1, LENGTHPREY };
 
class Prey : public HasName, public LivesOnAreas {
public:
  Prey(CommentStream& infile, const IntVector& areas, const char* givenname,
  const TimeClass* const TimeInfo, Keeper* const keeper);
  Prey(CommentStream& infile, const char* givenname, const IntVector& areas,
  const TimeClass* const TimeInfo, Keeper* const keeper);
  virtual ~Prey() = 0;
  virtual void Sum(const AgeBandMatrix& Alkeys, int area) {};
  virtual void Sum(const PopInfoVector& NumberInArea, int area) {};
  virtual void Reset(const TimeClass* const TimeInfo);
  virtual void Print(ofstream& outfile) const;
  virtual void setCI(const LengthGroupDivision* const GivenLDiv);
  virtual void Subtract(AgeBandMatrix& Alkeys, int area);
  void addBiomassConsumption(int area, const DoubleVector& predcons);
  void addNumbersConsumption(int area, const DoubleVector& predcons);
  virtual void checkConsumption(int area, const TimeClass* const TimeInfo);
  double getBiomass(int area, int length) const { return biomass[this->areaNum(area)][length]; };
  double getNumber(int area, int length) const { return preynumber[this->areaNum(area)][length].N; };
  double getTotalBiomass(int area) const { return total[this->areaNum(area)]; };
  int isOverConsumption(int area);
  const DoubleVector& getRatio(int area) const { return ratio[this->areaNum(area)]; };
  const DoubleVector& getUseRatio(int area) const { return useratio[this->areaNum(area)]; };
  const DoubleVector& getConsumption(int area) const { return consumption[this->areaNum(area)]; };
  const DoubleVector& getOverConsumption(int area) const { return overconsumption[this->areaNum(area)]; };
  double getTotalOverConsumption(int area) const;
  const LengthGroupDivision* getLengthGroupDiv() const { return LgrpDiv; };
  const PopInfoVector& getConsumptionPopInfo(int area) const { return preynumber[this->areaNum(area)]; };
  double getEnergy() const { return energy; };
  int isPreyArea(int area);
  PreyType getType() const { return type; };
protected:
  ConversionIndex* CI;
  LengthGroupDivision* LgrpDiv;
  PopInfoMatrix preynumber;
  ModelVariable energy;
  DoubleMatrix biomass;
  DoubleVector total;
  DoubleMatrix ratio;
  DoubleMatrix useratio;
  DoubleMatrix consumption;
  DoubleMatrix cons;
  DoubleMatrix overconsumption;
  IntVector isoverconsumption;
  DoubleMatrix consratio;
  PreyType type;
};
 
#endif