it is now enforced that indices of a metric tensor are of the same dimension

diff --git a/check/exam_archive.cpp b/check/exam_archive.cpp
index fca415811e42cca1bf608ea01a38d90f98bfac85..f51bb46ea5f3178471c85869a602135badbec2ec 100644 (file)
--- a/check/exam_archive.cpp
+++ b/check/exam_archive.cpp
@@ -40,15 +40,15 @@ unsigned exam_archive(void)
            * epsilon_tensor(idx(fail(), 3), idx(0, 3), idx(y/2, 3))
          + lorentz_g(
              varidx(lst(x, -11*y, acos(2*x).series(x==3-5*I, 3)) * color_ONE()
-               * metric_tensor(spinidx(0, 5, false, true), varidx(2, 4)), zeta(3)),
+               * metric_tensor(varidx(log(cos(128.0/(x*y))), 5), varidx(2, 5)), zeta(3)),
              varidx(diag_matrix(lst(-1, Euler, atan(x/y==-15*I/17)))
-               * delta_tensor(idx(x, 2), idx(wild(7), 3)), log(cos(128.0/(x*y))), true),
+               * delta_tensor(idx(x, 2), idx(wild(7), 3)), zeta(3), true),
              true
            )
          + dirac_gamma(varidx(mu, dim)) * dirac_gamma(varidx(mu, 4-dim, true))
            * color_T(idx(x, 8), 1) * color_h(idx(x, 8), idx(y, 8), idx(2, 8))
            * indexed(x, sy_anti(), idx(2*y+1, x), varidx(-mu, 5))
-         - 2.4275 * spinor_metric(spinidx(x), spinidx(y))
+         - 2.4275 * spinor_metric(spinidx(0, 2, false, true), spinidx(y))
          + abs(x).series(x == y, 4);
 
        archive ar;

diff --git a/ginac/tensor.cpp b/ginac/tensor.cpp
index 0ef65cc8ae66e963874e717384b3376dafacba4f..587a8170c38422c44b5719d5ca6bde2a0c9571a9 100644 (file)
--- a/ginac/tensor.cpp
+++ b/ginac/tensor.cpp
@@ -566,6 +566,9 @@ ex metric_tensor(const ex & i1, const ex & i2)
 {
        if (!is_ex_of_type(i1, varidx) || !is_ex_of_type(i2, varidx))
                throw(std::invalid_argument("indices of metric tensor must be of type varidx"));
+       ex dim = ex_to<idx>(i1).get_dim();
+       if (!dim.is_equal(ex_to<idx>(i2).get_dim()))
+               throw(std::invalid_argument("all indices of metric tensor must have the same dimension"));
 
        return indexed(tensmetric(), sy_symm(), i1, i2);
 }
@@ -574,6 +577,9 @@ ex lorentz_g(const ex & i1, const ex & i2, bool pos_sig)
 {
        if (!is_ex_of_type(i1, varidx) || !is_ex_of_type(i2, varidx))
                throw(std::invalid_argument("indices of metric tensor must be of type varidx"));
+       ex dim = ex_to<idx>(i1).get_dim();
+       if (!dim.is_equal(ex_to<idx>(i2).get_dim()))
+               throw(std::invalid_argument("all indices of metric tensor must have the same dimension"));
 
        return indexed(minkmetric(pos_sig), sy_symm(), i1, i2);
 }