Lines Matching refs:f
73 PetscInt coordDim, Nf, *Nc, f, spDim, d, v, tp; in DMProjectPoint_Func_Private() local
84 for (f = 0, v = 0, tp = 0; f < Nf; ++f) { in DMProjectPoint_Func_Private()
85 void *const ctx = ctxs ? ctxs[f] : NULL; in DMProjectPoint_Func_Private()
88 if (!sp[f]) continue; in DMProjectPoint_Func_Private()
89 PetscCall(PetscDSGetCohesive(ds, f, &cohesive)); in DMProjectPoint_Func_Private()
90 PetscCall(PetscDualSpaceGetDimension(sp[f], &spDim)); in DMProjectPoint_Func_Private()
91 if (funcs[f]) { in DMProjectPoint_Func_Private()
92 if (isFE[f]) { in DMProjectPoint_Func_Private()
100 PetscCall(PetscDualSpaceGetDM(sp[f], &rdm)); in DMProjectPoint_Func_Private()
101 PetscCall(PetscDualSpaceGetAllData(sp[f], &allPoints, NULL)); in DMProjectPoint_Func_Private()
103 PetscCall(DMGetWorkArray(rdm, numPoints * Nc[f], MPIU_SCALAR, &pointEval)); in DMProjectPoint_Func_Private()
105 PetscCall(PetscArrayzero(pointEval, numPoints * Nc[f])); in DMProjectPoint_Func_Private()
143 PetscCall((*funcs[f])(coordDim, time, v0, Nc[f], &pointEval[Nc[f] * q], ctx)); in DMProjectPoint_Func_Private()
146 PetscCall(PetscDualSpacePullback(sp[f], fegeom, numPoints, Nc[f], pointEval)); in DMProjectPoint_Func_Private()
147 PetscCall(PetscDualSpaceApplyAll(sp[f], pointEval, &values[v])); in DMProjectPoint_Func_Private()
149 PetscCall(DMRestoreWorkArray(rdm, numPoints * Nc[f], MPIU_SCALAR, &pointEval)); in DMProjectPoint_Func_Private()
155 …; d < spDim; ++d, ++v) PetscCall(PetscDualSpaceApplyFVM(sp[f], d, time, fvgeom, Nc[f], funcs[f], c… in DMProjectPoint_Func_Private()
209 PetscInt numConstants, Nf, NfIn, NfAux = 0, f, spDim, d, v, inp, tp = 0; in DMProjectPoint_Field_Private() local
288 for (f = 0, v = 0; f < Nf; ++f) { in DMProjectPoint_Field_Private()
296 if (!sp[f]) continue; in DMProjectPoint_Field_Private()
297 PetscCall(PetscDSGetCohesive(ds, f, &cohesive)); in DMProjectPoint_Field_Private()
298 PetscCall(PetscDualSpaceGetDimension(sp[f], &spDim)); in DMProjectPoint_Field_Private()
299 if (!funcs[f]) { in DMProjectPoint_Field_Private()
307 PetscCall(PetscDualSpaceGetDM(sp[f], &dm)); in DMProjectPoint_Field_Private()
308 PetscCall(PetscDualSpaceGetSymmetries(sp[f], &perms, NULL)); in DMProjectPoint_Field_Private()
309 PetscCall(PetscDualSpaceGetAllData(sp[f], &allPoints, NULL)); in DMProjectPoint_Field_Private()
311 PetscCall(DMGetWorkArray(dm, numPoints * Nc[f], MPIU_SCALAR, &pointEval)); in DMProjectPoint_Field_Private()
341 …(*funcs[f])(dE, NfIn, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, a_t, a_x, time, fegeom.v,… in DMProjectPoint_Field_Private()
343 PetscCall(PetscDualSpaceApplyAll(sp[f], pointEval, &values[v])); in DMProjectPoint_Field_Private()
344 PetscCall(DMRestoreWorkArray(dm, numPoints * Nc[f], MPIU_SCALAR, &pointEval)); in DMProjectPoint_Field_Private()
368 PetscInt numConstants, Nf, NfIn, NfAux = 0, f, spDim, d, v, inp, tp = 0; in DMProjectPoint_BdField_Private() local
459 for (f = 0, v = 0; f < Nf; ++f) { in DMProjectPoint_BdField_Private()
467 if (!sp[f]) continue; in DMProjectPoint_BdField_Private()
468 PetscCall(PetscDSGetCohesive(ds, f, &cohesive)); in DMProjectPoint_BdField_Private()
469 PetscCall(PetscDualSpaceGetDimension(sp[f], &spDim)); in DMProjectPoint_BdField_Private()
470 if (!funcs[f]) { in DMProjectPoint_BdField_Private()
477 PetscCall(PetscDualSpaceGetDM(sp[f], &dm)); in DMProjectPoint_BdField_Private()
478 PetscCall(PetscDualSpaceGetAllData(sp[f], &allPoints, NULL)); in DMProjectPoint_BdField_Private()
480 PetscCall(DMGetWorkArray(dm, numPoints * Nc[f], MPIU_SCALAR, &pointEval)); in DMProjectPoint_BdField_Private()
517 …funcs[f])(dE, NfIn, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, a_t, a_x, time, fegeom.v, f… in DMProjectPoint_BdField_Private()
519 PetscCall(PetscDualSpaceApplyAll(sp[f], pointEval, &values[v])); in DMProjectPoint_BdField_Private()
520 PetscCall(DMRestoreWorkArray(dm, numPoints * Nc[f], MPIU_SCALAR, &pointEval)); in DMProjectPoint_BdField_Private()
563 PetscInt f, numPoints; in PetscDualSpaceGetAllPointsUnion() local
571 for (f = 0; f < Nf; ++f) { in PetscDualSpaceGetAllPointsUnion()
572 if (funcs[f]) { in PetscDualSpaceGetAllPointsUnion()
576 PetscCall(PetscDualSpaceGetAllData(sp[f], &fAllPoints, NULL)); in PetscDualSpaceGetAllPointsUnion()
583 for (f = 0; f < Nf; ++f) { in PetscDualSpaceGetAllPointsUnion()
584 if (funcs[f]) { in PetscDualSpaceGetAllPointsUnion()
589 PetscCall(PetscDualSpaceGetAllData(sp[f], &fAllPoints, NULL)); in PetscDualSpaceGetAllPointsUnion()
724 … = 0, minHeight, maxHeight, minHeightIn, minHeightAux, h, regionNum, Nf, NfIn, NfAux = 0, NfTot, f; in DMProjectLocal_Generic_Plex() local
855 for (f = 0; f < Nf; ++f) { in DMProjectLocal_Generic_Plex()
858 PetscCall(PetscDSGetDiscType_Internal(ds, f, &disctype)); in DMProjectLocal_Generic_Plex()
862 isFE[f] = PETSC_TRUE; in DMProjectLocal_Generic_Plex()
864 PetscCall(PetscDSGetDiscretization(ds, f, (PetscObject *)&fe)); in DMProjectLocal_Generic_Plex()
865 PetscCall(PetscFEGetDualSpace(fe, &cellsp[f])); in DMProjectLocal_Generic_Plex()
869 isFE[f] = PETSC_FALSE; in DMProjectLocal_Generic_Plex()
871 PetscCall(PetscDSGetDiscretization(ds, f, (PetscObject *)&fv)); in DMProjectLocal_Generic_Plex()
872 PetscCall(PetscFVGetDualSpace(fv, &cellsp[f])); in DMProjectLocal_Generic_Plex()
874 isFE[f] = PETSC_FALSE; in DMProjectLocal_Generic_Plex()
875 cellsp[f] = NULL; in DMProjectLocal_Generic_Plex()
878 for (f = 0; f < NfIn; ++f) { in DMProjectLocal_Generic_Plex()
881 PetscCall(PetscDSGetDiscType_Internal(dsIn, f, &disctype)); in DMProjectLocal_Generic_Plex()
885 PetscCall(PetscDSGetDiscretization(dsIn, f, (PetscObject *)&fe)); in DMProjectLocal_Generic_Plex()
886 PetscCall(PetscFEGetDualSpace(fe, &cellspIn[f])); in DMProjectLocal_Generic_Plex()
890 PetscCall(PetscDSGetDiscretization(dsIn, f, (PetscObject *)&fv)); in DMProjectLocal_Generic_Plex()
891 PetscCall(PetscFVGetDualSpace(fv, &cellspIn[f])); in DMProjectLocal_Generic_Plex()
893 cellspIn[f] = NULL; in DMProjectLocal_Generic_Plex()
896 for (f = 0; f < Nf; ++f) { in DMProjectLocal_Generic_Plex()
898 sp[f] = cellsp[f]; in DMProjectLocal_Generic_Plex()
899 } else PetscCall(PetscDualSpaceGetHeightSubspace(cellsp[f], htInc, &sp[f])); in DMProjectLocal_Generic_Plex()
911 for (f = 0; f < NfIn; ++f) { in DMProjectLocal_Generic_Plex()
913 spIn[f] = cellspIn[f]; in DMProjectLocal_Generic_Plex()
914 } else PetscCall(PetscDualSpaceGetHeightSubspace(cellspIn[f], htIncIn, &spIn[f])); in DMProjectLocal_Generic_Plex()
916 PetscCall(PetscDSGetDiscType_Internal(dsIn, f, &disctype)); in DMProjectLocal_Generic_Plex()
918 PetscCall(PetscDSGetDiscretization(dsIn, f, (PetscObject *)&fem)); in DMProjectLocal_Generic_Plex()
922 PetscCall(PetscDSGetJetDegree(dsIn, f, &k)); in DMProjectLocal_Generic_Plex()
923 PetscCall(PetscFECreateTabulation(subfem, 1, numPoints, points, k, &T[f])); in DMProjectLocal_Generic_Plex()
925 for (f = 0; f < NfAux; ++f) { in DMProjectLocal_Generic_Plex()
926 PetscCall(PetscDSGetDiscType_Internal(dsAux, f, &disctype)); in DMProjectLocal_Generic_Plex()
928 PetscCall(PetscDSGetDiscretization(dsAux, f, (PetscObject *)&fem)); in DMProjectLocal_Generic_Plex()
932 PetscCall(PetscDSGetJetDegree(dsAux, f, &k)); in DMProjectLocal_Generic_Plex()
933 PetscCall(PetscFECreateTabulation(subfem, 1, numPoints, points, k, &TAux[f])); in DMProjectLocal_Generic_Plex()
951 for (f = 0; f < Nf; ++f) PetscCall(PetscDualSpaceGetHeightSubspace(cellsp[f], hEff, &sp[f])); in DMProjectLocal_Generic_Plex()
962 for (f = 0; f < Nf; ++f) { in DMProjectLocal_Generic_Plex()
965 if (!sp[f]) continue; in DMProjectLocal_Generic_Plex()
966 PetscCall(PetscDSGetCohesive(ds, f, &cohesive)); in DMProjectLocal_Generic_Plex()
967 PetscCall(PetscDualSpaceGetDimension(sp[f], &spDim)); in DMProjectLocal_Generic_Plex()
984 for (f = 0; f < NfIn; ++f) { in DMProjectLocal_Generic_Plex()
987 if (!spIn[f]) continue; in DMProjectLocal_Generic_Plex()
988 PetscCall(PetscDSGetCohesive(dsIn, f, &cohesive)); in DMProjectLocal_Generic_Plex()
989 PetscCall(PetscDualSpaceGetDimension(spIn[f], &spDim)); in DMProjectLocal_Generic_Plex()
1007 for (f = 0; f < NfTot; ++f) fieldActive[f] = PETSC_FALSE; in DMProjectLocal_Generic_Plex()
1008 … for (f = 0; f < Nf; ++f) fieldActive[fields[f]] = (funcs[f] && sp[f]) ? PETSC_TRUE : PETSC_FALSE; in DMProjectLocal_Generic_Plex()
1100 for (f = 0; f < NfIn; ++f) PetscCall(PetscTabulationDestroy(&T[f])); in DMProjectLocal_Generic_Plex()
1101 for (f = 0; f < NfAux; ++f) PetscCall(PetscTabulationDestroy(&TAux[f])); in DMProjectLocal_Generic_Plex()