improved the CLN C-source output (integers are printed as integers, and
[ginac.git] / check / exam_indexed.cpp
index a8c7cbd698f6da73a6673bb32e60cdfd44c387c9..23362a83e79b70c0a4e4b9ddcb6504a04698750c 100644 (file)
@@ -148,6 +148,11 @@ static unsigned epsilon_check(void)
        return result;
 }
 
+DECLARE_FUNCTION_2P(symm_fcn)
+REGISTER_FUNCTION(symm_fcn, set_symmetry(sy_symm(0, 1)));
+DECLARE_FUNCTION_2P(anti_fcn)
+REGISTER_FUNCTION(anti_fcn, set_symmetry(sy_anti(0, 1)));
+
 static unsigned symmetry_check(void)
 {
        // check symmetric/antisymmetric objects
@@ -155,7 +160,7 @@ static unsigned symmetry_check(void)
        unsigned result = 0;
 
        idx i(symbol("i"), 3), j(symbol("j"), 3), k(symbol("k"), 3), l(symbol("l"), 3);
-       symbol A("A"), B("B");
+       symbol A("A"), B("B"), C("C");
        ex e;
 
        result += check_equal(indexed(A, sy_symm(), i, j), indexed(A, sy_symm(), j, i));
@@ -187,6 +192,13 @@ static unsigned symmetry_check(void)
        result += check_equal(symmetrize(e), 0);
        result += check_equal(antisymmetrize(e), e);
 
+       e = (indexed(A, sy_anti(), i, j, k, l) * (indexed(B, j) * indexed(C, k) + indexed(B, k) * indexed(C, j)) + indexed(B, i, l)).expand();
+       result += check_equal_simplify(e, indexed(B, i, l));
+
+       result += check_equal(symm_fcn(0, 1) + symm_fcn(1, 0), 2*symm_fcn(0, 1));
+       result += check_equal(anti_fcn(0, 1) + anti_fcn(1, 0), 0);
+       result += check_equal(anti_fcn(0, 0), 0);
+
        return result;
 }