#include <limits>
#include <vector>

#include "current_info.hpp"
#include "graph.hpp"
#include <cholmod.h>

#ifdef FRACTURE_LONGINT

#define CHOL_F(x) cholmod_l_##x
#define CHOL_INT long int

#else

#define CHOL_F(x) cholmod_##x
#define CHOL_INT int

#endif

class problem {
private:
  const graph& G;
  unsigned axis;
  cholmod_dense* b;
  cholmod_factor* factor;
  cholmod_sparse* voltcurmat;
  cholmod_common* c;

public:
  current_info sol;

  problem(const graph&, unsigned, cholmod_common*);
  problem(const graph&, unsigned, cholmod_sparse*, cholmod_common*);
  problem(const problem&);
  ~problem();
  void solve(std::vector<bool>& fuses);
  void break_edge(unsigned, bool unbreak = false);
};