include/geometry/simulator.hpp
author Dmitriy Morozov <dmitriy@mrzv.org>
Thu, 07 Sep 2017 17:35:28 -0700
changeset 290 21b56c6c9512
parent 272 29306411272b
permissions -rw-r--r--
Add Vanessa Robins' periodic alpha shapes code

#include "utilities/log.h"
#include "utilities/counter.h"

#ifdef LOGGING
static rlog::RLogChannel* rlSimulator =             DEF_CHANNEL("geometry/simulator", rlog::Log_Debug);

#endif // LOGGING

#ifdef COUNTERS
static Counter*  cSimulatorProcess =                GetCounter("simulator/process");
#endif // COUNTERS


template<class FuncKernel_, template<class Event> class EventComparison_>
template<class Event_>
typename Simulator<FuncKernel_, EventComparison_>::Key
Simulator<FuncKernel_, EventComparison_>::
add(const Event_& e)
{
    Event* ee = new Event_(e);
    return queue_.push(ee);
}

template<class FuncKernel_, template<class Event> class EventComparison_>
template<class Event_>
typename Simulator<FuncKernel_, EventComparison_>::Key
Simulator<FuncKernel_, EventComparison_>::
add(const Function& f, const Event_& e)
{
    Event* ee = new Event_(e);
    rLog(rlSimulator, "Solving: %s", tostring(f).c_str());
    int sign = FunctionKernel::sign_at_negative_infinity(f);        // going to be sign after current time
    rLog(rlSimulator, "Sign at -infinity: %i", sign);
    if (sign != 0)
    {
        FunctionKernel::solve(f, ee->root_stack());
        rLog(rlSimulator, "Got solution with root stack size: %i", ee->root_stack().size());
    }

    while (!ee->root_stack().empty() && ee->root_stack().top() < current_time())
    {
        // rLog(rlSimulator, "Popping expired root: %f", ee->root_stack().top());
        ee->root_stack().pop();
        sign *= -1;
    }

    if (sign == -1)
    {
        rLog(rlSimulator, "Popping the root because of negative sign (degeneracy)");
        // rLog(rlSimulator, "Popping the root because of negative sign (degeneracy): %f", ee->root_stack().top());
        // rLog(rlSimulator, "  Current time: %f", current_time());
        // AssertMsg(ee->root_stack().top() == current_time(),
                 // "If sign is negative, we must be in the degenerate case");
        ee->root_stack().pop();
    }

    if (ee->root_stack().empty())
        rLog(rlSimulator, "Pushing event with empty root stack");
    else
    {
        rLog(rlSimulator, "Root stack size: %i", ee->root_stack().size());
        rLog(rlSimulator, "Pushing: %s", tostring(ee->root_stack().top()).c_str());
    }
    Key k = queue_.push(ee);
    return k;
}

//template<class FuncKernel_, template<class Event> class EventComparison_>
//void
//Simulator<FuncKernel_, EventComparison_>::
//update(Key k, const Function& f)
//{
//    Event* ee = *k;
//    Event* e = new Event;
//    e->root_stack() = RootStack();                              // no clear() in std::stack
//    FunctionKernel::solve(f, e->root_stack());
//    while (!e->root_stack().empty() && e->root_stack().top() < current_time())
//        e->root_stack().pop();
//    queue_.replace(k, e);
//    delete ee;
//}

template<class FuncKernel_, template<class Event> class EventComparison_>
void
Simulator<FuncKernel_, EventComparison_>::
process()
{
    Count(cSimulatorProcess);
    if (reached_infinity()) return;
    rLog(rlSimulator, "Queue size: %i", queue_.size());
    Key top = queue_.top();
    Event* e = *top;
    rLog(rlSimulator, "Processing event: %s", intostring(*e).c_str());

    current_ = e->root_stack().top(); e->root_stack().pop();
    queue_.demoted(top);

    // Get the top element out of the queue, put it back depending on what process() says
    if (!(e->process(this)))            { queue_.remove(top);  delete e; }

    ++count_;
}

template<class FuncKernel_, template<class Event> class EventComparison_>
typename Simulator<FuncKernel_, EventComparison_>::Time
Simulator<FuncKernel_, EventComparison_>::
audit_time() const
{
    const_Key top = queue_.top();
    Event* e = *top;

    if (e->root_stack().empty()) return current_ + 1;
    else return FunctionKernel::between(e->root_stack().top(), current_);
}

template<class FuncKernel_, template<class Event> class EventComparison_>
std::ostream&
Simulator<FuncKernel_, EventComparison_>::
operator<<(std::ostream& out) const
{
    out << "Simulator: " << std::endl;
    return queue_.print(out, "  ");
}

template<class FuncKernel_, template<class Event> class EventComparison_>
std::ostream&
operator<<(std::ostream& out, const Simulator<FuncKernel_, EventComparison_>& s)
{
    return s.operator<<(out);
}

template<class FuncKernel_, template<class Event> class EventComparison_>
std::ostream&
operator<<(std::ostream& out, const typename Simulator<FuncKernel_, EventComparison_>::Event& e)
{
    return e.operator<<(out);
}