Lines Matching full:rho
29 PetscScalar rho, c; member
75 PetscCall(DMDASetFieldName(da, 0, "rho")); in main()
141 PetscScalar rho, c, rhoxx, cxx, cx, rhox, kcxrhox; in IFunction() local
175 f[0].rho = udot[0].rho; /* u[0].rho - 0.0; */ in IFunction()
181 f[Mx - 1].rho = udot[Mx - 1].rho; /* u[Mx-1].rho - 1.0; */ in IFunction()
190 rho = u[i].rho; in IFunction()
191 rhoxx = (-2.0 * rho + u[i - 1].rho + u[i + 1].rho) * sx; in IFunction()
196 rhox = .5 * (u[i + 1].rho - u[i - 1].rho) / hx; in IFunction()
198 kcxrhox = appctx->kappa * (cxx * rho + cx * rhox); in IFunction()
200 …rhox = appctx->kappa * ((u[i + 1].c - u[i].c) * u[i + 1].rho - (u[i].c - u[i - 1].c) * u[i].rho) *… in IFunction()
203 …rho = udot[i].rho - appctx->epsilon * rhoxx + kcxrhox - appctx->mu * PetscAbsScalar(rho) * (1.0 - … in IFunction()
204 …f[i].c = udot[i].c - appctx->delta * cxx + appctx->lambda * c + appctx->alpha * rho * c / (appct… in IFunction()
244 if (i < Mx - 1) u[i].rho = 0.0; in InitialConditions()
245 else u[i].rho = 1.0; in InitialConditions()