#include <iostream>
#include <vector>

using namespace std;

extern "C" {
#include <taucs.h>
}

int main()
{
	taucs_logfile("stdout");
	cout << 10 << endl;
	taucs_ccs_matrix *A = taucs_ccs_read_ijv("ConstraintMatrix.mmf", TAUCS_DOUBLE | TAUCS_SYMMETRIC | TAUCS_LOWER);
	cout << A->n << " " << A->m << endl;
  int* perm;
  int* invperm;
	taucs_ccs_matrix*  Aod;
	void* F;	
	cout << 20 << endl;
	taucs_ccs_order(A, &perm, &invperm, "metis");
	cout << 30 << endl;
  Aod = taucs_ccs_permute_symmetrically(A, perm, invperm);
	cout << 40 << endl;
	F = taucs_ccs_factor_llt_mf(Aod);	
	cout << 50 << endl;
	F = taucs_ccs_factor_ldlt(Aod);	
	cout << 60 << endl;
}