| /petsc/src/snes/utils/ |
| H A D | convest.c | 27 PetscErrorCode PetscConvEstDestroy(PetscConvEst *ce) in PetscConvEstDestroy() argument
30 if (!*ce) PetscFunctionReturn(PETSC_SUCCESS); in PetscConvEstDestroy()
31 PetscValidHeaderSpecific(*ce, PETSC_OBJECT_CLASSID, 1); in PetscConvEstDestroy()
32 if (--((PetscObject)*ce)->refct > 0) { in PetscConvEstDestroy()
33 *ce = NULL; in PetscConvEstDestroy()
36 PetscCall(PetscFree3((*ce)->initGuess, (*ce)->exactSol, (*ce)->ctxs)); in PetscConvEstDestroy()
37 PetscCall(PetscFree2((*ce)->dofs, (*ce)->errors)); in PetscConvEstDestroy()
38 PetscCall(PetscHeaderDestroy(ce)); in PetscConvEstDestroy()
54 PetscErrorCode PetscConvEstSetFromOptions(PetscConvEst ce) in PetscConvEstSetFromOptions() argument
57 …PetscOptionsBegin(PetscObjectComm((PetscObject)ce), "", "Convergence Estimator Options", "PetscCon… in PetscConvEstSetFromOptions()
[all …]
|
| /petsc/src/ts/utils/ |
| H A D | tsconvest.c | 7 static PetscErrorCode PetscConvEstSetTS_Private(PetscConvEst ce, PetscObject solver) in PetscConvEstSetTS_Private() argument
12 PetscCall(PetscObjectGetClassId(ce->solver, &id)); in PetscConvEstSetTS_Private()
13 …PetscCheck(id == TS_CLASSID, PetscObjectComm((PetscObject)ce), PETSC_ERR_ARG_WRONG, "Solver was no… in PetscConvEstSetTS_Private()
14 PetscCall(TSGetDM((TS)ce->solver, &ce->idm)); in PetscConvEstSetTS_Private()
18 static PetscErrorCode PetscConvEstInitGuessTS_Private(PetscConvEst ce, PetscInt r, DM dm, Vec u) in PetscConvEstInitGuessTS_Private() argument
21 PetscCall(TSComputeInitialCondition((TS)ce->solver, u)); in PetscConvEstInitGuessTS_Private()
25 static PetscErrorCode PetscConvEstComputeErrorTS_Private(PetscConvEst ce, PetscInt r, DM dm, Vec u,… in PetscConvEstComputeErrorTS_Private() argument
27 TS ts = (TS)ce->solver; in PetscConvEstComputeErrorTS_Private()
39 for (f = 1; f < ce->Nf; ++f) errors[f] = errors[0]; in PetscConvEstComputeErrorTS_Private()
45 PetscCall(DMComputeL2FieldDiff(dm, t, ce->exactSol, ce->ctxs, u, errors)); in PetscConvEstComputeErrorTS_Private()
[all …]
|
| /petsc/src/dm/impls/da/ |
| H A D | dageometry.c | 39 …e DMGetLocalBoundingBox_DA(DM da, PetscReal lmin[], PetscReal lmax[], PetscInt cs[], PetscInt ce[]) in DMGetLocalBoundingBox_DA() argument
72 if (ce) ce[d] = lmax[d]; in DMGetLocalBoundingBox_DA()
122 if (ce) { in DMGetLocalBoundingBox_DA()
123 ce[0] = xe; in DMGetLocalBoundingBox_DA()
124 if (dim > 1) ce[1] = ye; in DMGetLocalBoundingBox_DA()
125 if (dim > 2) ce[2] = ze; in DMGetLocalBoundingBox_DA()
133 PetscInt cs[2], ce[2]; in private_DMDALocatePointsIS_2D_Regular() local
142 PetscCall(DMGetLocalBoundingBox_DA(dmregular, gmin_l, gmax_l, cs, ce)); in private_DMDALocatePointsIS_2D_Regular()
145 xe = ce[0]; in private_DMDALocatePointsIS_2D_Regular()
146 ye = ce[1]; in private_DMDALocatePointsIS_2D_Regular()
[all …]
|
| /petsc/src/mat/graphops/color/utils/ |
| H A D | bipartite.c | 7 PetscInt i, j, rs, re, cs, ce, cn; in MatColoringCreateBipartiteGraph() local
17 PetscCall(MatGetOwnershipRangeColumn(m, &cs, &ce)); in MatColoringCreateBipartiteGraph()
18 cn = ce - cs; in MatColoringCreateBipartiteGraph()
|
| /petsc/src/tao/constrained/tutorials/ |
| H A D | ex1.c | 41 Vec ce, ci, bl, bu, Xseq; member
86 … PetscCall(TaoSetEqualityConstraintsRoutine(tao, user.ce, FormEqualityConstraints, (void *)&user)); in main()
202 PetscCall(VecCreate(PETSC_COMM_WORLD, &user->ce)); /* a 1x1 vec for equality constraints */ in InitializeProblem()
203 PetscCall(VecSetSizes(user->ce, neloc, user->ne)); in InitializeProblem()
204 PetscCall(VecSetFromOptions(user->ce)); in InitializeProblem()
205 PetscCall(VecSetUp(user->ce)); in InitializeProblem()
241 if (!user->noeqflag) PetscCall(VecDestroy(&user->ce)); in DestroyProblem()
387 PetscScalar ce; in FormEqualityConstraints() local
403 ce = x[0] * x[0] + x[1] - 2.0; in FormEqualityConstraints()
404 PetscCall(VecSetValue(CE, 0, ce, INSERT_VALUES)); in FormEqualityConstraints()
|
| H A D | maros.c | 240 PetscErrorCode FormEqualityConstraints(Tao tao, Vec x, Vec ce, PetscCtx ctx) in FormEqualityConstraints() argument
245 PetscCall(MatMult(user->Aeq, x, ce)); in FormEqualityConstraints()
246 PetscCall(VecAXPY(ce, -1.0, user->beq)); in FormEqualityConstraints()
|
| /petsc/src/tao/interface/ |
| H A D | taosolver_bounds.c | 329 …rCode TaoSetEqualityConstraintsRoutine(Tao tao, Vec ce, PetscErrorCode (*func)(Tao tao, Vec x, Vec… in TaoSetEqualityConstraintsRoutine() argument
333 if (ce) PetscValidHeaderSpecific(ce, VEC_CLASSID, 2); in TaoSetEqualityConstraintsRoutine()
334 PetscCall(PetscObjectReference((PetscObject)ce)); in TaoSetEqualityConstraintsRoutine()
337 tao->constraints_equality = ce; in TaoSetEqualityConstraintsRoutine()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex8.c | 179 CECtx *ce; in CECreate() local
182 PetscCall(PetscMalloc(sizeof(CECtx), &ce)); in CECreate()
183 p->data = (void *)ce; in CECreate()
193 ce->lambda = 10; in CECreate()
196 …, "Parameter controlling stiffness: xdot + lambda*(x - cos(t))", "", ce->lambda, &ce->lambda, NULL… in CECreate()
|
| /petsc/src/tao/constrained/impls/ipm/ |
| H A D | pdipm.h | 35 Vec ce, ci; /* equality and inequality constraints */ member
|
| H A D | pdipm.c | 66 PetscCall(VecGetArrayWrite(pdipm->ce, &carr)); in TaoPDIPMUpdateConstraints()
84 PetscCall(VecRestoreArrayWrite(pdipm->ce, &carr)); in TaoPDIPMUpdateConstraints()
504 PetscCall(VecGetArrayRead(pdipm->ce, &carr)); in TaoSNESFunction_PDIPM()
506 PetscCall(VecRestoreArrayRead(pdipm->ce, &carr)); in TaoSNESFunction_PDIPM()
599 PetscCall(VecPlaceArray(pdipm->ce, Farr + pdipm->off_lambdae)); in TaoSNESFunction_PDIPM_residual()
600 PetscCall(VecNorm(pdipm->ce, NORM_2, &cnorm[0])); in TaoSNESFunction_PDIPM_residual()
601 PetscCall(VecResetArray(pdipm->ce)); in TaoSNESFunction_PDIPM_residual()
901 PetscCall(VecCreate(comm, &pdipm->ce)); in TaoSetup_PDIPM()
902 PetscCall(VecSetSizes(pdipm->ce, pdipm->nce, pdipm->Nce)); in TaoSetup_PDIPM()
903 PetscCall(VecSetFromOptions(pdipm->ce)); in TaoSetup_PDIPM()
[all …]
|
| /petsc/src/ksp/pc/impls/gamg/ |
| H A D | classical.c | 190 PetscInt fn, cn, fs, fe, cs, ce, i, j, ncols, col, row_f, row_c, cmax = 0, idx, noff; in PCGAMGProlongator_Classical_Direct() local
237 PetscCall(VecGetOwnershipRange(C, &cs, &ce)); in PCGAMGProlongator_Classical_Direct()
562 PetscInt fs, fe, cs, ce, nl, i, j, k, li, lni, ci, ncols, maxcols, fn, cn, cid; in PCGAMGProlongator_Classical_Standard() local
608 PetscCall(VecGetOwnershipRange(cv, &cs, &ce)); in PCGAMGProlongator_Classical_Standard()
658 if (lni >= cs && lni < ce) { in PCGAMGProlongator_Classical_Standard()
671 if (lni >= cs && lni < ce) { in PCGAMGProlongator_Classical_Standard()
|
| /petsc/src/binding/petsc4py/test/ |
| H A D | test_mat_aij.py | 283 cs, ce = self.A.getOwnershipRangeColumn()
285 cols = N.array(range(cs, ce), dtype=PETSc.IntType)
309 cs, ce = self.A.getOwnershipRangeColumn()
311 cols = N.array(range(cs, ce), dtype=PETSc.IntType)
|
| H A D | test_vec.py | 248 cs, ce = A.getOwnershipRangeColumn()
256 for i in range(cs, ce):
|
| /petsc/src/ksp/ksp/interface/ |
| H A D | itfunc.c | 2614 PetscErrorCode KSPComputeConvergenceRate(KSP ksp, PetscReal *cr, PetscReal *rRsq, PetscReal *ce, Pe… in KSPComputeConvergenceRate() argument
2644 if (ce || eRsq) { in KSPComputeConvergenceRate()
2647 if (ce) *ce = 0.0; in KSPComputeConvergenceRate()
2663 if (ce) *ce = PetscExpReal(slope); in KSPComputeConvergenceRate()
|
| /petsc/src/binding/petsc4py/src/petsc4py/PETSc/ |
| H A D | DMSwarm.pyx | 857 cdef PetscInt ce = asInt(e)
861 CHKERR(DMSwarmSortGetPointsPerCell(self.dm, ce, &cnpoints, &cpidlist))
881 cdef PetscInt ce = asInt(e)
883 CHKERR(DMSwarmSortGetNumberOfPointsPerCell(self.dm, ce, &npoints))
|
| /petsc/doc/tutorials/ |
| H A D | in_person.md | 66 …[Slides](http://www.tacc.utexas.edu/c/document_library/get_file?uuid=802633ce-f165-407f-bec8-f0efc…
|
| /petsc/src/mat/impls/hypre/ |
| H A D | mhypre.c | 1618 PetscInt i, rs, re, cs, ce, bs; local
1627 ce = A->cmap->rend;
1629 PetscCallHYPRE(HYPRE_IJMatrixCreate(hA->comm, rs, re - 1, cs, ce - 1, &hA->ij));
1635 …ce - cs, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Inconsistent local cols: IJMatrix [%" PetscHYPRE_BigInt…
|
| /petsc/doc/manual/ |
| H A D | tao.md | 2094 &ce(x) = 0 \\
2099 Here, $ce(x)$ is set of equality constraints that include
2108 L_{\mu}(x,\lambda_{ce},\lambda_{ci},z) = f(x) + \lambda_{ce}^Tce(x) - \lambda_{ci}^T(ci(x) - z) - \…
2111 where, $\lambda_{ce}$ and $\lambda_{ci}$ are the Lagrangian
2117 \nabla L_{\mu}(x,\lambda_{ce},\lambda_{ci},z) =
2119 \nabla f(x) + \nabla ce(x)^T\lambda_{ce} - \nabla ci(x)^T \lambda_{ci} \\
2120 ce(x) \\
|
| /petsc/src/dm/interface/ |
| H A D | dmcoordinates.c | 881 …calBoundingBox_Coordinates(DM dm, PetscReal lmin[], PetscReal lmax[], PetscInt cs[], PetscInt ce[]) in DMGetLocalBoundingBox_Coordinates() argument
|
| /petsc/src/mat/impls/baij/mpi/ |
| H A D | mpibaij.c | 3230 PetscInt *rowidxs, *colidxs, rs, cs, ce; in MatLoad_MPIBAIJ_Binary() local
3263 PetscCall(PetscLayoutGetRange(mat->cmap, &cs, &ce)); in MatLoad_MPIBAIJ_Binary()
3298 if (col >= cs && col < ce) { in MatLoad_MPIBAIJ_Binary()
|
| /petsc/doc/ |
| H A D | petsc.bib | 4982 url = {http://www.ce.washington.edu/em03/proceedings/papers/477.pdf},
8111 url = {http://www.ce.berkeley.edu/projects/feap/},
10612 title = {{I}ce {S}heet {S}ystem {M}odel},
|