Function test_Exp_OP_MUL

tests/unit/test_algebraic_expression.c:650–693

}

void test_Exp_OP_MUL() {
	// Exp = A * I
	RG_Matrix A;
	RG_Matrix I;
	RG_Matrix res;

	// A
	// 1 1
	// 0 0
	RG_Matrix_new(&A, GrB_BOOL, 2, 2);
	RG_Matrix_setElement_BOOL(A, 0, 0);
	RG_Matrix_setElement_BOOL(A, 0, 1);
	RG_Matrix_wait(A, true); // force flush

	// I
	// 1 0
	// 0 1
	RG_Matrix_new(&I, GrB_BOOL, 2, 2);
	RG_Matrix_setElement_BOOL(I, 0, 0);
	RG_Matrix_setElement_BOOL(I, 1, 1);

	rax *matrices = raxNew();
	raxInsert(matrices, (unsigned char *)"A", strlen("A"), A, NULL);
	raxInsert(matrices, (unsigned char *)"I", strlen("I"), I, NULL);
	AlgebraicExpression *exp = AlgebraicExpression_FromString("A*I", matrices);

	// Matrix used for intermidate computations of AlgebraicExpression_Eval
	// but also contains the result of expression evaluation.
	RG_Matrix_new(&res, GrB_BOOL, 2, 2);
	AlgebraicExpression_Eval(exp, res);

	// Using the A matrix described above,
	// A * I = A.
	GrB_Matrix expected;
	RG_Matrix_export(&expected, A);
	TEST_ASSERT(_compare_matrices(expected, res));

	raxFree(matrices);
	RG_Matrix_free(&A);
	RG_Matrix_free(&I);
	RG_Matrix_free(&res);
	GrB_Matrix_free(&expected);
	AlgebraicExpression_Free(exp);
}

void test_Exp_OP_ADD_Transpose() {
	// Exp = A + Transpose(A)