#include "fracture.h"

int main() {
	cholmod_common c;
	CHOL_F(start)(&c);

	unsigned int width = 64;
	unsigned int n = pow(width / 2 + 1, 2) + pow((width + 1) / 2, 2);

	graph_t *network = ini_square_network(width, true, false, &c);
	net_t *instance = net_create(network, 1e-14, 0.001, 0, true, &c);
	net_fracture(instance, &c, 1e-10);

	double *corr = get_corr(instance, NULL, &c);

	for (int i = 0; i < 2 * width; i++) {
		printf("%f ", corr[i]);
	}
	printf("\n");

	net_free(instance, &c);
	graph_free(network, &c);

	CHOL_F(finish)(&c);

	return 0;
}