Function _compare_matrices

tests/unit/test_algebraic_expression.c:223–275

}

bool _compare_matrices(GrB_Matrix expected, RG_Matrix actual) {
	GrB_Matrix a = expected;
	GrB_Matrix b = NULL;
	RG_Matrix_export(&b, actual);

	GrB_Index acols, arows, avals;
	GrB_Index bcols, brows, bvals;

	GrB_Matrix_ncols(&acols, a);
	GrB_Matrix_nrows(&arows, a);
	GrB_Matrix_nvals(&avals, a);
	GrB_Matrix_ncols(&bcols, b);
	GrB_Matrix_nrows(&brows, b);
	GrB_Matrix_nvals(&bvals, b);

	if(acols != bcols || arows != brows || avals != bvals) {
		printf("acols: %lu bcols: %lu\n", acols, bcols);
		printf("arows: %lu brows: %lu\n", arows, brows);
		printf("avals: %lu bvals: %lu\n", avals, bvals);

		GrB_Matrix_free(&b);
		return false;
	}

	GrB_Index aI[avals];     // array for returning row indices of tuples
	GrB_Index aJ[avals];     // array for returning col indices of tuples
	bool aX[avals];          // array for returning values of tuples
	GrB_Index bI[bvals];     // array for returning row indices of tuples
	GrB_Index bJ[bvals];     // array for returning col indices of tuples
	bool bX[bvals];          // array for returning values of tuples

	GrB_Matrix_wait(a, GrB_MATERIALIZE);
	GrB_Matrix_wait(b, GrB_MATERIALIZE);

	GrB_Matrix_extractTuples_BOOL(aI, aJ, aX, &avals, a);
	GrB_Matrix_extractTuples_BOOL(bI, bJ, bX, &bvals, b);

	for(int i = 0; i < avals; i++) {
		if(aI[i] != bI[i] || aJ[i] != bJ[i] || aX[i] != bX[i]) {
			printf("Matrix A \n");
			_print_matrix(a);
			printf("\n\n");
			printf("Matrix B \n");
			_print_matrix(b);
			printf("\n\n");
			GrB_Matrix_free(&b);
			return false;
		}
	}

	GrB_Matrix_free(&b);
	return true;
}

void _compare_algebraic_operand(AlgebraicExpression *a, AlgebraicExpression *b) {
	TEST_ASSERT(a->type == AL_OPERAND);