Home
last modified time | relevance | path

Searched refs:coloring (Results 1 – 25 of 45) sorted by relevance

12

/petsc/src/mat/impls/aij/mpi/
H A Dfdmpiaij.c31 PetscErrorCode MatFDColoringApply_BAIJ(Mat J, MatFDColoring coloring, Vec x1, void *sctx) in MatFDColoringApply_BAIJ() argument
33 …rorCode (*f)(void *, Vec, Vec, void *) = (PetscErrorCode (*)(void *, Vec, Vec, void *))coloring->f; in MatFDColoringApply_BAIJ()
35 PetscScalar dx = 0.0, *w3_array, *dy_i, *dy = coloring->dy; in MatFDColoringApply_BAIJ()
38 PetscReal epsilon = coloring->error_rel, umin = coloring->umin, unorm; in MatFDColoringApply_BAIJ()
39 Vec w1 = coloring->w1, w2 = coloring->w2, w3, vscale = coloring->vscale; in MatFDColoringApply_BAIJ()
40 void *fctx = coloring->fctx; in MatFDColoringApply_BAIJ()
41 PetscInt ctype = coloring->ctype, nxloc, nrows_k; in MatFDColoringApply_BAIJ()
43 MatEntry *Jentry = coloring->matentry; in MatFDColoringApply_BAIJ()
44 MatEntry2 *Jentry2 = coloring->matentry2; in MatFDColoringApply_BAIJ()
45 …const PetscInt ncolors = coloring->ncolors, *ncolumns = coloring->ncolumns, *nrows = coloring-… in MatFDColoringApply_BAIJ()
[all …]
/petsc/src/mat/tests/
H A Dex24.c30 MatColoring coloring; in main() local
48 PetscCall(MatColoringCreate(J, &coloring)); in main()
49 PetscCall(MatColoringSetType(coloring, MATCOLORINGGREEDY)); in main()
50 PetscCall(MatColoringSetFromOptions(coloring)); in main()
51 PetscCall(MatColoringApply(coloring, &iscoloring)); in main()
56 PetscCall(MatColoringDestroy(&coloring)); in main()
H A Dex199.c12 ISColoring coloring; in main() local
28 PetscCall(MatColoringApply(ctx, &coloring)); in main()
29 PetscCall(MatColoringTest(ctx, coloring)); in main()
32 PetscCall(MatISColoringTest(C, coloring)); in main()
36 PetscCall(ISColoringDestroy(&coloring)); in main()
/petsc/src/mat/matfd/
H A Dfdmatrix.c545 PetscErrorCode MatFDColoringGetPerturbedColumns(MatFDColoring coloring, PetscInt *n, const PetscInt… in MatFDColoringGetPerturbedColumns() argument
548 if (coloring->currentcolor >= 0) { in MatFDColoringGetPerturbedColumns()
549 *n = coloring->ncolumns[coloring->currentcolor]; in MatFDColoringGetPerturbedColumns()
550 *cols = coloring->columns[coloring->currentcolor]; in MatFDColoringGetPerturbedColumns()
579 PetscErrorCode MatFDColoringApply(Mat J, MatFDColoring coloring, Vec x1, void *sctx) in MatFDColoringApply() argument
585 PetscValidHeaderSpecific(coloring, MAT_FDCOLORING_CLASSID, 2); in MatFDColoringApply()
587 PetscCall(PetscObjectCompareId((PetscObject)J, coloring->matid, &eq)); in MatFDColoringApply()
589 …PetscCheck(coloring->f, PetscObjectComm((PetscObject)J), PETSC_ERR_ARG_WRONGSTATE, "Must call MatF… in MatFDColoringApply()
590 …PetscCheck(coloring->setupcalled, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONGSTATE, "Must call MatFDColor… in MatFDColoringApply()
593 PetscCall(PetscLogEventBegin(MAT_FDColoringApply, coloring, J, x1, 0)); in MatFDColoringApply()
[all …]
/petsc/src/snes/impls/gs/
H A Dgssecant.c9 ISColoring coloring = gs->coloring; in SNESComputeNGSDefaultSecant() local
35 if (!coloring) { in SNESComputeNGSDefaultSecant()
39 PetscCall(DMCreateColoring(dm, IS_COLORING_GLOBAL, &coloring)); in SNESComputeNGSDefaultSecant()
45 PetscCall(MatColoringApply(mc, &coloring)); in SNESComputeNGSDefaultSecant()
48 gs->coloring = coloring; in SNESComputeNGSDefaultSecant()
50 PetscCall(ISColoringGetIS(coloring, PETSC_USE_POINTER, &ncolors, &coloris)); in SNESComputeNGSDefaultSecant()
122 PetscCall(ISColoringRestoreIS(coloring, PETSC_USE_POINTER, &coloris)); in SNESComputeNGSDefaultSecant()
H A Dgsimpl.h14 ISColoring coloring; member
/petsc/src/mat/graphops/color/impls/minpack/
H A Dcolor.c29 PetscInt *list, *work, clique, *seq, *coloring, n; in MatColoringApply_SL() local
66 PetscCall(PetscMalloc1(n, &coloring)); in MatColoringApply_SL()
67 PetscCall(MINPACKseq(&n, cja, cia, rja, ria, list, coloring, &ncolors, work)); in MatColoringApply_SL()
77 ISColoringValue *s = (ISColoringValue *)coloring; in MatColoringApply_SL()
78 for (i = 0; i < n; i++) s[i] = (ISColoringValue)(coloring[i] - 1); in MatColoringApply_SL()
129 PetscInt *list, *work, *seq, *coloring, n; in MatColoringApply_LF() local
167 PetscCall(PetscMalloc1(n, &coloring)); in MatColoringApply_LF()
168 PetscCall(MINPACKseq(&n, cja, cia, rja, ria, list, coloring, &ncolors, work)); in MatColoringApply_LF()
179 ISColoringValue *s = (ISColoringValue *)coloring; in MatColoringApply_LF()
180 for (i = 0; i < n; i++) s[i] = (ISColoringValue)(coloring[i] - 1); in MatColoringApply_LF()
[all …]
/petsc/src/mat/graphops/color/interface/
H A Dmatcoloring.c350 PetscErrorCode MatColoringApply(MatColoring mc, ISColoring *coloring) in MatColoringApply() argument
359 PetscAssertPointer(coloring, 2); in MatColoringApply()
361 PetscUseTypeMethod(mc, apply, coloring); in MatColoringApply()
365 if (mc->valid) PetscCall(MatColoringTest(mc, *coloring)); in MatColoringApply()
366 if (mc->valid_iscoloring) PetscCall(MatISColoringTest(mc->mat, *coloring)); in MatColoringApply()
374 PetscCall(ISColoringGetIS(*coloring, PETSC_USE_POINTER, &ncolors, NULL)); in MatColoringApply()
377 if (nc <= 1000) PetscCall(ISColoringView(*coloring, viewer)); in MatColoringApply()
/petsc/src/tao/interface/
H A Dfdiff.c166 MatFDColoring coloring = (MatFDColoring)ctx; in TaoDefaultComputeHessianColor() local
169 PetscValidHeaderSpecific(coloring, MAT_FDCOLORING_CLASSID, 5); in TaoDefaultComputeHessianColor()
171 PetscCall(MatFDColoringApply(B, coloring, V, ctx)); in TaoDefaultComputeHessianColor()
/petsc/src/mat/impls/aij/seq/
H A Dmatrart.c29 MatColoring coloring; in MatRARtSymbolic_SeqAIJ_SeqAIJ_colorrart() local
55 PetscCall(MatColoringCreate(C, &coloring)); in MatRARtSymbolic_SeqAIJ_SeqAIJ_colorrart()
56 PetscCall(MatColoringSetDistance(coloring, 2)); in MatRARtSymbolic_SeqAIJ_SeqAIJ_colorrart()
57 PetscCall(MatColoringSetType(coloring, MATCOLORINGSL)); in MatRARtSymbolic_SeqAIJ_SeqAIJ_colorrart()
58 PetscCall(MatColoringSetFromOptions(coloring)); in MatRARtSymbolic_SeqAIJ_SeqAIJ_colorrart()
59 PetscCall(MatColoringApply(coloring, &iscoloring)); in MatRARtSymbolic_SeqAIJ_SeqAIJ_colorrart()
60 PetscCall(MatColoringDestroy(&coloring)); in MatRARtSymbolic_SeqAIJ_SeqAIJ_colorrart()
H A Dmatmatmult.c1249 MatColoring coloring; in MatMatTransposeMultSymbolic_SeqAIJ_SeqAIJ() local
1256 PetscCall(MatColoringCreate(C, &coloring)); in MatMatTransposeMultSymbolic_SeqAIJ_SeqAIJ()
1257 PetscCall(MatColoringSetDistance(coloring, 2)); in MatMatTransposeMultSymbolic_SeqAIJ_SeqAIJ()
1258 PetscCall(MatColoringSetType(coloring, MATCOLORINGSL)); in MatMatTransposeMultSymbolic_SeqAIJ_SeqAIJ()
1259 PetscCall(MatColoringSetFromOptions(coloring)); in MatMatTransposeMultSymbolic_SeqAIJ_SeqAIJ()
1260 PetscCall(MatColoringApply(coloring, &iscoloring)); in MatMatTransposeMultSymbolic_SeqAIJ_SeqAIJ()
1261 PetscCall(MatColoringDestroy(&coloring)); in MatMatTransposeMultSymbolic_SeqAIJ_SeqAIJ()
1752 PetscErrorCode MatTransColoringApplySpToDen_SeqAIJ(MatTransposeColoring coloring, Mat B, Mat Btdens… in MatTransColoringApplySpToDen_SeqAIJ() argument
1759 …PetscInt *columns = coloring->columns, *colorforcol = coloring->colorforcol, ncolors = colorin… in MatTransColoringApplySpToDen_SeqAIJ()
1765 ncolumns = coloring->ncolumns[k]; in MatTransColoringApplySpToDen_SeqAIJ()
/petsc/src/snes/utils/
H A Ddmlocalsnes.c154 ISColoring coloring; in SNESComputeJacobian_DMLocal() local
156 PetscCall(DMCreateColoring(dm, dm->coloringtype, &coloring)); in SNESComputeJacobian_DMLocal()
157 PetscCall(MatFDColoringCreate(B, coloring, &fdcoloring)); in SNESComputeJacobian_DMLocal()
158 PetscCall(ISColoringDestroy(&coloring)); in SNESComputeJacobian_DMLocal()
168 PetscCall(MatFDColoringSetUp(B, coloring, fdcoloring)); in SNESComputeJacobian_DMLocal()
H A Ddmdasnes.c175 ISColoring coloring; in SNESComputeJacobian_DMDA() local
177 PetscCall(DMCreateColoring(dm, dm->coloringtype, &coloring)); in SNESComputeJacobian_DMDA()
178 PetscCall(MatFDColoringCreate(B, coloring, &fdcoloring)); in SNESComputeJacobian_DMDA()
188 PetscCall(MatFDColoringSetUp(B, coloring, fdcoloring)); in SNESComputeJacobian_DMDA()
189 PetscCall(ISColoringDestroy(&coloring)); in SNESComputeJacobian_DMDA()
/petsc/src/ts/utils/
H A Ddmlocalts.c155 ISColoring coloring; in TSComputeIJacobian_DMLocal() local
157 PetscCall(DMCreateColoring(dm, dm->coloringtype, &coloring)); in TSComputeIJacobian_DMLocal()
158 PetscCall(MatFDColoringCreate(B, coloring, &fdcoloring)); in TSComputeIJacobian_DMLocal()
159 PetscCall(ISColoringDestroy(&coloring)); in TSComputeIJacobian_DMLocal()
169 PetscCall(MatFDColoringSetUp(B, coloring, fdcoloring)); in TSComputeIJacobian_DMLocal()
/petsc/src/mat/ftn-mod/
H A Dpetscmat.h902 subroutine MatFDColoringRestorePerturbedColumns(coloring, n, cols, ierr)
8 MatFDColoring coloring
/petsc/src/snes/tutorials/output/
H A Dex19_ngs_fd.out18 Use finite difference secant approximation with coloring with h = 1e-08
/petsc/src/vec/is/is/utils/
H A Discoloring.c7 PetscErrorCode ISColoringReference(ISColoring coloring) in ISColoringReference() argument
10 coloring->refct++; in ISColoringReference()
32 PetscErrorCode ISColoringSetType(ISColoring coloring, ISColoringType type) in ISColoringSetType() argument
35 coloring->ctype = type; in ISColoringSetType()
54 PetscErrorCode ISColoringGetType(ISColoring coloring, ISColoringType *type) in ISColoringGetType() argument
57 *type = coloring->ctype; in ISColoringGetType()
/petsc/src/mat/graphops/color/utils/
H A Dvalid.c6 PETSC_EXTERN PetscErrorCode MatColoringTest(MatColoring mc, ISColoring coloring) in MatColoringTest() argument
26 PetscCall(ISColoringGetIS(coloring, PETSC_USE_POINTER, &ncolors, &colors)); in MatColoringTest()
/petsc/src/ts/tests/output/
H A Dex21_1_bdf_ngmres_fas_ngs.out104 Use finite difference secant approximation with coloring with h = 1.49012e-08
112 Use finite difference secant approximation with coloring with h = 1.49012e-08
H A Dex21_2_bdf_ngmres_fas_ngs.out132 Use finite difference secant approximation with coloring with h = 1.49012e-08
140 Use finite difference secant approximation with coloring with h = 1.49012e-08
/petsc/doc/changes/
H A D212.md18 - Allow coloring of periodic grids
/petsc/src/dm/impls/da/
H A Dfdda.c187 PetscErrorCode DMCreateColoring_DA(DM da, ISColoringType ctype, ISColoring *coloring) in DMCreateColoring_DA() argument
245 if (dim == 1) PetscCall(DMCreateColoring_DA_1d_MPIAIJ(da, ctype, coloring)); in DMCreateColoring_DA()
246 else if (dim == 2) PetscCall(DMCreateColoring_DA_2d_MPIAIJ(da, ctype, coloring)); in DMCreateColoring_DA()
247 else if (dim == 3) PetscCall(DMCreateColoring_DA_3d_MPIAIJ(da, ctype, coloring)); in DMCreateColoring_DA()
259 PetscErrorCode DMCreateColoring_DA_2d_MPIAIJ(DM da, ISColoringType ctype, ISColoring *coloring) in DMCreateColoring_DA_2d_MPIAIJ() argument
283 PetscCall(DMCreateColoring_DA_2d_5pt_MPIAIJ(da, ctype, coloring)); in DMCreateColoring_DA_2d_MPIAIJ()
297 *coloring = dd->localcoloring; in DMCreateColoring_DA_2d_MPIAIJ()
316 *coloring = dd->ghostedcoloring; in DMCreateColoring_DA_2d_MPIAIJ()
319 PetscCall(ISColoringReference(*coloring)); in DMCreateColoring_DA_2d_MPIAIJ()
323 PetscErrorCode DMCreateColoring_DA_3d_MPIAIJ(DM da, ISColoringType ctype, ISColoring *coloring) in DMCreateColoring_DA_3d_MPIAIJ() argument
[all …]
/petsc/src/snes/tutorials/
H A Dex14.c65 …PetscBool matrix_free = PETSC_FALSE, coloring = PETSC_FALSE, coloring_ds = PETSC_FALSE, local_… in main() local
121 PetscCall(PetscOptionsGetBool(NULL, NULL, "-fdcoloring", &coloring, NULL)); in main()
127 if (coloring) { in main()
/petsc/src/dm/impls/redundant/
H A Ddmredundant.c175 static PetscErrorCode DMCreateColoring_Redundant(DM dm, ISColoringType ctype, ISColoring *coloring) in DMCreateColoring_Redundant() argument
194 …loringCreate(PetscObjectComm((PetscObject)dm), red->N, nloc, colors, PETSC_OWN_POINTER, coloring)); in DMCreateColoring_Redundant()
195 PetscCall(ISColoringSetType(*coloring, ctype)); in DMCreateColoring_Redundant()
/petsc/src/ts/tutorials/advection-diffusion-reaction/output/
H A Dex4_1.out85 Jacobian is built using finite differences with coloring

12