Grab a node once instead of an index every iteration.
[libfirm] / kaps.c
1 #include "adt/array.h"
2
3 #include "kaps.h"
4 #include "matrix.h"
5 #include "pbqp_edge.h"
6 #include "pbqp_edge_t.h"
7 #include "pbqp_node.h"
8 #include "pbqp_node_t.h"
9 #include "vector.h"
10
11 pbqp_node *get_node(pbqp *pbqp, unsigned index)
12 {
13         return pbqp->nodes[index];
14 }
15
16 pbqp_edge *get_edge(pbqp *pbqp, unsigned src_index, unsigned tgt_index)
17 {
18         int i;
19         int len;
20
21         if (tgt_index < src_index) {
22                 unsigned tmp = src_index;
23                 src_index    = tgt_index;
24                 tgt_index    = tmp;
25         }
26
27         pbqp_node *src_node = get_node(pbqp, src_index);
28         pbqp_node *tgt_node = get_node(pbqp, tgt_index);
29         assert(src_node);
30         assert(tgt_node);
31
32         len = ARR_LEN(src_node->edges);
33
34         for (i = 0; i < len; ++i) {
35                 pbqp_edge *cur_edge = src_node->edges[i];
36                 if (cur_edge->tgt == tgt_node) {
37                         return cur_edge;
38                 }
39         }
40
41         return NULL;
42 }
43
44 pbqp *alloc_pbqp(unsigned number_nodes)
45 {
46         pbqp* pbqp = xmalloc(sizeof(*pbqp));
47
48         obstack_init(&pbqp->obstack);
49
50         pbqp->solution = 0;
51         pbqp->num_nodes = number_nodes;
52         pbqp->dump_file = NULL;
53         pbqp->nodes = obstack_alloc(&pbqp->obstack, number_nodes
54                         * sizeof(*pbqp->nodes));
55         memset(pbqp->nodes, 0, number_nodes * sizeof(*pbqp->nodes));
56
57         return pbqp;
58 }
59
60 void free_pbqp(pbqp *pbqp)
61 {
62         obstack_free(&pbqp->obstack, NULL);
63         xfree(pbqp);
64 }
65
66 void add_node_costs(pbqp *pbqp, unsigned node_index, vector *costs)
67 {
68         pbqp_node *node = get_node(pbqp, node_index);
69
70         if (node == NULL) {
71                 node = alloc_node(pbqp, node_index, costs);
72                 pbqp->nodes[node_index] = node;
73         } else {
74                 vector_add(node->costs, costs);
75         }
76 }
77
78 void add_edge_costs(pbqp *pbqp, unsigned src_index, unsigned tgt_index,
79                 pbqp_matrix *costs)
80 {
81         pbqp_edge *edge = get_edge(pbqp, src_index, tgt_index);
82
83         if (edge == NULL) {
84                 edge = alloc_edge(pbqp, src_index, tgt_index, costs);
85         } else {
86                 pbqp_matrix_add(edge->costs, costs);
87         }
88 }
89
90 num get_node_solution(pbqp *pbqp, unsigned node_index)
91 {
92         pbqp_node *node = get_node(pbqp, node_index);
93         assert(node);
94
95         return node->solution;
96 }
97
98 num get_solution(pbqp *pbqp)
99 {
100         return pbqp->solution;
101 }
102
103 void set_dumpfile(pbqp *pbqp, FILE *f)
104 {
105         assert(pbqp);
106         pbqp->dump_file = f;
107 }