8 vector *vector_alloc(pbqp *pbqp, unsigned length)
11 vector *vec = obstack_alloc(&pbqp->obstack, sizeof(*vec) + sizeof(*vec->entries) * length);
15 memset(vec->entries, 0, sizeof(*vec->entries) * length);
20 vector *vector_copy(pbqp *pbqp, vector *v)
22 unsigned len = v->len;
23 vector *copy = obstack_copy(&pbqp->obstack, v, sizeof(*copy) + sizeof(*copy->entries) * len);
29 void vector_add(vector *sum, vector *summand)
36 assert(sum->len == summand->len);
40 for (i = 0; i < len; ++i) {
41 if (sum->entries[i].data == INF_COSTS) continue;
43 if (summand->entries[i].data == INF_COSTS) {
44 sum->entries[i].data = INF_COSTS;
46 sum->entries[i].data += summand->entries[i].data;
51 void vector_set(vector *vec, unsigned index, num value)
53 assert(index < vec->len);
54 vec->entries[index].data = value;
58 void vector_set_description(vector *vec, unsigned index, char *name)
60 assert(index < vec->len);
61 vec->entries[index].name = name;
65 void vector_add_value(vector *vec, num value)
74 for (index = 0; index < len; ++index) {
75 if (vec->entries[index].data == INF_COSTS) continue;
77 if (value == INF_COSTS) {
78 vec->entries[index].data = INF_COSTS;
80 vec->entries[index].data += value;
85 void vector_add_matrix_col(vector *vec, pbqp_matrix *mat, unsigned col_index)
92 assert(vec->len == mat->rows);
93 assert(col_index < mat->cols);
97 for (index = 0; index < len; ++index) {
98 if (vec->entries[index].data == INF_COSTS) continue;
100 if (mat->entries[index * mat->cols + col_index] == INF_COSTS) {
101 vec->entries[index].data = INF_COSTS;
103 vec->entries[index].data += mat->entries[index * mat->cols + col_index];
108 void vector_add_matrix_row(vector *vec, pbqp_matrix *mat, unsigned row_index)
115 assert(vec->len == mat->cols);
116 assert(row_index < mat->rows);
120 for (index = 0; index < len; ++index) {
121 if (vec->entries[index].data == INF_COSTS) continue;
123 if (mat->entries[row_index * mat->cols + index] == INF_COSTS) {
124 vec->entries[index].data = INF_COSTS;
126 vec->entries[index].data += mat->entries[row_index * mat->cols + index];
131 num vector_get_min(vector *vec)
142 for (index = 0; index < len; ++index) {
143 num elem = vec->entries[index].data;
153 unsigned vector_get_min_index(vector *vec)
157 unsigned min_index = 0;
165 for (index = 0; index < len; ++index) {
166 num elem = vec->entries[index].data;