ecosystem.h

Go to the documentation of this file.

#ifndef ecosystem_h
#define ecosystem_h
 
#include "areatime.h"
#include "keeper.h"
#include "errorhandler.h"
#include "tagptrvector.h"
#include "baseclassptrvector.h"
#include "stockptrvector.h"
#include "fleetptrvector.h"
#include "otherfoodptrvector.h"
#include "printerptrvector.h"
#include "likelihoodptrvector.h"
#include "optinfoptrvector.h"
#include "stock.h"
#include "fleet.h"
#include "otherfood.h"
#include "printer.h"
#include "maininfo.h"
#include "printinfo.h"
#include "optinfo.h"
#include "gadget.h"
 
class Ecosystem {
public:
  Ecosystem() {};
  Ecosystem(const MainInfo& main);
  ~Ecosystem();
  void readMain(CommentStream& infile, const MainInfo& main);
  void readLikelihood(CommentStream& infile);
  void readPrinters(CommentStream& infile);
  void readFleet(CommentStream& infile);
  void readTagging(CommentStream& infile);
  void readOtherFood(CommentStream& infile);
  void readStock(CommentStream& infile);
  void readOptimisation(CommentStream& infile);
  void writeStatus(const char* filename) const;
  void writeInitialInformation(const char* const filename);
  void writeValues();
  void writeParams(const char* const filename, int prec) const;
  void writeBestValues() { keeper->writeBestValues(); };
  void writeOptValues();
  void Optimise();
  void Initialise();
  void Update(const StochasticData* const Stochastic) const {keeper->Update(Stochastic); };
  void Reset();
  void storeVariables(double likvalue, const DoubleVector& point) { keeper->storeVariables(likvalue, point); };
  void resetVariables() { keeper->resetVariables(); };
  void scaleVariables() { keeper->scaleVariables(); };
  void getOptInitialValues(DoubleVector& val) const { keeper->getOptInitialValues(val); };
  void getOptScaledValues(DoubleVector& val) const { keeper->getOptScaledValues(val); };
  void getOptLowerBounds(DoubleVector& lbs) const { keeper->getOptLowerBounds(lbs); };
  void getOptUpperBounds(DoubleVector& ubs) const { keeper->getOptUpperBounds(ubs); };
  void checkBounds() const { keeper->checkBounds(likevec); };
  int numVariables() const { return keeper->numVariables(); };
  int numOptVariables() const { return keeper->numOptVariables(); };
  void Simulate(int print);
  void updatePredationOneArea(int area = 0);
  void updatePopulationOneArea(int area = 0);
  void updateAgesOneArea(int area = 0);
  double SimulateAndUpdate(const DoubleVector& x);
  double getLikelihood() const { return likelihood; };
  int getFuncEval() const { return funceval; };
  int getCurrentTime() const { return TimeInfo->getTime(); };
  int getCurrentYear() const { return TimeInfo->getYear(); };
  int getCurrentStep() const { return TimeInfo->getStep(); };
  int numTotalSteps() const { return TimeInfo->numTotalSteps(); };
  LikelihoodPtrVector& getModelLikelihoodVector() { return likevec; };
  PrinterPtrVector& getModelPrinterVector() { return printvec; };
  FleetPtrVector& getModelFleetVector() { return fleetvec; };
  StockPtrVector& getModelStockVector() { return stockvec; };
  volatile int interrupted;
protected:
  double likelihood;
  int funceval;
  BaseClassPtrVector basevec;
  LikelihoodPtrVector likevec;
  PrinterPtrVector printvec;
  OptInfoPtrVector optvec;
  TimeClass* TimeInfo;
  AreaClass* Area;
  Keeper* keeper;
  StockPtrVector stockvec;
  TagPtrVector tagvec;
  OtherFoodPtrVector otherfoodvec;
  FleetPtrVector fleetvec;
  PrintInfo printinfo;
  int printcount;
  DoubleVector initialval;
  DoubleVector currentval;
  IntVector optflag;
};
 
#endif