Function GfFillLoop

benchmarks/GfFill.cpp:46–62

// ===== Fill Green function Tensor Loop

static void GfFillLoop(benchmark::State &state) {
  long n_iw = state.range(0);

  auto iw_mesh_large = mesh::imfreq{1.0, Fermion, 2 * n_iw};
  auto M             = gf<mesh::prod<imfreq, imfreq>, scalar_valued>{{iw_mesh_large, iw_mesh_large}};

  auto iw_mesh = mesh::imfreq{1.0, Fermion, n_iw};
  auto iW_mesh = mesh::imfreq{1.0, Boson, n_iw};
  auto M4      = gf<prod<imfreq, imfreq, imfreq>, scalar_valued>{{iW_mesh, iw_mesh, iw_mesh}};

  for (auto _ : state)
    for (auto iW : iW_mesh)
      for (auto iw1 : iw_mesh)
        for (auto iw2 : iw_mesh) M4[iW, iw1, iw2] += M[iw1, iW + iw1] * M[iW + iw2, iw2];

  state.SetBytesProcessed(int64_t(state.iterations()) * 8 * n_iw * n_iw * n_iw * sizeof(scalar_t));
}
BENCHMARK(GfFillLoop)->RangeMultiplier(2)->Range(1 << 6, 1 << 8); //->Complexity();

// ===== Fill Green function Tensor Clefs