| /petsc/src/snes/tutorials/ |
| H A D | ex70.c | 97 PetscErrorCode StokesSetupPC(Stokes *s, KSP ksp) in StokesSetupPC() argument
105 PetscCall(PCFieldSplitSetIS(pc, "0", s->isg[0])); in StokesSetupPC()
106 PetscCall(PCFieldSplitSetIS(pc, "1", s->isg[1])); in StokesSetupPC()
107 if (s->userPC) PetscCall(PCFieldSplitSetSchurPre(pc, PC_FIELDSPLIT_SCHUR_PRE_USER, s->myS)); in StokesSetupPC()
108 if (s->userKSP) { in StokesSetupPC()
111 PetscCall(KSPSetOperators(subksp[1], s->myS, s->myS)); in StokesSetupPC()
117 PetscErrorCode StokesWriteSolution(Stokes *s) in StokesWriteSolution() argument
128 PetscCall(VecGetArrayRead(s->x, &array)); in StokesWriteSolution()
131 for (j = 0; j < s->ny; j++) { in StokesWriteSolution()
132 for (i = 0; i < s->nx; i++) { in StokesWriteSolution()
[all …]
|
| /petsc/src/dm/impls/stag/tests/ |
| H A D | ex17.c | 37 DMStagStencil s; in main() local
38 s.c = 0; in main()
39 s.loc = DMSTAG_ELEMENT; in main()
40 PetscCall(DMStagCreateISFromStencils(dm, 1, &s, &is)); in main()
47 DMStagStencil s; in main() local
48 s.c = 0; in main()
49 s.loc = DMSTAG_RIGHT; in main()
50 PetscCall(DMStagCreateISFromStencils(dm, 1, &s, &is)); in main()
57 DMStagStencil s[2]; in main() local
58 s[0].c = 0; in main()
[all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijsolvnat14.c | 13 PetscScalar s[14]; in MatSolve_SeqBAIJ_14_NaturalOrdering() local
69 s[0] = x[idt]; in MatSolve_SeqBAIJ_14_NaturalOrdering()
70 s[1] = x[1 + idt]; in MatSolve_SeqBAIJ_14_NaturalOrdering()
71 s[2] = x[2 + idt]; in MatSolve_SeqBAIJ_14_NaturalOrdering()
72 s[3] = x[3 + idt]; in MatSolve_SeqBAIJ_14_NaturalOrdering()
73 s[4] = x[4 + idt]; in MatSolve_SeqBAIJ_14_NaturalOrdering()
74 s[5] = x[5 + idt]; in MatSolve_SeqBAIJ_14_NaturalOrdering()
75 s[6] = x[6 + idt]; in MatSolve_SeqBAIJ_14_NaturalOrdering()
76 s[7] = x[7 + idt]; in MatSolve_SeqBAIJ_14_NaturalOrdering()
77 s[8] = x[8 + idt]; in MatSolve_SeqBAIJ_14_NaturalOrdering()
[all …]
|
| H A D | baijsolvnat11.c | 13 PetscScalar s[11]; in MatSolve_SeqBAIJ_11_NaturalOrdering() local
63 s[0] = x[idt]; in MatSolve_SeqBAIJ_11_NaturalOrdering()
64 s[1] = x[1 + idt]; in MatSolve_SeqBAIJ_11_NaturalOrdering()
65 s[2] = x[2 + idt]; in MatSolve_SeqBAIJ_11_NaturalOrdering()
66 s[3] = x[3 + idt]; in MatSolve_SeqBAIJ_11_NaturalOrdering()
67 s[4] = x[4 + idt]; in MatSolve_SeqBAIJ_11_NaturalOrdering()
68 s[5] = x[5 + idt]; in MatSolve_SeqBAIJ_11_NaturalOrdering()
69 s[6] = x[6 + idt]; in MatSolve_SeqBAIJ_11_NaturalOrdering()
70 s[7] = x[7 + idt]; in MatSolve_SeqBAIJ_11_NaturalOrdering()
71 s[8] = x[8 + idt]; in MatSolve_SeqBAIJ_11_NaturalOrdering()
[all …]
|
| /petsc/src/mat/impls/baij/seq/ftn-kernels/ |
| H A D | fsolvebaij.F90 | 21 PetscScalar :: s(0:3)
40 s(0:3) = b(idx + 0:idx + 3)
44 …s(0) = s(0) - (a(ax + 0)*x(jdx + 0) + a(ax + 4)*x(jdx + 1) + a(ax + 8)*x(jdx + 2) + a(ax + 12)*x(j…
45 …s(1) = s(1) - (a(ax + 1)*x(jdx + 0) + a(ax + 5)*x(jdx + 1) + a(ax + 9)*x(jdx + 2) + a(ax + 13)*x(j…
46 …s(2) = s(2) - (a(ax + 2)*x(jdx + 0) + a(ax + 6)*x(jdx + 1) + a(ax + 10)*x(jdx + 2) + a(ax + 14)*x(…
47 …s(3) = s(3) - (a(ax + 3)*x(jdx + 0) + a(ax + 7)*x(jdx + 1) + a(ax + 11)*x(jdx + 2) + a(ax + 15)*x(…
50 x(idx + 0:idx + 3) = s(0:3)
60 s(0:3) = x(idx + 0:idx + 3)
63 …s(0) = s(0) - (a(ax + 0)*x(jdx + 0) + a(ax + 4)*x(jdx + 1) + a(ax + 8)*x(jdx + 2) + a(ax + 12)*x(j…
64 …s(1) = s(1) - (a(ax + 1)*x(jdx + 0) + a(ax + 5)*x(jdx + 1) + a(ax + 9)*x(jdx + 2) + a(ax + 13)*x(j…
[all …]
|
| /petsc/src/vec/is/section/interface/ |
| H A D | section.c | 37 PetscErrorCode PetscSectionCreate(MPI_Comm comm, PetscSection *s) in PetscSectionCreate() argument
40 PetscAssertPointer(s, 2); in PetscSectionCreate()
43 …PetscCall(PetscHeaderCreate(*s, PETSC_SECTION_CLASSID, "PetscSection", "Section", "IS", comm, Pets… in PetscSectionCreate()
44 (*s)->pStart = -1; in PetscSectionCreate()
45 (*s)->pEnd = -1; in PetscSectionCreate()
46 (*s)->perm = NULL; in PetscSectionCreate()
47 (*s)->pointMajor = PETSC_TRUE; in PetscSectionCreate()
48 (*s)->includesConstraints = PETSC_TRUE; in PetscSectionCreate()
49 (*s)->atlasDof = NULL; in PetscSectionCreate()
50 (*s)->atlasOff = NULL; in PetscSectionCreate()
[all …]
|
| /petsc/src/sys/utils/ |
| H A D | segbuffer.c | 25 struct _PetscSegBufferLink *newlink, *s; in PetscSegBufferAlloc_Private() local
28 s = seg->head; in PetscSegBufferAlloc_Private()
30 …alloc = PetscMax(s->used + count, PetscMin(1000000 / ((PetscCount)seg->unitbytes) + 1, s->alloc + … in PetscSegBufferAlloc_Private()
34 newlink->tailused = s->used + s->tailused; in PetscSegBufferAlloc_Private()
35 newlink->tail = s; in PetscSegBufferAlloc_Private()
92 struct _PetscSegBufferLink *s; in PetscSegBufferGet() local
95 s = seg->head; in PetscSegBufferGet()
96 if (PetscUnlikely(s->used + count > s->alloc)) PetscCall(PetscSegBufferAlloc_Private(seg, count)); in PetscSegBufferGet()
97 s = seg->head; in PetscSegBufferGet()
98 *(char **)buf = &s->u.array[s->used * seg->unitbytes]; in PetscSegBufferGet()
[all …]
|
| H A D | psleep.c | 28 PetscErrorCode PetscSleep(PetscReal s) in PetscSleep() argument
31 if (s < 0) getc(stdin); in PetscSleep()
35 else if (s < 1) { in PetscSleep()
38 rq.tv_nsec = (long)(s * 1e9); in PetscSleep()
43 else if (s < 1) usleep((unsigned int)(s * 1e6)); in PetscSleep()
48 sleep((int)s); in PetscSleep()
50 else _sleep((int)(s * 1000)); in PetscSleep()
52 else _sleep((int)s); in PetscSleep()
|
| /petsc/src/dm/tutorials/ |
| H A D | ex26.c | 12 int s, n = 15; in main() local
28 for (s = 0; s < n / 2; s++) { in main()
29 zeros[2 * s + tick] = -gsl_sf_hermite_zero(n, s + 1); in main()
30 zeros[2 * s + 1 + tick] = gsl_sf_hermite_zero(n, s + 1); in main()
36 for (s = 0; s < n; s++) { in main()
37 h = gsl_sf_hermite(n - 1, (double)zeros[s]); in main()
38 weights[s] = scale / (h * h); in main()
51 …for (s = 0; s < n; s++) g += weights[s] * PetscPowRealInt(sqrt(2) * sigma * zeros[s] + mu, moment); in main()
|
| /petsc/config/BuildSystem/config/ |
| H A D | framework.py | 459 …lines = [s for s in lines if s.find("icc: command line remark #10148: option '-i-dynamic' not supp…
460 …lines = [s for s in lines if s.find("[: unexpected operator") < 0] # Deals with error in mpiicc a…
461 lines = [s for s in lines if s.find(': remark #10441:') < 0]
462 lines = [s for s in lines if s.find("'linker' input unused") < 0]
464 lines = [s for s in lines if not s.startswith('cc_r:')]
466 lines = [s for s in lines if s.find('license.dat') < 0]
468 lines = [s for s in lines if s.find('INFO: catamount target') < 0]
469 lines = [s for s in lines if s.find('INFO: linux target') < 0]
471 lines = [s for s in lines if s.find('Encountered 0 errors') < 0]
473 …lines = [s for s in lines if s.find('No supported cpu target is set, CRAY_CPU_TARGET=x86-64 will b…
[all …]
|
| /petsc/src/ts/impls/multirate/ |
| H A D | mprk.c | 32 PetscInt s; /* Number of stages */ member
69 static PetscErrorCode TSMPRKGenerateTableau2(PetscInt ratio, PetscInt s, const PetscReal Abase[], c… in TSMPRKGenerateTableau2() argument
76 for (i = 0; i < s; i++) in TSMPRKGenerateTableau2()
77 for (j = 0; j < s; j++) { in TSMPRKGenerateTableau2()
78 A1[(k * s + i) * ratio * s + k * s + j] = Abase[i * s + j]; in TSMPRKGenerateTableau2()
79 A2[(k * s + i) * ratio * s + k * s + j] = Abase[i * s + j] / ratio; in TSMPRKGenerateTableau2()
83 for (i = 0; i < s; i++) in TSMPRKGenerateTableau2()
84 for (j = 0; j < s; j++) A2[(k * s + i) * ratio * s + l * s + j] = bbase[j] / ratio; in TSMPRKGenerateTableau2()
85 for (j = 0; j < s; j++) { in TSMPRKGenerateTableau2()
86 b1[k * s + j] = bbase[j] / ratio; in TSMPRKGenerateTableau2()
[all …]
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexfluent.c | 117 static PetscErrorCode DMPlexCreateFluent_ReadSection(PetscViewer viewer, FluentSection *s) in DMPlexCreateFluent_ReadSection() argument
126 snum = sscanf(buffer, "%d", &s->index); in DMPlexCreateFluent_ReadSection()
129 s->index = -1; in DMPlexCreateFluent_ReadSection()
133 if (s->index == 0) { /* Comment */ in DMPlexCreateFluent_ReadSection()
136 } else if (s->index == 2) { /* Dimension */ in DMPlexCreateFluent_ReadSection()
138 snum = sscanf(buffer, "%d", &s->nd); in DMPlexCreateFluent_ReadSection()
141 } else if (s->index == 10 || s->index == 2010) { /* Vertices */ in DMPlexCreateFluent_ReadSection()
145 snum = sscanf(buffer, "(%x %x %x %d %d)", &s->zoneID, &s->first, &s->last, &s->type, &s->nd); in DMPlexCreateFluent_ReadSection()
147 if (s->zoneID > 0) { in DMPlexCreateFluent_ReadSection()
148 PetscInt numCoords = s->last - s->first + 1; in DMPlexCreateFluent_ReadSection()
[all …]
|
| /petsc/src/vec/vec/tests/ |
| H A D | ex45.c | 20 Vec v, s; in main() local
42 PetscCall(VecCreate(PETSC_COMM_WORLD, &s)); in main()
43 PetscCall(VecSetSizes(s, PETSC_DECIDE, n / 2)); in main()
44 PetscCall(VecSetBlockSize(s, 2)); in main()
45 PetscCall(VecSetFromOptions(s)); in main()
62 PetscCall(VecStrideSubSetGather(v, PETSC_DETERMINE, vidx, NULL, s, INSERT_VALUES)); in main()
63 PetscCall(VecView(s, PETSC_VIEWER_STDOUT_WORLD)); in main()
64 PetscCall(VecScale(s, 100.0)); in main()
66 PetscCall(VecStrideSubSetScatter(s, PETSC_DETERMINE, NULL, vidx, v, ADD_VALUES)); in main()
73 PetscCall(VecStrideSubSetGather(v, 2, vidx, sidx, s, INSERT_VALUES)); in main()
[all …]
|
| /petsc/src/vec/is/section/interface/hdf5/ |
| H A D | sectionhdf5.c | 7 static PetscErrorCode PetscSectionView_HDF5_SingleField(PetscSection s, PetscViewer viewer) in PetscSectionView_HDF5_SingleField() argument
16 PetscCall(PetscObjectGetComm((PetscObject)s, &comm)); in PetscSectionView_HDF5_SingleField()
17 PetscCall(PetscSectionGetChart(s, &pStart, &pEnd)); in PetscSectionView_HDF5_SingleField()
18 hasConstraints = (s->bc) ? PETSC_TRUE : PETSC_FALSE; in PetscSectionView_HDF5_SingleField()
21 PetscCall(PetscSectionGetDof(s, p, &dof)); in PetscSectionView_HDF5_SingleField()
24 PetscCall(PetscSectionGetConstraintDof(s, p, &cdof)); in PetscSectionView_HDF5_SingleField()
38 PetscCall(PetscSectionGetDof(s, p, &dof)); in PetscSectionView_HDF5_SingleField()
41 PetscCall(PetscSectionGetOffset(s, p, &offs[n])); in PetscSectionView_HDF5_SingleField()
45 PetscCall(PetscSectionGetConstraintDof(s, p, &cdof)); in PetscSectionView_HDF5_SingleField()
46 PetscCall(PetscSectionGetConstraintIndices(s, p, &cpinds)); in PetscSectionView_HDF5_SingleField()
[all …]
|
| /petsc/src/vec/vec/utils/ |
| H A D | vsection.c | 21 PetscErrorCode PetscSectionVecView(PetscSection s, Vec v, PetscViewer viewer) in PetscSectionVecView() argument
28 PetscValidHeaderSpecific(s, PETSC_SECTION_CLASSID, 1); in PetscSectionVecView()
37 if (s->numFields) { in PetscSectionVecView()
38 …PetscCall(PetscViewerASCIIPrintf(viewer, "%s with %" PetscInt_FMT " fields\n", name, s->numFields)… in PetscSectionVecView()
39 for (f = 0; f < s->numFields; ++f) { in PetscSectionVecView()
40 …viewer, " field %" PetscInt_FMT " with %" PetscInt_FMT " components\n", f, s->numFieldComponents[… in PetscSectionVecView()
42 PetscCall(PetscSectionArrayView_ASCII_Internal(s->field[f], array, PETSC_SCALAR, viewer)); in PetscSectionVecView()
48 PetscCall(PetscSectionArrayView_ASCII_Internal(s, array, PETSC_SCALAR, viewer)); in PetscSectionVecView()
72 PetscErrorCode VecGetValuesSection(Vec v, PetscSection s, PetscInt point, PetscScalar *values[]) in VecGetValuesSection() argument
75 const PetscInt p = point - s->pStart; in VecGetValuesSection()
[all …]
|
| /petsc/src/dm/impls/da/ |
| H A D | fdda.c | 261 PetscInt xs, ys, nx, ny, i, j, ii, gxs, gys, gnx, gny, m, n, M, N, dim, s, k, nc, col; in DMCreateColoring_DA_2d_MPIAIJ() local
275 PetscCall(DMDAGetInfo(da, &dim, &m, &n, NULL, &M, &N, NULL, &nc, &s, &bx, &by, NULL, &st)); in DMCreateColoring_DA_2d_MPIAIJ()
276 col = 2 * s + 1; in DMCreateColoring_DA_2d_MPIAIJ()
282 if (st == DMDA_STENCIL_STAR && s == 1) { in DMCreateColoring_DA_2d_MPIAIJ()
325 …PetscInt xs, ys, nx, ny, i, j, gxs, gys, gnx, gny, m, n, p, dim, s, k, nc, col, zs, gzs, i… in DMCreateColoring_DA_3d_MPIAIJ() local
338 PetscCall(DMDAGetInfo(da, &dim, &m, &n, &p, &M, &N, &P, &nc, &s, &bx, &by, &bz, &st)); in DMCreateColoring_DA_3d_MPIAIJ()
339 col = 2 * s + 1; in DMCreateColoring_DA_3d_MPIAIJ()
390 PetscInt xs, nx, i, i1, gxs, gnx, l, m, M, dim, s, nc, col; in DMCreateColoring_DA_1d_MPIAIJ() local
402 PetscCall(DMDAGetInfo(da, &dim, &m, NULL, NULL, &M, NULL, NULL, &nc, &s, &bx, NULL, NULL, NULL)); in DMCreateColoring_DA_1d_MPIAIJ()
403 col = 2 * s + 1; in DMCreateColoring_DA_1d_MPIAIJ()
[all …]
|
| /petsc/src/sys/dll/ |
| H A D | dl.c | 128 char libname[PETSC_MAX_PATH_LEN], par2[PETSC_MAX_PATH_LEN], *s; in PetscDLLibraryOpen() local
153 PetscCall(PetscStrrstr(libname, ".gz", &s)); in PetscDLLibraryOpen()
154 if (s && s[3] == 0) s[0] = 0; in PetscDLLibraryOpen()
155 PetscCall(PetscStrrstr(libname, ".a", &s)); in PetscDLLibraryOpen()
156 if (s && s[2] == 0) s[0] = 0; in PetscDLLibraryOpen()
158 PetscCall(PetscStrrstr(libname, suffix, &s)); in PetscDLLibraryOpen()
159 if (s) s[0] = 0; in PetscDLLibraryOpen()
171 for (s = basename; *s; s++) in PetscDLLibraryOpen()
172 if (*s == '-') *s = '_'; in PetscDLLibraryOpen()
217 char *symbol = NULL, *s = NULL; in PetscDLLibrarySym() local
[all …]
|
| /petsc/src/dm/impls/plex/tests/output/ |
| H A D | ex9_correctness_1.out | 1 [0] Cones: 64 Total time: 1.669e-06s Average time per cone: 2.608e-08s
2 [0] Closures: 64 Total time: 1.693e-05s Average time per cone: 2.645e-07s
3 [0] VecClosures: 64 Total time: 4.506e-05s Average time per vector closure: 7.041e-07s
4 [0] VecClosures with Index: 64 Total time: 4.292e-06s Average time per vector closure: 6.706e-08s
|
| H A D | ex9_correctness_2.out | 1 [0] Cones: 4096 Total time: 1.502e-05s Average time per cone: 3.667e-09s
2 [0] Closures: 4096 Total time: 6.642e-04s Average time per cone: 1.622e-07s
3 [0] VecClosures: 4096 Total time: 1.391e-03s Average time per vector closure: 3.396e-07s
4 [0] VecClosures with Index: 4096 Total time: 2.031e-04s Average time per vector closure: 4.959e-08s
|
| H A D | ex9_correctness_0.out | 1 [0] Cones: 4 Total time: 1.192e-06s Average time per cone: 2.980e-07s
2 [0] Closures: 4 Total time: 7.868e-06s Average time per cone: 1.967e-06s
3 [0] VecClosures: 4 Total time: 3.338e-06s Average time per vector closure: 8.345e-07s
4 [0] VecClosures with Index: 4 Total time: 7.153e-07s Average time per vector closure: 1.788e-07s
|
| /petsc/src/dm/dt/dualspace/impls/simple/ |
| H A D | dspacesimple.c | 6 PetscDualSpace_Simple *s = (PetscDualSpace_Simple *)sp->data; in PetscDualSpaceSetUp_Simple() local
16 PetscCall(PetscSectionSetDof(section, pStart, s->dim)); in PetscDualSpaceSetUp_Simple()
24 PetscDualSpace_Simple *s = (PetscDualSpace_Simple *)sp->data; in PetscDualSpaceDestroy_Simple() local
27 PetscCall(PetscFree(s->numDof)); in PetscDualSpaceDestroy_Simple()
28 PetscCall(PetscFree(s)); in PetscDualSpaceDestroy_Simple()
52 PetscDualSpace_Simple *s = (PetscDualSpace_Simple *)sp->data; in PetscDualSpaceSimpleSetDimension_Simple() local
57 for (f = 0; f < s->dim; ++f) PetscCall(PetscQuadratureDestroy(&sp->functional[f])); in PetscDualSpaceSimpleSetDimension_Simple()
59 s->dim = dim; in PetscDualSpaceSimpleSetDimension_Simple()
60 PetscCall(PetscCalloc1(s->dim, &sp->functional)); in PetscDualSpaceSimpleSetDimension_Simple()
61 PetscCall(PetscFree(s->numDof)); in PetscDualSpaceSimpleSetDimension_Simple()
[all …]
|
| /petsc/src/dm/impls/plex/tutorials/ |
| H A D | ex7.c | 7 PetscSection s; in SetupSection() local
12 PetscCall(PetscSectionCreate(PetscObjectComm((PetscObject)dm), &s)); in SetupSection()
13 PetscCall(PetscSectionSetNumFields(s, 1)); in SetupSection()
14 PetscCall(PetscSectionSetFieldComponents(s, 0, 1)); in SetupSection()
15 PetscCall(PetscSectionSetChart(s, vStart, vEnd)); in SetupSection()
17 PetscCall(PetscSectionSetDof(s, v, 1)); in SetupSection()
18 PetscCall(PetscSectionSetFieldDof(s, v, 0, 1)); in SetupSection()
20 PetscCall(PetscSectionSetUp(s)); in SetupSection()
21 PetscCall(DMSetLocalSection(dm, s)); in SetupSection()
22 PetscCall(PetscSectionDestroy(&s)); in SetupSection()
|
| /petsc/doc/ |
| H A D | build_c2html.py | 37 dirs[:] = [d for d in dirs if not any([s for s in SKIPDIRS if d.startswith(s)])]
41 …d([os.path.join(loc,root,f+'.html') for f in files if any([s for s in SUFFIXES if f.endswith(s)])])
65 if any([s for s in SUFFIXES if f.endswith(s)]):
70 s = line.find('\\n')
71 line = line[l + 15:s]
77 if any([f for f in files if any([s for s in SUFFIXES if f.endswith(s)])]):
80 if any([s for s in SUFFIXES if f.endswith(s)]):
|
| /petsc/src/dm/dt/fe/impls/composite/ |
| H A D | fecomposite.c | 26 PetscInt dim, pdim, spdim, j, s; in PetscFESetUp_Composite() local
44 for (s = 0; s < cmp->numSubelements; ++s) { in PetscFESetUp_Composite()
49 PetscCall(DMPlexGetTransitiveClosure(K, s, PETSC_TRUE, &closureSize, &closure)); in PetscFESetUp_Composite()
56 for (k = 0; k < dof; k++) cmp->embedding[s * spdim + sd++] = off + k; in PetscFESetUp_Composite()
58 PetscCall(DMPlexRestoreTransitiveClosure(K, s, PETSC_TRUE, &closureSize, &closure)); in PetscFESetUp_Composite()
59 …ement %" PetscInt_FMT " has %" PetscInt_FMT " dual basis vectors != %" PetscInt_FMT, s, sd, spdim); in PetscFESetUp_Composite()
70 for (s = 0; s < cmp->numSubelements; ++s) { in PetscFESetUp_Composite()
77 PetscCall(PetscDualSpaceGetFunctional(fem->dualSpace, cmp->embedding[s * spdim + j], &f)); in PetscFESetUp_Composite()
83 invVscalar[(s * spdim + j) * spdim + k] = 0.0; in PetscFESetUp_Composite()
84 …for (q = 0; q < Nq; ++q) invVscalar[(s * spdim + j) * spdim + k] += Bf[q * spdim + k] * weights[q]; in PetscFESetUp_Composite()
[all …]
|
| /petsc/src/dm/impls/plex/tests/ |
| H A D | ex47.c | 385 PetscSection s; in main() local
437 PetscSection s; in main() local
439 PetscCall(PetscSectionCreate(PETSC_COMM_WORLD, &s)); in main()
440 PetscCall(PetscSectionSetNumFields(s, 1)); in main()
441 PetscCall(PetscSectionSetFieldComponents(s, 0, 1)); in main()
442 PetscCall(PetscSectionSetChart(s, vStart, vEnd)); in main()
453 PetscCall(PetscSectionSetDof(s, c, nbElemVertex)); in main()
454 PetscCall(PetscSectionSetFieldDof(s, c, 0, nbElemVertex)); in main()
456 PetscCall(PetscSectionSetUp(s)); in main()
458 PetscCall(PetscSectionView(s, PETSC_VIEWER_STDOUT_WORLD)); in main()
[all …]
|