Searched refs:u_div (Results 1 – 1 of 1) sorted by relevance
810 PetscReal u_div = 0.0; in f1_conduct_v() local814 for (c = 0; c < Nc; ++c) u_div += PetscRealPart(u_x[uOff_x[VEL] + c * dim + c]); in f1_conduct_v()820 f1[c * dim + c] -= 2.0 * coef / 3.0 * u_div; in f1_conduct_v()