| /petsc/src/ts/tutorials/multirate/ |
| H A D | ex8.c | 225 PetscScalar *x, *f, *slope; in FVRHSFunction_3WaySplit() local
243 …PetscCall(DMDAGetArray(da, PETSC_TRUE, &slope)); /* contains ghost points … in FVRHSFunction_3WaySplit()
282 slope[i * dof + j] = tmp; in FVRHSFunction_3WaySplit()
293 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2; in FVRHSFunction_3WaySplit()
294 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxs / 2; in FVRHSFunction_3WaySplit()
305 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2; in FVRHSFunction_3WaySplit()
306 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxm / 2; in FVRHSFunction_3WaySplit()
317 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxm / 2; in FVRHSFunction_3WaySplit()
318 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxs / 2; in FVRHSFunction_3WaySplit()
329 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxm / 2; in FVRHSFunction_3WaySplit()
[all …]
|
| H A D | ex6.c | 210 PetscScalar *x, *f, *slope; in FVRHSFunction_2WaySplit() local
227 …PetscCall(DMDAGetArray(da, PETSC_TRUE, &slope)); /* contains ghost points … in FVRHSFunction_2WaySplit()
265 slope[i * dof + j] = tmp; in FVRHSFunction_2WaySplit()
276 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2; in FVRHSFunction_2WaySplit()
277 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxs / 2; in FVRHSFunction_2WaySplit()
288 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2; in FVRHSFunction_2WaySplit()
289 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxf / 2; in FVRHSFunction_2WaySplit()
300 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxf / 2; in FVRHSFunction_2WaySplit()
301 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxf / 2; in FVRHSFunction_2WaySplit()
312 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxf / 2; in FVRHSFunction_2WaySplit()
[all …]
|
| H A D | ex5.c | 192 PetscScalar *x, *f, *slope; in FVRHSFunctionslow() local
206 PetscCall(DMDAGetArray(da, PETSC_TRUE, &slope)); in FVRHSFunctionslow()
245 slope[i * dof + j] = tmp; in FVRHSFunctionslow()
257 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hx / 2; in FVRHSFunctionslow()
258 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hx / 2; in FVRHSFunctionslow()
273 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hx / 2; in FVRHSFunctionslow()
274 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hx / 2; in FVRHSFunctionslow()
285 PetscCall(DMDARestoreArray(da, PETSC_TRUE, &slope)); in FVRHSFunctionslow()
297 PetscScalar *x, *f, *slope; in FVRHSFunctionfast() local
311 PetscCall(DMDAGetArray(da, PETSC_TRUE, &slope)); in FVRHSFunctionfast()
[all …]
|
| H A D | ex4.c | 554 PetscScalar *x, *f, *slope; in FVRHSFunction_2WaySplit() local
571 …PetscCall(DMDAGetArray(da, PETSC_TRUE, &slope)); /* contains ghost points … in FVRHSFunction_2WaySplit()
634 slope[i * dof + j] = tmp; in FVRHSFunction_2WaySplit()
645 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2; in FVRHSFunction_2WaySplit()
646 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxs / 2; in FVRHSFunction_2WaySplit()
657 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2; in FVRHSFunction_2WaySplit()
658 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxf / 2; in FVRHSFunction_2WaySplit()
669 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxf / 2; in FVRHSFunction_2WaySplit()
670 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxf / 2; in FVRHSFunction_2WaySplit()
681 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxf / 2; in FVRHSFunction_2WaySplit()
[all …]
|
| H A D | finitevolume1d.c | 528 PetscScalar *x, *f, *slope; in FVRHSFunction() local
542 PetscCall(DMDAGetArray(da, PETSC_TRUE, &slope)); /* contains ghost points */ in FVRHSFunction()
580 slope[i * dof + j] = tmp; in FVRHSFunction()
590 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hx / 2; in FVRHSFunction()
591 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hx / 2; in FVRHSFunction()
604 PetscCall(DMDARestoreArray(da, PETSC_TRUE, &slope)); in FVRHSFunction()
|
| /petsc/src/sys/utils/ |
| H A D | mathfit.c | 20 …gression(PetscInt n, const PetscReal x[], const PetscReal y[], PetscReal *slope, PetscReal *interc… in PetscLinearRegression() argument
30 PetscAssertPointer(slope, 4); in PetscLinearRegression()
69 *slope = beta[1]; in PetscLinearRegression()
|
| /petsc/src/ts/utils/ |
| H A D | tsconvest.c | 54 PetscReal *dt, *x, *y, slope, intercept; in PetscConvEstGetConvRateTS_Temporal_Private() local
96 PetscCall(PetscLinearRegression(Nr + 1, x, y, &slope, &intercept)); in PetscConvEstGetConvRateTS_Temporal_Private()
98 alpha[f] = slope; in PetscConvEstGetConvRateTS_Temporal_Private()
122 PetscReal *x, *y, slope, intercept; in PetscConvEstGetConvRateTS_Spatial_Private() local
244 PetscCall(PetscLinearRegression(Nr + 1, x, y, &slope, &intercept)); in PetscConvEstGetConvRateTS_Spatial_Private()
246 alpha[f] = -slope * dim; in PetscConvEstGetConvRateTS_Spatial_Private()
|
| H A D | dmplexts.c | 366 PetscReal dt, shift, slope, intercept; in DMTSCheckJacobian() local
454 PetscCall(PetscLinearRegression(Nv, hs, es, &slope, &intercept)); in DMTSCheckJacobian()
458 if (convRate) *convRate = slope; in DMTSCheckJacobian()
461 …PetscAbsReal(2 - slope) <= tol, comm, PETSC_ERR_ARG_WRONG, "Taylor approximation convergence rate … in DMTSCheckJacobian()
463 …n) PetscCall(PetscPrintf(comm, "Taylor approximation converging at order %3.2f\n", (double)slope)); in DMTSCheckJacobian()
|
| /petsc/src/snes/utils/ |
| H A D | dmplexsnes.c | 902 PetscReal slope, intercept; in DMSNESCheckJacobian() local
984 PetscCall(PetscLinearRegression(Nv, hs, es, &slope, &intercept)); in DMSNESCheckJacobian()
988 …PetscAbsReal(2 - slope) <= tol, comm, PETSC_ERR_ARG_WRONG, "Taylor approximation convergence rate … in DMSNESCheckJacobian()
991 if (convRate) *convRate = slope; in DMSNESCheckJacobian()
993 …n) PetscCall(PetscPrintf(comm, "Taylor approximation converging at order %3.2f\n", (double)slope)); in DMSNESCheckJacobian()
|
| H A D | convest.c | 338 PetscReal *x, *y, slope, intercept; in PetscConvEstGetConvRateSNES_Private() local
439 PetscCall(PetscLinearRegression(Nr + 1, x, y, &slope, &intercept)); in PetscConvEstGetConvRateSNES_Private()
441 alpha[f] = -slope * dim; in PetscConvEstGetConvRateSNES_Private()
|
| /petsc/src/ksp/ksp/interface/ |
| H A D | itfunc.c | 2617 PetscReal *x, *y, slope, intercept, mean = 0.0, var = 0.0, res = 0.0; in KSPComputeConvergenceRate() local
2634 PetscCall(PetscLinearRegression(n, x, y, &slope, &intercept)); in KSPComputeConvergenceRate()
2636 res += PetscSqr(y[k] - (slope * x[k] + intercept)); in KSPComputeConvergenceRate()
2640 if (cr) *cr = PetscExpReal(slope); in KSPComputeConvergenceRate()
2657 PetscCall(PetscLinearRegression(n, x, y, &slope, &intercept)); in KSPComputeConvergenceRate()
2659 res += PetscSqr(y[k] - (slope * x[k] + intercept)); in KSPComputeConvergenceRate()
2663 if (ce) *ce = PetscExpReal(slope); in KSPComputeConvergenceRate()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex9.c | 1165 PetscScalar *x, *f, *slope; in FVRHSFunction() local
1181 PetscCall(DMDAGetArray(da, PETSC_TRUE, &slope)); in FVRHSFunction()
1219 slope[i * dof + j] = tmp; in FVRHSFunction()
1229 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hx / 2; in FVRHSFunction()
1230 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hx / 2; in FVRHSFunction()
1245 PetscCall(DMDARestoreArray(da, PETSC_TRUE, &slope)); in FVRHSFunction()
|
| /petsc/doc/miscellaneous/ |
| H A D | applications_publications.md | 73 deformation rate and a steeper initial slope whereas on the other case tectonic
|
| /petsc/doc/community/meetings/2025/ |
| H A D | index.md | 214 …nite element method and finite volume method, as well as preconditioners, slope limiters, and more…
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexgeometry.c | 509 PetscReal slope = (x - x_j) * (y_i - y_j) - (x_i - x_j) * (y - y_j); in DMPlexLocatePoint_Quad_2D_Linear_Internal() local
510 if (slope == 0) { in DMPlexLocatePoint_Quad_2D_Linear_Internal()
515 if ((slope < 0) != (y_i < y_j)) ++crossings; in DMPlexLocatePoint_Quad_2D_Linear_Internal()
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | testcase3D.cas | 1849 (vof-to-dpm/lump-charact/calc-vof-slope-from-node-values? #t)
1962 (dpm/dpm-vof/close-packing/static-gpress-slope 60.)
1972 (dpm/dpm-vof/static-pileup/bott-wall-max-slope-angle-in-deg 60.)
1985 (dpm/dpm-vof/static-pileup/frict-drag/slope-at-onset 1000.)
2116 (dpm/film-splash-angle-slope 0.266)
2948 (melt/slope-limit-fac 1.)
2949 (melt/slope-cutoff-fac 0.01)
3167 (pdf/auto-refine/slope 0.25)
3196 (prepdf/strained-fla-auto-refine/slope 0.1)
3223 (prepdf/ckcfd-prmx-fla-auto-refine/slope 0.1)
[all …]
|
| /petsc/doc/tutorials/physics/ |
| H A D | guide_to_stokes.md | 187 … of refined meshes, calculates the error on each mesh, and determines the slope on a logarithmic s…
|
| /petsc/doc/ |
| H A D | petsc.bib | 20325 title = {Analysis of slope limiters on irregular grids},
|