Home
last modified time | relevance | path

Searched refs:ic (Results 1 – 25 of 36) sorted by relevance

12

/petsc/src/mat/graphops/color/impls/minpack/
H A Dido.c13 PetscInt jcol = 0, ncomp = 0, ic, ip, jp, ir, maxinc, numinc, numord, maxlst, numwgt, numlst; in MINPACKido() local
105 ic = iwa4[jp]; in MINPACKido()
107 iwa2[ic] = 0; in MINPACKido()
108 iwa3[ic] = iwa1[0]; in MINPACKido()
109 if (iwa1[0] > 0) iwa2[iwa1[0]] = ic; in MINPACKido()
110 iwa1[0] = ic; in MINPACKido()
190 ic = indcol[ip]; in MINPACKido()
195 if (iwa4[ic] < numord) { in MINPACKido()
196 iwa4[ic] = numord; in MINPACKido()
200 numinc = list[ic]; in MINPACKido()
[all …]
H A Dslo.c11 PetscInt jcol, ic, ip, jp, ir, mindeg, numdeg, numord; in MINPACKslo() local
165 ic = indcol[ip]; in MINPACKslo()
170 if (iwa4[ic] > numord) { in MINPACKslo()
171 iwa4[ic] = numord; in MINPACKslo()
175 numdeg = list[ic]; in MINPACKslo()
176 --list[ic]; in MINPACKslo()
179 i__4 = list[ic]; in MINPACKslo()
184 if (!iwa2[ic]) iwa1[numdeg] = iwa3[ic]; in MINPACKslo()
185 else iwa3[iwa2[ic]] = iwa3[ic]; in MINPACKslo()
187 if (iwa3[ic] > 0) iwa2[iwa3[ic]] = iwa2[ic]; in MINPACKslo()
[all …]
H A Ddegr.c11 PetscInt jcol, ic, ip, jp, ir; in MINPACKdegr() local
88 ic = indcol[ip]; in MINPACKdegr()
94 if (iwa[ic] < jcol) { in MINPACKdegr()
95 iwa[ic] = jcol; in MINPACKdegr()
96 ++ndeg[ic]; in MINPACKdegr()
H A Dseq.c11 PetscInt jcol, j, ic, ip, jp, ir; in MINPACKseq() local
101 ic = indcol[ip]; in MINPACKseq()
106 iwa[ngrp[ic]] = j; in MINPACKseq()
/petsc/src/mat/impls/baij/seq/
H A Dbaijsolvtran5.c10 PetscInt i, nz, idx, idt, ii, ic, ir, oidx; in MatSolveTranspose_SeqBAIJ_5_inplace() local
28 ic = 5 * c[i]; in MatSolveTranspose_SeqBAIJ_5_inplace()
29 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_5_inplace()
30 t[ii + 1] = b[ic + 1]; in MatSolveTranspose_SeqBAIJ_5_inplace()
31 t[ii + 2] = b[ic + 2]; in MatSolveTranspose_SeqBAIJ_5_inplace()
32 t[ii + 3] = b[ic + 3]; in MatSolveTranspose_SeqBAIJ_5_inplace()
33 t[ii + 4] = b[ic + 4]; in MatSolveTranspose_SeqBAIJ_5_inplace()
120 PetscInt nz, idx, idt, j, i, oidx, ii, ic, ir; in MatSolveTranspose_SeqBAIJ_5() local
139 ic = bs * c[i]; in MatSolveTranspose_SeqBAIJ_5()
140 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_5()
[all …]
H A Dbaijsolvtran7.c10 PetscInt i, nz, idx, idt, ii, ic, ir, oidx; in MatSolveTranspose_SeqBAIJ_7_inplace() local
28 ic = 7 * c[i]; in MatSolveTranspose_SeqBAIJ_7_inplace()
29 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_7_inplace()
30 t[ii + 1] = b[ic + 1]; in MatSolveTranspose_SeqBAIJ_7_inplace()
31 t[ii + 2] = b[ic + 2]; in MatSolveTranspose_SeqBAIJ_7_inplace()
32 t[ii + 3] = b[ic + 3]; in MatSolveTranspose_SeqBAIJ_7_inplace()
33 t[ii + 4] = b[ic + 4]; in MatSolveTranspose_SeqBAIJ_7_inplace()
34 t[ii + 5] = b[ic + 5]; in MatSolveTranspose_SeqBAIJ_7_inplace()
35 t[ii + 6] = b[ic + 6]; in MatSolveTranspose_SeqBAIJ_7_inplace()
135 PetscInt nz, idx, idt, j, i, oidx, ii, ic, ir; in MatSolveTranspose_SeqBAIJ_7() local
[all …]
H A Dbaijsolvtran4.c10 PetscInt i, nz, idx, idt, ii, ic, ir, oidx; in MatSolveTranspose_SeqBAIJ_4_inplace() local
28 ic = 4 * c[i]; in MatSolveTranspose_SeqBAIJ_4_inplace()
29 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_4_inplace()
30 t[ii + 1] = b[ic + 1]; in MatSolveTranspose_SeqBAIJ_4_inplace()
31 t[ii + 2] = b[ic + 2]; in MatSolveTranspose_SeqBAIJ_4_inplace()
32 t[ii + 3] = b[ic + 3]; in MatSolveTranspose_SeqBAIJ_4_inplace()
112 PetscInt nz, idx, idt, j, i, oidx, ii, ic, ir; in MatSolveTranspose_SeqBAIJ_4() local
131 ic = bs * c[i]; in MatSolveTranspose_SeqBAIJ_4()
132 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_4()
133 t[ii + 1] = b[ic + 1]; in MatSolveTranspose_SeqBAIJ_4()
[all …]
H A Dbaijsolvtran6.c10 PetscInt i, nz, idx, idt, ii, ic, ir, oidx; in MatSolveTranspose_SeqBAIJ_6_inplace() local
28 ic = 6 * c[i]; in MatSolveTranspose_SeqBAIJ_6_inplace()
29 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_6_inplace()
30 t[ii + 1] = b[ic + 1]; in MatSolveTranspose_SeqBAIJ_6_inplace()
31 t[ii + 2] = b[ic + 2]; in MatSolveTranspose_SeqBAIJ_6_inplace()
32 t[ii + 3] = b[ic + 3]; in MatSolveTranspose_SeqBAIJ_6_inplace()
33 t[ii + 4] = b[ic + 4]; in MatSolveTranspose_SeqBAIJ_6_inplace()
34 t[ii + 5] = b[ic + 5]; in MatSolveTranspose_SeqBAIJ_6_inplace()
128 PetscInt nz, idx, idt, j, i, oidx, ii, ic, ir; in MatSolveTranspose_SeqBAIJ_6() local
147 ic = bs * c[i]; in MatSolveTranspose_SeqBAIJ_6()
[all …]
H A Dbaijsolvtran3.c10 PetscInt i, nz, idx, idt, ii, ic, ir, oidx; in MatSolveTranspose_SeqBAIJ_3_inplace() local
28 ic = 3 * c[i]; in MatSolveTranspose_SeqBAIJ_3_inplace()
29 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_3_inplace()
30 t[ii + 1] = b[ic + 1]; in MatSolveTranspose_SeqBAIJ_3_inplace()
31 t[ii + 2] = b[ic + 2]; in MatSolveTranspose_SeqBAIJ_3_inplace()
104 PetscInt nz, idx, idt, j, i, oidx, ii, ic, ir; in MatSolveTranspose_SeqBAIJ_3() local
123 ic = bs * c[i]; in MatSolveTranspose_SeqBAIJ_3()
124 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_3()
125 t[ii + 1] = b[ic + 1]; in MatSolveTranspose_SeqBAIJ_3()
126 t[ii + 2] = b[ic + 2]; in MatSolveTranspose_SeqBAIJ_3()
H A Dbaijsolvtran2.c10 PetscInt i, nz, idx, idt, ii, ic, ir, oidx; in MatSolveTranspose_SeqBAIJ_2_inplace() local
28 ic = 2 * c[i]; in MatSolveTranspose_SeqBAIJ_2_inplace()
29 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_2_inplace()
30 t[ii + 1] = b[ic + 1]; in MatSolveTranspose_SeqBAIJ_2_inplace()
96 PetscInt nz, idx, idt, j, i, oidx, ii, ic, ir; in MatSolveTranspose_SeqBAIJ_2() local
115 ic = bs * c[i]; in MatSolveTranspose_SeqBAIJ_2()
116 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_2()
117 t[ii + 1] = b[ic + 1]; in MatSolveTranspose_SeqBAIJ_2()
H A Dbaijfact4.c11 const PetscInt *r, *ic; in MatILUFactorNumeric_SeqBAIJ_N_inplace() local
25 PetscCall(ISGetIndices(isicol, &ic)); in MatILUFactorNumeric_SeqBAIJ_N_inplace()
40 … for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ic[ajtmpold[j]], v + bs2 * j, bs2)); in MatILUFactorNumeric_SeqBAIJ_N_inplace()
77 PetscCall(ISRestoreIndices(isicol, &ic)); in MatILUFactorNumeric_SeqBAIJ_N_inplace()
H A Dbaijfact3.c156 const PetscInt *r, *ic; in MatLUFactorSymbolic_SeqBAIJ() local
177 PetscCall(ISGetIndices(isicol, &ic)); in MatLUFactorSymbolic_SeqBAIJ()
201 PetscCall(PetscLLAddPerm(nnz, ajtmp, ic, n, &nlnk, lnk, lnkbt)); in MatLUFactorSymbolic_SeqBAIJ()
243 PetscCall(ISRestoreIndices(isicol, &ic)); in MatLUFactorSymbolic_SeqBAIJ()
H A Dbaijfact2.c119 const PetscInt *r, *ic; in MatLUFactorNumeric_SeqBAIJ_N() local
130 PetscCall(ISGetIndices(isicol, &ic)); in MatLUFactorNumeric_SeqBAIJ_N()
154 for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ic[ajtmp[j]], v + bs2 * j, bs2)); in MatLUFactorNumeric_SeqBAIJ_N()
203 PetscCall(ISRestoreIndices(isicol, &ic)); in MatLUFactorNumeric_SeqBAIJ_N()
315 const PetscInt *r, *ic; in MatILUFactorSymbolic_SeqBAIJ() local
372 PetscCall(ISGetIndices(isicol, &ic)); in MatILUFactorSymbolic_SeqBAIJ()
400 PetscCall(PetscIncompleteLLInit(nnz, cols, n, ic, &nlnk, lnk, lnk_lvl, lnkbt)); in MatILUFactorSymbolic_SeqBAIJ()
458 PetscCall(ISRestoreIndices(isicol, &ic)); in MatILUFactorSymbolic_SeqBAIJ()
H A Dbaijfact13.c14 const PetscInt *r, *ic; in MatILUFactorNumeric_SeqBAIJ_3_inplace() local
32 PetscCall(ISGetIndices(isicol, &ic)); in MatILUFactorNumeric_SeqBAIJ_3_inplace()
49 x = rtmp + 9 * ic[ajtmpold[j]]; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
150 PetscCall(ISRestoreIndices(isicol, &ic)); in MatILUFactorNumeric_SeqBAIJ_3_inplace()
172 const PetscInt *r, *ic; in MatLUFactorNumeric_SeqBAIJ_3() local
183 PetscCall(ISGetIndices(isicol, &ic)); in MatLUFactorNumeric_SeqBAIJ_3()
206 for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ic[ajtmp[j]], v + bs2 * j, bs2)); in MatLUFactorNumeric_SeqBAIJ_3()
261 PetscCall(ISRestoreIndices(isicol, &ic)); in MatLUFactorNumeric_SeqBAIJ_3()
H A Dbaijfact.c12 const PetscInt *r, *ic; in MatLUFactorNumeric_SeqBAIJ_2() local
22 PetscCall(ISGetIndices(isicol, &ic)); in MatLUFactorNumeric_SeqBAIJ_2()
45 …for (PetscInt j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ic[ajtmp[j]], v + bs2 * j, … in MatLUFactorNumeric_SeqBAIJ_2()
106 PetscCall(ISRestoreIndices(isicol, &ic)); in MatLUFactorNumeric_SeqBAIJ_2()
228 const PetscInt *r, *ic; in MatILUFactorNumeric_SeqBAIJ_2_inplace() local
246 PetscCall(ISGetIndices(isicol, &ic)); in MatILUFactorNumeric_SeqBAIJ_2_inplace()
262 x = rtmp + 4 * ic[ajtmpold[j]]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
324 PetscCall(ISRestoreIndices(isicol, &ic)); in MatILUFactorNumeric_SeqBAIJ_2_inplace()
454 const PetscInt *r, *ic, *ics; in MatLUFactorNumeric_SeqBAIJ_1() local
489 PetscCall(ISGetIndices(isicol, &ic)); in MatLUFactorNumeric_SeqBAIJ_1()
[all …]
/petsc/src/dm/impls/stag/
H A Dstagmulti.c454 PetscInt ir, ic[2]; in DMStagPopulateInterpolation1d_Internal() local
486 PetscCall(DMStagStencilToIndexLocal(dmc, dim, count, colc, ic)); in DMStagPopulateInterpolation1d_Internal()
487 PetscCall(MatSetValuesLocal(A, 1, &ir, count, ic, weight, INSERT_VALUES)); in DMStagPopulateInterpolation1d_Internal()
501 PetscCall(DMStagStencilToIndexLocal(dmc, dim, 1, colc, ic)); in DMStagPopulateInterpolation1d_Internal()
502 PetscCall(MatSetValuesLocal(A, 1, &ir, 1, ic, weight, INSERT_VALUES)); in DMStagPopulateInterpolation1d_Internal()
513 PetscInt ir, ic[4]; in DMStagPopulateInterpolation2d_Internal() local
552 PetscCall(DMStagStencilToIndexLocal(dmc, dim, count, colc, ic)); in DMStagPopulateInterpolation2d_Internal()
553 PetscCall(MatSetValuesLocal(A, 1, &ir, count, ic, weight, INSERT_VALUES)); in DMStagPopulateInterpolation2d_Internal()
586 PetscCall(DMStagStencilToIndexLocal(dmc, dim, count, colc, ic)); in DMStagPopulateInterpolation2d_Internal()
587 PetscCall(MatSetValuesLocal(A, 1, &ir, count, ic, weight, INSERT_VALUES)); in DMStagPopulateInterpolation2d_Internal()
[all …]
H A Dstagstencil.c308 PetscInt *ir, *ic; in DMStagMatGetValuesStencil() local
314 PetscCall(PetscMalloc2(nRow, &ir, nCol, &ic)); in DMStagMatGetValuesStencil()
316 PetscCall(DMStagStencilToIndexLocal(dm, dim, nCol, posCol, ic)); in DMStagMatGetValuesStencil()
317 PetscCall(MatGetValuesLocal(mat, nRow, ir, nCol, ic, val)); in DMStagMatGetValuesStencil()
318 PetscCall(PetscFree2(ir, ic)); in DMStagMatGetValuesStencil()
346 PetscInt *ir, *ic; in DMStagMatSetValuesStencil() local
351 PetscCall(PetscMalloc2(nRow, &ir, nCol, &ic)); in DMStagMatSetValuesStencil()
353 PetscCall(DMStagStencilToIndexLocal(dm, dm->dim, nCol, posCol, ic)); in DMStagMatSetValuesStencil()
354 PetscCall(MatSetValuesLocal(mat, nRow, ir, nCol, ic, val, insertMode)); in DMStagMatSetValuesStencil()
355 PetscCall(PetscFree2(ir, ic)); in DMStagMatSetValuesStencil()
/petsc/src/mat/impls/baij/mpi/ftn-custom/
H A Dzmpibaijf.c10 PETSC_EXTERN void matmpibaijgetseqbaij_(Mat *A, Mat *Ad, Mat *Ao, PetscInt *ic, size_t *iic, PetscE… in matmpibaijgetseqbaij_() argument
15 *iic = PetscIntAddressToFortran(ic, (PetscInt *)i); in matmpibaijgetseqbaij_()
/petsc/src/ksp/ksp/tutorials/
H A Dex56.c27 PetscInt m, nn, M, Istart, Iend, i, j, k, ii, jj, kk, ic, ne = 4, id; in main() local
102 for (i = Ni0, ic = 0; i < Ni1; i++) { in main()
109 for (jj = 0; jj < 3; jj++, ic++) { in main()
110 d_nnz[ic] = 3 * (27 - osz[nbc]); in main()
111 o_nnz[ic] = 3 * osz[nbc]; in main()
116 …PetscCheck(ic == m, PETSC_COMM_SELF, PETSC_ERR_PLIB, "ic %" PetscInt_FMT " does not equal m %" Pet… in main()
178 for (i = Ni0, ic = 0, ii = 0; i < Ni1; i++, ii++) { in main()
180 for (k = Nk0, kk = 0; k < Nk1; k++, kk++, ic++) { in main()
182 x = coords[3 * ic] = h * (PetscReal)i; in main()
183 y = coords[3 * ic + 1] = h * (PetscReal)j; in main()
[all …]
/petsc/src/ts/tutorials/autodiff/
H A Dex16opt_ic.cxx126 Vec ic, r; in main() local
230 PetscCall(MatCreateVecs(user.A, &ic, NULL)); in main()
231 PetscCall(VecGetArray(ic, &x_ptr)); in main()
234 PetscCall(VecRestoreArray(ic, &x_ptr)); in main()
236 PetscCall(TaoSetSolution(tao, ic)); in main()
265 PetscCall(VecDestroy(&ic)); in main()
/petsc/src/tao/unconstrained/tutorials/
H A Dspectraladjointassimilation.c66 Vec ic; member
165 PetscCall(VecDuplicate(u, &appctx.dat.ic)); in main()
243 PetscCall(InitialConditions(appctx.dat.ic, &appctx)); in main()
255 PetscCall(TaoSetSolution(tao, appctx.dat.ic)); in main()
269 PetscCall(VecDestroy(&appctx.dat.ic)); in main()
580 PetscErrorCode FormFunctionGradient(Tao tao, Vec ic, PetscReal *f, Vec G, PetscCtx ctx) in FormFunctionGradient() argument
589 PetscCall(VecCopy(ic, appctx->dat.curr_sol)); in FormFunctionGradient()
622 PetscCall(VecDuplicate(appctx->dat.ic, &temp)); in MonitorError()
623 PetscCall(VecWAXPY(temp, -1.0, appctx->dat.ic, appctx->dat.true_solution)); in MonitorError()
H A Dburgers_spectral.c59 Vec ic; member
154 PetscCall(VecDuplicate(u, &appctx.dat.ic)); in main()
237 PetscCall(InitialConditions(appctx.dat.ic, &appctx)); in main()
248 PetscCall(TaoSetSolution(tao, appctx.dat.ic)); in main()
261 PetscCall(VecDestroy(&appctx.dat.ic)); in main()
582 PetscCall(VecWAXPY(temp, -1.0, appctx->dat.ic, appctx->dat.true_solution)); in FormFunctionGradient()
610 PetscCall(VecDuplicate(appctx->dat.ic, &temp)); in MonitorError()
611 PetscCall(VecWAXPY(temp, -1.0, appctx->dat.ic, appctx->dat.true_solution)); in MonitorError()
H A Dadjointreads.m16 %ic - initial condition (the one to be optimized)
26 % plot(xg,ic,'g-','LineWidth',2,'Markersize',12);
104 % plot(xg,ic,'ks-','Markersize',6,'LineWidth',2); drawnow;
/petsc/src/ksp/pc/impls/bjacobi/bjkokkos/
H A Dbjkokkos.kokkos.cxx46 …t *glb_Aaj, const PetscScalar *glb_Aaa, const PetscInt *r, const PetscInt *ic, const PetscInt star… in MatMult() argument
49 int rowa = ic[rowb]; in MatMult()
62 …t *glb_Aaj, const PetscScalar *glb_Aaa, const PetscInt *r, const PetscInt *ic, const PetscInt star… in MatMultTranspose() argument
67 int rowa = ic[rowb]; in MatMultTranspose()
88 …t *glb_Aaj, const PetscScalar *glb_Aaa, const PetscInt *r, const PetscInt *ic, PetscScalar *work_s… in BJSolve_TFQMR() argument
159 int rowa = ic[rowb]; in BJSolve_TFQMR()
201 static_cast<void>(MatMult(team, glb_Aai, glb_Aaj, glb_Aaa, r, ic, start, end, T, V)); in BJSolve_TFQMR()
223 static_cast<void>(MatMult(team, glb_Aai, glb_Aaj, glb_Aaa, r, ic, start, end, T, AUQ)); in BJSolve_TFQMR()
302 static_cast<void>(MatMult(team, glb_Aai, glb_Aaj, glb_Aaa, r, ic, start, end, T, V)); in BJSolve_TFQMR()
315 int rowa = ic[rowb]; in BJSolve_TFQMR()
[all …]
/petsc/src/dm/impls/moab/
H A Ddmmbmg.cxx224 for (unsigned ic = 0; ic < connp.size(); ++ic) { in DMCreateInterpolation_Moab() local
227 …if (found.find(connp[ic]) != found.end()) continue; /* make sure we don't doubl… in DMCreateInterpolation_Moab()
228 …if (dmbp->vghost->find(connp[ic]) != dmbp->vghost->end()) onz[ldof]++; /* update out-of-proc onz */ in DMCreateInterpolation_Moab()
230 found.insert(connp[ic]); in DMCreateInterpolation_Moab()

12