6 matrix *matrix_copy(pbqp *pbqp, matrix *m)
10 matrix *copy = obstack_alloc(&pbqp->obstack, sizeof(*copy));
14 len = m->rows * m->cols;
16 for (i = 0; i < len; ++i) {
17 copy->entries[i] = m->entries[i];
23 void matrix_add(matrix *sum, matrix *summand)
30 assert(sum->cols == summand->cols);
31 assert(sum->rows == summand->rows);
33 len = sum->rows * sum->cols;
35 for (i = 0; i < len; ++i) {
36 sum->entries[i] += summand->entries[i];