| /petsc/src/dm/impls/stag/tests/ |
| H A D | ex30.c | 66 PetscInt startx, starty, startz, N[3], nx, ny, nz, ex, ey, ez, d; in CreateMat() local 98 for (ey = starty; ey < starty + ny; ++ey) { in CreateMat() 105 row.j = ey; in CreateMat() 111 if (ey == N[1] - 1) { in CreateMat() 116 row.j = ey; in CreateMat() 127 row.j = ey; in CreateMat() 140 row.j = ey; in CreateMat() 152 row.j = ey; in CreateMat() 156 if (ey == 0) { in CreateMat() 160 col[0].j = ey; in CreateMat() [all …]
|
| H A D | ex8.c | 20 PetscInt startx, starty, startz, nx, ny, nz, ex, ey, ez, nExtrax, nExtray, nExtraz; in main() local 53 for (ey = starty; ey < starty + ny + nExtray; ++ey) { in main() 55 arrSol[ez][ey][ex][iux] = 2 * PRESSURE_CONST; in main() 56 arrRHS[ez][ey][ex][iux] = 0.0; in main() 57 arrSol[ez][ey][ex][iuy] = 2 * PRESSURE_CONST; in main() 58 arrRHS[ez][ey][ex][iuy] = 0.0; in main() 59 arrSol[ez][ey][ex][iuz] = 2 * PRESSURE_CONST; in main() 60 arrRHS[ez][ey][ex][iuz] = 0.0; in main() 61 if (ex < startx + nx && ey < starty + ny && ez < startz + nz) { in main() 62 arrSol[ez][ey][ex][ip] = PRESSURE_CONST; in main() [all …]
|
| H A D | ex11.c | 20 PetscInt startx, starty, nx, ny, ex, ey, nExtrax, nExtray; in main() local 51 for (ey = starty; ey < starty + ny + nExtray; ++ey) { in main() 53 arrSol[ey][ex][iux] = 2 * PRESSURE_CONST; in main() 54 arrRHS[ey][ex][iux] = 0.0; in main() 55 arrSol[ey][ex][iuy] = 2 * PRESSURE_CONST; in main() 56 arrRHS[ey][ex][iuy] = 0.0; in main() 57 if (ex < startx + nx && ey < starty + ny) { in main() 58 arrSol[ey][ex][ip] = PRESSURE_CONST; in main() 59 arrRHS[ey][ex][ip] = PRESSURE_CONST; in main() 120 …PetscInt startx, starty, nx, ny, nExtrax, nExtray, ex, ey, idxP, idxUx, idxUy, startGhostx, … in ApplyOperator() local [all …]
|
| H A D | ex23.c | 9 PetscInt ex, ey, ez, n[3], start[3], nExtra[3], iNext, iPrev, iCenter, d, round; in main() local 34 for (ey = start[1]; ey < start[1] + n[1]; ++ey) { in main() 35 cArrY[ey][iPrev] *= 1.1; in main() 36 cArrY[ey][iNext] = cArrY[ey][iPrev] + 0.1; in main() 37 cArrY[ey][iCenter] = 0.5 * (cArrY[ey][iPrev] + cArrY[ey][iNext]); in main() 51 for (ey = start[1]; ey < start[1] + n[1]; ++ey) { in main() 52 cArrY[ey][iPrev] = 2.0; in main() 53 cArrY[ey][iNext] = 2.0; in main() 54 cArrY[ey][iCenter] = 2.0; in main()
|
| H A D | ex15.c | 244 PetscInt ex, ey, startx, starty, nx, ny; in CreateSystem2d() local 258 for (ey = starty; ey < starty + ny; ++ey) { in CreateSystem2d() 266 row.j = ey; in CreateSystem2d() 273 if (ey == N[1] - 1) { in CreateSystem2d() 279 row.j = ey; in CreateSystem2d() 287 if (ey == 0) { in CreateSystem2d() 293 row.j = ey; in CreateSystem2d() 306 row.j = ey; in CreateSystem2d() 312 col[0].j = ey; in CreateSystem2d() 317 col[1].j = ey - 1; in CreateSystem2d() [all …]
|
| H A D | ex18.c | 252 PetscInt ex, ey, startx, starty, nx, ny; in CreateSystem() local 291 for (ey = starty; ey < starty + ny; ++ey) { in CreateSystem() 300 row.j = ey; in CreateSystem() 305 valRhs = bogusScale * uxRef(cArrX[ex][inext], cArrY[ey][icenter]); in CreateSystem() 309 if (ey == N[1] - 1) { in CreateSystem() 315 row.j = ey; in CreateSystem() 320 valRhs = uyRef(cArrX[ex][icenter], cArrY[ey][inext]); in CreateSystem() 325 if (ey == 0) { in CreateSystem() 331 row.j = ey; in CreateSystem() 336 valRhs = uyRef(cArrX[ex][icenter], cArrY[ey][iprev]); in CreateSystem() [all …]
|
| H A D | ex40.c | 358 for (PetscInt ey = start[1]; ey < start[1] + n[1] + n_extra[1]; ++ey) { in FormFunction2DNoCoupling() local 365 row.j = ey; in FormFunction2DNoCoupling() 378 row.j = ey; in FormFunction2DNoCoupling() 386 if (ey < N[1]) { in FormFunction2DNoCoupling() 392 row.j = ey; in FormFunction2DNoCoupling() 400 if (ex < N[0] && ey < N[1]) { in FormFunction2DNoCoupling() 406 row.j = ey; in FormFunction2DNoCoupling() 436 for (PetscInt ey = start[1]; ey < start[1] + n[1] + n_extra[1]; ++ey) { in FormJacobian2DNoCoupling() local 443 row.j = ey; in FormJacobian2DNoCoupling() 456 row.j = ey; in FormJacobian2DNoCoupling() [all …]
|
| /petsc/src/dm/impls/stag/tutorials/ |
| H A D | ex3.c | 170 PetscInt startx, starty, startz, N[3], nx, ny, nz, ex, ey, ez, d; in CreateSystem() local 204 for (ey = starty; ey < starty + ny; ++ey) { in CreateSystem() 212 row.j = ey; in CreateSystem() 217 …valRhs = uxRef(arrCoord[ez][ey][ex][icux_right[0]], arrCoord[ez][ey][ex][icux_right[1]], arrCoord[… in CreateSystem() 220 if (ey == N[1] - 1) { in CreateSystem() 226 row.j = ey; in CreateSystem() 231 …valRhs = uyRef(arrCoord[ez][ey][ex][icuy_up[0]], arrCoord[ez][ey][ex][icuy_up[1]], arrCoord[ez][ey… in CreateSystem() 240 row.j = ey; in CreateSystem() 245 …valRhs = uzRef(arrCoord[ez][ey][ex][icuz_front[0]], arrCoord[ez][ey][ex][icuz_front[1]], arrCoord[… in CreateSystem() 256 row.j = ey; in CreateSystem() [all …]
|
| H A D | ex6.c | 228 PetscInt N[3], ex, ey, ez, startx, starty, startz, nx, ny, nz, extrax, extray, extraz; in CreateLame() local 236 for (ey = starty; ey < starty + ny; ++ey) { in CreateLame() 241 pos.j = ey; in CreateLame() 246 pos.j = ey; in CreateLame() 253 for (ey = starty; ey < starty + ny + extray; ++ey) { in CreateLame() 258 pos.j = ey; in CreateLame() 267 for (ey = starty; ey < starty + ny; ++ey) { in CreateLame() 272 pos.j = ey; in CreateLame() 278 pos.j = ey; in CreateLame() 288 for (ey = starty; ey < starty + ny + extray; ++ey) { in CreateLame() [all …]
|
| H A D | ex2.c | 148 PetscInt ex, ey, startx, starty, nx, ny; in CreateSystem() local 180 …for (ey = starty; ey < starty + ny; ++ey) { /* With DMStag, always iterate x fastest, y second fas… in CreateSystem() 188 row.j = ey; in CreateSystem() 192 valRhs = uxRef(cArrX[ex][inext], cArrY[ey][icenter]); in CreateSystem() 195 if (ey == N[1] - 1) { in CreateSystem() 201 row.j = ey; in CreateSystem() 205 valRhs = uyRef(cArrX[ex][icenter], cArrY[ey][inext]); in CreateSystem() 209 if (ey == 0) { in CreateSystem() 215 row.j = ey; in CreateSystem() 219 valRhs = uyRef(cArrX[ex][icenter], cArrY[ey][iprev]); in CreateSystem() [all …]
|
| H A D | ex4.c | 576 PetscInt ex, ey, startx, starty, nx, ny; in CreateSystem2d() local 615 …for (ey = starty; ey < starty + ny; ++ey) { /* With DMStag, always iterate x fastest, y second fas… in CreateSystem2d() 619 const PetscBool bottom_boundary = (PetscBool)(ey == 0); in CreateSystem2d() 620 const PetscBool top_boundary = (PetscBool)(ey == N[1] - 1); in CreateSystem2d() 622 if (ey == N[1] - 1) { in CreateSystem2d() 629 row.j = ey; in CreateSystem2d() 636 if (ey == 0) { in CreateSystem2d() 643 row.j = ey; in CreateSystem2d() 660 row.j = ey; in CreateSystem2d() 666 rhoPoint[0].j = ey; in CreateSystem2d() [all …]
|
| /petsc/src/dm/impls/stag/ |
| H A D | stagda.c | 348 PetscInt ex, ey; in DMStagTransferCoordinatesToDMDA() local 356 for (ey = start[1]; ey < start[1] + n[1] + extraPoint[1]; ++ey) { in DMStagTransferCoordinatesToDMDA() 358 for (d = 0; d < 2; ++d) cArrDa[ey][ex][d] = cArrStag[ey][ex][slot + d]; in DMStagTransferCoordinatesToDMDA() 365 for (ey = start[1]; ey < start[1] + n[1] + extraPoint[1]; ++ey) { in DMStagTransferCoordinatesToDMDA() 367 cArrDa[ey][ex][0] = cArrX[ex][0]; in DMStagTransferCoordinatesToDMDA() 368 cArrDa[ey][ex][1] = cArrY[ey][0]; in DMStagTransferCoordinatesToDMDA() 375 PetscInt ex, ey, ez; in DMStagTransferCoordinatesToDMDA() local 384 for (ey = start[1]; ey < start[1] + n[1] + extraPoint[1]; ++ey) { in DMStagTransferCoordinatesToDMDA() 386 for (d = 0; d < 3; ++d) cArrDa[ez][ey][ex][d] = cArrStag[ez][ey][ex][slot + d]; in DMStagTransferCoordinatesToDMDA() 395 for (ey = start[1]; ey < start[1] + n[1] + extraPoint[1]; ++ey) { in DMStagTransferCoordinatesToDMDA() [all …]
|
| H A D | stag2d.c | 1345 for (PetscInt ey = start[1]; ey < start[1] + n[1] + n_extra[1]; ++ey) { in DMCreateMatrix_Stag_2D_AIJ_Assemble() local 1350 row_vertex[c].j = ey; in DMCreateMatrix_Stag_2D_AIJ_Assemble() 1357 row_face_down[c].j = ey; in DMCreateMatrix_Stag_2D_AIJ_Assemble() 1361 if (ey < N[1]) { in DMCreateMatrix_Stag_2D_AIJ_Assemble() 1364 row_face_left[c].j = ey; in DMCreateMatrix_Stag_2D_AIJ_Assemble() 1368 if (ex < N[0] && ey < N[1]) { in DMCreateMatrix_Stag_2D_AIJ_Assemble() 1371 row_element[c].j = ey; in DMCreateMatrix_Stag_2D_AIJ_Assemble() 1436 for (PetscInt ey = start[1]; ey < start[1] + n[1] + n_extra[1]; ++ey) { in DMCreateMatrix_Stag_2D_AIJ_Assemble() local 1440 row[i].j = ey; in DMCreateMatrix_Stag_2D_AIJ_Assemble() 1443 const PetscInt ey_offset = ey + offset_y; in DMCreateMatrix_Stag_2D_AIJ_Assemble()
|
| H A D | stag3d.c | 3463 for (PetscInt ey = start[1]; ey < start[1] + n[1] + n_extra[1]; ++ey) { in DMCreateMatrix_Stag_3D_AIJ_Assemble() local 3467 row_vertex[c].j = ey; in DMCreateMatrix_Stag_3D_AIJ_Assemble() 3475 row_edge_down_left[c].j = ey; in DMCreateMatrix_Stag_3D_AIJ_Assemble() 3481 if (ey < N[1]) { in DMCreateMatrix_Stag_3D_AIJ_Assemble() 3484 row_edge_back_left[c].j = ey; in DMCreateMatrix_Stag_3D_AIJ_Assemble() 3490 if (ey < N[0]) { in DMCreateMatrix_Stag_3D_AIJ_Assemble() 3493 row_edge_back_down[c].j = ey; in DMCreateMatrix_Stag_3D_AIJ_Assemble() 3499 if (ey < N[1] && ez < N[2]) { in DMCreateMatrix_Stag_3D_AIJ_Assemble() 3502 row_face_left[c].j = ey; in DMCreateMatrix_Stag_3D_AIJ_Assemble() 3511 row_face_down[c].j = ey; in DMCreateMatrix_Stag_3D_AIJ_Assemble() [all …]
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexcreate.c | 612 PetscInt e, ex, ey; in DMPlexCreateBoxSurfaceMesh_Tensor_2D_Internal() local 618 for (ey = 0; ey < edges[1]; ey++) { in DMPlexCreateBoxSurfaceMesh_Tensor_2D_Internal() 619 PetscInt edge = vx * edges[1] + ey + edges[0] * (edges[1] + 1); in DMPlexCreateBoxSurfaceMesh_Tensor_2D_Internal() 620 PetscInt vertex = ey * (edges[0] + 1) + vx + numEdges; in DMPlexCreateBoxSurfaceMesh_Tensor_2D_Internal() 629 if (ey == edges[1] - 1) { in DMPlexCreateBoxSurfaceMesh_Tensor_2D_Internal() 636 if (ey == edges[1] - 1) { in DMPlexCreateBoxSurfaceMesh_Tensor_2D_Internal() 1309 PetscInt c, e, ex, ey; in DMPlexCreateSquareMesh_Simplex_CrissCross() local 1317 for (ey = 0; ey < numYEdges; ++ey) { in DMPlexCreateSquareMesh_Simplex_CrissCross() 1320 const PetscInt eyp1 = (ey + 1) % numYVertices; in DMPlexCreateSquareMesh_Simplex_CrissCross() 1321 const PetscInt ez = firstZEdge + 4 * (ey * numXEdges + ex); in DMPlexCreateSquareMesh_Simplex_CrissCross() [all …]
|
| /petsc/doc/ |
| H A D | petsc.bib | 24220 S. Banerjee and S. Karmesin and S. Atlas and K. Keahey and W. F. Humphr ey}
|