| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex43.c | 89 PetscInt ei, ej, p, cum; in glvis_extract_eta() local
100 for (ei = sex; ei < sex + mx; ei++) { in glvis_extract_eta()
101 for (p = 0; p < GAUSS_POINTS; p++) array[cum++] = props[ej][ei].eta[p]; in glvis_extract_eta()
688 static PetscErrorCode GetElementCoords(DMDACoor2d **_coords, PetscInt ei, PetscInt ej, PetscScalar … in GetElementCoords() argument
692 el_coords[NSD * 0] = _coords[ej][ei].x; in GetElementCoords()
693 el_coords[NSD * 0 + 1] = _coords[ej][ei].y; in GetElementCoords()
694 el_coords[NSD * 1] = _coords[ej + 1][ei].x; in GetElementCoords()
695 el_coords[NSD * 1 + 1] = _coords[ej + 1][ei].y; in GetElementCoords()
696 el_coords[NSD * 2] = _coords[ej + 1][ei + 1].x; in GetElementCoords()
697 el_coords[NSD * 2 + 1] = _coords[ej + 1][ei + 1].y; in GetElementCoords()
[all …]
|
| H A D | ex42.c | 904 PetscInt ei, ej, ek; in AssembleA_Stokes() local
924 for (ei = sex; ei < sex + mx; ei++) { in AssembleA_Stokes()
926 PetscCall(GetElementCoords3D(_coords, ei, ej, ek, el_coords)); in AssembleA_Stokes()
928 PetscCall(CellPropertiesGetCell(cell_properties, ei, ej, ek, &props)); in AssembleA_Stokes()
947 PetscCall(DMDAGetElementEqnums3D_up(u_eqn, p_eqn, ei, ej, ek)); in AssembleA_Stokes()
990 PetscInt ei, ej, ek; in AssembleA_PCStokes() local
1011 for (ei = sex; ei < sex + mx; ei++) { in AssembleA_PCStokes()
1013 PetscCall(GetElementCoords3D(_coords, ei, ej, ek, el_coords)); in AssembleA_PCStokes()
1015 PetscCall(CellPropertiesGetCell(cell_properties, ei, ej, ek, &props)); in AssembleA_PCStokes()
1032 PetscCall(DMDAGetElementEqnums3D_up(u_eqn, p_eqn, ei, ej, ek)); in AssembleA_PCStokes()
[all …]
|
| H A D | ex49.c | 551 static PetscErrorCode GetElementCoords(DMDACoor2d **_coords, PetscInt ei, PetscInt ej, PetscScalar … in GetElementCoords() argument
555 el_coords[NSD * 0 + 0] = _coords[ej][ei].x; in GetElementCoords()
556 el_coords[NSD * 0 + 1] = _coords[ej][ei].y; in GetElementCoords()
557 el_coords[NSD * 1 + 0] = _coords[ej + 1][ei].x; in GetElementCoords()
558 el_coords[NSD * 1 + 1] = _coords[ej + 1][ei].y; in GetElementCoords()
559 el_coords[NSD * 2 + 0] = _coords[ej + 1][ei + 1].x; in GetElementCoords()
560 el_coords[NSD * 2 + 1] = _coords[ej + 1][ei + 1].y; in GetElementCoords()
561 el_coords[NSD * 3 + 0] = _coords[ej][ei + 1].x; in GetElementCoords()
562 el_coords[NSD * 3 + 1] = _coords[ej][ei + 1].y; in GetElementCoords()
573 PetscInt ei, ej; in AssembleA_Elasticity() local
[all …]
|
| /petsc/src/tao/util/ |
| H A D | tao_util.c | 252 PetscReal ai, bi, ci, di, ei; in MatDFischer() local
345 ei = Fischer(PetscRealPart(u[i]) - PetscRealPart(x[i]), -PetscRealPart(f[i])); in MatDFischer()
346 ai = fischnorm(PetscRealPart(x[i]) - PetscRealPart(l[i]), ei); in MatDFischer()
349 bi = ei / ai - 1.0; in MatDFischer()
398 PetscReal ai, bi, ci, di, ei, fi; in MatDSFischer() local
448 ei = SFischer(PetscRealPart(u[i]) - PetscRealPart(x[i]), -PetscRealPart(f[i]), mu); in MatDSFischer()
449 ai = fischsnorm(PetscRealPart(x[i]) - PetscRealPart(l[i]), ei, mu); in MatDSFischer()
452 bi = ei / ai - 1.0; in MatDSFischer()
453 ei = 2.0 * mu / ei; in MatDSFischer()
458 dm[i] = ei + bi * fi; in MatDSFischer()
|
| /petsc/src/sys/logging/handler/impls/perfstubs/ |
| H A D | logperfstubs.c | 74 PetscLogEvent ei; in PetscLogHandlerPSUpdateEvents() local
76 PetscCall(PetscMPIIntCast(i, &ei)); in PetscLogHandlerPSUpdateEvents()
77 PetscCall(PetscLogStateEventGetInfo(state, ei, &event_info)); in PetscLogHandlerPSUpdateEvents()
|
| /petsc/src/snes/tutorials/ |
| H A D | ex16.c | 562 PetscInt ii, jj, kk, ll, ei, ej, ek, el; in ApplyBCsElement() local
570 for (ei = 0; ei < NB; ei++) { in ApplyBCsElement()
572 …if (OnBoundary(i + ii, j + jj, k + kk, mx, my, mz) || OnBoundary(i + ei, j + ej, k + ek, mx, my, m… in ApplyBCsElement()
573 PetscInt teidx = el + 3 * (ei + ej * NB + ek * NB * NB); in ApplyBCsElement()
|
| /petsc/src/ts/utils/dmplexlandau/ |
| H A D | plexland.c | 165 PetscInt cStart, cEnd, ei; in LandauFormJacobian_Internal() local
170 for (ei = cStart; ei < cEnd; ++ei) { in LandauFormJacobian_Internal()
172 PetscCall(DMPlexVecGetClosure(ctx->plex[grid], section[grid], locX2, ei, NULL, &coef)); in LandauFormJacobian_Internal()
174 … PetscCall(DMPlexVecRestoreClosure(ctx->plex[grid], section[grid], locX2, ei, NULL, &coef)); in LandauFormJacobian_Internal()
|