Lines Matching refs:dim
10 PetscInt dim; member
19 …PetscHashCombine(PetscHashCombine(PetscHashCombine(PetscHashInt((key).dim), PetscHashInt((key).ord…
22 …(((k1).dim == (k2).dim) ? ((k1).order == (k2).order) ? ((k1).formDegree == (k2).formDegree) ? ((k1…
26 static PetscErrorCode ExpectedNumDofs_Total(PetscInt dim, PetscInt order, PetscInt formDegree, Pets…
27 static PetscErrorCode ExpectedNumDofs_Interior(PetscInt dim, PetscInt order, PetscInt formDegree, P…
29 static PetscErrorCode ExpectedNumDofs_Total(PetscInt dim, PetscInt order, PetscInt formDegree, Pets… in ExpectedNumDofs_Total() argument
40 PetscCall(PetscDTBinomialInt(order + dim, order + formDegree, &rnchooserk)); in ExpectedNumDofs_Total()
47 PetscCall(PetscDTBinomialInt(order + dim, order + formDegree, &rnchooserk)); in ExpectedNumDofs_Total()
56 PetscCall(PetscDTBinomialInt(dim, formDegree, &nchoosek)); in ExpectedNumDofs_Total()
58 rp1pownmk = PetscPowInt(order + 1, dim - formDegree); in ExpectedNumDofs_Total()
64 PetscCall(PetscDTBinomialInt(dim, formDegree, &nchoosek)); in ExpectedNumDofs_Total()
65 rp1pown = PetscPowInt(order + 1, dim); in ExpectedNumDofs_Total()
74 if (formDegree < dim) { in ExpectedNumDofs_Total()
75 PetscCall(ExpectedNumDofs_Total(dim - 1, order, formDegree, trimmed, 0, 1, &tracek)); in ExpectedNumDofs_Total()
79 PetscCall(ExpectedNumDofs_Total(dim - 1, order, formDegree - 1, trimmed, 0, 1, &tracekm1)); in ExpectedNumDofs_Total()
87 static PetscErrorCode ExpectedNumDofs_Interior(PetscInt dim, PetscInt order, PetscInt formDegree, P… in ExpectedNumDofs_Interior() argument
95 if (order + formDegree > dim) { in ExpectedNumDofs_Interior()
96 PetscInt eorder = order + formDegree - dim - 1; in ExpectedNumDofs_Interior()
97 PetscInt eformDegree = dim - formDegree; in ExpectedNumDofs_Interior()
101 PetscCall(PetscDTBinomialInt(eorder + dim, eorder + eformDegree, &rnchooserk)); in ExpectedNumDofs_Interior()
109 if (order + formDegree > dim) { in ExpectedNumDofs_Interior()
110 PetscInt eorder = order + formDegree - dim; in ExpectedNumDofs_Interior()
111 PetscInt eformDegree = dim - formDegree; in ExpectedNumDofs_Interior()
115 PetscCall(PetscDTBinomialInt(eorder + dim, eorder + eformDegree, &rnchooserk)); in ExpectedNumDofs_Interior()
123 if (dim < 2) { in ExpectedNumDofs_Interior()
124 PetscCall(ExpectedNumDofs_Interior(dim, order, formDegree, trimmed, 0, nCopies, nDofs)); in ExpectedNumDofs_Interior()
131 if (formDegree < dim) { in ExpectedNumDofs_Interior()
132 …PetscCall(ExpectedNumDofs_Interior(dim - 1, order, formDegree, trimmed, dim > 2 ? 1 : 0, 1, &trace… in ExpectedNumDofs_Interior()
136 …PetscCall(ExpectedNumDofs_Interior(dim - 1, order, formDegree - 1, trimmed, dim > 2 ? 1 : 0, 1, &t… in ExpectedNumDofs_Interior()
147 if (formDegree < dim) { in ExpectedNumDofs_Interior()
148 PetscCall(ExpectedNumDofs_Interior(dim - 1, order, formDegree, trimmed, 0, 1, &tracek)); in ExpectedNumDofs_Interior()
152 PetscCall(ExpectedNumDofs_Interior(dim - 1, order, formDegree - 1, trimmed, 0, 1, &tracekm1)); in ExpectedNumDofs_Interior()
160 PetscErrorCode testLagrange(PetscHashLag lagTable, DM K, PetscInt dim, PetscInt order, PetscInt for… in testLagrange() argument
171 PetscCall(PetscDTBinomialInt(dim, PetscAbsInt(formDegree), &Nk)); in testLagrange()
181 …etscInt_FMT ", formDegree %" PetscInt_FMT ", nCopies %" PetscInt_FMT "\n", dim, order, (PetscInt)t… in testLagrange()
182 PetscCall(ExpectedNumDofs_Total(dim, order, formDegree, trimmed, tensorCell, nCopies, &exspdim)); in testLagrange()
183 if (continuous && dim > 0 && order > 0) { in testLagrange()
184 …PetscCall(ExpectedNumDofs_Interior(dim, order, formDegree, trimmed, tensorCell, nCopies, &exspintd… in testLagrange()
193 key.dim = dim; in testLagrange()
273 if (dim <= 2 && spintdim) { in testLagrange()
301 PetscInt dim; in main() local
311 dim = 3; in main()
316 PetscCall(PetscOptionsRangeInt("-dim", "The spatial dimension", "ex1.c", dim, &dim, NULL, 0, 3)); in main()
324 …PetscCall(DMPlexCreateReferenceCell(PETSC_COMM_SELF, DMPolytopeTypeSimpleShape(dim, (PetscBool)(te… in main()
329 ordermax = tensorCell == 2 ? 4 : tensorCell == 1 ? 3 : dim + 2; in main()
333 for (formDegree = PetscMin(0, -dim + 1); formDegree <= dim; formDegree++) { in main()
336 …for (nCopies = 1; nCopies <= 3; nCopies++) PetscCall(testLagrange(lagTable, dm, dim, order, formDe… in main()