Searched refs:divp (Results 1 – 1 of 1) sorted by relevance
1612 PetscReal divp = 0.; in f0_momentum() local1618 divp += u_x[uOff_x[1] + e * dim + e]; in f0_momentum()1620 f0[d] += divp * u[uOff[1] + d] / n; // (div p) p / n in f0_momentum()1635 PetscReal divp = 0.; in f0_energy() local1641 divp += u_x[uOff_x[1] + d * dim + d]; in f0_energy()1643 f0[0] += 4. * divp * pr / n; // 4 div p pr / n in f0_energy()