Function test_Exp_OP_ADD_Transpose

tests/unit/test_algebraic_expression.c:695–746

}

void test_Exp_OP_ADD_Transpose() {
	// Exp = A + Transpose(A)
	RG_Matrix res;
	GrB_Matrix expected;
	/* We must use matrices that we've added to the graph, or else
	 * the transpose optimization will erroneously use the adjacency matrix
	 * (since there is no schema associated with the operands).

	 * The 'visit' matrix represented by V has the form:
	 * 0 0 1 1
	 * 0 0 1 0
	 * 0 0 0 0
	 * 0 0 0 0
	 *
	 * Its transpose, TV, has the form:
	 * 0 0 0 0
	 * 0 0 0 0
	 * 1 1 0 0
	 * 1 0 0 0
	 *
	 * The expected result of the addition is:
	 * 0 0 1 1
	 * 0 0 1 0
	 * 1 1 0 0
	 * 1 0 0 0
	 */

	GraphContext *gc = QueryCtx_GetGraphCtx();
	Graph *g = gc->g;
	GrB_Index n = Graph_RequiredMatrixDim(g);

	GrB_Matrix_new(&expected, GrB_BOOL, n, n);
	GrB_Matrix_setElement_BOOL(expected, true, 0, 2);
	GrB_Matrix_setElement_BOOL(expected, true, 0, 3);
	GrB_Matrix_setElement_BOOL(expected, true, 1, 2);
	GrB_Matrix_setElement_BOOL(expected, true, 2, 0);
	GrB_Matrix_setElement_BOOL(expected, true, 2, 1);
	GrB_Matrix_setElement_BOOL(expected, true, 3, 0);
	// Matrix used for intermidate computations of AlgebraicExpression_Eval
	// but also contains the result of expression evaluation.
	RG_Matrix_new(&res, GrB_BOOL, n, n);
	AlgebraicExpression *exp = AlgebraicExpression_FromString("V+tV", _matrices);
	AlgebraicExpression_Eval(exp, res);

	// Using the A matrix described above,
	// A + Transpose(A) = B.
	TEST_ASSERT(_compare_matrices(expected, res));

	RG_Matrix_free(&res);
	GrB_Matrix_free(&expected);
	AlgebraicExpression_Free(exp);
}

void test_Exp_OP_MUL_Transpose() {
	// Exp = Transpose(A) * A