Home
last modified time | relevance | path

Searched refs:ys (Results 1 – 25 of 120) sorted by relevance

12345

/petsc/src/tao/unconstrained/tutorials/
H A Dminsurf2.c166 PetscInt xs, xm, gxs, gxm, ys, ym, gys, gym; in FormFunction() local
177 PetscCall(DMDAGetCorners(user->dm, &xs, &ys, NULL, &xm, &ym, NULL)); in FormFunction()
188 for (j = ys; j < ys + ym; j++) { in FormFunction()
193 xl = user->left[j - ys + 1]; in FormFunction()
205 xr = user->right[j - ys + 1]; in FormFunction()
235 for (j = ys; j < ys + ym; j++) { in FormFunction()
236 d3 = (user->left[j - ys + 1] - user->left[j - ys + 2]) * rhy; in FormFunction()
237 d2 = (user->left[j - ys + 1] - x[j][0]) * rhx; in FormFunction()
241 if (ys == 0) { /* bottom side */ in FormFunction()
249 for (j = ys; j < ys + ym; j++) { in FormFunction()
[all …]
H A Deptorsion2.c159 PetscInt xs, ys, xm, ym, gxm, gym, gxs, gys, xe, ye; in FormInitialGuess() local
164 PetscCall(DMDAGetCorners(user->dm, &xs, &ys, NULL, &xm, &ym, NULL)); in FormInitialGuess()
169 ye = ys + ym; in FormInitialGuess()
170 for (j = ys; j < ye; j++) { /* for (j=0; j<my; j++) */ in FormInitialGuess()
201 PetscInt xs, ys, xm, ym, gxm, gym, gxs, gys; in FormFunctionGradient() local
228 PetscCall(DMDAGetCorners(user->dm, &xs, &ys, NULL, &xm, &ym, NULL)); in FormFunctionGradient()
233 ye = ys + ym; in FormFunctionGradient()
236 if (ys == 0) ysm = ys - 1; in FormFunctionGradient()
237 else ysm = ys; in FormFunctionGradient()
276 for (j = ys; j < yep; j++) { /* for (j=0; j<=my; j++) */ in FormFunctionGradient()
[all …]
H A Deptorsion2f.F9057 PetscInt xs, ys, xm, ym
66 PetscCall(DMDAGetCorners(dm, xs, ys, PETSC_NULL_INTEGER, xm, ym, PETSC_NULL_INTEGER, ierr))
71 ye = ys + ym
72 do j = ys, ye - 1
123 PetscInt xs, ys, xm, ym
145 PetscCall(DMDAGetCorners(dm, xs, ys, PETSC_NULL_INTEGER, xm, ym, PETSC_NULL_INTEGER, ierr))
153 ye = ys + ym
159 if (ys == 0) then
160 ysm = ys - 1
162 ysm = ys
[all …]
/petsc/src/dm/impls/da/
H A Ddagetelem.c41 PetscInt j, ys, ye, Ys, Ye; in DMDAGetElements_2D() local
65 PetscCall(DMDAGetCorners(dm, &xs, &ys, NULL, &xe, &ye, NULL)); in DMDAGetElements_2D()
70 ye += ys; in DMDAGetElements_2D()
72 if (ys != Ys) ys -= 1; in DMDAGetElements_2D()
73 da->ne = ns * (xe - xs - 1) * (ye - ys - 1); in DMDAGetElements_2D()
75 for (j = ys; j < ye - 1; j++) { in DMDAGetElements_2D()
90 corners[0] = (xs - Xs) + (ys - Ys) * (Xe - Xs); in DMDAGetElements_2D()
91 corners[1] = (xe - 1 - Xs) + (ys - Ys) * (Xe - Xs); in DMDAGetElements_2D()
106 PetscInt j, ys, ye, Ys, Ye; in DMDAGetElements_3D() local
131 PetscCall(DMDAGetCorners(dm, &xs, &ys, &zs, &xe, &ye, &ze)); in DMDAGetElements_3D()
[all …]
H A Ddagetarray.c45 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in DMDAVecGetArray() local
51 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in DMDAVecGetArray()
62 gys = ys; in DMDAVecGetArray()
94 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in DMDAVecRestoreArray() local
100 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in DMDAVecRestoreArray()
111 gys = ys; in DMDAVecRestoreArray()
169 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in DMDAVecGetArrayWrite() local
179 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in DMDAVecGetArrayWrite()
190 gys = ys; in DMDAVecGetArrayWrite()
221 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in DMDAVecRestoreArrayWrite() local
[all …]
H A Ddageometry.c42 PetscInt ys, ye, Ys, Ye; in DMGetLocalBoundingBox_DA() local
52 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xe, &ye, &ze)); in DMGetLocalBoundingBox_DA()
58 ye += ys; in DMGetLocalBoundingBox_DA()
64 if (ys != Ys && Ys >= 0) ys -= 1; in DMGetLocalBoundingBox_DA()
83 c0 = (xs - Xs) + (ys - Ys) * (Xe - Xs); in DMGetLocalBoundingBox_DA()
87 c0 = (xs - Xs) + (ys - Ys) * (Xe - Xs) + (zs - Zs) * (Xe - Xs) * (Ye - Ys); in DMGetLocalBoundingBox_DA()
119 if (dim > 1) cs[1] = ys; in DMGetLocalBoundingBox_DA()
135 PetscInt ys, ye, mylocal; in private_DMDALocatePointsIS_2D_Regular() local
144 ys = cs[1]; in private_DMDALocatePointsIS_2D_Regular()
150 mylocal = ye - ys - 1; in private_DMDALocatePointsIS_2D_Regular()
[all …]
/petsc/src/ksp/pc/impls/kaczmarz/
H A Dkaczmarz.c18 PetscInt xs, xe, ys, ye, ncols, i, j; in PCApply_Kaczmarz() local
28 PetscCall(MatGetOwnershipRangeColumn(pc->pmat, &ys, &ye)); in PCApply_Kaczmarz()
38 if (cols[j] >= ys && cols[j] < ye) r -= yarray[cols[j] - ys] * vals[j]; in PCApply_Kaczmarz()
43 if (cols[j] >= ys && cols[j] < ye) yarray[cols[j] - ys] += vals[j] * lambda * r / anrm; in PCApply_Kaczmarz()
54 if (cols[j] >= ys && cols[j] < ye) r -= yarray[cols[j] - ys] * vals[j]; in PCApply_Kaczmarz()
59 if (cols[j] >= ys && cols[j] < ye) yarray[cols[j] - ys] += vals[j] * lambda * r / anrm; in PCApply_Kaczmarz()
/petsc/src/tao/bound/tutorials/
H A Dplate2f.F9068 PetscInt ys, ym
93 PetscCall(DMDAGetCorners(dm, xs, ys, PETSC_NULL_INTEGER, xm, ym, PETSC_NULL_INTEGER, ierr))
112 do j = ys, ys + ym - 1
124 xl = left_v(1 + j - ys + 1)
125 xlt = left_v(1 + j - ys + 2)
138 xr = right_v(1 + j - ys + 1)
139 xrb = right_v(1 + j - ys)
206 do j = ys, ys + ym - 1
207 d3 = (left_v(1 + j - ys + 1) - left_v(1 + j - ys + 2))*rhy
208 d2 = (left_v(1 + j - ys + 1) - x_v(1 + (j - gys)*gxm))*rhx
[all …]
H A Dplate2.c202 PetscInt xs, xm, gxs, gxm, ys, ym, gys, gym; in FormFunctionGradient() local
214 PetscCall(DMDAGetCorners(user->dm, &xs, &ys, NULL, &xm, &ym, NULL)); in FormFunctionGradient()
233 for (j = ys; j < ys + ym; j++) { in FormFunctionGradient()
241 xl = left[j - ys + 1]; in FormFunctionGradient()
242 xlt = left[j - ys + 2]; in FormFunctionGradient()
255 xr = right[j - ys + 1]; in FormFunctionGradient()
256 xrb = right[j - ys]; in FormFunctionGradient()
318 for (j = ys; j < ys + ym; j++) { in FormFunctionGradient()
319 d3 = (left[j - ys + 1] - left[j - ys + 2]) * rhy; in FormFunctionGradient()
320 d2 = (left[j - ys + 1] - x[(j - gys) * gxm]) * rhx; in FormFunctionGradient()
[all …]
/petsc/src/dm/tests/
H A Dex10k.kokkos.cxx30 PetscInt M = 5, N = 7, xm, ym, xs, ys; in main() local
52 PetscCall(DMDAGetCorners(da, &xs, &ys, 0, &xm, &ym, 0)); in main()
60 for (PetscInt j = ys; j < ys + ym; j++) { /* run on host */ in main()
70 …"init 1", MDRangePolicy<Rank<3, Iterate::Right, Iterate::Right>>({ys, xs, 0}, {ys + ym, xs + xm, d… in main()
82 for (PetscInt j = ys; j < ys + ym; j++) { in main()
94 …"stencil 1", MDRangePolicy<Rank<3, Iterate::Right, Iterate::Right>>({ys, xs, 0}, {ys + ym, xs + xm… in main()
125 PetscCall(DMDAGetCorners(da, &xs, &ys, 0, &xm, &ym, 0)); in main()
133 for (PetscInt j = ys; j < ys + ym; j++) { in main()
146 …"init 2", MDRangePolicy<Rank<2, Iterate::Right, Iterate::Right>>({ys, xs}, {ys + ym, xs + xm}), KO… in main()
160 for (PetscInt j = ys; j < ys + ym; j++) { in main()
[all …]
H A Dex2k.kokkos.cxx36 PetscInt xm, ym, zm, xs, ys, zs, gxm, gym, gzm, gxs, gys, gzs; in Update1() local
39 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in Update1()
44 for (j = ys; j < ys + ym; j++) { in Update1()
59 PetscInt xm, ym, zm, xs, ys, zs, gxm, gym, gzm, gxs, gys, gzs; in Update2() local
62 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in Update2()
65 #define Y2(k, j, i) y2[(k - zs) * ym * xm + (j - ys) * xm + (i - xs)] in Update2()
69 for (j = ys; j < ys + ym; j++) { in Update2()
84 PetscInt xm, ym, zm, xs, ys, zs, gxm, gym, gzm, gxs, gys, gzs; in main() local
111 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in main()
146 …HostExecutionSpace, Rank<3, Iterate::Right, Iterate::Right>>({zs, ys, xs}, {zs + zm, ys + ym, xs +… in main()
/petsc/src/dm/tutorials/
H A Dex11f90.F9014 PetscInt ys, yl
53 PetscCallA(DMDAGetCorners(ada, xs, ys, PETSC_NULL_INTEGER, xl, yl, PETSC_NULL_INTEGER, ierr))
56 do j = ys, ys + yl - 1
76 PetscCallA(DMDAGetCorners(ada, xs, ys, zs, xl, yl, zl, ierr))
79 do j = ys, ys + yl - 1
126 PetscCallA(DMDAGetCorners(ada, xs, ys, PETSC_NULL_INTEGER, xl, yl, PETSC_NULL_INTEGER, ierr))
129 do j = ys, ys + yl - 1
145 PetscCallA(DMDAGetCorners(ada, xs, ys, zs, xl, yl, zl, ierr))
148 do j = ys, ys + yl - 1
/petsc/src/vec/vec/tests/
H A Dex2k.c24 Vec x, *ys; in main() local
73 PetscCall(VecDuplicateVecs(x, maxN, &ys)); in main()
75 for (i = 0; i < maxN; i++) PetscCall(VecSetRandom(ys[i], rnd)); in main()
81 if (testMDot) PetscCall(VecMDot(x, N, ys, vals)); in main()
82 else if (testMAXPY) PetscCall(VecMAXPY(x, N, vals, ys)); in main()
90 if (testMDot) PetscCall(VecMDot(x, N, ys, vals)); in main()
91 else if (testMAXPY) PetscCall(VecMAXPY(x, N, vals, ys)); in main()
115 PetscCall(VecDestroyVecs(maxN, &ys)); in main()
/petsc/src/snes/tutorials/
H A Dex55k.kokkos.cxx41 …PetscInt xs = info->xs, ys = info->ys, xm = info->xm, ym = info->ym, mx = info->mx, my = info->my; in FormFunctionLocalVec() local
61 …alVec", MDRangePolicy<Rank<2, Iterate::Right, Iterate::Right>>(exec, {ys, xs}, {ys + ym, xs + xm})… in FormFunctionLocalVec()
118 …PetscInt xs = info->xs, ys = info->ys, xm = info->xm, ym = info->ym, mx = info->mx, my = info->my; in FormObjectiveLocalVec() local
140 …alVec", MDRangePolicy<Rank<2, Iterate::Right, Iterate::Right>>(exec, {ys, xs}, {ys + ym, xs + xm}), in FormObjectiveLocalVec()
176 …PetscInt xs = info->xs, ys = info->ys, xm = info->xm, ym = info->ym, mx = info->mx, my = info-… in FormJacobianLocalVec() local
190 hx = xm > 1 ? PetscRealPart(coords[ys][xs + 1].x) - PetscRealPart(coords[ys][xs].x) : 1.0; in FormJacobianLocalVec()
191 hy = ym > 1 ? PetscRealPart(coords[ys + 1][xs].y) - PetscRealPart(coords[ys][xs].y) : 1.0; in FormJacobianLocalVec()
208 for (j = ys; j < ys + ym; j++) { in FormJacobianLocalVec()
264 …alVec", MDRangePolicy<Rank<2, Iterate::Right, Iterate::Right>>(exec, {ys, xs}, {ys + ym, xs + xm})… in FormJacobianLocalVec()
265 PetscInt p = ((j - ys) * xm + (i - xs)) * 5; in FormJacobianLocalVec()
H A Dex16.c601 PetscInt xs = info->xs, ys = info->ys, zs = info->zs; in FormJacobianLocal() local
619 yes = ys; in FormJacobianLocal()
622 yee = ys + ym; in FormJacobianLocal()
625 if (ys > 0) yes = ys - 1; in FormJacobianLocal()
628 if (ys + ym == my) yee = ys + ym - 1; in FormJacobianLocal()
646 …if (i + ii >= xs && i + ii < xm + xs && j + jj >= ys && j + jj < ys + ym && k + kk >= zs && k + kk… in FormJacobianLocal()
677 for (j = ys; j < ys + ym; j++) { in FormJacobianLocal()
710 PetscInt xs = info->xs, ys = info->ys, zs = info->zs; in FormFunctionLocal() local
723 PetscCall(DMDAGetCorners(info->da, &xs, &ys, &zs, &xm, &ym, &zm)); in FormFunctionLocal()
727 for (j = ys; j < ys + ym; j++) { in FormFunctionLocal()
[all …]
H A Dex5f.F9040 PetscInt ys, ye, ym, gys, gye, gym
114 PetscScalar x(xs:xe, ys:ye)
129 do j = ys, ye
165 PetscScalar f(xs:xe, ys:ye)
176 ys = info%YS + 1
177 ye = ys + info%YM - 1
191 do j = ys, ye
287 do j = ys, ye
444 call DMDAGetCorners(da, xs, ys, PETSC_NULL_INTEGER, xm, ym, PETSC_NULL_INTEGER, ierr)
453 ys = ys + 1
[all …]
/petsc/src/dm/impls/da/ftn-custom/
H A Dzda1f90.c68 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in dmdavecgetarray1_() local
71 *ierr = DMDAGetCorners(*da, &xs, &ys, &zs, &xm, &ym, &zm); in dmdavecgetarray1_()
86 gys = ys; in dmdavecgetarray1_()
109 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in dmdavecgetarray2_() local
112 *ierr = DMDAGetCorners(*da, &xs, &ys, &zs, &xm, &ym, &zm); in dmdavecgetarray2_()
127 gys = ys; in dmdavecgetarray2_()
156 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in dmdavecgetarray3_() local
159 *ierr = DMDAGetCorners(*da, &xs, &ys, &zs, &xm, &ym, &zm); in dmdavecgetarray3_()
174 gys = ys; in dmdavecgetarray3_()
205 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof, zero = 0; in dmdavecgetarray4_() local
[all …]
/petsc/src/ksp/ksp/tutorials/
H A Dex66.c66 PetscInt i, j, M, N, xm, ym, xs, ys; in ComputeRHS() local
81 PetscCall(DMDAGetCorners(da, &xs, &ys, 0, &xm, &ym, 0)); /* Fine grid */ in ComputeRHS()
83 for (j = ys; j < ys + ym; j++) { in ComputeRHS()
100 PetscInt i, j, M, N, xm, ym, xs, ys; in ComputeJacobian() local
113 PetscCall(DMDAGetCorners(da, &xs, &ys, 0, &xm, &ym, 0)); in ComputeJacobian()
114 for (j = ys; j < ys + ym; j++) { in ComputeJacobian()
H A Dex14f.F9092 PetscInt xs, ys, xm
113 PetscCall(DMDAGetCorners(da, xs, ys, PETSC_NULL_INTEGER, xm, ym, PETSC_NULL_INTEGER, ierr))
117 do j = ys, ys + ym - 1
120 row = i - xs + (j - ys)*xm + 1
149 PetscInt i, j, row, xs, ys, xm, ym, gxs
181 PetscCall(DMDAGetCorners(da, xs, ys, PETSC_NULL_INTEGER, xm, ym, PETSC_NULL_INTEGER, ierr))
186 do j = ys, ys + ym - 1
232 PetscInt xs, ys, xm, ym
269 PetscCall(DMDAGetCorners(da, xs, ys, PETSC_NULL_INTEGER, xm, ym, PETSC_NULL_INTEGER, ierr))
288 do j = ys, ys + ym - 1
/petsc/src/ts/tutorials/autodiff/
H A Dadr_ex5adj_mf.cxx184 PetscInt i, j, xs, ys, xm, ym, Mx, My; in InitialConditions() local
202 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in InitialConditions()
207 for (j = ys; j < ys + ym; j++) { in InitialConditions()
228 PetscInt i, j, Mx, My, xs, ys, xm, ym; in InitializeLambda() local
236 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in InitializeLambda()
238 if (xs <= i && i < xs + xm && ys <= j && j < ys + ym) { in InitializeLambda()
251 PetscInt i, j, xs, ys, xm, ym; in IFunctionLocalPassive() local
264 ys = info->ys; in IFunctionLocalPassive()
268 for (j = ys; j < ys + ym; j++) { in IFunctionLocalPassive()
291 PetscInt i, j, xs, ys, xm, ym, gxs, gys, gxm, gym; in IFunctionActive() local
[all …]
/petsc/src/binding/petsc4py/demo/legacy/wrap-swig/
H A DBratu3D.c24 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormInitGuess() local
55 PetscCall(DMDAGetCorners(da,&xs,&ys,&zs,&xm,&ym,&zm)); in FormInitGuess()
62 for (j=ys; j<ys+ym; j++) { in FormInitGuess()
84 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormFunction() local
120 PetscCall(DMDAGetCorners(da,&xs,&ys,&zs,&xm,&ym,&zm)); in FormFunction()
126 for (j=ys; j<ys+ym; j++) { in FormFunction()
161 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormJacobian() local
197 PetscCall(DMDAGetCorners(da,&xs,&ys,&zs,&xm,&ym,&zm)); in FormJacobian()
215 for (j=ys; j<ys+ym; j++) { in FormJacobian()
/petsc/src/binding/petsc4py/demo/legacy/wrap-cython/
H A DBratu3Dimpl.c24 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormInitGuess() local
56 PetscCall(DMDAGetCorners(da,&xs,&ys,&zs,&xm,&ym,&zm)); in FormInitGuess()
63 for (j=ys; j<ys+ym; j++) { in FormInitGuess()
85 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormFunction() local
124 PetscCall(DMDAGetCorners(da,&xs,&ys,&zs,&xm,&ym,&zm)); in FormFunction()
130 for (j=ys; j<ys+ym; j++) { in FormFunction()
165 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormJacobian() local
204 PetscCall(DMDAGetCorners(da,&xs,&ys,&zs,&xm,&ym,&zm)); in FormJacobian()
222 for (j=ys; j<ys+ym; j++) { in FormJacobian()
/petsc/src/mat/tests/
H A Dex197.c8 Vec x, y, ys; in main() local
42 PetscCall(VecDuplicate(y, &ys)); in main()
55 PetscCall(MatMultHermitianTransposeAdd(A, x, y, ys)); in main()
56 PetscCall(VecView(ys, PETSC_VIEWER_STDOUT_WORLD)); in main()
71 PetscCall(VecDestroy(&ys)); in main()
/petsc/src/ts/tutorials/advection-diffusion-reaction/
H A Dex5opt_ic.c82 PetscInt i, j, xs, ys, xm, ym, Mx, My; in InitialConditions() local
94 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in InitialConditions()
97 for (j = ys; j < ys + ym; j++) { in InitialConditions()
114 PetscInt i, j, xs, ys, xm, ym, Mx, My; in PerturbedInitialConditions() local
126 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in PerturbedInitialConditions()
129 for (j = ys; j < ys + ym; j++) { in PerturbedInitialConditions()
146 PetscInt i, j, xs, ys, xm, ym, Mx, My; in PerturbedInitialConditions2() local
158 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in PerturbedInitialConditions2()
161 for (j = ys; j < ys + ym; j++) { in PerturbedInitialConditions2()
179 PetscInt i, j, xs, ys, xm, ym, Mx, My; in PerturbedInitialConditions3() local
[all …]
H A Dex5adj.c22 PetscInt i, j, Mx, My, xs, ys, xm, ym; in InitializeLambda() local
30 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in InitializeLambda()
32 if (xs <= i && i < xs + xm && ys <= j && j < ys + ym) { in InitializeLambda()
155 PetscInt i, j, xs, ys, xm, ym, Mx, My; in InitialConditions() local
173 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in InitialConditions()
178 for (j = ys; j < ys + ym; j++) { in InitialConditions()

12345