Lines Matching refs:hs
173 PetscReal hf, hs, cfl_idt = 0; in FVRHSFunction() local
182 hs = (ctx->xmax - ctx->xmin) / 2.0 * (ctx->hratio + 1.0) / Mx; in FVRHSFunction()
214 cfl_idt = PetscMax(cfl_idt, PetscAbsScalar(maxspeed / hs)); in FVRHSFunction()
216 for (j = 0; j < dof; j++) f[(i - 1) * dof + j] -= ctx->flux[j] / hs; in FVRHSFunction()
219 for (j = 0; j < dof; j++) f[i * dof + j] += ctx->flux[j] / hs; 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()
232 for (j = 0; j < dof; j++) f[(i - 1) * dof + j] -= ctx->flux[j] / hs; 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()
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()
283 for (j = 0; j < dof; j++) f[i * dof + j] += ctx->flux[j] / hs; 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()
296 for (j = 0; j < dof; j++) f[(i - 1) * dof + j] -= ctx->flux[j] / hs; in FVRHSFunction()
299 for (j = 0; j < dof; j++) f[i * dof + j] += ctx->flux[j] / hs; in FVRHSFunction()
322 PetscReal hf, hs; in FVRHSFunctionslow() local
331 hs = (ctx->xmax - ctx->xmin) / 2.0 * (ctx->hratio + 1.0) / Mx; in FVRHSFunctionslow()
364 for (j = 0; j < dof; j++) f[(islow - 1) * dof + j] -= ctx->flux[j] / hs; in FVRHSFunctionslow()
367 for (j = 0; j < dof; j++) f[islow * dof + j] += ctx->flux[j] / hs; 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()
381 for (j = 0; j < dof; j++) f[(islow - 1) * dof + j] -= ctx->flux[j] / hs; 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()
394 for (j = 0; j < dof; j++) f[islow * dof + j] += ctx->flux[j] / 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()
408 for (j = 0; j < dof; j++) f[(islow - 1) * dof + j] -= ctx->flux[j] / hs; in FVRHSFunctionslow()
411 for (j = 0; j < dof; j++) f[islow * dof + j] += ctx->flux[j] / hs; in FVRHSFunctionslow()
425 for (j = 0; j < dof; j++) f[(islow - 1) * dof + j] -= ctx->flux[j] / hs; in FVRHSFunctionslow()
428 for (j = 0; j < dof; j++) f[islow * dof + j] += ctx->flux[j] / hs; in FVRHSFunctionslow()
444 PetscReal hf, hs; in FVRHSFunctionfast() local
453 hs = (ctx->xmax - ctx->xmin) / 2.0 * (ctx->hratio + 1.0) / 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()
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()
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()
559 const PetscReal hs = (ctx->xmax - ctx->xmin) / 2.0 * (ctx->hratio + 1.0) / Mx; in FVSample() local
564 if (i * hs + 0.5 * hs < (ctx->xmax - ctx->xmin) * 0.25) { in FVSample()
565 xi = ctx->xmin + 0.5 * hs + i * hs; in FVSample()
569 xj = xi + hs * (j - N / 2) / (PetscReal)N; in FVSample()
583 …i = ctx->xmin + (ctx->xmax - ctx->xmin) * 0.75 + 0.5 * hs + (i - count_slow / 2 - count_fast) * hs; in FVSample()
587 xj = xi + hs * (j - N / 2) / (PetscReal)N; in FVSample()
642 const PetscReal hs = (ctx->xmax - ctx->xmin) / 2.0 * (ctx->hratio + 1.0) / Mx; in SolutionErrorNorms() local
649 …if (i < count_slow / 2 || i > count_slow / 2 + count_fast - 1) *nrm1 += hs * PetscAbs(ptr_X[i] - p… in SolutionErrorNorms()
768 const PetscReal hs = (ctx.xmax - ctx.xmin) / 2.0 / count_slow; in main() local
781 mass_initial = mass_initial + hs * ptr_X0[i * dof + k]; in main()
782 mass_final = mass_final + hs * ptr_X[i * dof + k]; in main()
818 …if (i < count_slow / 2 || i > count_slow / 2 + count_fast - 1) nrm1 = nrm1 + hs * PetscAbs(ptr_X[i… in main()