474   PetscReal  moments_0[5], moments_1a[5], moments_1b[5], dt_init;  in go()  local
533   for (PetscInt i = 0; i < 5; i++) moments_0[i] = moments_1a[i] = moments_1b[i] = 0;  in go()
637                       moments_0[0] += w;                   // not thread safe  in go()
638                       moments_0[1] += w * ctx->v_0 * x[1]; // z-momentum  in go()
639                       moments_0[2] += w * 0.5 * ctx->v_0 * ctx->v_0 * v2;  in go()
670             const PetscReal N_inv = 1 / moments_0[0];  in go()
677                   moments_0[3] -= ww * PetscLogReal(ww);  in go()
679                 } else moments_0[4] -= w;  in go()
999       energy_error_rel = PetscAbsReal(moments_1b[2] - moments_0[2]) / moments_0[2];  in go()
1001 …g %%\n", (double)moments_0[0], (double)moments_0[1], (double)moments_0[2], (double)moments_0[3], 1…  in go()
1002 …a[1], (double)moments_1a[2], (double)moments_1a[3], 100 * (double)(moments_1a[4] / moments_0[0])));  in go()
1003 …b[1], (double)moments_1b[2], (double)moments_1b[3], 100 * (double)(moments_1b[4] / moments_0[0])));  in go()
1004 …ntropy generation = %e\n", (double)(moments_1a[3] - moments_0[3]), (double)(moments_1b[3] - moment…  in go()
1005 …ning = %e, Landau = %e (%g %d)\n", (double)(moments_1a[2] - moments_0[2]) / (double)moments_0[2], …  in go()