Home
last modified time | relevance | path

Searched refs:Nk (Results 1 – 14 of 14) sorted by relevance

/petsc/src/dm/dt/interface/
H A Ddtaltv.c103 PetscInt Nk, Nf; in PetscDTAltVApply() local
109 PetscCall(PetscDTBinomialInt(N, k, &Nk)); in PetscDTAltVApply()
111 for (i = 0; i < Nk; i++) { in PetscDTAltVApply()
268 PetscInt Nk; in PetscDTAltVWedgeMatrix() local
272 PetscCall(PetscDTBinomialInt(N, k, &Nk)); in PetscDTAltVWedgeMatrix()
276 for (i = 0; i < Njk * Nk; i++) awedgeMat[i] = 0.; in PetscDTAltVWedgeMatrix()
290 awedgeMat[i * Nk + kInd] += jkOdd ? -a[jInd] : a[jInd]; in PetscDTAltVWedgeMatrix()
325 PetscInt i, j, Nk, Mk; in PetscDTAltVPullback() local
332 PetscCall(PetscDTBinomialInt(N, PetscAbsInt(k), &Nk)); in PetscDTAltVPullback()
336 for (i = 0; i < Nk; i++) { in PetscDTAltVPullback()
[all …]
H A Ddt.c1100 PetscInt Nk, Ndeg; in PetscDTPKDEvalJet() local
1106 PetscCall(PetscDTBinomialInt(dim + k, k, &Nk)); in PetscDTPKDEvalJet()
1133 for (i = 0; i < npoints; i++) p[degidx * Nk * npoints + i] = 1.; in PetscDTPKDEvalJet()
1134 for (i = 0; i < (Nk - 1) * npoints; i++) p[(degidx * Nk + 1) * npoints + i] = 0.; in PetscDTPKDEvalJet()
1175 for (kidx = 0; kidx < Nk; kidx++) { in PetscDTPKDEvalJet()
1180 …p[(degidx * Nk + kidx) * npoints + pt] = (cnm1 * thetanm1 + cnm1x * thetanm1x) * p[(m1idx * Nk + k… in PetscDTPKDEvalJet()
1181 …if (m2idx >= 0) p[(degidx * Nk + kidx) * npoints + pt] -= cnm2 * thetanm2 * p[(m2idx * Nk + kidx) … in PetscDTPKDEvalJet()
1196 …p[(degidx * Nk + kidx) * npoints + pt] += mplty * 0.5 * (cnm1x - cnm1) * p[(m1idx * Nk + km1idx) *… in PetscDTPKDEvalJet()
1198 …p[(degidx * Nk + kidx) * npoints + pt] += mplty * cnm2 * thetanm1 * p[(m2idx * Nk + km1idx) * npoi… in PetscDTPKDEvalJet()
1210 …p[(degidx * Nk + kidx) * npoints + pt] -= 0.5 * factor * cnm2 * p[(m2idx * Nk + km2idx) * npoints … in PetscDTPKDEvalJet()
[all …]
H A Ddtweakform.c956 PetscInt Nf = wf->Nf, Nk, k; in PetscWeakFormViewTable_Ascii() local
959 PetscCall(PetscHMapFormGetSize(map, &Nk)); in PetscWeakFormViewTable_Ascii()
960 if (Nk) { in PetscWeakFormViewTable_Ascii()
968 PetscCall(PetscMalloc6(Nk, &keys, Nk, &names, Nk, &values, Nk, &idx1, Nk, &idx2, Nk, &idx)); in PetscWeakFormViewTable_Ascii()
973 for (k = 0; k < Nk; ++k) { in PetscWeakFormViewTable_Ascii()
977 PetscCall(PetscSortIntWithPermutation(Nk, values, idx1)); in PetscWeakFormViewTable_Ascii()
979 for (k = 0; k < Nk; ++k) { in PetscWeakFormViewTable_Ascii()
984 PetscCall(PetscSortStrWithPermutation(Nk, names, idx2)); in PetscWeakFormViewTable_Ascii()
985 for (k = 0; k < Nk; ++k) { in PetscWeakFormViewTable_Ascii()
993 for (k = 0; k < Nk; ++k) { in PetscWeakFormViewTable_Ascii()
[all …]
/petsc/src/dm/dt/tests/
H A Dex7.c8 PetscInt Nk, Mk, i, j, l; in CheckPullback() local
17 PetscCall(PetscDTBinomialInt(N, k, &Nk)); in CheckPullback()
26 PetscCall(PetscMalloc2(Nk, &Lstarw, M * k, &Lx)); in CheckPullback()
27 PetscCall(PetscMalloc2(Nk * Mk, &Lstar, Nk, &Lstarwcheck)); in CheckPullback()
33 PetscCall(PetscMalloc1(Nk, &sLsw)); in CheckPullback()
50 for (i = 0; i < Nk; i++) { in CheckPullback()
67 if (Nk * Mk > 0) PetscCall(PetscRealView(Nk * Mk, Lstar, viewer)); in CheckPullback()
72 if (Nk > 0) PetscCall(PetscRealView(Nk, Lstarw, viewer)); in CheckPullback()
133 PetscInt j, Nk, M; in main() local
137 PetscCall(PetscDTBinomialInt(N, k, &Nk)); in main()
[all …]
H A Dex8.c18 PetscInt j, k, Nk, kchk; in main() local
20 PetscCall(PetscDTBinomialInt(d + n, d, &Nk)); in main()
21 for (k = 0; k < Nk; k++) { in main()
38 PetscCall(PetscDTIndexToGradedOrder(d, Nk - 1 - k, gtup)); in main()
40 …Nk - 1 - k, PETSC_COMM_SELF, PETSC_ERR_PLIB, "PetscDTGradedOrderToIndex, d = %" PetscInt_FMT ", n … in main()
41 … d = %" PetscInt_FMT ", n = %" PetscInt_FMT ", k = %" PetscInt_FMT " incorrect", d, n, Nk - 1 - k); in main()
H A Dex13.c11 PetscInt Nk; // jet size in constructTabulationAndMass() local
17 PetscCall(PetscDTBinomialInt(dim + jetDegree, dim, &Nk)); in constructTabulationAndMass()
18 PetscCall(PetscMalloc1(Nbpt * Nf * Nk * npoints, &p_trimmed)); in constructTabulationAndMass()
29 const PetscReal *p_i = &p_trimmed[(i * Nf + f) * Nk * npoints]; in constructTabulationAndMass()
30 const PetscReal *p_j = &p_trimmed[(j * Nf + f) * Nk * npoints]; in constructTabulationAndMass()
39 *_Nk = Nk; in constructTabulationAndMass()
52 PetscInt Nk; // jet size in test() local
74 …nAndMass(dim, deg, form, jetDegree, npoints, points, weights, &Nbpt, &Nf, &Nk, &p_trimmed, &M_trim… in test()
77 PetscCall(PetscMalloc1(Nbp * Nk * npoints, &p_scalar)); in test()
108 const PetscReal *p_i = &p_scalar[i * Nk * npoints]; in test()
[all …]
H A Dex9.c35 PetscInt Np, Nk, i, j, l, d, npoints; in testDerivativesLegendre() local
54 PetscCall(PetscDTBinomialInt(dim + k, dim, &Nk)); in testDerivativesLegendre()
112 PetscCall(PetscMalloc3(Nk * Np, &pkd_jet_basis, Nk, &lgndre_jet, Nk, &pkd_jet)); in testDerivativesLegendre()
115 for (i = 0; i < Nk; i++) { in testDerivativesLegendre()
117 for (j = 0; j < Np; j++) val += pkd_coeffs[j] * pkd_jet_basis[j * Nk + i]; in testDerivativesLegendre()
128 for (j = 0; j < Nk; j++) lgndre_jet[j] = 0.; in testDerivativesLegendre()
133 for (j = 0; j < Nk; j++) { in testDerivativesLegendre()
144 for (i = 0; i < Nk; i++) { in testDerivativesLegendre()
/petsc/src/dm/dt/dualspace/impls/lagrange/
H A Ddspacelagrange.c558 PetscInt depth, pdepth, Nk, pNk; in PetscLagNodeIndicesPushForward() local
632 PetscCall(PetscDTBinomialInt(dim, PetscAbsInt(formDegree), &Nk)); in PetscLagNodeIndicesPushForward()
634 PetscCall(DMGetWorkArray(dm, pNk * Nk, MPIU_REAL, &Jstar)); in PetscLagNodeIndicesPushForward()
637 for (i = 0; i < Nk; i++) { in PetscLagNodeIndicesPushForward()
639 for (j = 0; j < pNk; j++) val += nodeVec[n * pNk + j] * Jstar[j * Nk + i]; in PetscLagNodeIndicesPushForward()
640 pfNodeVec[n * Nk + i] = val; in PetscLagNodeIndicesPushForward()
643 PetscCall(DMRestoreWorkArray(dm, pNk * Nk, MPIU_REAL, &Jstar)); in PetscLagNodeIndicesPushForward()
655 PetscInt Nk, NkT, NkF; in PetscLagNodeIndicesTensor() local
666 PetscCall(PetscDTBinomialInt(dim, PetscAbsInt(formDegree), &Nk)); in PetscLagNodeIndicesTensor()
673 ni->nodeVecDim = Nk; in PetscLagNodeIndicesTensor()
[all …]
/petsc/include/
H A Dpetscdt.h458 PetscInt Nk; in PetscDTEnumSubset() local
461 PetscCall(PetscDTBinomialInt(n, k, &Nk)); in PetscDTEnumSubset()
463 PetscInt Nminuskminus = (Nk * (k - l)) / (n - i); in PetscDTEnumSubset()
464 PetscInt Nminusk = Nk - Nminuskminus; in PetscDTEnumSubset()
468 Nk = Nminuskminus; in PetscDTEnumSubset()
471 Nk = Nminusk; in PetscDTEnumSubset()
498 PetscInt j = 0, Nk; in PetscDTSubsetIndex() local
502 PetscCall(PetscDTBinomialInt(n, k, &Nk)); in PetscDTSubsetIndex()
504 PetscInt Nminuskminus = (Nk * (k - l)) / (n - i); in PetscDTSubsetIndex()
505 PetscInt Nminusk = Nk - Nminuskminus; in PetscDTSubsetIndex()
[all …]
/petsc/src/snes/utils/
H A Ddmplexsnes.c359 PetscInt Nm = 2, m, Nk = 0, k, kp, off = 0; in DMPlexSNESComputeResidualDS() local
366 Nk += Nkm; in DMPlexSNESComputeResidualDS()
368 PetscCall(PetscMalloc1(Nk, &reskeys)); in DMPlexSNESComputeResidualDS()
370 …PetscCheck(off == Nk, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Number of keys %" PetscInt_FMT " should… in DMPlexSNESComputeResidualDS()
371 PetscCall(PetscFormKeySort(Nk, reskeys)); in DMPlexSNESComputeResidualDS()
372 for (k = 0, kp = 1; kp < Nk; ++kp) { in DMPlexSNESComputeResidualDS()
378 Nk = k; in DMPlexSNESComputeResidualDS()
381 for (k = 0; k < Nk; ++k) { in DMPlexSNESComputeResidualDS()
474 PetscInt Nm = 4, m, Nk = 0, k, kp, off = 0; in DMSNESComputeJacobianAction() local
481 Nk += Nkm; in DMSNESComputeJacobianAction()
[all …]
/petsc/src/dm/dt/dualspace/impls/lagrange/tests/
H A Dex1.c164 PetscInt Nk; in testLagrange() local
171 PetscCall(PetscDTBinomialInt(dim, PetscAbsInt(formDegree), &Nk)); in testLagrange()
177 PetscCall(PetscDualSpaceSetNumComponents(sp, nCopies * Nk)); in testLagrange()
225 PetscCheck(nodeVecDim == Nk, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Incorrect nodeVecDim"); in testLagrange()
/petsc/src/dm/impls/plex/tutorials/
H A Dex6.c65 static PetscErrorCode LoadData3D(DM dm, PetscInt Ni, PetscInt Nj, PetscInt Nk, PetscInt clSize, Vec… in LoadData3D() argument
72 for (k = 0; k < Nk; ++k) { in LoadData3D()
/petsc/src/dm/dt/dualspace/impls/sum/
H A Ddualspacesum.c488 PetscInt k, Nc, Nk, Nknew, Ncnew, Ns; in PetscDualSpaceCreateFacetSubspace_Sum() local
520 PetscCall(PetscDTBinomialInt(dim, PetscAbsInt(k), &Nk)); in PetscDualSpaceCreateFacetSubspace_Sum()
521 Ncopies = Nc / Nk; in PetscDualSpaceCreateFacetSubspace_Sum()
538 subNcopies = subNc / Nk; in PetscDualSpaceCreateFacetSubspace_Sum()
/petsc/src/dm/dt/dualspace/interface/
H A Ddualspace.c1887 PetscInt dim, v, c, Nk; in PetscDualSpaceTransform() local
1897 PetscCall(PetscDTBinomialInt(dim, PetscAbsInt(dsp->k), &Nk)); in PetscDualSpaceTransform()
1901 switch (Nk) { in PetscDualSpaceTransform()
1912 …)dsp), PETSC_ERR_ARG_OUTOFRANGE, "Unsupported form size %" PetscInt_FMT " for transformation", Nk); in PetscDualSpaceTransform()