Add timer for residual and omega

Author: filikat

highs/ipm/hipo/ipm/Info.h

@@ -46,6 +46,8 @@ struct Info {
   double matrix_structure_time{};
   double factor_time{};
   double solve_time{};
+  double residual_time{};
+  double omega_time{};
 
   // Counters
   Int factor_number{};

highs/ipm/hipo/ipm/Refine.cpp

@@ -5,13 +5,20 @@
 namespace hipo {
 
 void Solver::refine(NewtonDir& delta) {
+  Clock clock;
+
   NewtonDir correction(m_, n_);
   NewtonDir temp(m_, n_);
 
   // compute the residuals of the linear system in it_->ires
+  clock.start();
   it_->residuals6x6(delta);
+  info_.residual_time += clock.stop();
 
+  clock.start();
   double omega = computeOmega(delta);
+  info_.omega_time += clock.stop();
+
   double old_omega{};
 
   for (Int iter = 0; iter < kMaxIterRefine; ++iter) {
@@ -23,10 +30,15 @@ void Solver::refine(NewtonDir& delta) {
     temp = delta;
     temp.add(correction);
 
+    clock.start();
     it_->residuals6x6(temp);
+    info_.residual_time += clock.stop();
 
     old_omega = omega;
+
+    clock.start();
     omega = computeOmega(temp);
+    info_.omega_time += clock.stop();
 
     if (omega < old_omega) {
       delta = temp;

highs/ipm/hipo/ipm/Solver.cpp

@@ -1209,6 +1209,10 @@ void Solver::printSummary() const {
                << fix(info_.factor_time, 0, 2) << '\n';
     log_stream << textline("Solve time:") << fix(info_.solve_time, 0, 2)
                << '\n';
+    log_stream << textline("Residual time:") << fix(info_.residual_time, 0, 2)
+               << '\n';
+    log_stream << textline("Omega time:") << fix(info_.omega_time, 0, 2)
+               << '\n';
     log_stream << textline("Factorisations:") << integer(info_.factor_number)
                << '\n';
     log_stream << textline("Solves:") << integer(info_.solve_number) << '\n';