Lines Matching refs:nsubspaces
504 PetscErrorCode PCPatchSetDiscretisationInfo(PC pc, PetscInt nsubspaces, DM *dms, PetscInt *bs, Pets… in PCPatchSetDiscretisationInfo() argument
516 PetscCall(PetscMalloc1(nsubspaces, &sfs)); in PCPatchSetDiscretisationInfo()
517 PetscCall(PetscMalloc1(nsubspaces, &patch->dofSection)); in PCPatchSetDiscretisationInfo()
518 PetscCall(PetscMalloc1(nsubspaces, &patch->bs)); in PCPatchSetDiscretisationInfo()
519 PetscCall(PetscMalloc1(nsubspaces, &patch->nodesPerCell)); in PCPatchSetDiscretisationInfo()
520 PetscCall(PetscMalloc1(nsubspaces, &patch->cellNodeMap)); in PCPatchSetDiscretisationInfo()
521 PetscCall(PetscMalloc1(nsubspaces + 1, &patch->subspaceOffsets)); in PCPatchSetDiscretisationInfo()
523 patch->nsubspaces = nsubspaces; in PCPatchSetDiscretisationInfo()
525 for (i = 0; i < nsubspaces; ++i) { in PCPatchSetDiscretisationInfo()
536 PetscCall(PCPatchCreateDefaultSF_Private(pc, nsubspaces, sfs, patch->bs)); in PCPatchSetDiscretisationInfo()
539 patch->subspaceOffsets[nsubspaces] = subspaceOffsets[nsubspaces]; in PCPatchSetDiscretisationInfo()
555 PetscCall(DMGetNumFields(dm, &patch->nsubspaces)); in PCPatchSetDiscretisationInfoCombined()
556 PetscCall(PetscCalloc1(patch->nsubspaces, &patch->dofSection)); in PCPatchSetDiscretisationInfoCombined()
557 PetscCall(PetscMalloc1(patch->nsubspaces, &patch->bs)); in PCPatchSetDiscretisationInfoCombined()
558 PetscCall(PetscMalloc1(patch->nsubspaces, &patch->nodesPerCell)); in PCPatchSetDiscretisationInfoCombined()
559 PetscCall(PetscMalloc1(patch->nsubspaces, &patch->cellNodeMap)); in PCPatchSetDiscretisationInfoCombined()
560 PetscCall(PetscCalloc1(patch->nsubspaces + 1, &patch->subspaceOffsets)); in PCPatchSetDiscretisationInfoCombined()
563 …Call(PetscSectionGetStorageSize(patch->dofSection[0], &patch->subspaceOffsets[patch->nsubspaces])); in PCPatchSetDiscretisationInfoCombined()
565 for (i = 0; i < patch->nsubspaces; ++i) { in PCPatchSetDiscretisationInfoCombined()
812 for (g = 0; g < patch->nsubspaces; ++g) { in PCPatchGetGlobalDofs()
819 for (g = 0; g < patch->nsubspaces; ++g) { in PCPatchGetGlobalDofs()
846 for (k = 0; k < patch->nsubspaces; ++k) { in PCPatchGetPointDofs()
1196 PetscInt Nf = patch->nsubspaces; in PCPatchCreateCellPatchDiscretisationInfo()
1366 for (k = 0; k < patch->nsubspaces; ++k) { in PCPatchCreateCellPatchDiscretisationInfo()
1485 for (k = 0; k < patch->nsubspaces; ++k) { in PCPatchCreateCellPatchDiscretisationInfo()
1581 if (patch->nsubspaces > 1) { /* for nsubspaces = 1, data we need is already in dofsArray */ in PCPatchCreateCellPatchDiscretisationInfo()
1587 for (k = 0; k < patch->nsubspaces; ++k) { in PCPatchCreateCellPatchDiscretisationInfo()
1619 if (1 == patch->nsubspaces) { in PCPatchCreateCellPatchDiscretisationInfo()
1635 PetscCall(PetscSectionSetNumFields(patch->patchSection, patch->nsubspaces)); in PCPatchCreateCellPatchDiscretisationInfo()
1639 …nsubspaces, PetscObjectComm((PetscObject)pc), PETSC_ERR_ARG_WRONG, "Mismatch between number of sec… in PCPatchCreateCellPatchDiscretisationInfo()
1647 for (f = 0; f < patch->nsubspaces; ++f) { in PCPatchCreateCellPatchDiscretisationInfo()
1656 for (f = 0; f < patch->nsubspaces; ++f) { in PCPatchCreateCellPatchDiscretisationInfo()
1662 for (f = 0; f < patch->nsubspaces; ++f) { in PCPatchCreateCellPatchDiscretisationInfo()
1897 PetscInt Nf = patch->nsubspaces, Np, poff, p, f; in PCPatchComputeFunction_DMPlex_Private()
1917 …if (patch->nsubspaces == 1) PetscCall(PetscSectionSetOffset(patch->patchSection, point, oarray[(po… in PCPatchComputeFunction_DMPlex_Private()
1979 PetscInt Nf = patch->nsubspaces, Np, poff, p, f; in PCPatchComputeOperator_DMPlex_Private()
1998 …if (patch->nsubspaces == 1) PetscCall(PetscSectionSetOffset(patch->patchSection, point, oarray[(po… in PCPatchComputeOperator_DMPlex_Private()
2473 if (!patch->nsubspaces) { in PCSetUp_PATCH()
2539 localSize = patch->subspaceOffsets[patch->nsubspaces]; in PCSetUp_PATCH()
2905 for (i = 0; i < patch->nsubspaces; i++) PetscCall(PetscSectionDestroy(&patch->dofSection[i])); in PCReset_PATCH()
2910 for (i = 0; i < patch->nsubspaces; i++) PetscCall(PetscFree(patch->cellNodeMap[i])); in PCReset_PATCH()
2954 patch->nsubspaces = 0; in PCReset_PATCH()
3047 if (patch->nsubspaces <= 0) { in PCSetFromOptions_PATCH()
3050 nfields = patch->nsubspaces; in PCSetFromOptions_PATCH()