Lines Matching refs:hf
173 PetscReal hf, hs, cfl_idt = 0; in FVRHSFunction() local
183 hf = (ctx->xmax - ctx->xmin) / 2.0 * (1.0 + 1.0 / ctx->hratio) / Mx; in FVRHSFunction()
223 alpha[0] = hs * hf / (hs + hs) / (hs + hs + hf); in FVRHSFunction()
224 gamma[0] = hs * (hs + hs) / (hs + hf) / (hs + hs + hf); in FVRHSFunction()
235 for (j = 0; j < dof; j++) f[i * dof + j] += ctx->flux[j] / hf; in FVRHSFunction()
239 alpha[0] = hf * hf / (hs + hf) / (hs + hf + hf); in FVRHSFunction()
240 gamma[0] = hf * (hs + hf) / (hf + hf) / (hs + hf + hf); in FVRHSFunction()
248 for (j = 0; j < dof; j++) f[(i - 1) * dof + j] -= ctx->flux[j] / hf; in FVRHSFunction()
251 for (j = 0; j < dof; j++) f[i * dof + j] += ctx->flux[j] / hf; in FVRHSFunction()
264 for (j = 0; j < dof; j++) f[(i - 1) * dof + j] -= ctx->flux[j] / hf; in FVRHSFunction()
267 for (j = 0; j < dof; j++) f[i * dof + j] += ctx->flux[j] / hf; in FVRHSFunction()
271 alpha[0] = hf * hs / (hf + hf) / (hf + hf + hs); in FVRHSFunction()
272 gamma[0] = hf * (hf + hf) / (hf + hs) / (hf + hf + hs); in FVRHSFunction()
280 for (j = 0; j < dof; j++) f[(i - 1) * dof + j] -= ctx->flux[j] / hf; in FVRHSFunction()
287 alpha[0] = hs * hs / (hf + hs) / (hf + hs + hs); in FVRHSFunction()
288 gamma[0] = hs * (hf + hs) / (hs + hs) / (hf + hs + hs); in FVRHSFunction()
322 PetscReal hf, hs; in FVRHSFunctionslow() local
332 hf = (ctx->xmax - ctx->xmin) / 2.0 * (1.0 + 1.0 / ctx->hratio) / Mx; in FVRHSFunctionslow()
372 alpha[0] = hs * hf / (hs + hs) / (hs + hs + hf); in FVRHSFunctionslow()
373 gamma[0] = hs * (hs + hs) / (hs + hf) / (hs + hs + hf); in FVRHSFunctionslow()
385 alpha[0] = hf * hs / (hf + hf) / (hf + hf + hs); in FVRHSFunctionslow()
386 gamma[0] = hf * (hf + hf) / (hf + hs) / (hf + hf + hs); in FVRHSFunctionslow()
399 alpha[0] = hs * hs / (hf + hs) / (hf + hs + hs); in FVRHSFunctionslow()
400 gamma[0] = hs * (hf + hs) / (hs + hs) / (hf + hs + hs); in FVRHSFunctionslow()
444 PetscReal hf, hs; in FVRHSFunctionfast() local
454 hf = (ctx->xmax - ctx->xmin) / 2.0 * (1.0 + 1.0 / ctx->hratio) / Mx; in FVRHSFunctionfast()
477 alpha[0] = hs * hf / (hs + hs) / (hs + hs + hf); in FVRHSFunctionfast()
478 gamma[0] = hs * (hs + hs) / (hs + hf) / (hs + hs + hf); in FVRHSFunctionfast()
486 for (j = 0; j < dof; j++) f[ifast * dof + j] += ctx->flux[j] / hf; in FVRHSFunctionfast()
491 alpha[0] = hf * hf / (hs + hf) / (hs + hf + hf); in FVRHSFunctionfast()
492 gamma[0] = hf * (hs + hf) / (hf + hf) / (hs + hf + hf); in FVRHSFunctionfast()
500 for (j = 0; j < dof; j++) f[(ifast - 1) * dof + j] -= ctx->flux[j] / hf; in FVRHSFunctionfast()
503 for (j = 0; j < dof; j++) f[ifast * dof + j] += ctx->flux[j] / hf; in FVRHSFunctionfast()
517 for (j = 0; j < dof; j++) f[(ifast - 1) * dof + j] -= ctx->flux[j] / hf; in FVRHSFunctionfast()
520 for (j = 0; j < dof; j++) f[ifast * dof + j] += ctx->flux[j] / hf; in FVRHSFunctionfast()
525 alpha[0] = hf * hs / (hf + hf) / (hf + hf + hs); in FVRHSFunctionfast()
526 gamma[0] = hf * (hf + hf) / (hf + hs) / (hf + hf + hs); in FVRHSFunctionfast()
534 for (j = 0; j < dof; j++) f[(ifast - 1) * dof + j] -= ctx->flux[j] / hf; in FVRHSFunctionfast()
560 const PetscReal hf = (ctx->xmax - ctx->xmin) / 2.0 * (1.0 + 1.0 / ctx->hratio) / Mx; in FVSample() local
573 …} else if ((ctx->xmax - ctx->xmin) * 0.25 + (i - count_slow / 2) * hf + 0.5 * hf < (ctx->xmax - ct… in FVSample()
574 xi = ctx->xmin + (ctx->xmax - ctx->xmin) * 0.25 + 0.5 * hf + (i - count_slow / 2) * hf; in FVSample()
578 xj = xi + hf * (j - N / 2) / (PetscReal)N; in FVSample()
643 const PetscReal hf = (ctx->xmax - ctx->xmin) / 2.0 * (1.0 + 1.0 / ctx->hratio) / Mx; in SolutionErrorNorms() local
650 else *nrm1 += hf * PetscAbs(ptr_X[i] - ptr_Y[i]); in SolutionErrorNorms()
769 const PetscReal hf = (ctx.xmax - ctx.xmin) / 2.0 / count_fast; in main() local
786 mass_initial = mass_initial + hf * ptr_X0[i * dof + k]; in main()
787 mass_final = mass_final + hf * ptr_X[i * dof + k]; in main()
819 else nrm1 = nrm1 + hf * PetscAbs(ptr_X[i] - ptr_XR[i]); in main()