| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex43.c | 89 PetscInt ei, ej, p, cum; in glvis_extract_eta() local
99 for (ej = sey; ej < sey + my; ej++) { 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
923 for (ej = sey; ej < sey + my; ej++) { 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
1010 for (ej = sey; ej < sey + my; ej++) { 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/snes/tutorials/ |
| H A D | ex16.c | 461 void FormElementJacobian(Field *ex, CoordField *ec, Field *ef, Field *eq, PetscScalar *ej, AppCtx *… in FormElementJacobian() argument
470 if (ej) in FormElementJacobian()
471 for (i = 0; i < NPB * NPB; i++) ej[i] = 0.; in FormElementJacobian()
523 if (ej) { in FormElementJacobian()
545 …ej[teidx + NPB * tridx] += scl * (lgrad[0] * dFS[3 * ll + 0] + lgrad[1] * dFS[3 * ll + 1] + lgrad[… in FormElementJacobian()
562 PetscInt ii, jj, kk, ll, ei, ej, ek, el; in ApplyBCsElement() local
569 for (ej = 0; ej < NB; ej++) { 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()
596 PetscScalar ej[NPB * NPB]; in FormJacobianLocal() local
[all …]
|
| /petsc/src/ts/utils/dmplexlandau/ |
| H A D | plexland.c | 1516 for (PetscInt ej = cStart, eidx = 0; ej < cEnd; ++ej, ++eidx, glb_elem_idx++) { in CreateStaticData() local
1525 …PetscCall(DMPlexGetClosureIndices(ctx->plex[grid], section[grid], globsection[grid], ej, PETSC_TRU… in CreateStaticData()
1590 …exRestoreClosureIndices(ctx->plex[grid], section[grid], globsection[grid], ej, PETSC_TRUE, &numind… in CreateStaticData()
1604 for (PetscInt ej = 0; ej < maps[grid].num_reduced; ++ej) { in CreateStaticData() local
1606 maps[grid].c_maps[ej][q].scale = pointMaps[ej][q].scale; in CreateStaticData()
1607 maps[grid].c_maps[ej][q].gid = pointMaps[ej][q].gid; in CreateStaticData()
1626 for (PetscInt ej = 0; ej < numCells[grid]; ++ej, glb_elem_idx++) { in CreateStaticData() local
1628 …const LandauIdx *const Idxs = &maps[grid].gIdx[ej][0][0]; // just use field-0 maps, Th… in CreateStaticData()
1653 for (PetscInt ej = 0; ej < numCells[grid]; ++ej, glb_elem_idx++) { in CreateStaticData() local
1657 const LandauIdx *const Idxs = &maps[grid].gIdx[ej][fieldA][0]; in CreateStaticData()
[all …]
|
| /petsc/src/ts/utils/dmplexlandau/kokkos/ |
| H A D | landau.kokkos.cxx | 804 for (PetscInt ej = cStart; ej < cEnd; ++ej) { in LandauKokkosJacobian() local
805 const PetscScalar *elMat = &h_elem_mats(b_id, grid, ej - cStart, 0); in LandauKokkosJacobian()
806 … PetscCall(DMPlexMatSetClosure(plex[grid], section, globalSection, B, ej, elMat, ADD_VALUES)); in LandauKokkosJacobian()
807 if (grid == 0 && ej == -1) { in LandauKokkosJacobian()
|