DECLARE_EXPORT void SolverMRP::solve(const Flow* fl, void* v) // @todo implement search mode
{
// Note: This method is only called for consuming flows and for the leading
// flow of an alternate group. See SolverMRP::checkOperation
SolverMRPdata* data = static_cast<SolverMRPdata*>(v);
if (fl->hasAlternates())
{
// CASE I: It is an alternate flow.
// We ask each alternate flow in order of priority till we find a flow
// that has a non-zero reply.
// 1) collect a list of alternates
list<const Flow*> thealternates;
const Flow *x = fl->hasAlternates() ? fl : fl->getAlternate();
for (Operation::flowlist::const_iterator i = fl->getOperation()->getFlows().begin();
i != fl->getOperation()->getFlows().end(); ++i)
if ((i->getAlternate() == x || &*i == x)
&& i->getEffective().within(data->state->q_flowplan->getDate()))
thealternates.push_front(&*i);
// 2) Sort the list
thealternates.sort(sortFlow);
// 3) Control the planning mode
bool originalPlanningMode = data->constrainedPlanning;
data->constrainedPlanning = true;
const Flow *firstAlternate = NULL;
double firstQuantity = 0.0;
// Remember the top constraint
bool originalLogConstraints = data->logConstraints;
//Problem* topConstraint = data->planningDemand->getConstraints().top();
// 4) Loop through the alternates till we find a non-zero reply
Date min_next_date(Date::infiniteFuture);
double ask_qty;
FlowPlan *flplan = data->state->q_flowplan;
for (list<const Flow*>::const_iterator i = thealternates.begin();
i != thealternates.end();)
{
const Flow *curflow = *i;
data->state->q_flowplan = flplan; // because q_flowplan can change
// 4a) Switch to this flow
if (data->state->q_flowplan->getFlow() != curflow)
data->state->q_flowplan->setFlow(curflow);
// 4b) Call the Python user exit if there is one
if (userexit_flow)
{
PythonObject result = userexit_flow.call(data->state->q_flowplan, PythonObject(data->constrainedPlanning));
if (!result.getBool())
{
// Return value is false, alternate rejected
if (data->getSolver()->getLogLevel()>1)
logger << indent(curflow->getOperation()->getLevel())
<< " User exit disallows consumption from '"
<< (*i)->getBuffer()->getName() << "'" << endl;
// Move to the next alternate
if (++i != thealternates.end() && data->getSolver()->getLogLevel()>1)
logger << indent(curflow->getOperation()->getLevel()) << " Alternate flow switches from '"
<< curflow->getBuffer()->getName() << "' to '"
<< (*i)->getBuffer()->getName() << "'" << endl;
continue;
}
}
// Remember the first alternate
if (!firstAlternate)
{
firstAlternate = *i;
firstQuantity = data->state->q_flowplan->getQuantity();
}
// Constraint tracking
if (*i != firstAlternate)
// Only enabled on first alternate
data->logConstraints = false;
else
// Keep track of constraints, if enabled
data->logConstraints = originalLogConstraints;
// 4c) Ask the buffer
data->state->q_qty = ask_qty = - data->state->q_flowplan->getQuantity();
data->state->q_date = data->state->q_flowplan->getDate();
CommandManager::Bookmark* topcommand = data->setBookmark();
curflow->getBuffer()->solve(*this,data);
// 4d) A positive reply: exit the loop
if (data->state->a_qty > ROUNDING_ERROR)
{
// Update the opplan, which is required to (1) update the flowplans
// and to (2) take care of lot sizing constraints of this operation.
if (data->state->a_qty < ask_qty - ROUNDING_ERROR)
{
flplan->setQuantity(-data->state->a_qty, true);
data->state->a_qty = -flplan->getQuantity();
}
if (data->state->a_qty > ROUNDING_ERROR)
{
data->constrainedPlanning = originalPlanningMode;
data->logConstraints = originalLogConstraints;
return;
}
}
// 4e) Undo the plan on the alternate
data->rollback(topcommand);
// 4f) Prepare for the next alternate
if (data->state->a_date < min_next_date)
min_next_date = data->state->a_date;
if (++i != thealternates.end() && data->getSolver()->getLogLevel()>1)
logger << indent(curflow->getOperation()->getLevel()) << " Alternate flow switches from '"
<< curflow->getBuffer()->getName() << "' to '"
<< (*i)->getBuffer()->getName() << "'" << endl;
}
// 5) No reply found, all alternates are infeasible
if (!originalPlanningMode)
{
assert(firstAlternate);
// Unconstrained plan: Plan on the primary alternate
// Switch to this flow
if (flplan->getFlow() != firstAlternate)
flplan->setFlow(firstAlternate);
// Message
if (data->getSolver()->getLogLevel()>1)
logger << indent(fl->getOperation()->getLevel())
<< " Alternate flow plans unconstrained on alternate '"
<< firstAlternate->getBuffer()->getName() << "'" << endl;
// Plan unconstrained
data->constrainedPlanning = false;
data->state->q_flowplan = flplan; // because q_flowplan can change
flplan->setQuantity(firstQuantity, true);
data->state->q_qty = ask_qty = - flplan->getQuantity();
data->state->q_date = flplan->getDate();
firstAlternate->getBuffer()->solve(*this,data);
data->state->a_qty = -flplan->getQuantity();
// Restore original planning mode
data->constrainedPlanning = originalPlanningMode;
}
else
{
// Constrained plan: Return 0
data->state->a_date = min_next_date;
data->state->a_qty = 0;
if (data->getSolver()->getLogLevel()>1)
logger << indent(fl->getOperation()->getLevel()) <<
" Alternate flow doesn't find supply on any alternate : "
<< data->state->a_qty << " " << data->state->a_date << endl;
}
}
else
{
// CASE II: Not an alternate flow.
// In this case, this method is passing control on to the buffer.
data->state->q_qty = - data->state->q_flowplan->getQuantity();
data->state->q_date = data->state->q_flowplan->getDate();
if (data->state->q_qty != 0.0)
{
fl->getBuffer()->solve(*this,data);
if (data->state->a_date > fl->getEffective().getEnd())
{
// The reply date must be less than the effectivity end date: after
// that date the flow in question won't consume any material any more.
if (data->getSolver()->getLogLevel()>1
&& data->state->a_qty < ROUNDING_ERROR)
logger << indent(fl->getBuffer()->getLevel()) << " Buffer '"
<< fl->getBuffer()->getName() << "' answer date is adjusted to "
<< fl->getEffective().getEnd()
<< " because of a date effective flow" << endl;
data->state->a_date = fl->getEffective().getEnd();
}
}
else
{
// It's a zero quantity flowplan.
// E.g. because it is not effective.
data->state->a_date = data->state->q_date;
data->state->a_qty = 0.0;
}
}
}
void OperationItemSupplier::trimExcess(bool zero_or_minimum) const
{
// This method can only trim operations not loading a resource
if (getLoads().begin() != getLoads().end())
return;
for (Operation::flowlist::const_iterator fliter = getFlows().begin();
fliter != getFlows().end(); ++fliter)
{
if (fliter->getQuantity() <= 0)
// Strange, shouldn't really happen
continue;
FlowPlan* candidate = nullptr;
double curmin = 0;
double oh = 0;
double excess_min = DBL_MAX;
for (Buffer::flowplanlist::const_iterator flplniter = fliter->getBuffer()->getFlowPlans().begin();
flplniter != fliter->getBuffer()->getFlowPlans().end();
++flplniter)
{
// For any operationplan we get the onhand when its successor
// replenishment arrives. If that onhand is higher than the minimum
// onhand value we can resize it.
// This is only valid in unconstrained plans and when there are
// no upstream activities.
if (flplniter->getEventType() == 3 && zero_or_minimum)
curmin = flplniter->getMin();
else if (flplniter->getEventType() == 1)
{
const FlowPlan* flpln = static_cast<const FlowPlan*>(&*flplniter);
if (oh - curmin < excess_min)
{
excess_min = oh - curmin;
if (excess_min < 0)
excess_min = 0;
}
if (flpln->getQuantity() > 0 && !flpln->getOperationPlan()->getLocked() && (!candidate || candidate->getDate() != flpln->getDate()))
{
if (candidate
&& excess_min > ROUNDING_ERROR
&& candidate->getQuantity() > excess_min + ROUNDING_ERROR
&& candidate->getQuantity() > getSizeMinimum() + ROUNDING_ERROR
)
{
// This candidate can now be resized
candidate->setQuantity(candidate->getQuantity() - excess_min, false);
candidate = nullptr;
}
else if (flpln->getOperation() == this)
candidate = const_cast<FlowPlan*>(flpln);
else
candidate = nullptr;
excess_min = DBL_MAX;
}
}
oh = flplniter->getOnhand();
}
if (candidate
&& excess_min > ROUNDING_ERROR
&& candidate->getQuantity() > excess_min + ROUNDING_ERROR
&& candidate->getQuantity() > getSizeMinimum() + ROUNDING_ERROR
)
// Resize the last candidate at the end of the horizon
candidate->setQuantity(candidate->getQuantity() - excess_min, false);
}
}
void OperatorDelete::solve(const Buffer* b, void* v)
{
if (getLogLevel()>1)
logger << "Scanning " << b << " for excess" << endl;
Buffer::flowplanlist::const_iterator fiter = b->getFlowPlans().rbegin();
Buffer::flowplanlist::const_iterator fend = b->getFlowPlans().end();
if (fiter == fend)
return; // There isn't a single flowplan in the buffer
double excess = fiter->getOnhand() - fiter->getMin();
// Find the earliest occurence of the excess
fiter = b->getFlowPlans().begin();
while (excess > ROUNDING_ERROR && fiter != fend)
{
if (fiter->getQuantity() <= 0)
{
// Not a producer
++fiter;
continue;
}
FlowPlan* fp = nullptr;
if (fiter->getEventType() == 1)
fp = const_cast<FlowPlan*>(static_cast<const FlowPlan*>(&*fiter));
double cur_excess = b->getFlowPlans().getExcess(&*fiter);
if (!fp || fp->getOperationPlan()->getLocked() || cur_excess < ROUNDING_ERROR)
{
// No excess producer, or it's locked
++fiter;
continue;
}
assert(fp);
// Increment the iterator here, because it can get invalidated later on
while (
fiter != fend
&& fiter->getEventType() == 1
&& static_cast<const FlowPlan*>(&*fiter)->getOperationPlan()->getTopOwner()==fp->getOperationPlan()->getTopOwner()
)
++fiter;
if (cur_excess >= fp->getQuantity() - ROUNDING_ERROR)
{
// The complete operationplan is excess.
// Reduce the excess
excess -= fp->getQuantity();
// Add upstream buffers to the stack
pushBuffers(fp->getOperationPlan(), true);
// Log message
if (getLogLevel()>0)
logger << "Removing excess operationplan: '"
<< fp->getOperationPlan()->getOperation()
<< "' " << fp->getOperationPlan()->getDates()
<< " " << fp->getOperationPlan()->getQuantity()
<< endl;
// Delete operationplan
if (cmds)
cmds->add(new CommandDeleteOperationPlan(fp->getOperationPlan()));
else
delete fp->getOperationPlan();
}
else
{
// Reduce the operationplan
double newsize = fp->setQuantity(fp->getQuantity() - cur_excess, false, false);
if (newsize == fp->getQuantity())
// No resizing is feasible
continue;
// Add upstream buffers to the stack
pushBuffers(fp->getOperationPlan(), true);
// Reduce the excess
excess -= fp->getQuantity() - newsize;
if (getLogLevel()>0)
logger << "Resizing excess operationplan to " << newsize << ": '"
<< fp->getOperationPlan()->getOperation()
<< "' " << fp->getOperationPlan()->getDates()
<< " " << fp->getOperationPlan()->getQuantity()
<< endl;
// Resize operationplan
if (cmds)
cmds->add(new CommandMoveOperationPlan(
fp->getOperationPlan(), Date::infinitePast,
fp->getOperationPlan()->getDates().getEnd(), newsize)
);
else
fp->getOperationPlan()->setQuantity(newsize);
}
}
}