| /petsc/src/ts/tutorials/ |
| H A D | ex22f_mf.F90 | 53 PetscInt xs, xe, gxs, gxe, dof, gdof
107 gdof = i2*(gxe - gxs + 1)
108 PetscCallA(MatCreateShell(PETSC_COMM_WORLD, dof, dof, gdof, gdof, shell_shift, A, ierr))
149 gdof = i2*(gxe - gxs + 1)
150 call SaveSolutionToDisk(da, X, gdof, xs, xe)
486 subroutine SaveSolutionToDisk(da, X, gdof, xs, xe) argument
493 PetscInt gdof, i
495 PetscScalar data2(2, xs:xe), data(gdof)
501 data2 = reshape(xx(gdof:gdof), (/two, xe - xs + 1/))
502 data = reshape(data2, (/gdof/))
[all …]
|
| /petsc/src/vec/vec/utils/ |
| H A D | vsection.c | 203 PetscInt gdof, fdof = 0; in PetscSectionGetField_Internal() local
205 PetscCall(PetscSectionGetDof(sectionGlobal, p, &gdof)); in PetscSectionGetField_Internal()
206 if (gdof > 0) PetscCall(PetscSectionGetFieldDof(section, p, field, &fdof)); in PetscSectionGetField_Internal()
211 PetscInt gdof, goff; in PetscSectionGetField_Internal() local
213 PetscCall(PetscSectionGetDof(sectionGlobal, p, &gdof)); in PetscSectionGetField_Internal()
214 if (gdof > 0) { in PetscSectionGetField_Internal()
|
| /petsc/src/vec/is/sf/utils/ |
| H A D | sfutils.c | 154 PetscInt gdof, gcdof; in PetscSFSetGraphSection() local
156 PetscCall(PetscSectionGetDof(globalSection, p, &gdof)); in PetscSFSetGraphSection()
158 …gdof < 0 ? -(gdof + 1) : gdof), PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Point %" PetscInt_FMT … in PetscSFSetGraphSection()
159 nleaves += gdof < 0 ? -(gdof + 1) - gcdof : gdof - gcdof; in PetscSFSetGraphSection()
165 PetscInt dof, cdof, off, gdof, gcdof, goff, gsize, d, c; in PetscSFSetGraphSection() local
171 PetscCall(PetscSectionGetDof(globalSection, p, &gdof)); in PetscSFSetGraphSection()
174 if (!gdof) continue; /* Censored point */ in PetscSFSetGraphSection()
175 gsize = gdof < 0 ? -(gdof + 1) - gcdof : gdof - gcdof; in PetscSFSetGraphSection()
188 if (gdof < 0) { in PetscSFSetGraphSection()
196 …cInt_FMT " mapped to invalid process %d (%" PetscInt_FMT ", %" PetscInt_FMT ")", p, r, gdof, goff); in PetscSFSetGraphSection()
|
| /petsc/src/dm/interface/ |
| H A D | dmi.c | 111 PetscInt gdof, pSubSize = 0; in PetscSectionSelectFields_Private() local
113 PetscCall(PetscSectionGetDof(gs, point, &gdof)); in PetscSectionSelectFields_Private()
114 if (gdof > 0) { in PetscSectionSelectFields_Private()
168 PetscInt gdof, goff; in PetscSectionSelectFields_Private() local
170 PetscCall(PetscSectionGetDof(gs, point, &gdof)); in PetscSectionSelectFields_Private()
171 if (gdof > 0) { in PetscSectionSelectFields_Private()
485 PetscInt gdof, gcdof, gtdof, d; in DMCreateSectionSuperDM() local
487 PetscCall(PetscSectionGetDof(sectionGlobals[i], p, &gdof)); in DMCreateSectionSuperDM()
489 gtdof = gdof - gcdof; in DMCreateSectionSuperDM()
490 if (gdof > 0 && gtdof) { in DMCreateSectionSuperDM()
|
| H A D | dm.c | 1980 PetscInt gdof; in DMCreateFieldIS() local
1982 PetscCall(PetscSectionGetDof(sectionGlobal, p, &gdof)); in DMCreateFieldIS()
1983 if (gdof > 0) { in DMCreateFieldIS()
2003 PetscInt gdof, goff; in DMCreateFieldIS() local
2005 PetscCall(PetscSectionGetDof(sectionGlobal, p, &gdof)); in DMCreateFieldIS()
2006 if (gdof > 0) { in DMCreateFieldIS()
3165 PetscInt dof, gdof, cdof, gcdof, off, goff, d, e; in DMLocalToGlobalBegin() local
3168 PetscCall(PetscSectionGetDof(gs, p, &gdof)); in DMLocalToGlobalBegin()
3174 if (!gdof || goff < 0) continue; in DMLocalToGlobalBegin()
3175 …gdof, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Inconsistent sizes, p: %" PetscInt_FMT " dof: %" PetscI… in DMLocalToGlobalBegin()
[all …]
|
| /petsc/src/dm/impls/plex/hdf5/ |
| H A D | plexhdf5.c | 589 PetscInt gdof, fdof = 0, val; in VecView_Plex_Local_HDF5_Internal() local
591 PetscCall(PetscSectionGetDof(sectionGlobal, p, &gdof)); in VecView_Plex_Local_HDF5_Internal()
592 if (gdof > 0) PetscCall(PetscSectionGetFieldDof(section, p, f, &fdof)); in VecView_Plex_Local_HDF5_Internal()
605 PetscInt gdof, goff, val; in VecView_Plex_Local_HDF5_Internal() local
607 PetscCall(PetscSectionGetDof(sectionGlobal, p, &gdof)); in VecView_Plex_Local_HDF5_Internal()
608 if (gdof > 0) { in VecView_Plex_Local_HDF5_Internal()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexfem.c | 660 PetscInt fdof, gdof, fpdof, gpdof; in DMPlexBasisTransformFieldTensor_Internal() local
668 PetscCall(PetscSectionGetFieldDof(ts, pg, g, &gdof)); in DMPlexBasisTransformFieldTensor_Internal()
681 switch (gdof) { in DMPlexBasisTransformFieldTensor_Internal()
698 switch (gdof) { in DMPlexBasisTransformFieldTensor_Internal()
743 PetscInt Nf, f, g, Np, cpf, cpg, fdof, gdof, r, c = 0; in DMPlexBasisTransformPointTensor_Internal() local
754 PetscCall(PetscSectionGetFieldDof(s, points[cpg], g, &gdof)); in DMPlexBasisTransformPointTensor_Internal()
756 c += gdof; in DMPlexBasisTransformPointTensor_Internal()
|
| H A D | plex.c | 10882 PetscInt gdof, pSubSize = 0; in DMCreateSubDomainDM_Plex() local
10884 PetscCall(PetscSectionGetDof(sectionGlobal, p, &gdof)); in DMCreateSubDomainDM_Plex()
10885 if (gdof > 0) { in DMCreateSubDomainDM_Plex()
10915 PetscInt gdof, goff; in DMCreateSubDomainDM_Plex() local
10917 PetscCall(PetscSectionGetDof(subsectionGlobal, p, &gdof)); in DMCreateSubDomainDM_Plex()
10918 if (gdof > 0) { in DMCreateSubDomainDM_Plex()
|