]> www.ginac.de Git - ginac.git/commitdiff
- Disabled three badly written checks that fail with the new tinfo mechanism.
authorJens Vollinga <vollinga@thep.physik.uni-mainz.de>
Thu, 19 Jan 2006 19:37:16 +0000 (19:37 +0000)
committerJens Vollinga <vollinga@thep.physik.uni-mainz.de>
Thu, 19 Jan 2006 19:37:16 +0000 (19:37 +0000)
check/exam_clifford.cpp
check/exam_normalization.cpp

index c0650a478bc43f7bc119f8655ca88bf82a27e8da..53b57906557727df751f9bf6b3146f69b4514f9e 100644 (file)
@@ -236,7 +236,9 @@ static unsigned clifford_check3()
        e = dirac_gamma(mu, 0) * dirac_gamma(mu.toggle_variance(), 1) * dirac_gamma(nu, 0) * dirac_gamma(nu.toggle_variance(), 1);
        result += check_equal_simplify(dirac_trace(e, 0), 4 * dim * dirac_ONE(1));
        result += check_equal_simplify(dirac_trace(e, 1), 4 * dim * dirac_ONE(0));
        e = dirac_gamma(mu, 0) * dirac_gamma(mu.toggle_variance(), 1) * dirac_gamma(nu, 0) * dirac_gamma(nu.toggle_variance(), 1);
        result += check_equal_simplify(dirac_trace(e, 0), 4 * dim * dirac_ONE(1));
        result += check_equal_simplify(dirac_trace(e, 1), 4 * dim * dirac_ONE(0));
-       result += check_equal_simplify(dirac_trace(e, 2), canonicalize_clifford(e)); // e will be canonicalized by the calculation of the trace
+       // Fails with new tinfo mechanism because the order of gamme matrices with different rl depends on luck. 
+       // TODO: better check.
+       //result += check_equal_simplify(dirac_trace(e, 2), canonicalize_clifford(e)); // e will be canonicalized by the calculation of the trace
        result += check_equal_simplify(dirac_trace(e, lst(0, 1)), 16 * dim);
 
        return result;
        result += check_equal_simplify(dirac_trace(e, lst(0, 1)), 16 * dim);
 
        return result;
index c0a2e232735287e8bdc33cf017565e35998db54b..b5eaf707cee577d42d3f6f63e2d9fb372c987675 100644 (file)
@@ -80,18 +80,24 @@ static unsigned exam_normal2()
        d = (x + y) * (w + z);
        result += check_normal(e, d);
        
        d = (x + y) * (w + z);
        result += check_normal(e, d);
        
-       e = (pow(x, 2) - pow(y, 2)) / pow(x-y, 3);
-       d = (x + y) / pow(x - y, 2);
-       result += check_normal(e, d);
+       // Fails stochastically with the new tinfo mechanism, because
+       // sometimes the equivalent answer ... / pow(y - x, 2) is calculated.
+       // TODO: make check for both cases.
+//     e = (pow(x, 2) - pow(y, 2)) / pow(x-y, 3);
+//     d = (x + y) / pow(x - y, 2);
+//     result += check_normal(e, d);
        
        e = (pow(x, -1) + x) / (pow(x , 2) * 2 + 2);
        d = pow(x * 2, -1);
        result += check_normal(e, d);
        
        
        e = (pow(x, -1) + x) / (pow(x , 2) * 2 + 2);
        d = pow(x * 2, -1);
        result += check_normal(e, d);
        
+       // Fails stochastically with the new tinfo mechanism, because
+       // sometimes the equivalent answer ... / pow(y - x, 2) is calculated.
+       // TODO: make check for both cases.
        // Fraction cancellation with rational coefficients
        // Fraction cancellation with rational coefficients
-       e = (pow(x, 2) - pow(y, 2)) / pow(x/2 - y/2, 3);
-       d = (8 * x + 8 * y) / pow(x - y, 2);
-       result += check_normal(e, d);
+//     e = (pow(x, 2) - pow(y, 2)) / pow(x/2 - y/2, 3);
+//     d = (8 * x + 8 * y) / pow(x - y, 2);
+//     result += check_normal(e, d);
        
        // Fraction cancellation with rational coefficients
        e = z/5 * (x/7 + y/10) / (x/14 + y/20);
        
        // Fraction cancellation with rational coefficients
        e = z/5 * (x/7 + y/10) / (x/14 + y/20);