IASolver.cpp

Go to the documentation of this file.

// IASolver.cpp
// Interval Assignment for Meshkit
//
#include "meshkit/IASolver.hpp"
#include "meshkit/IPData.hpp"
// #include "meshkit/IAMilp.hpp"
#include "meshkit/IASolverRelaxed.hpp"
#include "meshkit/IASolverInt.hpp"

#include <cstdlib>
#include <stdio.h>
#include <math.h>
#include <limits.h>

// timing tests
#include <time.h> //zzyk
#include <iostream>

namespace MeshKit 
{
    
IASolver::IASolver(IAData *ia_data_ptr, IASolution *ia_solution_ptr)
  : IASolverToolInt( ia_data_ptr, ia_solution_ptr ),
   debugging(true) 
  // debugging(false) 
  {}

IASolver::~IASolver() {}


bool IASolver::solve_relaxed()
{
  IASolverRelaxed relaxed(ia_data(), ia_solution(), !debugging); 
  bool succeeded=relaxed.solve();
  return succeeded;
}  


/* milp solution is too slow
bool IASolver::solve_int()
{
  IAMilp milp(this,this);
  bool succeeded=milp.solve();
  return succeeded;
}  
*/

bool IASolver::solve_int()
{
  // set relaxed solution to this's solution (from solve_relaxed)
  IASolverInt solver_int(ia_data(), ia_solution(), !debugging);
  bool succeeded=solver_int.solve();
  return succeeded;
}  

bool IASolver::solve_even()
{
  bool succeeded=true;
  return succeeded;
}  


bool IASolver::solve()
{

  if (debugging)
  {
    printf("Solving subproblem %p\n", this);
    print_problem();
  }
  
  // todo: subdivide problem into independent sub-problems for speed
  //zzyk
  clock_t t = clock();
  bool relaxed_succeeded = solve_relaxed();
  print_solution();

  t = clock()-t;
  float seconds = ((float)t)/CLOCKS_PER_SEC;
  std::cout << "relaxed solution " << ((relaxed_succeeded) ? "succeeded" : "failed") <<  " took " << seconds << " seconds" << std::endl;

  t = clock();
  bool int_succeeded = relaxed_succeeded && solve_int();
  t = clock()-t;
  seconds = ((float)t)/CLOCKS_PER_SEC;
  std::cout << "integer solution " << ((int_succeeded) ? "succeeded" : "failed") << " took " << seconds << " seconds" << std::endl;

  bool even_succeeded = int_succeeded && solve_even();
  
  if (debugging)
  {
    printf("==========IA summary:\n"); 
    if (relaxed_succeeded)
    {
      printf(" relaxed succeeded\n");
      if (int_succeeded)
      {
        printf(" integer succeeded\n");
        if (even_succeeded)
        {
          printf(" even succeeded\n");
        }
        else
        {
          printf(" even failed\n");        
        }
      }
      else
      {
        printf(" integer failed\n");        
      }
    }
    else
    {
      printf(" relaxed failed\n");
    }  
    print_solution();
  }
          
  return even_succeeded;
}


} // namespace MeshKit