test_T4_007.cpp

Go to the documentation of this file.

#include "../../../src/FTensor.hpp"
#include "../test_for_zero.hpp"
#include <iostream>
using namespace FTensor;
using namespace std;
 
void test_T4_007(const Tensor4<double, 1, 2, 3, 4> &t4,
                 const Tensor3<double, 2, 3, 4> &t3_2) {
  Index<'i', 1> i;
  Index<'j', 2> j;
  Index<'k', 3> k;
  Index<'l', 4> l;
  Index<'n', 3> n;
 
  Number<0> N0;
  Number<1> N1;
  Number<2> N2;
 
  Tensor3<double, 2, 3, 4> t32;
  t32(j, k, l) = 0.;
 
  // Yield tensor 4 setting 4th slot
  {
    Tensor4<double, 2, 3, 4, 3> t_4;
    t_4(j, k, l, N0) = t3_2(j, k, l);
    t_4(j, k, l, N1) = t32(j, k, l);
    t_4(j, k, l, N2) = t32(j, k, l);
    for (int ii = 0; ii != 2; ++ii)
      for (int jj = 0; jj != 3; ++jj)
        for (int kk = 0; kk != 4; ++kk) {
          test_for_zero(t_4(ii, jj, kk, 0) - t3_2(ii, jj, kk),
                        "T4(i,j,k,N0)(" + to_string(ii) + "," + to_string(jj) +
                            "," + to_string(kk) + ")");
          for (int ll : {1, 2}) {
            test_for_zero(t_4(ii, jj, kk, ll), "T4(i,j,k,ll)(" + to_string(ii) +
                                                   "," + to_string(jj) + "," +
                                                   to_string(kk) + "," +
                                                   to_string(ll) + ")");
          }
        }
  }
  {
    Tensor4<double, 2, 3, 4, 3> t_4;
    t_4(j, k, l, N0) = t32(j, k, l);
    t_4(j, k, l, N1) = t3_2(j, k, l);
    t_4(j, k, l, N2) = t32(j, k, l);
    for (int ii = 0; ii != 2; ++ii)
      for (int jj = 0; jj != 3; ++jj)
        for (int kk = 0; kk != 4; ++kk) {
          test_for_zero(t_4(ii, jj, kk, 1) - t3_2(ii, jj, kk),
                        "T4(i,j,k,N1)(" + to_string(ii) + "," + to_string(jj) +
                            "," + to_string(kk) + ")");
          for (int ll : {0, 2}) {
            test_for_zero(t_4(ii, jj, kk, ll), "T4(i,j,k,ll)(" + to_string(ii) +
                                                   "," + to_string(jj) + "," +
                                                   to_string(kk) + "," +
                                                   to_string(ll) + ")");
          }
        }
  }
  {
    Tensor4<double, 2, 3, 4, 3> t_4;
    t_4(j, k, l, N0) = t32(j, k, l);
    t_4(j, k, l, N1) = t32(j, k, l);
    t_4(j, k, l, N2) = t3_2(j, k, l);
    for (int ii = 0; ii != 2; ++ii)
      for (int jj = 0; jj != 3; ++jj)
        for (int kk = 0; kk != 4; ++kk) {
          test_for_zero(t_4(ii, jj, kk, 2) - t3_2(ii, jj, kk),
                        "T4(i,j,k,N2)(" + to_string(ii) + "," + to_string(jj) +
                            "," + to_string(kk) + ")");
          for (int ll : {0, 1}) {
            test_for_zero(t_4(ii, jj, kk, ll), "T4(i,j,k,ll)(" + to_string(ii) +
                                                   "," + to_string(jj) + "," +
                                                   to_string(kk) + "," +
                                                   to_string(ll) + ")");
          }
        }
  }
 
  // Yield tensor 4 setting 3th slot
  {
    Tensor4<double, 2, 3, 3, 4> t_4;
    t_4(j, k, N0, l) = t3_2(j, k, l);
    t_4(j, k, N1, l) = t32(j, k, l);
    t_4(j, k, N2, l) = t32(j, k, l);
    for (int ii = 0; ii != 2; ++ii)
      for (int jj = 0; jj != 3; ++jj)
        for (int kk = 0; kk != 4; ++kk) {
          test_for_zero(t_4(ii, jj, 0, kk) - t3_2(ii, jj, kk),
                        "T4(i,j,N0,k)(" + to_string(ii) + "," + to_string(jj) +
                            "," + to_string(kk) + ")");
          for (int ll : {1, 2}) {
            test_for_zero(t_4(ii, jj, ll, kk), "T4(i,j,l, k)(" + to_string(ii) +
                                                   "," + to_string(jj) + "," +
                                                   to_string(ll) + "," +
                                                   to_string(kk) + ")");
          }
        }
  }
  {
    Tensor4<double, 2, 3, 3, 4> t_4;
    t_4(j, k, N0, l) = t32(j, k, l);
    t_4(j, k, N1, l) = t3_2(j, k, l);
    t_4(j, k, N2, l) = t32(j, k, l);
    for (int ii = 0; ii != 2; ++ii)
      for (int jj = 0; jj != 3; ++jj)
        for (int kk = 0; kk != 4; ++kk) {
          test_for_zero(t_4(ii, jj, 1, kk) - t3_2(ii, jj, kk),
                        "T4(i,j,N1,k)(" + to_string(ii) + "," + to_string(jj) +
                            "," + to_string(kk) + ")");
          for (int ll : {0, 2}) {
            test_for_zero(t_4(ii, jj, ll, kk), "T4(i,j,k,ll)(" + to_string(ii) +
                                                   "," + to_string(jj) + "," +
                                                   to_string(ll) + "," +
                                                   to_string(kk) + ")");
          }
        }
  }
  {
    Tensor4<double, 2, 3, 3, 4> t_4;
    t_4(j, k, N0, l) = t32(j, k, l);
    t_4(j, k, N1, l) = t32(j, k, l);
    t_4(j, k, N2, l) = t3_2(j, k, l);
    for (int ii = 0; ii != 2; ++ii)
      for (int jj = 0; jj != 3; ++jj)
        for (int kk = 0; kk != 4; ++kk) {
          test_for_zero(t_4(ii, jj, 2, kk) - t3_2(ii, jj, kk),
                        "T4(i,j,N2,k)(" + to_string(ii) + "," + to_string(jj) +
                            "," + to_string(kk) + ")");
          for (int ll : {0, 1}) {
            test_for_zero(t_4(ii, jj, ll, kk), "T4(i,j,k,ll)(" + to_string(ii) +
                                                   "," + to_string(jj) + "," +
                                                   to_string(ll) + "," +
                                                   to_string(kk) + ")");
          }
        }
  }
 
  //  Tensor4 to a Tensor4, yielding a Ddg.
  {
    Index<'i', 3> i;
    Index<'j', 3> j;
    Index<'k', 2> k;
    Index<'l', 2> l;
    Ddg<double, 3, 2> t_ddg;
    Tensor4<double, 3, 3, 2, 2> t_4;
    for (int ii = 0; ii != 3; ++ii)
      for (int jj = 0; jj != 3; ++jj)
        for (int kk = 0; kk != 2; ++kk)
          for (int ll = 0; ll != 2; ++ll) {
            t_4(ii, jj, kk, ll) = 1 + ii + 10. * jj + 100 * kk + 1000 * ll;
          }
    t_ddg(i, j, k, l) = t_4(i, j, k, l) || t_4(j, i, l, k);
    for (int ii = 0; ii != 3; ++ii)
      for (int jj = 0; jj != 3; ++jj)
        for (int kk = 0; kk != 2; ++kk)
          for (int ll = 0; ll != 2; ++ll) {
            test_for_zero(t_ddg(ii, jj, ll, kk) - t_4(ii, jj, kk, ll) -
                              t_4(jj, ii, ll, kk),
                          "T4(i,j,k,l)||T4(j,i,l,k)(" + to_string(ii) + "," +
                              to_string(jj) + "," + to_string(ll) + "," +
                              to_string(kk) + ")");
          }
  }
 
  //  Tensor4 times tensor 3 yields tensor 3
  {
    Tensor4<double, 1, 2, 3, 4> t_4;
    Tensor2<double, 1, 2> t_2_1;
    Tensor2<double, 3, 4> t_2_2;
    for (int ii = 0; ii != 1; ++ii)
      for (int jj = 0; jj != 2; ++jj) {
        t_2_1(ii, jj) = 1 + ii + 10. * jj;
      }
    for (int kk = 0; kk != 3; ++kk)
      for (int ll = 0; ll != 4; ++ll) {
        t_2_2(kk, ll) = 1 + 100 * kk + 1000 * ll;
      }
    t_4(i, j, k, l) = t_2_1(i, j) * t_2_2(k, l);
    Tensor3<double, 3, 4, 3> t_3_1;
    for (int kk = 0; kk != 3; ++kk)
      for (int ll = 0; ll != 4; ++ll)
        for (int mm = 0; mm != 3; ++mm) {
          t_3_1(kk, ll, mm) = 1 + 100 * kk + 1000 * ll + 10000 * mm;
        }
    Tensor3<double, 1, 2, 3> t_3_2;
    Index<'m', 3> m;
    t_3_2(i, j, m) = t_4(i, j, k, l) * t_3_1(k, l, m) -
                     t_2_1(i, j) * (t_3_1(k, l, m) * t_2_2(k, l));
    for (int ii = 0; ii != 1; ++ii)
      for (int jj = 0; jj != 2; ++jj)
        for (int mm = 0; mm != 3; ++mm) {
          test_for_zero(t_3_2(ii, jj, mm),
                        "T4(i,j,k,l)*T3(k,l,m)(" + to_string(ii) + "," +
                            to_string(jj) + "," + to_string(mm) + ")");
        }
  }
 
  //  Tensor4 times tensor 3 yields tensor 3
  {
    Tensor4<double, 1, 2, 3, 4> t_4;
    Tensor2<double, 1, 4> t_2_1;
    Tensor2<double, 2, 3> t_2_2;
    for (int ii = 0; ii != 1; ++ii)
      for (int ll = 0; ll != 4; ++ll) {
        t_2_1(ii, ll) = 1 + ii + 10. * ll;
      }
    for (int jj = 0; jj != 2; ++jj)
      for (int kk = 0; kk != 3; ++kk) {
        t_2_2(jj, kk) = 1 + 100 * jj + 1000 * kk;
      }
    t_4(i, j, k, l) = t_2_1(i, l) * t_2_2(j, k);
    Tensor3<double, 3, 3, 2> t_3_1;
    for (int mm = 0; mm != 3; ++mm)
      for (int kk = 0; kk != 3; ++kk)
        for (int jj = 0; jj != 2; ++jj) {
          t_3_1(mm, kk, jj) = 1 + 100 * mm + 1000 * kk + 10000 * jj;
        }
    Tensor3<double, 1, 4, 3> t_3_2;
    Index<'m', 3> m;
    t_3_2(i, l, m) = t_4(i, j, k, l) * t_3_1(m, k, j) -
                     t_2_1(i, l) * (t_3_1(m, k, j) * t_2_2(j, k));
    for (int ii = 0; ii != 1; ++ii)
      for (int ll = 0; ll != 4; ++ll)
        for (int mm = 0; mm != 3; ++mm) {
          test_for_zero(t_3_2(ii, ll, mm),
                        "T4(i,j,k,l)*T3(m,k,j)(" + to_string(ii) + "," +
                            to_string(ll) + "," + to_string(mm) + ")");
        }
  }
 
}