#include "fracture.h"

int main(int argc, char *argv[]) {
	int opt;

	// defining variables to be (potentially) set by command line flags
	unsigned int N, L, filename_len;
	double beta, inf, cutoff;
	bound_t boundary;

	filename_len = 100;

	N = 100;
	L = 16;
	beta = .3;
	inf = 1e10;
	cutoff = 1e-10;
	boundary = FREE_BOUND;

	int boundary_int;
	char boundc2 = 'f';

	while ((opt = getopt(argc, argv, "n:L:b:B:dcoNsCrt")) != -1) {
		switch (opt) {
			case 'n':
				N = atoi(optarg);
				break;
			case 'L':
				L = atoi(optarg);
				break;
			case 'b':
				beta = atof(optarg);
				break;
			case 'B':
				boundary_int = atoi(optarg);
				switch (boundary_int) {
					case 0:
						boundary = FREE_BOUND;
						boundc2 = 'f';
						break;
					case 1:
						boundary = CYLINDER_BOUND;
						boundc2 = 'c';
						break;
					default:
						printf("boundary specifier must be 0 (FREE_BOUND) or 1 (CYLINDER_BOUND).\n");
				}
				break;
			default: /* '?' */
				exit(EXIT_FAILURE);
		}
	}

	char *break_filename = (char *)malloc(filename_len * sizeof(char));
	snprintf(break_filename, filename_len, "breaks_v_vc_%c_%u_%g.txt", boundc2, L, beta);
	FILE *break_out = fopen(break_filename, "a");
	free(break_filename);


	// start cholmod
	cholmod_common c;
	CHOL_F(start)(&c);

	(&c)->supernodal = CHOLMOD_SIMPLICIAL;


	graph_t *network = ini_voronoi_network(L, false, boundary, genfunc_hyperuniform, &c);
	net_t *perm_voltage_instance = create_instance(network, inf, true, true, &c);
	net_t *perm_current_instance = create_instance(network, inf, false, true, &c);
	double *fuse_thres = gen_fuse_thres(network->ne, network->edge_coords, beta, beta_scaling_flat);
	net_t *voltage_instance = copy_instance(perm_voltage_instance, &c);
	net_t *current_instance = copy_instance(perm_current_instance, &c);
	data_t *breaking_data_voltage = fracture_network(voltage_instance, fuse_thres, &c, cutoff);
	data_t *breaking_data_current = fracture_network(current_instance, fuse_thres, &c, cutoff);
	free_instance(voltage_instance, &c);
	free_instance(current_instance, &c);
	free_instance(perm_voltage_instance, &c);
	free_instance(perm_current_instance, &c);
	free(fuse_thres);

	FILE *net_out = fopen("network.txt", "w");
	for (unsigned int j = 0; j < network->nv; j++) {
		fprintf(net_out, "%f %f ", network->vert_coords[2 * j],
						network->vert_coords[2 * j + 1]);
	}
	fprintf(net_out, "\n");
	for (unsigned int j = 0; j < network->ne; j++) {
		fprintf(net_out, "%u %u ", network->ev[2 * j],
						network->ev[2 * j + 1]);
	}
	fprintf(net_out, "\n");
	for (unsigned int j = 0; j < network->dnv; j++) {
		fprintf(net_out, "%f %f ", network->dual_vert_coords[2 * j],
						network->dual_vert_coords[2 * j + 1]);
	}
	fprintf(net_out, "\n");
	for (unsigned int j = 0; j < network->ne; j++) {
		fprintf(net_out, "%u %u ", network->dev[2 * j],
						network->dev[2 * j + 1]);
	}

	free_net(network, &c);

	for (unsigned int j = 0; j < breaking_data_voltage->num_broken; j++) {
		fprintf(break_out, "%u %g %g ", breaking_data_voltage->break_list[j],
						breaking_data_voltage->extern_field[j], breaking_data_voltage->conductivity[j]);
	}
	fprintf(break_out, "\n");
	for (unsigned int j = 0; j < breaking_data_current->num_broken; j++) {
		fprintf(break_out, "%u %g %g ", breaking_data_current->break_list[j],
						breaking_data_current->extern_field[j], breaking_data_current->conductivity[j]);
	}
	fprintf(break_out, "\n");

	free_break_data(breaking_data_voltage);
	free_break_data(breaking_data_current);
	fclose(break_out);

	CHOL_F(finish)(&c);

	return 0;
}