/beegfs/hackathon/users/eoseret/qaas_runs_test/178-663-0018/intel/TeaLeaf/build/TeaLeaf/src/omp/solver_methods.cpp: 100 - 105
--------------------------------------------------------------------------------

100:   #pragma omp parallel for
101: #endif
102:   for (int jj = halo_depth; jj < y - halo_depth; ++jj) {
103:     for (int kk = halo_depth; kk < x - halo_depth; ++kk) {
104:       const int index = kk + jj * x;
105:       energy[index] = u[index] / density[index];
