isl/dist/isl_mat.c

536 	__isl_take isl_mat *right)
541 	if (!left || !right)
544 	isl_assert(left->ctx, left->n_row == right->n_row, goto error);
547 	isl_assert(left->ctx, right->n_col >= 1, goto error);
552 	    isl_seq_first_non_zero(right->row[0]+1, right->n_col-1) == -1,
555 	sum = isl_mat_alloc(left->ctx, left->n_row, left->n_col + right->n_col - 1);
558 	isl_int_lcm(sum->row[0][0], left->row[0][0], right->row[0][0]);
560 	isl_int_divexact(right->row[0][0], sum->row[0][0], right->row[0][0]);
566 				right->row[0][0], right->row[i][0]);
570 				right->row[i]+1, right->row[0][0],
571 				right->n_col-1);
575 	isl_int_divexact(right->row[0][0], sum->row[0][0], right->row[0][0]);
577 	isl_mat_free(right);
581 	isl_mat_free(right);
977 	__isl_keep isl_mat **right, unsigned i, unsigned j)
980 	*right = isl_mat_swap_rows(*right, i, j);
982 	if (!*left || !*right)
988 	__isl_keep isl_mat *left, __isl_keep isl_mat *right, unsigned row)
991 	isl_seq_neg(right->row[row], right->row[row], right->n_col);
994 static void inv_subtract(__isl_keep isl_mat *left, __isl_keep isl_mat *right,
1002 	isl_seq_combine(right->row[i], right->ctx->one, right->row[i],
1003 			m, right->row[row], right->n_col);
1006 /* Compute inv(left)*right
1009 	__isl_take isl_mat *right)
1014 	if (!left || !right)
1018 	isl_assert(left->ctx, left->n_row == right->n_row, goto error);
1022 		return right;
1026 	right = isl_mat_cow(right);
1027 	if (!left || !right)
1042 			if (inv_exchange(&left, &right, pivot, row) < 0)
1045 			inv_oppose(left, right, row);
1052 			inv_subtract(left, right, row, first, a);
1054 				if (inv_exchange(&left, &right, row, first) < 0)
1071 			isl_seq_combine(right->row[i], a, right->row[i],
1072 					b, right->row[row], right->n_col);
1088 		isl_seq_scale(right->row[row], right->row[row],
1089 				left->row[row][row], right->n_col);
1094 	return right;
1097 	isl_mat_free(right);
1272 	__isl_take isl_mat *right)
1277 	if (!left || !right)
1279 	isl_assert(left->ctx, left->n_col == right->n_row, goto error);
1280 	prod = isl_mat_alloc(left->ctx, left->n_row, right->n_col);
1287 		isl_mat_free(right);
1293 				    left->row[i][0], right->row[0][j]);
1299 					    left->row[i][k], right->row[k][j]);
1303 	isl_mat_free(right);
1307 	isl_mat_free(right);