+int process::run()
+{
+ bool one_more = true;
+
+ INFO1("Initial load: %g", load);
+ VERB0("Starting...");
+ iter = 0;
+ while (one_more) {
+ ++iter;
+
+ if (opt::log_rate && iter % opt::log_rate == 0) {
+ if (opt::bookkeeping)
+ INFO3("(%u) current load: %g ; expected: %g",
+ iter, load, expected_load);
+ else
+ INFO2("(%u) current load: %g",
+ iter, load);
+ }
+
+ compute();
+
+ bool close_received = !receive(false);
+
+ if (opt::bookkeeping)
+ expected_load -= load_balance(expected_load);
+ else
+ load -= load_balance(load);
+
+ send();
+ comm.flush(false);
+
+ if (opt::exit_on_close && close_received)
+ one_more = false;
+ if (opt::maxiter && iter >= opt::maxiter)
+ one_more = false;
+ }
+ VERB0("Going to finalize...");
+ finalize();
+
+ /* Open Questions :
+ * - definition of load on heterogeneous hosts ?
+ * - how to detect convergence ?
+ * - how to manage link failures ?
+ */
+
+ VERB0("Done.");
+ if (opt::bookkeeping)
+ INFO4("Final load after %d iteration%s: %g ; expected: %g",
+ iter, ESSE(iter), load, expected_load);
+ else
+ INFO3("Final load after %d iteration%s: %g", iter, ESSE(iter), load);
+ return 0;
+}
+
+void process::compute()
+{
+ // fixme: shall we do something special when duration is 0 ?
+ double duration = opt::comp_cost(load);
+ m_task_t task = MSG_task_create("computation", duration, 0.0, NULL);
+ DEBUG2("compute %g flop%s.", duration, ESSE(duration));
+ MSG_task_execute(task);
+ MSG_task_destroy(task);
+}
+
+double process::load_balance(double /*my_load*/)
+{
+ return 0.0;
+}
+
+void process::send1_no_bookkeeping(neighbor& nb)
+{
+ comm.send(nb.get_ctrl_mbox(), new message(message::INFO, load));
+ double load_to_send = nb.get_to_send();
+ if (load_to_send > 0.0) {
+ comm.send(nb.get_data_mbox(), new message(message::LOAD, load_to_send));
+ nb.set_to_send(0.0);
+ }
+}
+
+void process::send1_bookkeeping(neighbor& nb)
+{
+ comm.send(nb.get_ctrl_mbox(), new message(message::INFO, expected_load));
+ double load_to_send;
+ double new_debt;
+ double debt_to_send = nb.get_to_send();
+ if (debt_to_send > 0.0) {
+ comm.send(nb.get_ctrl_mbox(),
+ new message(message::CREDIT, debt_to_send));
+ nb.set_to_send(0.0);
+ new_debt = nb.get_debt() + debt_to_send;
+ } else {
+ new_debt = nb.get_debt();
+ }
+ if (load <= new_debt) {
+ load_to_send = load;
+ nb.set_debt(new_debt - load_to_send);
+ load = 0.0;
+ } else {
+ load_to_send = new_debt;
+ nb.set_debt(0.0);
+ load -= load_to_send;
+ }
+ if (load_to_send > 0.0)
+ comm.send(nb.get_data_mbox(), new message(message::LOAD, load_to_send));
+}
+
+void process::send()
+{
+ // fixme: shall we send data at all iterations?
+ if (opt::bookkeeping) {
+ std::for_each(neigh.begin(), neigh.end(),
+ std::tr1::bind(&process::send1_bookkeeping,
+ this, std::tr1::placeholders::_1));
+ } else {
+ std::for_each(neigh.begin(), neigh.end(),
+ std::tr1::bind(&process::send1_no_bookkeeping,
+ this, std::tr1::placeholders::_1));
+ }
+}
+
+// Returns false if a CLOSE message was received.
+bool process::receive(bool wait_for_close)
+{
+ bool result = true;
+ message* msg;
+ m_host_t from;
+ while ((ctrl_close_pending ||
+ data_close_pending) && comm.recv(msg, from, wait_for_close)) {
+ DEBUG2("received %s from %s",
+ msg->to_string().c_str(), MSG_host_get_name(from));
+ switch (msg->get_type()) {
+ case message::INFO: {
+ neighbor* n = rev_neigh[from];
+ n->set_load(msg->get_amount());
+ break;
+ }
+ case message::CREDIT:
+ expected_load += msg->get_amount();
+ break;
+ case message::LOAD:
+ load += msg->get_amount();
+ break;
+ case message::CTRL_CLOSE:
+ if (--ctrl_close_pending == 1)
+ comm.next_close_on_ctrl_is_last();
+ DEBUG1("ctrl_close_pending = %d", ctrl_close_pending);
+ result = false;
+ break;
+ case message::DATA_CLOSE:
+ if (--data_close_pending == 1)
+ comm.next_close_on_data_is_last();
+ DEBUG1("data_close_pending = %d", data_close_pending);
+ result = false;
+ break;
+ }
+ delete msg;
+ }
+ return result;
+}
+
+void process::finalize1(neighbor& nb)
+{
+ comm.send(nb.get_ctrl_mbox(), new message(message::CTRL_CLOSE, 0.0));
+ comm.send(nb.get_data_mbox(), new message(message::DATA_CLOSE, 0.0));
+}
+
+void process::finalize()
+{
+ DEBUG2("send CLOSE to %d neighbor%s.",
+ (int )neigh.size(), ESSE(neigh.size()));
+ std::for_each(neigh.begin(), neigh.end(),
+ std::tr1::bind(&process::finalize1,
+ this, std::tr1::placeholders::_1));
+
+ DEBUG2("wait for CLOSE from %d neighbor%s.",
+ (int )neigh.size(), ESSE(neigh.size()));
+ receive(true);
+
+ comm.flush(true);
+}
+