| /petsc/src/dm/dt/fe/impls/composite/ |
| H A D | fecomposite.c | 26 PetscInt dim, pdim, spdim, j, s; in PetscFESetUp_Composite() local
40 PetscCall(PetscSpaceGetDimension(fem->basisSpace, &spdim)); in PetscFESetUp_Composite()
41 PetscCall(PetscMalloc1(cmp->numSubelements * spdim, &cmp->embedding)); in PetscFESetUp_Composite()
56 for (k = 0; k < dof; k++) cmp->embedding[s * spdim + sd++] = off + k; in PetscFESetUp_Composite()
59 …spdim, PetscObjectComm((PetscObject)fem), PETSC_ERR_PLIB, "Subelement %" PetscInt_FMT " has %" Pet… in PetscFESetUp_Composite()
63 PetscCall(PetscMalloc1(cmp->numSubelements * spdim * spdim, &fem->invV)); in PetscFESetUp_Composite()
64 PetscCall(PetscMalloc2(spdim, &pivots, spdim, &work)); in PetscFESetUp_Composite()
66 PetscCall(PetscMalloc1(cmp->numSubelements * spdim * spdim, &invVscalar)); in PetscFESetUp_Composite()
71 for (j = 0; j < spdim; ++j) { in PetscFESetUp_Composite()
77 PetscCall(PetscDualSpaceGetFunctional(fem->dualSpace, cmp->embedding[s * spdim + j], &f)); in PetscFESetUp_Composite()
[all …]
|
| /petsc/src/dm/dt/space/impls/tensor/ |
| H A D | spacetensor.c | 330 PetscInt sNv, sNc, spdim; in PetscSpaceEvaluate_Tensor() local
335 PetscCall(PetscSpaceGetDimension(tens->tensspaces[s], &spdim)); in PetscSpaceEvaluate_Tensor()
336 …spdim), PETSC_COMM_SELF, PETSC_ERR_PLIB, "Bad tensor loop: Nv %" PetscInt_FMT ", Ns %" PetscInt_FM… in PetscSpaceEvaluate_Tensor()
337 …", s %" PetscInt_FMT ", cstep %" PetscInt_FMT ", sNc %" PetscInt_FMT, Nv, Ns, Nc, s, cstep, spdim); in PetscSpaceEvaluate_Tensor()
338 vskip = pdim / (vstep * spdim); in PetscSpaceEvaluate_Tensor()
346 for (PetscInt si = 0; si < spdim; si++) { in PetscSpaceEvaluate_Tensor()
348 PetscInt i = (k * spdim + si) * vstep + j; in PetscSpaceEvaluate_Tensor()
355 …for (PetscInt p = 0; p < npoints; p++) B[(pdim * p + i) * Nc + c] *= sB[(spdim * p + si) * sNc + s… in PetscSpaceEvaluate_Tensor()
365 for (PetscInt si = 0; si < spdim; si++) { in PetscSpaceEvaluate_Tensor()
367 PetscInt i = (k * spdim + si) * vstep + j; in PetscSpaceEvaluate_Tensor()
[all …]
|
| /petsc/src/dm/dt/dualspace/impls/lagrange/tests/ |
| H A D | ex1.c | 168 PetscInt spdim, spintdim, exspdim, exspintdim; in testLagrange() local
189 PetscCall(PetscDualSpaceGetDimension(sp, &spdim)); in testLagrange()
191 …scCheck(spdim == exspdim, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Expected dual space dimension %" Petsc… in testLagrange()
226 PetscCheck(nNodes == spdim, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Incorrect nNodes"); in testLagrange()
235 for (i = 0; i < spdim; i++) { in testLagrange()
|
| /petsc/src/dm/dt/dualspace/impls/refined/ |
| H A D | dualspacerefined.c | 66 PetscInt cStart, cEnd, c, spdim; in PetscDualSpaceSetUp_Refined() local
123 PetscCall(PetscDualSpaceGetDimension(sp, &spdim)); in PetscDualSpaceSetUp_Refined()
124 PetscCall(PetscMalloc1(spdim, &sp->functional)); in PetscDualSpaceSetUp_Refined()
|
| /petsc/src/dm/dt/space/impls/sum/ |
| H A D | spacesum.c | 447 PetscInt sNv, spdim, sNc, p; in PetscSpaceEvaluate_Sum() local
451 PetscCall(PetscSpaceGetDimension(sum->sumspaces[s], &spdim)); in PetscSpaceEvaluate_Sum()
452 …PetscCheck(offset + spdim <= pdimfull, PetscObjectComm((PetscObject)sp), PETSC_ERR_ARG_OUTOFRANGE,… in PetscSpaceEvaluate_Sum()
458 for (j = 0; j < spdim; ++j) { in PetscSpaceEvaluate_Sum()
467 sBInd = (p * spdim + j) * sNc + c; in PetscSpaceEvaluate_Sum()
495 offset += spdim; in PetscSpaceEvaluate_Sum()
|
| /petsc/src/dm/dt/dualspace/interface/ |
| H A D | dualspace.c | 456 s->spdim = -1; in PetscDualSpaceCreate()
694 if (sp->spdim < 0) { in PetscDualSpaceGetDimension()
698 if (section) PetscCall(PetscSectionGetStorageSize(section, &sp->spdim)); in PetscDualSpaceGetDimension()
699 else sp->spdim = 0; in PetscDualSpaceGetDimension()
701 *dim = sp->spdim; in PetscDualSpaceGetDimension()
1295 PetscInt spdim; in PetscDualSpaceCreateAllDataDefault() local
1306 PetscCall(PetscDualSpaceGetDimension(sp, &spdim)); in PetscDualSpaceCreateAllDataDefault()
1307 if (!spdim) { in PetscDualSpaceCreateAllDataDefault()
1311 nrows = spdim; in PetscDualSpaceCreateAllDataDefault()
1315 for (f = 1; f < spdim; f++) { in PetscDualSpaceCreateAllDataDefault()
[all …]
|
| /petsc/src/dm/impls/plex/ |
| H A D | plextree.c | 1143 PetscInt spdim; in DMPlexComputeAnchorMatrix_Tree_Direct() local
1158 PetscCall(DMGetDimension(dm, &spdim)); in DMPlexComputeAnchorMatrix_Tree_Direct()
1159 …PetscCall(PetscMalloc6(spdim, &v0, spdim, &v0parent, spdim, &vtmp, spdim * spdim, &J, spdim * spdi… in DMPlexComputeAnchorMatrix_Tree_Direct()
1194 PetscCall(PetscSpaceSetNumVariables(bspace, spdim)); in DMPlexComputeAnchorMatrix_Tree_Direct()
1200 for (i = 0, maxDof = 0; i <= spdim; i++) maxDof = PetscMax(maxDof, numDof[i]); in DMPlexComputeAnchorMatrix_Tree_Direct()
1219 …scCall(PetscMalloc7(fSize, &sizes, nPoints * Nc, &weights, spdim * nPoints, &pointsRef, spdim * nP… in DMPlexComputeAnchorMatrix_Tree_Direct()
1230 PetscCall(PetscArraycpy(pointsRef + spdim * offset, p, spdim * qPoints)); in DMPlexComputeAnchorMatrix_Tree_Direct()
1266 CoordinatesRefToReal(spdim, spdim, xi0, v0, J, &pointsRef[i * spdim], vtmp); in DMPlexComputeAnchorMatrix_Tree_Direct()
1267 CoordinatesRealToRef(spdim, spdim, xi0, v0parent, invJparent, vtmp, &pointsReal[i * spdim]); in DMPlexComputeAnchorMatrix_Tree_Direct()
|
| /petsc/src/dm/dt/dualspace/impls/lagrange/ |
| H A D | dspacelagrange.c | 1697 PetscInt nNodes, spdim; in PetscDualSpaceComputeFunctionalsFromAllData() local
1710 PetscCall(PetscSectionGetStorageSize(section, &spdim)); in PetscDualSpaceComputeFunctionalsFromAllData()
1711 PetscCheck(spdim == nDofs, PETSC_COMM_SELF, PETSC_ERR_PLIB, "incompatible all matrix size"); in PetscDualSpaceComputeFunctionalsFromAllData()
2162 PetscInt spdim, f; in PetscDualSpaceSetUp_Lagrange() local
2170 PetscCall(PetscDualSpaceGetDimension(spcont, &spdim)); in PetscDualSpaceSetUp_Lagrange()
2171 sp->spdim = sp->spintdim = spdim; in PetscDualSpaceSetUp_Lagrange()
2172 PetscCall(PetscSectionSetDof(section, 0, spdim)); in PetscDualSpaceSetUp_Lagrange()
2174 PetscCall(PetscMalloc1(spdim, &sp->functional)); in PetscDualSpaceSetUp_Lagrange()
2175 for (f = 0; f < spdim; f++) { in PetscDualSpaceSetUp_Lagrange()
2448 PetscCall(PetscSectionGetStorageSize(section, &sp->spdim)); in PetscDualSpaceSetUp_Lagrange()
|
| /petsc/include/petsc/private/ |
| H A D | petscfeimpl.h | 120 PetscInt spdim; /* The dual-space dimension */ member
|