problem-counters.cpp

Go to the documentation of this file.

#include <alpaqa/problem/problem-counters.hpp>
 
#include <iomanip>
#include <iostream>
 
namespace alpaqa {
 
namespace {
struct CountResult {
    unsigned count;
    std::chrono::nanoseconds time;
};
std::ostream &operator<<(std::ostream &os, const CountResult &t) {
    auto sec = [](auto t) { return std::chrono::duration<double>(t).count(); };
    os << std::setw(8);
    if (t.count > 0) {
        os << t.count << "  (";
        auto old  = os.flags();
        auto prec = os.precision(3);
        os << std::scientific << std::setw(9) << 1e6 * sec(t.time) << " µs, "
           << std::setw(9) << 1e6 * sec(t.time) / static_cast<double>(t.count)
           << " µs/call)\r\n";
        os.precision(prec);
        os.flags(old);
    } else {
        os << '-' << "\r\n";
    }
    return os;
}
} // namespace
 
std::ostream &operator<<(std::ostream &os, const EvalCounter &c) {
    os << "        proj_diff_g:" //
       << CountResult{c.proj_diff_g, c.time.proj_diff_g};
    os << "   proj_multipliers:" //
       << CountResult{c.proj_multipliers, c.time.proj_multipliers};
    os << "     prox_grad_step:" //
       << CountResult{c.prox_grad_step, c.time.prox_grad_step};
    os << "                  f:" //
       << CountResult{c.f, c.time.f};
    os << "             grad_f:" //
       << CountResult{c.grad_f, c.time.grad_f};
    os << "           f_grad_f:" //
       << CountResult{c.f_grad_f, c.time.f_grad_f};
    os << "                f_g:" //
       << CountResult{c.f_g, c.time.f_g};
    os << " grad_f_grad_g_prod:" //
       << CountResult{c.grad_f_grad_g_prod, c.time.grad_f_grad_g_prod};
    os << "                  g:" //
       << CountResult{c.g, c.time.g};
    os << "        grad_g_prod:" //
       << CountResult{c.grad_g_prod, c.time.grad_g_prod};
    os << "            grad_gi:" //
       << CountResult{c.grad_gi, c.time.grad_gi};
    os << "              jac_g:" //
       << CountResult{c.jac_g, c.time.jac_g};
    os << "             grad_L:" //
       << CountResult{c.grad_L, c.time.grad_L};
    os << "        hess_L_prod:" //
       << CountResult{c.hess_L_prod, c.time.hess_L_prod};
    os << "             hess_L:" //
       << CountResult{c.hess_L, c.time.hess_L};
    os << "        hess_ψ_prod:" //
       << CountResult{c.hess_ψ_prod, c.time.hess_ψ_prod};
    os << "             hess_ψ:" //
       << CountResult{c.hess_ψ, c.time.hess_ψ};
    os << "                  ψ:" //
       << CountResult{c.ψ, c.time.ψ};
    os << "             grad_ψ:" //
       << CountResult{c.grad_ψ, c.time.grad_ψ};
    os << "           ψ_grad_ψ:" //
       << CountResult{c.ψ_grad_ψ, c.time.ψ_grad_ψ};
    os << "+ -----------------:\n" //
       << "              total:" //
       << CountResult{
              c.proj_diff_g + c.proj_multipliers + c.prox_grad_step + c.f +
                  c.grad_f + c.f_grad_f + c.f_g + c.grad_f_grad_g_prod + c.g +
                  c.grad_g_prod + c.grad_gi + c.jac_g + c.grad_L +
                  c.hess_L_prod + c.hess_L + c.hess_ψ_prod + c.hess_ψ + c.ψ +
                  c.grad_ψ + c.ψ_grad_ψ,
              c.time.proj_diff_g + c.time.proj_multipliers +
                  c.time.prox_grad_step + c.time.f + c.time.grad_f +
                  c.time.f_grad_f + c.time.f_g + c.time.grad_f_grad_g_prod +
                  c.time.g + c.time.grad_g_prod + c.time.grad_gi +
                  c.time.jac_g + c.time.grad_L + c.time.hess_L_prod +
                  c.time.hess_L + c.time.hess_ψ_prod + c.time.hess_ψ +
                  c.time.ψ + c.time.grad_ψ + c.time.ψ_grad_ψ};
    return os;
}
 
} // namespace alpaqa