Lines Matching refs:moments
278 …globMpArray[], Mat g_Mass[], Vec t_fhat[][EX30_MAX_NUM_THRDS], PetscReal moments[], Vec globXArray… in gridToParticles_private() argument
315 moments[0] += w; in gridToParticles_private()
316 moments[1] += w * ctx->v_0 * coords[p * dim + 1]; // z-momentum in gridToParticles_private()
317 moments[2] += w * 0.5 * ctx->v_0 * ctx->v_0 * v2; in gridToParticles_private()
322 const PetscReal N_inv = 1 / moments[0]; in gridToParticles_private()
323 …fo(grid_dm[0], "gridToParticles_private [%" PetscInt_FMT "], n = %g\n", v_id, (double)moments[0])); in gridToParticles_private()
335 moments[3] -= ww * PetscLogReal(ww); in gridToParticles_private()
337 } else moments[4] -= w; // keep track of density that is lost in gridToParticles_private()
375 PetscReal moments[5], e_grid[LANDAU_MAX_GRIDS]; in PostStep() local
382 for (PetscInt i = 0; i < 5; i++) moments[i] = 0; in PostStep()
410 moments[0] += w; in PostStep()
411 moments[1] += w * ctx->v_0 * coords[pp * dim + 1]; // z-momentum in PostStep()
412 moments[2] += w * 0.5 * ctx->v_0 * ctx->v_0 * v2; in PostStep()
419 const PetscReal N_inv = 1 / moments[0]; in PostStep()
431 moments[3] -= ww * PetscLogReal(ww); in PostStep()
432 } else moments[4] -= w; in PostStep()
437 …int)n, (double)t, (double)moments[0], (double)moments[1], (double)moments[2], (double)moments[3], … in PostStep()