src/growth_region.cc

Go to the documentation of this file.

// Implements the GrowthRegion class
#include "growth_region.h"
 
namespace cycamore {
 
GrowthRegion::GrowthRegion(cyclus::Context* ctx)
    : cyclus::Region(ctx) {
  #if !CYCLUS_HAS_COIN
    throw cyclus::Error("Growth Region requires that Cyclus & Cycamore be compiled "
                        "with COIN support.");
  #endif
}
 
GrowthRegion::~GrowthRegion() {}
 
void GrowthRegion::AddCommodityDemand_(std::string commod,
                                       Demand& demand) {
 
 
  cyclus::toolkit::PiecewiseFunctionFactory pff;
  cyclus::toolkit::BasicFunctionFactory bff;
  bool continuous = false;
  int time;
  std::string type, params;
  Demand::const_iterator it;
  for (it = demand.begin(); it != demand.end(); it++) {
    time = it->first;
    type = it->second.first;
    params = it->second.second;
    pff.AddFunction(bff.GetFunctionPtr(type, params), time, continuous);
    continuous = true; // only the first entry is not continuous
  }
 
  // register the commodity and demand
  cyclus::toolkit::Commodity c(commod);
  sdmanager_.RegisterCommodity(c, pff.GetFunctionPtr());
}
 
void GrowthRegion::EnterNotify() {
  cyclus::Region::EnterNotify();
  std::set<cyclus::Agent*>::iterator ait;
  for (ait = cyclus::Agent::children().begin();
       ait != cyclus::Agent::children().end();
       ++ait) {
    Agent* a = *ait;
    Register_(a);
  }
 
  std::map<std::string, Demand>::iterator it;
  for (it = commodity_demand.begin(); it != commodity_demand.end(); ++it) {
    LOG(cyclus::LEV_INFO3, "greg") << "Adding demand for commodity "
                                   << it->first;
    AddCommodityDemand_(it->first, it->second);
  }
  
  InitializePosition();
}
 
void GrowthRegion::DecomNotify(Agent* a) {
  Unregister_(a);
}
 
void GrowthRegion::Register_(cyclus::Agent* agent) {
  using cyclus::toolkit::CommodityProducerManager;
  using cyclus::toolkit::Builder;
#if CYCLUS_HAS_COIN
  CommodityProducerManager* cpm_cast =
      dynamic_cast<CommodityProducerManager*>(agent);
  if (cpm_cast != NULL) {
    LOG(cyclus::LEV_INFO3, "greg") << "Registering agent "
                                   << agent->prototype() << agent->id()
                                   << " as a commodity producer manager.";
    sdmanager_.RegisterProducerManager(cpm_cast);
  }
 
  Builder* b_cast = dynamic_cast<Builder*>(agent);
  if (b_cast != NULL) {
    LOG(cyclus::LEV_INFO3, "greg") << "Registering agent "
                                   << agent->prototype() << agent->id()
                                   << " as a builder.";
    buildmanager_.Register(b_cast);
  }
#else
  throw cyclus::Error("Growth Region requires that Cyclus & Cycamore be compiled "
                      "with COIN support.");
#endif
}
 
void GrowthRegion::Unregister_(cyclus::Agent* agent) {
  using cyclus::toolkit::CommodityProducerManager;
  using cyclus::toolkit::Builder;
#if CYCLUS_HAS_COIN
  CommodityProducerManager* cpm_cast =
    dynamic_cast<CommodityProducerManager*>(agent);
  if (cpm_cast != NULL)
    sdmanager_.UnregisterProducerManager(cpm_cast);
 
  Builder* b_cast = dynamic_cast<Builder*>(agent);
  if (b_cast != NULL)
    buildmanager_.Unregister(b_cast);
#else
  throw cyclus::Error("Growth Region requires that Cyclus & Cycamore be compiled "
                      "with COIN support.");
#endif
}
 
void GrowthRegion::Tick() {
  double demand, supply, unmetdemand;
  cyclus::toolkit::Commodity commod;
  int time = context()->time();
  std::map<std::string, Demand>::iterator it;
  for (it = commodity_demand.begin(); it != commodity_demand.end(); ++it) {
    commod = cyclus::toolkit::Commodity(it->first);
    demand = sdmanager_.Demand(commod, time);
    supply = sdmanager_.Supply(commod);
    unmetdemand = demand - supply;
 
    LOG(cyclus::LEV_INFO3, "greg") << "GrowthRegion: " << prototype()
                                   << " at time: " << time
                                   << " has the following values regarding "
                                   << " commodity: " << commod.name();
    LOG(cyclus::LEV_INFO3, "greg") << "  * demand = " << demand;
    LOG(cyclus::LEV_INFO3, "greg") << "  * supply = " << supply;
    LOG(cyclus::LEV_INFO3, "greg") << "  * unmet demand = " << unmetdemand;
 
    if (unmetdemand > 0) {
      OrderBuilds(commod, unmetdemand);
    }
  }
  cyclus::Region::Tick();
}
 
void GrowthRegion::OrderBuilds(cyclus::toolkit::Commodity& commodity,
                               double unmetdemand) {
#if CYCLUS_HAS_COIN
  using std::vector;
  vector<cyclus::toolkit::BuildOrder> orders =
    buildmanager_.MakeBuildDecision(commodity, unmetdemand);
 
  LOG(cyclus::LEV_INFO3, "greg")
      << "The build orders have been determined. "
      << orders.size()
      << " different type(s) of prototypes will be built.";
 
  cyclus::toolkit::BuildOrder* order;
  cyclus::Institution* instcast;
  cyclus::Agent* agentcast;
  for (int i = 0; i < orders.size(); i++) {
    order = &orders.at(i);
    instcast = dynamic_cast<cyclus::Institution*>(order->builder);
    agentcast = dynamic_cast<cyclus::Agent*>(order->producer);
    if (!instcast || !agentcast) {
      throw cyclus::CastError("growth_region has tried to incorrectly "
                              "cast an already known entity.");
    }
 
    LOG(cyclus::LEV_INFO3, "greg")
        << "A build order for " << order->number
        << " prototype(s) of type "
        << dynamic_cast<cyclus::Agent*>(agentcast)->prototype()
        << " from builder " << instcast->prototype()
        << " is being placed.";
 
    for (int j = 0; j < order->number; j++) {
      LOG(cyclus::LEV_DEBUG2, "greg") << "Ordering build number: " << j + 1;
      context()->SchedBuild(instcast, agentcast->prototype());
    }
  }
#else
  throw cyclus::Error("Growth Region requires that Cyclus & Cycamore be compiled "
                      "with COIN support.");
#endif
}
 
extern "C" cyclus::Agent* ConstructGrowthRegion(cyclus::Context* ctx) {
  return new GrowthRegion(ctx);
}
 
}  // namespace cycamore