/* #include "igraph_visitor.h"
#include "igraph_memory.h"
#include "igraph_adjlist.h"
#include "igraph_interface.h"
#include "igraph_dqueue.h"
#include "igraph_stack.h"
#include "config.h"*/
int igraph_bfsr(const igraph_t *graph,
igraph_integer_t root,igraph_vector_t *roots,
igraph_neimode_t mode, igraph_bool_t unreachable,
const igraph_vector_t *restricted,
igraph_vector_t *order, igraph_vector_t *rank,
igraph_vector_t *father,
igraph_vector_t *pred, igraph_vector_t *succ,
igraph_vector_t *dist, igraph_bfshandler_t *callback,
void *extra) {
igraph_dqueue_t Q;
long int no_of_nodes=igraph_vcount(graph);
long int actroot=0;
igraph_vector_char_t added;
igraph_lazy_adjlist_t adjlist;
long int act_rank=0;
long int pred_vec=-1;
long int rootpos=0;
long int noroots= roots ? igraph_vector_size(roots) : 1;
if (!roots && (root < 0 || root >= no_of_nodes)) {
IGRAPH_ERROR("Invalid root vertex in BFS", IGRAPH_EINVAL);
}
if (roots) {
igraph_real_t min, max;
igraph_vector_minmax(roots, &min, &max);
if (min < 0 || max >= no_of_nodes) {
IGRAPH_ERROR("Invalid root vertex in BFS", IGRAPH_EINVAL);
}
}
if (restricted) {
igraph_real_t min, max;
igraph_vector_minmax(restricted, &min, &max);
if (min < 0 || max >= no_of_nodes) {
IGRAPH_ERROR("Invalid vertex id in restricted set", IGRAPH_EINVAL);
}
}
if (mode != IGRAPH_OUT && mode != IGRAPH_IN &&
mode != IGRAPH_ALL) {
IGRAPH_ERROR("Invalid mode argument", IGRAPH_EINVMODE);
}
if (!igraph_is_directed(graph)) { mode=IGRAPH_ALL; }
IGRAPH_CHECK(igraph_vector_char_init(&added, no_of_nodes));
IGRAPH_FINALLY(igraph_vector_char_destroy, &added);
IGRAPH_CHECK(igraph_dqueue_init(&Q, 100));
IGRAPH_FINALLY(igraph_dqueue_destroy, &Q);
IGRAPH_CHECK(igraph_lazy_adjlist_init(graph, &adjlist, mode, /*simplify=*/ 0));
IGRAPH_FINALLY(igraph_lazy_adjlist_destroy, &adjlist);
/* Mark the vertices that are not in the restricted set, as already
found. Special care must be taken for vertices that are not in
the restricted set, but are to be used as 'root' vertices. */
if (restricted) {
long int i, n=igraph_vector_size(restricted);
igraph_vector_char_fill(&added, 1);
for (i=0; i<n; i++) {
long int v=(long int) VECTOR(*restricted)[i];
VECTOR(added)[v]=0;
}
}
/* Resize result vectors, and fill them with IGRAPH_NAN */
# define VINIT(v) if (v) { \
igraph_vector_resize((v), no_of_nodes); \
igraph_vector_fill((v), IGRAPH_NAN); }
VINIT(order);
VINIT(rank);
VINIT(father);
VINIT(pred);
VINIT(succ);
VINIT(dist);
# undef VINIT
while (1) {
/* Get the next root vertex, if any */
if (roots && rootpos < noroots) {
/* We are still going through the 'roots' vector */
actroot=(long int) VECTOR(*roots)[rootpos++];
} else if (!roots && rootpos==0) {
/* We have a single root vertex given, and start now */
actroot=root;
rootpos++;
} else if (rootpos==noroots && unreachable) {
/* We finished the given root(s), but other vertices are also
tried as root */
actroot=0;
rootpos++;
} else if (unreachable && actroot+1 < no_of_nodes) {
/* We are already doing the other vertices, take the next one */
actroot++;
} else {
/* No more root nodes to do */
break;
}
/* OK, we have a new root, start BFS */
if (VECTOR(added)[actroot]) { continue; }
IGRAPH_CHECK(igraph_dqueue_push(&Q, actroot));
IGRAPH_CHECK(igraph_dqueue_push(&Q, 0));
VECTOR(added)[actroot] = 1;
if (father) { VECTOR(*father)[actroot] = -1; }
pred_vec=-1;
while (!igraph_dqueue_empty(&Q)) {
long int actvect=(long int) igraph_dqueue_pop(&Q);
long int actdist=(long int) igraph_dqueue_pop(&Q);
long int succ_vec;
igraph_vector_t *neis=igraph_lazy_adjlist_get(&adjlist,
(igraph_integer_t) actvect);
long int i, n=igraph_vector_size(neis);
if (pred) { VECTOR(*pred)[actvect] = pred_vec; }
if (rank) { VECTOR(*rank) [actvect] = act_rank; }
if (order) { VECTOR(*order)[act_rank++] = actvect; }
if (dist) { VECTOR(*dist)[actvect] = actdist; }
igraph_vector_shuffle(neis);
for (i=0; i<n; i++) {
long int nei=(long int) VECTOR(*neis)[i];
if (! VECTOR(added)[nei]) {
VECTOR(added)[nei] = 1;
IGRAPH_CHECK(igraph_dqueue_push(&Q, nei));
IGRAPH_CHECK(igraph_dqueue_push(&Q, actdist+1));
if (father) { VECTOR(*father)[nei] = actvect; }
}
}
succ_vec = igraph_dqueue_empty(&Q) ? -1L :
(long int) igraph_dqueue_head(&Q);
if (callback) {
igraph_bool_t terminate=
callback(graph, (igraph_integer_t) actvect, (igraph_integer_t)
pred_vec, (igraph_integer_t) succ_vec,
(igraph_integer_t) act_rank-1, (igraph_integer_t) actdist,
extra);
if (terminate) {
igraph_lazy_adjlist_destroy(&adjlist);
igraph_dqueue_destroy(&Q);
igraph_vector_char_destroy(&added);
IGRAPH_FINALLY_CLEAN(3);
return 0;
}
}
if (succ) { VECTOR(*succ)[actvect] = succ_vec; }
pred_vec=actvect;
} /* while Q !empty */
} /* for actroot < no_of_nodes */
igraph_lazy_adjlist_destroy(&adjlist);
igraph_dqueue_destroy(&Q);
igraph_vector_char_destroy(&added);
IGRAPH_FINALLY_CLEAN(3);
return 0;
}
int main() {
igraph_dqueue_t q;
int i;
/* igraph_dqueue_init, igraph_dqueue_destroy, igraph_dqueue_empty */
igraph_dqueue_init(&q, 5);
if (!igraph_dqueue_empty(&q)) {
return 1;
}
igraph_dqueue_destroy(&q);
/* igraph_dqueue_push, igraph_dqueue_pop */
igraph_dqueue_init(&q, 4);
igraph_dqueue_push(&q, 1);
igraph_dqueue_push(&q, 2);
igraph_dqueue_push(&q, 3);
igraph_dqueue_push(&q, 4);
if (igraph_dqueue_pop(&q) != 1) {
return 2;
}
if (igraph_dqueue_pop(&q) != 2) {
return 3;
}
if (igraph_dqueue_pop(&q) != 3) {
return 4;
}
if (igraph_dqueue_pop(&q) != 4) {
return 5;
}
igraph_dqueue_destroy(&q);
/* igraph_dqueue_clear, igraph_dqueue_size */
igraph_dqueue_init(&q, 0);
if (igraph_dqueue_size(&q) != 0) {
return 6;
}
igraph_dqueue_clear(&q);
if (igraph_dqueue_size(&q) != 0) {
return 7;
}
for (i=0; i<10; i++) {
igraph_dqueue_push(&q, i);
}
igraph_dqueue_clear(&q);
if (igraph_dqueue_size(&q) != 0) {
return 8;
}
igraph_dqueue_destroy(&q);
/* TODO: igraph_dqueue_full */
/* igraph_dqueue_head, igraph_dqueue_back, igraph_dqueue_pop_back */
igraph_dqueue_init(&q, 0);
for (i=0; i<10; i++) {
igraph_dqueue_push(&q, i);
}
for (i=0; i<10; i++) {
if (igraph_dqueue_head(&q) != 0) {
return 9;
}
if (igraph_dqueue_back(&q) != 9-i) {
return 10;
}
if (igraph_dqueue_pop_back(&q) != 9-i) {
return 11;
}
}
igraph_dqueue_destroy(&q);
/* print */
igraph_dqueue_init(&q, 4);
igraph_dqueue_push(&q, 1);
igraph_dqueue_push(&q, 2);
igraph_dqueue_push(&q, 3);
igraph_dqueue_push(&q, 4);
igraph_dqueue_pop(&q);
igraph_dqueue_pop(&q);
igraph_dqueue_push(&q, 5);
igraph_dqueue_push(&q, 6);
igraph_dqueue_print(&q);
igraph_dqueue_clear(&q);
igraph_dqueue_print(&q);
igraph_dqueue_destroy(&q);
if (IGRAPH_FINALLY_STACK_SIZE() != 0) return 12;
return 0;
}
