Home
last modified time | relevance | path

Searched refs:ind (Results 1 – 25 of 37) sorted by relevance

12

/petsc/src/tao/unconstrained/tutorials/
H A Deptorsion3.c331 PetscInt nx = user->mx, ny = user->my, ind, i, j, k; in FormGradient() local
357 ind = k; in FormGradient()
359 PetscCall(VecSetValues(G, 1, &ind, &val, ADD_VALUES)); in FormGradient()
362 ind = k + 1; in FormGradient()
364 PetscCall(VecSetValues(G, 1, &ind, &val, ADD_VALUES)); in FormGradient()
367 ind = k + nx; in FormGradient()
369 PetscCall(VecSetValues(G, 1, &ind, &val, ADD_VALUES)); in FormGradient()
387 ind = k - nx; in FormGradient()
389 PetscCall(VecSetValues(G, 1, &ind, &val, ADD_VALUES)); in FormGradient()
392 ind = k - 1; in FormGradient()
[all …]
H A Deptorsion1.c315 PetscInt nx = user->mx, ny = user->my, ind, i, j, k; in FormGradient() local
341 ind = k; in FormGradient()
343 PetscCall(VecSetValues(G, 1, &ind, &val, ADD_VALUES)); in FormGradient()
346 ind = k + 1; in FormGradient()
348 PetscCall(VecSetValues(G, 1, &ind, &val, ADD_VALUES)); in FormGradient()
351 ind = k + nx; in FormGradient()
353 PetscCall(VecSetValues(G, 1, &ind, &val, ADD_VALUES)); in FormGradient()
371 ind = k - nx; in FormGradient()
373 PetscCall(VecSetValues(G, 1, &ind, &val, ADD_VALUES)); in FormGradient()
376 ind = k - 1; in FormGradient()
[all …]
H A Drosenbrock1.c196 PetscInt i, ind[2]; in FormHessian() local
219 ind[0] = i; in FormHessian()
220 ind[1] = i + 1; in FormHessian()
221 PetscCall(MatSetValues(H, 2, ind, 2, ind, v[0], ADD_VALUES)); in FormHessian()
228 ind[0] = 2 * i; in FormHessian()
229 ind[1] = 2 * i + 1; in FormHessian()
230 PetscCall(MatSetValues(H, 2, ind, 2, ind, v[0], INSERT_VALUES)); in FormHessian()
H A Drosenbrock2.c187 PetscInt i, ind[2]; in FormHessian() local
210 ind[0] = i; in FormHessian()
211 ind[1] = i + 1; in FormHessian()
212 PetscCall(MatSetValues(H, 2, ind, 2, ind, v[0], ADD_VALUES)); in FormHessian()
219 ind[0] = 2 * i; in FormHessian()
220 ind[1] = 2 * i + 1; in FormHessian()
221 PetscCall(MatSetValues(H, 2, ind, 2, ind, v[0], INSERT_VALUES)); in FormHessian()
H A Drosenbrock3.c192 PetscInt i, ind[2]; in FormHessian() local
214 ind[0] = i; in FormHessian()
215 ind[1] = i + 1; in FormHessian()
216 PetscCall(MatSetValues(H, 2, ind, 2, ind, v[0], ADD_VALUES)); in FormHessian()
223 ind[0] = 2 * i; in FormHessian()
224 ind[1] = 2 * i + 1; in FormHessian()
225 PetscCall(MatSetValues(H, 2, ind, 2, ind, v[0], INSERT_VALUES)); in FormHessian()
H A Deptorsion2.c199 PetscInt i, j, k, ind; in FormFunctionGradient() local
256 ind = k; in FormFunctionGradient()
261 ind = k + 1; in FormFunctionGradient()
263 PetscCall(VecSetValuesLocal(G, 1, &ind, &val, ADD_VALUES)); in FormFunctionGradient()
266 ind = k + gxm; in FormFunctionGradient()
268 PetscCall(VecSetValuesLocal(G, 1, &ind, &val, ADD_VALUES)); in FormFunctionGradient()
288 ind = k - gxm; in FormFunctionGradient()
290 PetscCall(VecSetValuesLocal(G, 1, &ind, &val, ADD_VALUES)); in FormFunctionGradient()
293 ind = k - 1; in FormFunctionGradient()
295 PetscCall(VecSetValuesLocal(G, 1, &ind, &val, ADD_VALUES)); in FormFunctionGradient()
[all …]
H A Deptorsion2f.F90122 PetscInt xep, yep, i, j, k, ind
189 ind = k
194 ind = k + 1
196 PetscCall(VecSetValuesLocal(G, i1, [ind], [val], ADD_VALUES, ierr))
199 ind = k + gxm
201 PetscCall(VecSetValuesLocal(G, i1, [ind], [val], ADD_VALUES, ierr))
222 ind = k - gxm
224 PetscCall(VecSetValuesLocal(G, i1, [ind], [val], ADD_VALUES, ierr))
227 ind = k - 1
229 PetscCall(VecSetValuesLocal(G, i1, [ind], [val], ADD_VALUES, ierr))
[all …]
H A Drosenbrock1f.F90109 PetscInt i, nn, ind(0:1), i2
133 ind(0) = 2*i
134 ind(1) = 2*i + 1
135 PetscCall(MatSetValues(H, i2, ind, i2, ind, reshape(v, [i2*i2]), INSERT_VALUES, ierr))
H A Dburgers_spectral.c100 PetscInt i, xs, xm, ind, j, lenglob; in main() local
177 ind = i * (appctx.param.N - 1) + j; in main()
178 wrk_ptr1[ind] = x; in main()
179 wrk_ptr2[ind] = .5 * appctx.param.Le * appctx.SEMop.gll.weights[j]; in main()
180 if (j == 0) wrk_ptr2[ind] += .5 * appctx.param.Le * appctx.SEMop.gll.weights[j]; in main()
/petsc/src/snes/tests/
H A Dex2.c116 PetscInt N = 3, n = 0, dim, spaceDim, i, j, k, *ind, d; in CreatePoints_Grid() local
124 PetscCall(PetscCalloc1(spaceDim, &ind)); in CreatePoints_Grid()
140 ind[2] = k; in CreatePoints_Grid()
142 ind[1] = j; in CreatePoints_Grid()
144 ind[0] = i; in CreatePoints_Grid()
146 for (d = 0; d < spaceDim; ++d) (*pcoords)[n * spaceDim + d] = ind[d] * h[d]; in CreatePoints_Grid()
158 PetscCall(PetscFree(ind)); in CreatePoints_Grid()
166 PetscInt N = 3, n = 0, dim, spaceDim, i, j, k, *ind, d; in CreatePoints_GridReplicated() local
174 PetscCall(PetscCalloc1(spaceDim, &ind)); in CreatePoints_GridReplicated()
182 ind[2] = k; in CreatePoints_GridReplicated()
[all …]
/petsc/src/vec/is/is/tests/
H A Dex9.c63 PetscInt *ind; in main() local
66 PetscCall(PetscSFGetGraphLayout(sf, &lt, &nl, NULL, &ind)); in main()
70 …ind[i] == indices[i], PETSC_COMM_SELF, PETSC_ERR_PLIB, "PetscSFSetGraphLayout() gives global_roots… in main()
72 PetscCall(PetscFree(ind)); in main()
/petsc/src/dm/impls/plex/tests/
H A Dex34.c111 PetscInt *ind; in OrderHybridMesh() local
118 PetscCall(PetscMalloc1(pEnd - pStart, &ind)); in OrderHybridMesh()
119 for (p = 0; p < pEnd - pStart; ++p) ind[p] = p; in OrderHybridMesh()
133 if (coneSize == 6) ind[c] = off[1]++; in OrderHybridMesh()
134 else ind[c] = off[0]++; in OrderHybridMesh()
138 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, pEnd - pStart, ind, PETSC_OWN_POINTER, &perm)); in OrderHybridMesh()
/petsc/src/vec/is/utils/
H A Disio.c41 PetscInt *ind; in ISLoad_HDF5() local
48 PetscCall(PetscViewerHDF5Load(viewer, isname, is->map, H5T_NATIVE_LLONG, (void **)&ind)); in ISLoad_HDF5()
50 PetscCall(PetscViewerHDF5Load(viewer, isname, is->map, H5T_NATIVE_INT, (void **)&ind)); in ISLoad_HDF5()
52 PetscCall(ISGeneralSetIndices(is, is->map->n, ind, PETSC_OWN_POINTER)); in ISLoad_HDF5()
/petsc/src/tao/leastsquares/tutorials/matlab/more_wild_probs/
H A Djacobian.m32 for ind = 1:n % Do one coordinate direction at a time:
33 [g_fvec, fvec] = g_dfovec_1d(g_t, t, ind, m, n, g_x, x, nprob);
34 J(:,ind) = g_fvec;
/petsc/src/ksp/pc/impls/is/nn/
H A Dnn.c270 PetscMPIInt ind = 0; in PCNNCreateCoarseMatrix() local
272 PetscCallMPI(MPI_Waitany(n_neigh - 1, recv_request + 1, &ind, &stat)); in PCNNCreateCoarseMatrix()
273 ind++; in PCNNCreateCoarseMatrix()
275 mat[ind * n_neigh + 0] = 0.0; in PCNNCreateCoarseMatrix()
276 …for (k = 0; k < n_shared[ind]; k++) mat[ind * n_neigh + 0] += DZ_IN[ind][k] * pcis->work_N[shared[… in PCNNCreateCoarseMatrix()
/petsc/src/dm/impls/da/
H A Ddalocal.c364 PetscInt ind[3], d, off; in DMDASetVertexCoordinates() local
366 ind[0] = 0; in DMDASetVertexCoordinates()
367 ind[1] = 0; in DMDASetVertexCoordinates()
368 ind[2] = k + da->zs; in DMDASetVertexCoordinates()
370 ind[1] = j + da->ys; in DMDASetVertexCoordinates()
375 ind[0] = i + da->xs; in DMDASetVertexCoordinates()
376 for (d = 0; d < dim; ++d) coords[off + d] = h[d] * ind[d]; in DMDASetVertexCoordinates()
/petsc/doc/community/
H A Dindex.md1 (ind-contact)=
/petsc/src/dm/impls/stag/
H A Dstag2d.c134 PetscInt ind[2], start_ghost[2], n_ghost[2], s, c; in DMStagSetUniformCoordinatesExplicit_2d() local
160 for (ind[1] = start_ghost[1]; ind[1] < start_ghost[1] + n_ghost[1]; ++ind[1]) { in DMStagSetUniformCoordinatesExplicit_2d()
161 for (ind[0] = start_ghost[0]; ind[0] < start_ghost[0] + n_ghost[0]; ++ind[0]) { in DMStagSetUniformCoordinatesExplicit_2d()
164 …for (c = 0; c < 2; ++c) arr[ind[1]][ind[0]][idownleft + c] = min[c] + ((PetscReal)ind[c] + offs[c]… in DMStagSetUniformCoordinatesExplicit_2d()
168 …for (c = 0; c < 2; ++c) arr[ind[1]][ind[0]][idown + c] = min[c] + ((PetscReal)ind[c] + offs[c]) * … in DMStagSetUniformCoordinatesExplicit_2d()
172 …for (c = 0; c < 2; ++c) arr[ind[1]][ind[0]][ileft + c] = min[c] + ((PetscReal)ind[c] + offs[c]) * … in DMStagSetUniformCoordinatesExplicit_2d()
176 …for (c = 0; c < 2; ++c) arr[ind[1]][ind[0]][ielement + c] = min[c] + ((PetscReal)ind[c] + offs[c])… in DMStagSetUniformCoordinatesExplicit_2d()
H A Dstag1d.c115 for (PetscInt ind = start_ghost; ind < start_ghost + n_ghost; ++ind) { in DMStagSetUniformCoordinatesExplicit_1d() local
118 arr[ind][ileft] = min + ((PetscReal)ind + off) * h; in DMStagSetUniformCoordinatesExplicit_1d()
122 arr[ind][ielement] = min + ((PetscReal)ind + off) * h; in DMStagSetUniformCoordinatesExplicit_1d()
H A Dstag3d.c206 PetscInt ind[3], start_ghost[3], n_ghost[3], s, c; in DMStagSetUniformCoordinatesExplicit_3d() local
238 for (ind[2] = start_ghost[2]; ind[2] < start_ghost[2] + n_ghost[2]; ++ind[2]) { in DMStagSetUniformCoordinatesExplicit_3d()
239 for (ind[1] = start_ghost[1]; ind[1] < start_ghost[1] + n_ghost[1]; ++ind[1]) { in DMStagSetUniformCoordinatesExplicit_3d()
240 for (ind[0] = start_ghost[0]; ind[0] < start_ghost[0] + n_ghost[0]; ++ind[0]) { in DMStagSetUniformCoordinatesExplicit_3d()
243 …for (c = 0; c < 3; ++c) arr[ind[2]][ind[1]][ind[0]][ibackdownleft + c] = min[c] + ((PetscReal)ind[… in DMStagSetUniformCoordinatesExplicit_3d()
247 …for (c = 0; c < 3; ++c) arr[ind[2]][ind[1]][ind[0]][ibackdown + c] = min[c] + ((PetscReal)ind[c] +… in DMStagSetUniformCoordinatesExplicit_3d()
251 …for (c = 0; c < 3; ++c) arr[ind[2]][ind[1]][ind[0]][ibackleft + c] = min[c] + ((PetscReal)ind[c] +… in DMStagSetUniformCoordinatesExplicit_3d()
255 …for (c = 0; c < 3; ++c) arr[ind[2]][ind[1]][ind[0]][iback + c] = min[c] + ((PetscReal)ind[c] + off… in DMStagSetUniformCoordinatesExplicit_3d()
259 …for (c = 0; c < 3; ++c) arr[ind[2]][ind[1]][ind[0]][idownleft + c] = min[c] + ((PetscReal)ind[c] +… in DMStagSetUniformCoordinatesExplicit_3d()
263 …for (c = 0; c < 3; ++c) arr[ind[2]][ind[1]][ind[0]][idown + c] = min[c] + ((PetscReal)ind[c] + off… in DMStagSetUniformCoordinatesExplicit_3d()
[all …]
/petsc/src/vec/is/utils/hdf5/
H A Dhdf5io.c261 PetscInt *ind; in PetscViewerHDF5Load_Internal() local
265 PetscCall(PetscMalloc1(map->n, &ind)); in PetscViewerHDF5Load_Internal()
272 if (off >= map->rstart) ind[off - map->rstart] = h->cind[i * 3 + 2] + inc; in PetscViewerHDF5Load_Internal()
275 *newarr = ind; in PetscViewerHDF5Load_Internal()
/petsc/src/dm/dt/interface/
H A Ddtweakform.c173 …orm ht, DMLabel label, PetscInt value, PetscInt f, PetscInt part, PetscInt ind, PetscVoidFn **func) in PetscWeakFormGetIndexFunction_Private() argument
187 …PetscCheck(ind < chunk.size, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Index %" PetscInt_FMT " n… in PetscWeakFormGetIndexFunction_Private()
188 *func = ((PetscVoidFn **)&wf->funcs->array[chunk.start])[ind]; in PetscWeakFormGetIndexFunction_Private()
194 …Form ht, DMLabel label, PetscInt value, PetscInt f, PetscInt part, PetscInt ind, PetscVoidFn *func) in PetscWeakFormSetIndexFunction_Private() argument
207 PetscCall(PetscChunkBufferCreateChunk(wf->funcs, ind + 1, &chunk)); in PetscWeakFormSetIndexFunction_Private()
209 } else if (chunk.size <= ind) { in PetscWeakFormSetIndexFunction_Private()
210 PetscCall(PetscChunkBufferEnlargeChunk(wf->funcs, ind - chunk.size + 1, &chunk)); in PetscWeakFormSetIndexFunction_Private()
213 ((PetscVoidFn **)&wf->funcs->array[chunk.start])[ind] = func; in PetscWeakFormSetIndexFunction_Private()
217 …kForm wf, PetscHMapForm ht, DMLabel label, PetscInt value, PetscInt f, PetscInt part, PetscInt ind) in PetscWeakFormClearIndexFunction_Private() argument
230 } else if (!ind && chunk.size == 1) { in PetscWeakFormClearIndexFunction_Private()
[all …]
/petsc/src/vec/is/sf/utils/
H A Dsfutils.c454 PetscInt i, ind; in PetscSFCreateSectionSF() local
481 for (i = 0, ind = 0; i < numPoints; ++i) { in PetscSFCreateSectionSF()
491 for (d = 0; d < dof; ++d, ++ind) { in PetscSFCreateSectionSF()
492 localIndices[ind] = loff + d; in PetscSFCreateSectionSF()
493 remoteIndices[ind].rank = rank; in PetscSFCreateSectionSF()
494 remoteIndices[ind].index = remoteOffset + d; in PetscSFCreateSectionSF()
498 …PetscCheck(numIndices == ind, comm, PETSC_ERR_PLIB, "Inconsistency in indices, %" PetscInt_FMT " s… in PetscSFCreateSectionSF()
/petsc/src/ts/tutorials/
H A Dex50.c86 PetscInt i, xs, xm, ind, j, lenglob; in main() local
159 ind = i * (appctx.param.N - 1) + j; in main()
160 wrk_ptr1[ind] = x; in main()
161 wrk_ptr2[ind] = .5 * appctx.param.Le * appctx.SEMop.gll.weights[j]; in main()
162 if (j == 0) wrk_ptr2[ind] += .5 * appctx.param.Le * appctx.SEMop.gll.weights[j]; in main()
/petsc/src/dm/interface/
H A Ddmi.c124 const PetscInt *ind; in PetscSectionSelectFields_Private() local
129 PetscCall(PetscSectionGetFieldConstraintIndices(s, point, fields[f], &ind)); in PetscSectionSelectFields_Private()
188 const PetscInt *ind; in PetscSectionSelectFields_Private() local
191 PetscCall(PetscSectionGetFieldConstraintIndices(s, point, fields[f], &ind)); in PetscSectionSelectFields_Private()
196 if (ind[fcoff] == k) { in PetscSectionSelectFields_Private()

12