Home
last modified time | relevance | path

Searched refs:da (Results 1 – 25 of 298) sorted by relevance

12345678910>>...12

/petsc/src/dm/tutorials/
H A Dex19.c12 DM da; in main() local
18 PetscCall(DMDACreate(PETSC_COMM_WORLD, &da)); in main()
19 PetscCall(DMSetDimension(da, 1)); in main()
20 PetscCall(DMDASetSizes(da, 8, 1, 1)); in main()
21 PetscCall(DMDASetBoundaryType(da, DM_BOUNDARY_NONE, DM_BOUNDARY_NONE, DM_BOUNDARY_NONE)); in main()
22 PetscCall(DMDASetDof(da, 1)); in main()
23 PetscCall(DMDASetStencilWidth(da, 1)); in main()
24 PetscCall(DMDASetOverlap(da, 1, 1, 1)); in main()
25 PetscCall(DMSetFromOptions(da)); in main()
26 PetscCall(DMSetOptionsPrefix(da, "n1d_")); in main()
[all …]
H A Dex12.c18 DM da; in main() local
29 …reate2d(PETSC_COMM_WORLD, bx, by, stype, M, N, PETSC_DECIDE, PETSC_DECIDE, 1, 1, NULL, NULL, &da)); in main()
30 PetscCall(DMSetFromOptions(da)); in main()
31 PetscCall(DMSetUp(da)); in main()
32 PetscCall(DMGetGlobalVector(da, &global1)); in main()
33 PetscCall(DMGetGlobalVector(da, &global2)); in main()
34 PetscCall(DMRestoreGlobalVector(da, &global1)); in main()
35 PetscCall(DMRestoreGlobalVector(da, &global2)); in main()
36 PetscCall(DMGetGlobalVector(da, &global1)); in main()
37 PetscCall(DMGetGlobalVector(da, &global3)); in main()
[all …]
/petsc/src/dm/impls/da/
H A Ddacreate.c3 static PetscErrorCode DMSetFromOptions_DA(DM da, PetscOptionItems PetscOptionsObject) in DMSetFromOptions_DA() argument
5 DM_DA *dd = (DM_DA *)da->data; in DMSetFromOptions_DA()
6 PetscInt refine = 0, dim = da->dim, maxnlevels = 100, refx[100], refy[100], refz[100], n, i; in DMSetFromOptions_DA()
11 …PetscCheck(!da->setupcalled, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_WRONGSTATE, "Cannot c… in DMSetFromOptions_DA()
12 …PetscCheck(dd->M >= 0, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "Dimension must… in DMSetFromOptions_DA()
13 …PetscCheck(dd->N >= 0, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "Dimension must… in DMSetFromOptions_DA()
14 …PetscCheck(dd->P >= 0, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "Dimension must… in DMSetFromOptions_DA()
22 if (flg) PetscCall(DMDASetOverlap(da, dd->xol, dd->xol, dd->xol)); in DMSetFromOptions_DA()
28 if (flg) PetscCall(DMDASetNumLocalSubDomains(da, dd->Nsub)); in DMSetFromOptions_DA()
39 if (flg) PetscCall(DMDASetBoundaryType(da, bt, bt, bt)); in DMSetFromOptions_DA()
[all …]
H A Dda.c21 PetscErrorCode DMDASetSizes(DM da, PetscInt M, PetscInt N, PetscInt P) in DMDASetSizes() argument
23 DM_DA *dd = (DM_DA *)da->data; in DMDASetSizes()
26 PetscValidHeaderSpecificType(da, DM_CLASSID, 1, DMDA); in DMDASetSizes()
27 PetscValidLogicalCollectiveInt(da, M, 2); in DMDASetSizes()
28 PetscValidLogicalCollectiveInt(da, N, 3); in DMDASetSizes()
29 PetscValidLogicalCollectiveInt(da, P, 4); in DMDASetSizes()
30 …PetscCheck(!da->setupcalled, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_WRONGSTATE, "This fun… in DMDASetSizes()
31 …PetscCheck(M >= 0, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_SIZ, "Number of grid points in … in DMDASetSizes()
32 …PetscCheck(N >= 0, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_SIZ, "Number of grid points in … in DMDASetSizes()
33 …PetscCheck(P >= 0, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_SIZ, "Number of grid points in … in DMDASetSizes()
[all …]
H A Ddagtol.c7 PetscErrorCode DMGlobalToLocalBegin_DA(DM da, Vec g, InsertMode mode, Vec l) in DMGlobalToLocalBegin_DA() argument
9 DM_DA *dd = (DM_DA *)da->data; in DMGlobalToLocalBegin_DA()
12 PetscValidHeaderSpecific(da, DM_CLASSID, 1); in DMGlobalToLocalBegin_DA()
19 PetscErrorCode DMGlobalToLocalEnd_DA(DM da, Vec g, InsertMode mode, Vec l) in DMGlobalToLocalEnd_DA() argument
21 DM_DA *dd = (DM_DA *)da->data; in DMGlobalToLocalEnd_DA()
24 PetscValidHeaderSpecific(da, DM_CLASSID, 1); in DMGlobalToLocalEnd_DA()
31 PetscErrorCode DMLocalToGlobalBegin_DA(DM da, Vec l, InsertMode mode, Vec g) in DMLocalToGlobalBegin_DA() argument
33 DM_DA *dd = (DM_DA *)da->data; in DMLocalToGlobalBegin_DA()
36 PetscValidHeaderSpecific(da, DM_CLASSID, 1); in DMLocalToGlobalBegin_DA()
42 …DM_BOUNDARY_NONE || dd->s <= 0 || dd->m != 1, PetscObjectComm((PetscObject)da), PETSC_ERR_SUP, "Av… in DMLocalToGlobalBegin_DA()
[all …]
H A Ddareg.c3 PetscErrorCode DMSetUp_DA(DM da) in DMSetUp_DA() argument
5 DM_DA *dd = (DM_DA *)da->data; in DMSetUp_DA()
8 PetscValidHeaderSpecific(da, DM_CLASSID, 1); in DMSetUp_DA()
9 …PetscCheck(dd->w >= 1, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "Must have 1 or… in DMSetUp_DA()
10 …PetscCheck(dd->s >= 0, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "Stencil width … in DMSetUp_DA()
13 PetscCall(PetscCalloc1(da->dim, &dd->coordinatename)); in DMSetUp_DA()
14 if (da->dim == 1) { in DMSetUp_DA()
15 PetscCall(DMSetUp_DA_1D(da)); in DMSetUp_DA()
16 } else if (da->dim == 2) { in DMSetUp_DA()
17 PetscCall(DMSetUp_DA_2D(da)); in DMSetUp_DA()
[all …]
H A Ddagetelem.c5 DM_DA *da = (DM_DA *)dm->data; in DMDAGetElements_1D() local
10 if (!da->e) { in DMDAGetElements_1D()
13 …PetscCheck(da->s, PetscObjectComm((PetscObject)dm), PETSC_ERR_SUP, "Cannot get elements for DMDA w… in DMDAGetElements_1D()
19 da->ne = 1 * (xe - xs - 1); in DMDAGetElements_1D()
20 PetscCall(PetscMalloc1(1 + 2 * da->ne, &da->e)); in DMDAGetElements_1D()
22 da->e[cnt++] = (i - Xs); in DMDAGetElements_1D()
23 da->e[cnt++] = (i - Xs + 1); in DMDAGetElements_1D()
25 da->nen = 2; in DMDAGetElements_1D()
29 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, 2, corners, PETSC_COPY_VALUES, &da->ecorners)); in DMDAGetElements_1D()
31 *nel = da->ne; in DMDAGetElements_1D()
[all …]
H A Ddasub.c33 PetscErrorCode DMDAGetLogicalCoordinate(DM da, PetscScalar x, PetscScalar y, PetscScalar z, PetscIn… in DMDAGetLogicalCoordinate() argument
43 …PetscCheck(da->dim != 1, PetscObjectComm((PetscObject)da), PETSC_ERR_SUP, "Cannot get point from 1… in DMDAGetLogicalCoordinate()
44 …PetscCheck(da->dim != 3, PetscObjectComm((PetscObject)da), PETSC_ERR_SUP, "Cannot get point from 3… in DMDAGetLogicalCoordinate()
49 PetscCall(DMGetCoordinateDM(da, &dacoors)); in DMDAGetLogicalCoordinate()
51 PetscCall(DMGetCoordinates(da, &coors)); in DMDAGetLogicalCoordinate()
63 PetscCallMPI(MPIU_Allreduce(&D, &Dv, 1, MPIU_REAL, MPIU_MIN, PetscObjectComm((PetscObject)da))); in DMDAGetLogicalCoordinate()
71 PetscCallMPI(MPI_Comm_rank(PetscObjectComm((PetscObject)da), &rank)); in DMDAGetLogicalCoordinate()
74 PetscCallMPI(MPIU_Allreduce(&rank, &root, 1, MPI_INT, MPI_SUM, PetscObjectComm((PetscObject)da))); in DMDAGetLogicalCoordinate()
76 PetscCallMPI(MPI_Bcast(X, 1, MPIU_SCALAR, root, PetscObjectComm((PetscObject)da))); in DMDAGetLogicalCoordinate()
77 PetscCallMPI(MPI_Bcast(Y, 1, MPIU_SCALAR, root, PetscObjectComm((PetscObject)da))); in DMDAGetLogicalCoordinate()
[all …]
/petsc/src/binding/petsc4py/test/
H A Dtest_dmda.py16 self.da = PETSc.DMDA().create(
27 self.da = None
31 dim = self.da.getDim()
32 dof = self.da.getDof()
33 sizes = self.da.getSizes()
34 boundary = self.da.getBoundaryType()
35 stencil_type = self.da.getStencilType()
36 stencil_width = self.da.getStencilWidth()
45 dim = self.da.getDim()
46 ranges = self.da.getRanges()
[all …]
H A Dtest_dmstag.py16 self.da = PETSc.DMStag().create(
42 self.da = None
47 self.da.setCoordinateDMType('stag')
48 self.da.setUniformCoordinates(0, 1, 0, 1, 0, 1)
49 self.da.setUniformCoordinatesExplicit(0, 1, 0, 1, 0, 1)
50 cda = self.da.getCoordinateDM()
55 c = self.da.getCoordinatesLocal()
56 self.da.setCoordinatesLocal(c)
57 gc = self.da.getCoordinatesLocal()
61 c = self.da.getCoordinates()
[all …]
/petsc/src/dm/impls/da/kokkos/
H A Ddagetov.kokkos.cxx8 #define DMDA_VEC_GET_SHAPE(da, vec, xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, d… argument
10 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); \
11 PetscCall(DMDAGetGhostCorners(da, &gxs, &gys, &gzs, &gxm, &gym, &gzm)); \
12 …PetscCall(DMDAGetInfo(da, &dim, NULL, NULL, NULL, NULL, NULL, NULL, &dof, NULL, NULL, NULL, NULL, …
27 PetscErrorCode DMDAVecGetKokkosOffsetView_Private(DM da, Vec vec, PetscScalarKokkosOffsetView1DType… in DMDAVecGetKokkosOffsetView_Private() argument
33 PetscValidHeaderSpecificType(da, DM_CLASSID, 1, DMDA); in DMDAVecGetKokkosOffsetView_Private()
36 DMDA_VEC_GET_SHAPE(da, vec, xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof); in DMDAVecGetKokkosOffsetView_Private()
46 PetscErrorCode DMDAVecRestoreKokkosOffsetView_Private(DM da, Vec vec, PetscScalarKokkosOffsetView1D… in DMDAVecRestoreKokkosOffsetView_Private() argument
51 PetscValidHeaderSpecificType(da, DM_CLASSID, 1, DMDA); in DMDAVecRestoreKokkosOffsetView_Private()
61 PetscErrorCode DMDAVecGetKokkosOffsetView(DM da, Vec vec, ConstPetscScalarKokkosOffsetView1DType<Me… in DMDAVecGetKokkosOffsetView() argument
[all …]
/petsc/src/ts/tutorials/advection-diffusion-reaction/
H A Dex6.c40 DM da; in main() local
53 PetscCall(DMDACreate1d(PETSC_COMM_WORLD, DM_BOUNDARY_PERIODIC, 60, 1, 1, NULL, &da)); in main()
54 PetscCall(DMSetFromOptions(da)); in main()
55 PetscCall(DMSetUp(da)); in main()
58 PetscCall(DMCreateGlobalVector(da, &U)); in main()
62 PetscCall(TSSetDM(ts, da)); in main()
75 PetscCall(DMDAGetInfo(da, PETSC_IGNORE, &M, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)); in main()
96 PetscCall(DMDestroy(&da)); in main()
116 DM da; in InitialConditions() local
120 PetscCall(TSGetDM(ts, &da)); in InitialConditions()
[all …]
/petsc/src/binding/petsc4py/demo/legacy/wrap-f2py/
H A Drun_demo.py9 def __init__(self, da, lambda_): argument
10 assert da.getDim() == 2
11 self.da = da
16 da = self.da.fortran
18 ierr = Bratu2D.FormInitGuess(da, vec_X, self.lambda_)
23 da = self.da.fortran
26 ierr = Bratu2D.FormFunction(da, vec_X, vec_F, self.lambda_)
31 da = self.da.fortran
34 ierr = Bratu2D.FormJacobian(da, vec_X, mat_P, self.lambda_)
46 da = PETSc.DA().create([N, N], stencil_width=1) variable
[all …]
/petsc/src/dm/ftn-mod/
H A Dpetscdmda.h902 subroutine DMDARestoreNeighbors(da, ranks, ierr)
5 DM da
10 subroutine DMDARestoreOwnershipRanges(da, lx, ly, lz, ierr)
13 DM da
25 subroutine DMDAVecGetArray1(da, vec, array, ierr)
28 DM da
33 subroutine DMDAVecGetArray2(da, vec, array, ierr)
36 DM da
41 subroutine DMDAVecGetArray3(da, vec, array, ierr)
44 DM da
[all …]
/petsc/src/dm/impls/stag/tests/
H A Dex6.c22 DM da; in main() local
24 PetscCall(DMStagVecSplitToDMDA(dm, vec, DMSTAG_ELEMENT, -3, &da, &vecda)); in main()
25 PetscCall(DMDestroy(&da)); in main()
31 DM da; in main() local
33 PetscCall(DMStagVecSplitToDMDA(dm, vec, DMSTAG_ELEMENT, -5, &da, &vecda)); in main()
34 PetscCall(DMDestroy(&da)); in main()
40 DM da; in main() local
42 PetscCall(DMStagVecSplitToDMDA(dm, vec, DMSTAG_ELEMENT, -2, &da, &vecda)); in main()
43 PetscCall(DMDestroy(&da)); in main()
49 DM da; in main() local
[all …]
/petsc/src/dm/tests/
H A Dex47.c15 DM da; in test_3d() local
22 …_STENCIL_STAR, M, N, P, PETSC_DECIDE, PETSC_DECIDE, PETSC_DECIDE, dof, sw, NULL, NULL, NULL, &da)); in test_3d()
23 PetscCall(DMSetFromOptions(da)); in test_3d()
24 PetscCall(DMSetUp(da)); in test_3d()
26 PetscCall(DMDASetUniformCoordinates(da, 0.0, Lx, 0.0, Ly, 0.0, Lz)); in test_3d()
27 PetscCall(DMDAGetLocalInfo(da, &info)); in test_3d()
28 PetscCall(DMCreateGlobalVector(da, &v)); in test_3d()
29 PetscCall(DMDAVecGetArray(da, v, &va)); in test_3d()
40 PetscCall(DMDAVecRestoreArray(da, v, &va)); in test_3d()
45 PetscCall(DMDestroy(&da)); in test_3d()
[all …]
H A Dex42.c17 DM da; in test_3d() local
24 …_STENCIL_STAR, M, N, P, PETSC_DECIDE, PETSC_DECIDE, PETSC_DECIDE, dof, sw, NULL, NULL, NULL, &da)); in test_3d()
25 PetscCall(DMSetFromOptions(da)); in test_3d()
26 PetscCall(DMSetUp(da)); in test_3d()
28 PetscCall(DMDASetUniformCoordinates(da, 0.0, Lx, 0.0, Ly, 0.0, Lz)); in test_3d()
29 PetscCall(DMDAGetLocalInfo(da, &info)); in test_3d()
30 PetscCall(DMCreateGlobalVector(da, &v)); in test_3d()
31 PetscCall(DMDAVecGetArray(da, v, &va)); in test_3d()
42 PetscCall(DMDAVecRestoreArray(da, v, &va)); in test_3d()
47 PetscCall(DMDestroy(&da)); in test_3d()
[all …]
H A Dex48.c8 PetscErrorCode NameFields(DM da, PetscInt dof) in NameFields() argument
16 PetscCall(DMDASetFieldName(da, c, fieldname)); in NameFields()
29 DM da; in test_3d() local
36 …_STENCIL_STAR, M, N, P, PETSC_DECIDE, PETSC_DECIDE, PETSC_DECIDE, dof, sw, NULL, NULL, NULL, &da)); in test_3d()
37 PetscCall(DMSetFromOptions(da)); in test_3d()
38 PetscCall(DMSetUp(da)); in test_3d()
39 if (namefields) PetscCall(NameFields(da, dof)); in test_3d()
41 PetscCall(DMDASetUniformCoordinates(da, 0.0, Lx, 0.0, Ly, 0.0, Lz)); in test_3d()
42 PetscCall(DMDAGetLocalInfo(da, &info)); in test_3d()
43 PetscCall(DMCreateGlobalVector(da, &v)); in test_3d()
[all …]
H A Dex43.c16 PetscErrorCode PrintVecWithGhosts(DM da, Vec v) in PrintVecWithGhosts() argument
24 com = PetscObjectComm((PetscObject)da); in PrintVecWithGhosts()
26 PetscCall(DMDAGetLocalInfo(da, &info)); in PrintVecWithGhosts()
28 PetscCall(DMDAVecGetArray(da, v, &p)); in PrintVecWithGhosts()
33 PetscCall(DMDAVecRestoreArray(da, v, &p)); in PrintVecWithGhosts()
40 PetscErrorCode VecSetOwned(DM da, Vec v, PetscScalar value) in VecSetOwned() argument
45 PetscCall(DMDAGetCorners(da, &xs, &ys, 0, &xm, &ym, 0)); in VecSetOwned()
46 PetscCall(DMDAVecGetArray(da, v, &p)); in VecSetOwned()
50 PetscCall(DMDAVecRestoreArray(da, v, &p)); in VecSetOwned()
57 DM da; in main() local
[all …]
H A Dex10k.kokkos.cxx29 DM da; in main() local
49 …reate2d(PETSC_COMM_WORLD, bx, by, st, M, N, PETSC_DECIDE, PETSC_DECIDE, dof, sw, NULL, NULL, &da)); in main()
50 PetscCall(DMSetFromOptions(da)); in main()
51 PetscCall(DMSetUp(da)); in main()
52 PetscCall(DMDAGetCorners(da, &xs, &ys, 0, &xm, &ym, 0)); in main()
53 PetscCall(DMCreateGlobalVector(da, &g)); in main()
54 PetscCall(DMCreateLocalVector(da, &l)); in main()
55 PetscCall(DMCreateGlobalVector(da, &gg)); in main()
56 PetscCall(DMCreateLocalVector(da, &ll)); in main()
59 PetscCall(DMDAVecGetArrayDOFWrite(da, g, &garray)); in main()
[all …]
/petsc/src/binding/petsc4py/demo/legacy/poisson2d/
H A Dpoisson2d.py30 def __init__(self, da): argument
31 assert da.getDim() == 2
32 self.da = da
33 self.localX = da.createLocalVec()
36 b = self.da.getVecArray(B)
37 mx, my = self.da.getSizes()
39 (xs, xe), (ys, ye) = self.da.getRanges()
46 self.da.globalToLocal(X, self.localX)
47 x = self.da.getVecArray(self.localX)
48 y = self.da.getVecArray(Y)
[all …]
/petsc/src/snes/tutorials/
H A Dex3k.kokkos.cxx10 DM da; /* distributed array */ member
50 DM da = user->da; in CpuFunction() local
56 PetscCall(DMGetLocalVector(da, &xl)); in CpuFunction()
57 PetscCall(DMGlobalToLocal(da, x, INSERT_VALUES, xl)); in CpuFunction()
58 PetscCall(DMDAVecGetArray(da, xl, &X)); in CpuFunction()
59 PetscCall(DMDAVecGetArray(da, r, &R)); in CpuFunction()
60 PetscCall(DMDAVecGetArray(da, user->F, &F)); in CpuFunction()
62 PetscCall(DMDAGetCorners(da, &xs, NULL, NULL, &xm, NULL, NULL)); in CpuFunction()
63 …PetscCall(DMDAGetInfo(da, NULL, &M, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NU… in CpuFunction()
77 PetscCall(DMDAVecRestoreArray(da, xl, &X)); in CpuFunction()
[all …]
H A Dex78.c30 DM da; in main() local
52 PetscCall(DMDACreate1d(PETSC_COMM_WORLD, DM_BOUNDARY_PERIODIC, N, 1, 1, NULL, &da)); in main()
53 PetscCall(DMSetFromOptions(da)); in main()
54 PetscCall(DMSetUp(da)); in main()
60 PetscCall(DMCreateGlobalVector(da, &x)); in main()
71 PetscCall(SNESSetFunction(snes, r, FormFunction, da)); in main()
76 PetscCall(DMCreateMatrix(da, &J)); in main()
79 PetscCall(SNESSetJacobian(snes, J, J, FormJacobian, da)); in main()
89 PetscCall(DMDestroy(&da)); in main()
111 DM da = (DM)ctx; in FormFunction() local
[all …]
/petsc/src/ts/tutorials/phasefield/
H A Dbiharmonic2.c49 DM da; in main() local
82 PetscCall(DMDACreate1d(PETSC_COMM_WORLD, DM_BOUNDARY_PERIODIC, 10, 2, 2, NULL, &da)); in main()
83 PetscCall(DMSetFromOptions(da)); in main()
84 PetscCall(DMSetUp(da)); in main()
85 PetscCall(DMDASetFieldName(da, 0, "Biharmonic heat equation: w = -kappa*u_xx")); in main()
86 PetscCall(DMDASetFieldName(da, 1, "Biharmonic heat equation: u")); in main()
87 PetscCall(DMDAGetInfo(da, 0, &Mx, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)); in main()
94 PetscCall(DMCreateGlobalVector(da, &x)); in main()
101 PetscCall(TSSetDM(ts, da)); in main()
120 PetscCall(DMCreateColoring(da, IS_COLORING_GLOBAL, &iscoloring)); in main()
[all …]
/petsc/src/binding/petsc4py/demo/legacy/wrap-swig/
H A Drun_demo.py9 def __init__(self, da, lambda_): argument
10 assert da.getDim() == 3
11 self.da = da
17 Bratu3D.FormInitGuess(self.da, X, self.params)
21 Bratu3D.FormFunction(self.da, X, F, self.params)
25 Bratu3D.FormJacobian(self.da, X, P, self.params)
36 da = PETSc.DMDA().create([N, N, N], stencil_width=1) variable
37 app = App(da, lambda_)
40 F = da.createGlobalVec()
42 J = da.createMat()
[all …]

12345678910>>...12