int test_normalized_laplacian(igraph_bool_t dir) {
igraph_t g;
igraph_matrix_t m;
igraph_vector_t vec;
igraph_matrix_init(&m, 1, 1);
igraph_bool_t ok = 1;
/* Undirected graph, no loop or multiple edges */
igraph_ring(&g, 5, dir, 0, 1);
igraph_laplacian(&g, &m, 1);
ok = ok && check_laplacian(&g, &m);
/* Add some loop edges */
igraph_vector_init_real(&vec, 4, 1.0, 1.0, 2.0, 2.0);
igraph_add_edges(&g, &vec, 0);
igraph_vector_destroy(&vec);
igraph_laplacian(&g, &m, 1);
ok = ok && check_laplacian(&g, &m);
/* Duplicate some edges */
igraph_vector_init_real(&vec, 4, 1.0, 2.0, 3.0, 4.0);
igraph_add_edges(&g, &vec, 0);
igraph_vector_destroy(&vec);
igraph_laplacian(&g, &m, 1);
ok = ok && check_laplacian(&g, &m);
igraph_destroy(&g);
igraph_matrix_destroy(&m);
if (ok)
printf("OK\n");
return !ok;
}
int test_normalized_laplacian(igraph_vector_t *w, igraph_bool_t dir) {
igraph_t g;
igraph_matrix_t m, m2;
igraph_sparsemat_t sm;
igraph_vector_t vec, *weights = 0;
igraph_bool_t ok = 1;
igraph_matrix_init(&m, 1, 1);
igraph_sparsemat_init(&sm, 0, 0, 0);
if (w) {
weights = (igraph_vector_t*)calloc(1, sizeof(igraph_vector_t));
igraph_vector_copy(weights, w);
}
/* Undirected graph, no loop or multiple edges */
igraph_ring(&g, 5, dir, 0, 1);
igraph_laplacian(&g, &m, &sm, 1, weights);
igraph_matrix_init(&m2, 0, 0);
igraph_sparsemat_as_matrix(&m2, &sm);
if (!igraph_matrix_all_e_tol(&m, &m2, 0)) {
return 44;
}
igraph_matrix_destroy(&m2);
ok = ok && check_laplacian(&g, &m, weights);
/* Add some loop edges */
igraph_vector_init_real(&vec, 4, 1.0, 1.0, 2.0, 2.0);
igraph_add_edges(&g, &vec, 0);
igraph_vector_destroy(&vec);
if (weights) {
igraph_vector_push_back(weights, 2);
igraph_vector_push_back(weights, 2);
}
igraph_laplacian(&g, &m, &sm, 1, weights);
igraph_matrix_init(&m2, 0, 0);
igraph_sparsemat_as_matrix(&m2, &sm);
if (!igraph_matrix_all_e_tol(&m, &m2, 0)) {
return 45;
}
igraph_matrix_destroy(&m2);
ok = ok && check_laplacian(&g, &m, weights);
/* Duplicate some edges */
igraph_vector_init_real(&vec, 4, 1.0, 2.0, 3.0, 4.0);
igraph_add_edges(&g, &vec, 0);
igraph_vector_destroy(&vec);
if (weights) {
igraph_vector_push_back(weights, 3);
igraph_vector_push_back(weights, 3);
}
igraph_laplacian(&g, &m, &sm, 1, weights);
igraph_matrix_init(&m2, 0, 0);
igraph_sparsemat_as_matrix(&m2, &sm);
if (!igraph_matrix_all_e_tol(&m, &m2, 0)) {
return 46;
}
igraph_matrix_destroy(&m2);
ok = ok && check_laplacian(&g, &m, weights);
igraph_destroy(&g);
igraph_matrix_destroy(&m);
if (weights) {
igraph_vector_destroy(weights); free(weights);
}
if (ok)
printf("OK\n");
return !ok;
}
