#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 = create_instance(network, 1e-14, true, true, &c); double *fuse_thres = gen_fuse_thres(network->ne, network->edge_coords, 0.001, beta_scaling_flat); fracture_network(instance, fuse_thres, &c, 1e-10); double *corr = get_corr(instance, NULL, &c); for (int i = 0; i < 2 * width; i++) { printf("%f ", corr[i]); } printf("\n"); free_instance(instance, &c); free_net(network, &c); CHOL_F(finish)(&c); return 0; }