+
+xbt_edge_t* xbt_graph_spanning_tree_prim(xbt_graph_t g)
+{
+ int tree_size=0;
+ int tree_size_max=xbt_dynar_length(g->nodes)-1;
+ xbt_edge_t *tree = xbt_new0(xbt_edge_t,tree_size_max);
+ xbt_edge_t e,edge;
+ xbt_node_t node = NULL;
+ xbt_dynar_t edge_list = NULL;
+ xbt_heap_t heap = xbt_heap_new(10,NULL);
+ int cursor;
+
+ xbt_assert0(!(g->directed),
+ "Spanning trees do not make sense on directed graphs");
+
+ node = xbt_dynar_getfirst_as(g->nodes,xbt_node_t);
+ node->xbtdata = (void*) 1;
+ edge_list = node->out;
+ xbt_dynar_foreach(edge_list, cursor, e)
+ xbt_heap_push(heap,e, -(e->length));
+
+ while((edge=xbt_heap_pop(heap))) {
+ if((edge->src->xbtdata) && (edge->dst->xbtdata)) continue;
+ tree[tree_size++]=edge;
+ if(!(edge->src->xbtdata)) {
+ edge->src->xbtdata = (void*) 1;
+ edge_list = edge->src->out;
+ xbt_dynar_foreach(edge_list, cursor, e) {
+ xbt_heap_push(heap,e, -(e->length));
+ }
+ } else {
+ edge->dst->xbtdata = (void*) 1;
+ edge_list = edge->dst->out;
+ xbt_dynar_foreach(edge_list, cursor, e) {
+ xbt_heap_push(heap,e, -(e->length));
+ }
+ }
+ if(tree_size==tree_size_max) break;
+ }
+
+ xbt_heap_free(heap);
+
+ xbt_dynar_foreach(g->nodes, cursor, node) {
+ node->xbtdata = NULL;
+ }
+ return tree;
+}
+
+/********************* Import and Export ******************/