Home
last modified time | relevance | path

Searched refs:ci (Results 1 – 25 of 40) sorted by relevance

12

/petsc/src/tao/interface/
H A Dtaosolver_bounds.c366 …ode TaoGetEqualityConstraintsRoutine(Tao tao, Vec *ci, PetscErrorCode (**func)(Tao tao, Vec x, Vec… in TaoGetEqualityConstraintsRoutine() argument
370 if (ci) *ci = tao->constraints_equality; in TaoGetEqualityConstraintsRoutine()
397 …ode TaoSetInequalityConstraintsRoutine(Tao tao, Vec ci, PetscErrorCode (*func)(Tao tao, Vec x, Vec… in TaoSetInequalityConstraintsRoutine() argument
401 if (ci) PetscValidHeaderSpecific(ci, VEC_CLASSID, 2); in TaoSetInequalityConstraintsRoutine()
402 PetscCall(PetscObjectReference((PetscObject)ci)); in TaoSetInequalityConstraintsRoutine()
404 tao->constraints_inequality = ci; in TaoSetInequalityConstraintsRoutine()
434 …e TaoGetInequalityConstraintsRoutine(Tao tao, Vec *ci, PetscErrorCode (**func)(Tao tao, Vec x, Vec… in TaoGetInequalityConstraintsRoutine() argument
438 if (ci) *ci = tao->constraints_inequality; in TaoGetInequalityConstraintsRoutine()
/petsc/src/mat/impls/aij/seq/
H A Dmatptap.c69 PetscInt *ci, *cj, *ptadenserow, *ptasparserow, *ptaj, nspacedouble = 0; in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy() local
83 PetscCall(PetscMalloc1(pn + 1, &ci)); in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
84 ci[0] = 0; in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
143 ci[i + 1] = ci[i] + cnzi; in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
148 PetscCall(PetscMalloc1(ci[pn], &cj)); in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
153 PetscCall(PetscCalloc1(ci[pn], &ca)); in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
156 …PetscCall(MatSetSeqAIJWithArrays_private(PetscObjectComm((PetscObject)A), pn, pn, ci, cj, ca, ((Pe… in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
169 afill = (PetscReal)ci[pn] / (ai[am] + pi[pm] + 1.e-5); in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
178 if (ci[pn] != 0) { in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
194 PetscInt *ci = c->i, *cj = c->j, *cjj; in MatPtAPNumeric_SeqAIJ_SeqAIJ_SparseAxpy() local
[all …]
H A Dmatmatmult.c147 PetscInt *ai = a->i, *bi = b->i, *ci, *cj; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed() local
159 PetscCall(PetscMalloc1(am + 2, &ci)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
160 ci[0] = 0; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
204 ci[i + 1] = ci[i] + cnzi; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
210 PetscCall(PetscMalloc1(ci[am] + 1, &cj)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
215 …PetscCall(MatSetSeqAIJWithArrays_private(PetscObjectComm((PetscObject)A), am, bn, ci, cj, NULL, ((… in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
229 afill = (PetscReal)ci[am] / (ai[am] + bi[bm]) + 1.e-5; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
236 if (ci[am]) { in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
252 PetscInt *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j, *bjj, *ci = c->i, *cj = c->j; in MatMatMultNumeric_SeqAIJ_SeqAIJ_Sorted() local
264 PetscCall(PetscMalloc1(ci[cm] + 1, &ca)); in MatMatMultNumeric_SeqAIJ_SeqAIJ_Sorted()
[all …]
H A Dfdaij.c179 const PetscInt *is, *row, *ci, *cj; in MatFDColoringSetUp_SeqXAIJ() local
222 …PetscCall(MatGetColumnIJ_SeqBAIJ_Color(mat, 0, PETSC_FALSE, PETSC_FALSE, &ncols, &ci, &cj, &spidx,… in MatFDColoringSetUp_SeqXAIJ()
224 …PetscCall(MatGetColumnIJ_SeqSELL_Color(mat, 0, PETSC_FALSE, PETSC_FALSE, &ncols, &ci, &cj, &spidx,… in MatFDColoringSetUp_SeqXAIJ()
226 …PetscCall(MatGetColumnIJ_SeqAIJ_Color(mat, 0, PETSC_FALSE, PETSC_FALSE, &ncols, &ci, &cj, &spidx, … in MatFDColoringSetUp_SeqXAIJ()
247 tmp = ci[col]; in MatFDColoringSetUp_SeqXAIJ()
249 m = ci[col + 1] - tmp; in MatFDColoringSetUp_SeqXAIJ()
286 …PetscCall(MatRestoreColumnIJ_SeqBAIJ_Color(mat, 0, PETSC_FALSE, PETSC_FALSE, &ncols, &ci, &cj, &sp… in MatFDColoringSetUp_SeqXAIJ()
289 …PetscCall(MatRestoreColumnIJ_SeqSELL_Color(mat, 0, PETSC_FALSE, PETSC_FALSE, &ncols, &ci, &cj, &sp… in MatFDColoringSetUp_SeqXAIJ()
291 …PetscCall(MatRestoreColumnIJ_SeqAIJ_Color(mat, 0, PETSC_FALSE, PETSC_FALSE, &ncols, &ci, &cj, &spi… in MatFDColoringSetUp_SeqXAIJ()
/petsc/lib/petsc/bin/maint/
H A Drunjobs.py54 ci = fd.read() variable
60 if (ci.find(File_Header) >= 0):
63 saved_jobs = sorted(set((ci.splitlines()[1]).split()[2:]))
73 if ci.find('TEST_ARCH: arch-ci-'+arch) > -1:
96 for a in ci.split('\n\n'):
H A Dcheck-ci-settings.sh10 if git diff --exit-code HEAD..."${dest}" -- .gitlab-ci.yml lib/petsc/conf/rules; then
H A Dxclude5 petsc-dist/.gitlab-ci.yml
6 petsc-dist/.gitlab-alcf-ci.yml
H A Dbuilddist79 version_date_git=`git log -1 --pretty=format:%ci`
/petsc/src/mat/impls/aij/seq/cholmod/
H A Daijcholmod.c10 PetscInt m = A->rmap->n, i, j, k, nz, *ci, *cj; in MatWrapCholmod_seqaij() local
16 PetscCall(PetscMalloc2(m + 1, &ci, nz, &cj)); in MatWrapCholmod_seqaij()
23 ci[i] = k; in MatWrapCholmod_seqaij()
29 ci[i] = k; in MatWrapCholmod_seqaij()
39 C->p = ci; in MatWrapCholmod_seqaij()
/petsc/src/tao/constrained/tutorials/
H A Dex1.c41 Vec ce, ci, bl, bu, Xseq; member
89 …PetscCall(TaoSetInequalityConstraintsRoutine(tao, user.ci, FormInequalityConstraints, (void *)&use… in main()
208 PetscCall(VecCreate(PETSC_COMM_WORLD, &user->ci)); /* a 2x1 vec for inequality constraints */ in InitializeProblem()
209 PetscCall(VecSetSizes(user->ci, niloc, user->ni)); in InitializeProblem()
210 PetscCall(VecSetFromOptions(user->ci)); in InitializeProblem()
211 PetscCall(VecSetUp(user->ci)); in InitializeProblem()
242 PetscCall(VecDestroy(&user->ci)); in DestroyProblem()
354 PetscScalar ci; in FormInequalityConstraints() local
370 ci = x[0] * x[0] - x[1]; in FormInequalityConstraints()
371 PetscCall(VecSetValue(CI, 0, ci, INSERT_VALUES)); in FormInequalityConstraints()
[all …]
H A Dmaros.c231 PetscErrorCode FormInequalityConstraints(Tao tao, Vec x, Vec ci, PetscCtx ctx) in FormInequalityConstraints() argument
236 PetscCall(MatMult(user->Ain, x, ci)); in FormInequalityConstraints()
/petsc/src/ksp/pc/impls/gamg/
H A Dclassical.c562 PetscInt fs, fe, cs, ce, nl, i, j, k, li, lni, ci, ncols, maxcols, fn, cn, cid; in PCGAMGProlongator_Classical_Standard() local
645 ci = icol[j]; in PCGAMGProlongator_Classical_Standard()
647 PetscCall(MatGetRow(lA, ci, &ncols, &icol, NULL)); in PCGAMGProlongator_Classical_Standard()
651 PetscCall(MatRestoreRow(lA, ci, &ncols, &icol, NULL)); in PCGAMGProlongator_Classical_Standard()
665 ci = icol[j]; in PCGAMGProlongator_Classical_Standard()
667 PetscCall(MatGetRow(lA, ci, &ncols, &icol, NULL)); in PCGAMGProlongator_Classical_Standard()
679 PetscCall(MatRestoreRow(lA, ci, &ncols, &icol, NULL)); in PCGAMGProlongator_Classical_Standard()
713 ci = icol[j]; in PCGAMGProlongator_Classical_Standard()
716 PetscCall(MatGetRow(lA, ci, &ncols, &icol, &vcol)); in PCGAMGProlongator_Classical_Standard()
720 if (ci == icol[k]) jdiag = PetscRealPart(vcol[k]); in PCGAMGProlongator_Classical_Standard()
[all …]
/petsc/src/tao/util/
H A Dtao_util.c252 PetscReal ai, bi, ci, di, ei; in MatDFischer() local
327 ci = 1.0 / ai + 1.0; in MatDFischer()
334 ci = bi / ai + 1.0; in MatDFischer()
339 bi = ci + di * PetscRealPart(t2[i]); in MatDFischer()
353 da[i] = ai + bi * ci; in MatDFischer()
398 PetscReal ai, bi, ci, di, ei, fi; in MatDSFischer() local
444 ci = bi / ai + 1.0; in MatDSFischer()
456 da[i] = ai + bi * ci; in MatDSFischer()
/petsc/src/dm/partitioner/impls/simple/
H A Dpartsimple.c50 PetscInt Np = 1, Nr, np, nk, nj, ni, pk, pj, pi, ck, cj, ci, i; in PetscPartitionerPartition_Simple_Grid() local
98 for (ci = 0; ci < pcells[0]; ++ci) { in PetscPartitionerPartition_Simple_Grid()
99 …s[1] + ((nj * procs[1] + pj) * pcells[1] + cj)) * cells[0] + (ni * procs[0] + pi) * pcells[0] + ci; in PetscPartitionerPartition_Simple_Grid()
/petsc/src/tao/constrained/impls/ipm/
H A Dipm.c284 PetscCall(VecDuplicate(ipmP->s, &ipmP->ci)); in IPMInitializeBounds()
310 PetscCall(VecScatterCreate(tao->constraints_inequality, isuc, ipmP->ci, isc, &ipmP->ci_scat)); in IPMInitializeBounds()
338 PetscCall(VecScatterCreate(tao->XL, isx, ipmP->ci, isc, &ipmP->xl_scat)); in IPMInitializeBounds()
365 PetscCall(VecScatterCreate(tao->XU, isx, ipmP->ci, isc, &ipmP->xu_scat)); in IPMInitializeBounds()
464 PetscCall(VecDestroy(&ipmP->ci)); in TaoDestroy_IPM()
586 PetscCall(VecCopy(ipmP->ci, ipmP->rpi)); in IPMComputeKKT()
756 PetscCall(VecSet(ipmP->ci, 0.0)); in IPMUpdateAi()
760 …PetscCall(VecScatterBegin(ipmP->ci_scat, tao->constraints_inequality, ipmP->ci, INSERT_VALUES, SCA… in IPMUpdateAi()
761 …PetscCall(VecScatterEnd(ipmP->ci_scat, tao->constraints_inequality, ipmP->ci, INSERT_VALUES, SCATT… in IPMUpdateAi()
770 … PetscCall(VecScatterBegin(ipmP->xl_scat, ipmP->work, ipmP->ci, INSERT_VALUES, SCATTER_FORWARD)); in IPMUpdateAi()
[all …]
H A Dipm.h19 Vec ci; member
H A Dpdipm.h35 Vec ce, ci; /* equality and inequality constraints */ member
H A Dpdipm.c87 PetscCall(VecGetArrayWrite(pdipm->ci, &carr)); in TaoPDIPMUpdateConstraints()
127 PetscCall(VecRestoreArrayWrite(pdipm->ci, &carr)); in TaoPDIPMUpdateConstraints()
513 PetscCall(VecGetArrayRead(pdipm->ci, &carr)); in TaoSNESFunction_PDIPM()
520 PetscCall(VecRestoreArrayRead(pdipm->ci, &carr)); in TaoSNESFunction_PDIPM()
524 PetscCall(VecGetArrayRead(pdipm->ci, &carr)); in TaoSNESFunction_PDIPM()
531 PetscCall(VecRestoreArrayRead(pdipm->ci, &carr)); in TaoSNESFunction_PDIPM()
589 PetscCall(VecPlaceArray(pdipm->ci, Farr + pdipm->off_lambdai)); in TaoSNESFunction_PDIPM_residual()
590 PetscCall(VecNorm(pdipm->ci, NORM_2, &cnorm[1])); in TaoSNESFunction_PDIPM_residual()
591 PetscCall(VecResetArray(pdipm->ci)); in TaoSNESFunction_PDIPM_residual()
905 PetscCall(VecCreate(comm, &pdipm->ci)); in TaoSetup_PDIPM()
[all …]
/petsc/src/mat/impls/maij/
H A Dmaij.c592 const PetscInt *ci = c->i, *cj = c->j, *cjj; in MatPtAPNumeric_SeqAIJ_SeqMAIJ() local
603 PetscCall(PetscArrayzero(ca, ci[cm])); in MatPtAPNumeric_SeqAIJ_SeqMAIJ()
643 cjj = cj + ci[crow]; in MatPtAPNumeric_SeqAIJ_SeqMAIJ()
644 caj = ca + ci[crow]; in MatPtAPNumeric_SeqAIJ_SeqMAIJ()
675 PetscInt *ci, *cj, *ptadenserow, *ptasparserow, *denserow, *sparserow, *ptaj; in MatPtAPSymbolic_SeqAIJ_SeqMAIJ() local
688 PetscCall(PetscMalloc1(cn + 1, &ci)); in MatPtAPSymbolic_SeqAIJ_SeqMAIJ()
689 ci[0] = 0; in MatPtAPSymbolic_SeqAIJ_SeqMAIJ()
762 ci[i * ppdof + 1 + dof] = ci[i * ppdof + dof] + cnzi; in MatPtAPSymbolic_SeqAIJ_SeqMAIJ()
768 PetscCall(PetscMalloc1(ci[cn], &cj)); in MatPtAPSymbolic_SeqAIJ_SeqMAIJ()
773 PetscCall(PetscCalloc1(ci[cn], &ca)); in MatPtAPSymbolic_SeqAIJ_SeqMAIJ()
[all …]
/petsc/src/ksp/pc/impls/ml/
H A Dml.c176 PetscInt *ci, *cj, ncols; in MatConvert_MPIAIJ_ML() local
183 PetscCall(PetscMalloc1(1 + am, &ci)); in MatConvert_MPIAIJ_ML()
184 ci[0] = 0; in MatConvert_MPIAIJ_ML()
185 for (i = 0; i < am; i++) ci[i + 1] = ci[i] + (ai[i + 1] - ai[i]) + (bi[i + 1] - bi[i]); in MatConvert_MPIAIJ_ML()
186 PetscCall(PetscMalloc1(1 + ci[am], &cj)); in MatConvert_MPIAIJ_ML()
187 PetscCall(PetscMalloc1(1 + ci[am], &ca)); in MatConvert_MPIAIJ_ML()
205 PetscCheck(k == ci[am], PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "k: %d != ci[am]: %d", k, ci[am]); in MatConvert_MPIAIJ_ML()
209 PetscCall(MatCreateSeqAIJWithArrays(PETSC_COMM_SELF, am, an, ci, cj, ca, Aloc)); in MatConvert_MPIAIJ_ML()
220 ci = mat->i; in MatConvert_MPIAIJ_ML()
/petsc/src/mat/impls/aij/seq/spqr/
H A Daijspqr.c22 SuiteSparse_long *ci, *cj; /* SuiteSparse_long is the only choice for SPQR */ in MatWrapCholmod_SPQR_seqaij() local
58 PetscCall(PetscMalloc2(n + 1, &cj, nz, &ci)); in MatWrapCholmod_SPQR_seqaij()
67 ci[k] = ai[i]; in MatWrapCholmod_SPQR_seqaij()
88 C->i = ci; in MatWrapCholmod_SPQR_seqaij()
/petsc/doc/developers/contributing/
H A Dpipelines.md5 PETSc uses [GitLab Pipelines](https://docs.gitlab.com/ee/ci/pipelines/) for testing during continuo…
128 to your branch in the PETSc Discord channel `testing-ci-forum`.
/petsc/src/ksp/pc/impls/patch/
H A Dpcpatch.c71 PetscInt *closure = NULL, ci; in PCPatchConstruct_Vanka() local
76 for (ci = 0; ci < closureSize * 2; ci += 2) { in PCPatchConstruct_Vanka()
77 const PetscInt newpoint = closure[ci]; in PCPatchConstruct_Vanka()
147 PetscInt *closure = NULL, ci; in PCPatchConstruct_Pardecomp() local
151 for (ci = 0; ci < closureSize * 2; ci += 2) { in PCPatchConstruct_Pardecomp()
152 const PetscInt closstarp = closure[ci]; in PCPatchConstruct_Pardecomp()
743 PetscInt ignoredim, iStart = 0, iEnd = -1, starSize, closureSize, si, ci; in PCPatchCompleteCellPatch() local
767 for (ci = 0; ci < closureSize * 2; ci += 2) { in PCPatchCompleteCellPatch()
768 const PetscInt seenpoint = closure[ci]; in PCPatchCompleteCellPatch()
/petsc/src/mat/tests/
H A Dex23.c305 …PetscInt ri, ci, rr[3] = {0, 1, 0}, cr[4] = {1, 2, 0, 1}, rk[3] = {0, 2, 1}, ck[4] =… in main() local
312 for (ci = 0; ci < 2; ci++) { in main()
315 c = (PetscInt *)(ci == 0 ? cr : ck); in main()
324 …e (%" PetscInt_FMT " %" PetscInt_FMT ", bs %" PetscInt_FMT " %" PetscInt_FMT ")", ri, ci, rb, cb)); in main()
/petsc/src/mat/impls/aij/seq/kokkos/
H A Daijkok.kokkos.cxx716 MatRowMapKokkosView ai, bi, ci; in MatSeqAIJKokkosMergeMats() local
747 auto ci = MatRowMapKokkosView("i", ai.extent(0)); in MatSeqAIJKokkosMergeMats() local
757 ci(i) = coffset; in MatSeqAIJKokkosMergeMats()
758 if (i == m - 1) ci(m) = ai(m) + bi(m); in MatSeqAIJKokkosMergeMats()
771 PetscCallCXX(ckok = new Mat_SeqAIJKokkos(m, n, nnz, ci, cj, ca)); in MatSeqAIJKokkosMergeMats()
778 ci = ckok->i_dual.view_device(); in MatSeqAIJKokkosMergeMats()
785 if (k < alen) ca(ci(i) + k) = aa(ai(i) + k); in MatSeqAIJKokkosMergeMats()
786 else ca(ci(i) + k) = ba(bi(i) + k - alen); in MatSeqAIJKokkosMergeMats()

12