1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
|
#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;
}
|