log(-<realnumber>) now returns a real number

[ginac.git] / check / exam_indexed.cpp

diff --git a/check/exam_indexed.cpp b/check/exam_indexed.cpp
index 3d9b307db4f26fdcc33f8a6d245bb97fa41a973c..e20b56433b00ab8d7957f343ffaaa6f8a75733c9 100644 (file)
--- a/check/exam_indexed.cpp
+++ b/check/exam_indexed.cpp
@@ -375,6 +375,15 @@ static unsigned dummy_check()
        e = indexed(p, mu.toggle_variance(), mu) - indexed(p, nu, nu.toggle_variance());
        result += check_equal_simplify(e, 0);
 
+       // GiNaC 1.2.1 had a bug here because p.i*p.i -> (p.i)^2
+       e = indexed(p, i) * indexed(p, i) * indexed(p, j) + indexed(p, j);
+       ex fi = exprseq(e.get_free_indices());
+       if (!fi.is_equal(exprseq(j))) {
+               clog << "get_free_indices(" << e << ") erroneously returned "
+                    << fi << " instead of (.j)" << endl;
+               ++result;
+       }
+
        return result;
 }