Lines Matching refs:npoints
307 …PetscQuadrature q, PeOp PetscInt *dim, PeOp PetscInt *Nc, PeOp PetscInt *npoints, PeOp const Petsc… in PetscQuadratureGetData() argument
319 if (npoints) { in PetscQuadratureGetData()
320 PetscAssertPointer(npoints, 4); in PetscQuadratureGetData()
321 *npoints = q->numPoints; in PetscQuadratureGetData()
566 …QuadratureSetData(PetscQuadrature q, PetscInt dim, PetscInt Nc, PetscInt npoints, const PetscReal … in PetscQuadratureSetData() argument
572 if (npoints >= 0) q->numPoints = npoints; in PetscQuadratureSetData()
669 PetscInt dim, Nc, order, npoints, npointsRef, c, p, cp, d, e; in PetscQuadratureExpandComposite() local
679 PetscCall(PetscQuadratureGetData(q, &dim, &Nc, &npoints, &points, &weights)); in PetscQuadratureExpandComposite()
680 npointsRef = npoints * numSubelements; in PetscQuadratureExpandComposite()
684 for (p = 0; p < npoints; ++p) { in PetscQuadratureExpandComposite()
686 pointsRef[(c * npoints + p) * dim + d] = v0[c * dim + d]; in PetscQuadratureExpandComposite()
687 …for (e = 0; e < dim; ++e) pointsRef[(c * npoints + p) * dim + d] += jac[(c * dim + d) * dim + e] *… in PetscQuadratureExpandComposite()
690 …for (cp = 0; cp < Nc; ++cp) weightsRef[(c * npoints + p) * Nc + cp] = weights[p * Nc + cp] / numSu… in PetscQuadratureExpandComposite()
785 static PetscErrorCode PetscDTJacobiEval_Internal(PetscInt npoints, PetscReal a, PetscReal b, PetscI… in PetscDTJacobiEval_Internal() argument
797 for (i = 0; i < npoints; i++) in PetscDTJacobiEval_Internal()
801 for (i = 0; i < npoints; i++) { in PetscDTJacobiEval_Internal()
866 PetscErrorCode PetscDTJacobiEvalJet(PetscReal alpha, PetscReal beta, PetscInt npoints, const PetscR… in PetscDTJacobiEvalJet() argument
876 …<= k; i++) PetscCall(PetscDTJacobiEval_Internal(npoints, alpha, beta, i, points, 1, &zero, &p[i * … in PetscDTJacobiEvalJet()
880 PetscCall(PetscMalloc1((degree + 1) * npoints, &psingle)); in PetscDTJacobiEvalJet()
883 …PetscCall(PetscDTJacobiEval_Internal(npoints, alpha, beta, i, points, degree + 1, degrees, psingle… in PetscDTJacobiEvalJet()
885 … for (l = 0; l < npoints; l++) p[(j * (k + 1) + i) * npoints + l] = psingle[l * (degree + 1) + j]; in PetscDTJacobiEvalJet()
916 PetscErrorCode PetscDTJacobiEval(PetscInt npoints, PetscReal alpha, PetscReal beta, const PetscReal… in PetscDTJacobiEval() argument
921 if (!npoints || !ndegree) PetscFunctionReturn(PETSC_SUCCESS); in PetscDTJacobiEval()
922 …if (B) PetscCall(PetscDTJacobiEval_Internal(npoints, alpha, beta, 0, points, ndegree, degrees, B)); in PetscDTJacobiEval()
923 …if (D) PetscCall(PetscDTJacobiEval_Internal(npoints, alpha, beta, 1, points, ndegree, degrees, D)); in PetscDTJacobiEval()
924 …if (D2) PetscCall(PetscDTJacobiEval_Internal(npoints, alpha, beta, 2, points, ndegree, degrees, D2… in PetscDTJacobiEval()
948 PetscErrorCode PetscDTLegendreEval(PetscInt npoints, const PetscReal *points, PetscInt ndegree, con… in PetscDTLegendreEval() argument
951 PetscCall(PetscDTJacobiEval(npoints, 0., 0., points, ndegree, degrees, B, D, D2)); in PetscDTLegendreEval()
1097 PetscErrorCode PetscDTPKDEvalJet(PetscInt dim, PetscInt npoints, const PetscReal points[], PetscInt… in PetscDTPKDEvalJet() argument
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()
1163 for (pt = 0; pt < npoints; pt++) { 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()
1215 … p[(degidx * Nk + kidx) * npoints + pt] += mplty * cnm1x * p[(m1idx * Nk + km1idx) * npoints + pt]; in PetscDTPKDEvalJet()
1226 for (i = 0; i < Nk * npoints; i++) p[degidx * Nk * npoints + i] *= scale; in PetscDTPKDEvalJet()
1264 static PetscErrorCode PetscDTPTrimmedEvalJet_Internal(PetscInt dim, PetscInt npoints, const PetscRe… in PetscDTPTrimmedEvalJet_Internal() argument
1272 PetscCall(PetscDTPKDEvalJet(dim, npoints, points, degree, jetDegree, p)); in PetscDTPTrimmedEvalJet_Internal()
1276 PetscCall(PetscDTPKDEvalJet(dim, npoints, points, degree - 1, jetDegree, p)); in PetscDTPTrimmedEvalJet_Internal()
1285 PetscCall(PetscArrayzero(p, Nbpt * Nf * Nk * npoints)); in PetscDTPTrimmedEvalJet_Internal()
1290 PetscCall(PetscMalloc1(Nbpm1 * Nk * npoints, &p_scalar)); in PetscDTPTrimmedEvalJet_Internal()
1291 PetscCall(PetscDTPKDEvalJet(dim, npoints, points, degree - 1, jetDegree, p_scalar)); in PetscDTPTrimmedEvalJet_Internal()
1296 const PetscReal *src = &p_scalar[i * Nk * npoints]; in PetscDTPTrimmedEvalJet_Internal()
1298 PetscReal *dest = &p[(total++ * Nf + f) * Nk * npoints]; in PetscDTPTrimmedEvalJet_Internal()
1299 PetscCall(PetscArraycpy(dest, src, Nk * npoints)); in PetscDTPTrimmedEvalJet_Internal()
1321 const PetscReal *h_scalar = &p_scalar[(Nbpm1 - Nh) * Nk * npoints]; in PetscDTPTrimmedEvalJet_Internal()
1323 const PetscReal *h_s = &h_scalar[b * Nk * npoints]; in PetscDTPTrimmedEvalJet_Internal()
1331 PetscReal *p_f = &p[total++ * Nf * Nk * npoints]; in PetscDTPTrimmedEvalJet_Internal()
1342 PetscReal *p_i = &p_f[i * Nk * npoints]; in PetscDTPTrimmedEvalJet_Internal()
1344 const PetscReal *h_jet = &h_s[jet * npoints]; in PetscDTPTrimmedEvalJet_Internal()
1345 PetscReal *p_jet = &p_i[jet * npoints]; in PetscDTPTrimmedEvalJet_Internal()
1347 …for (PetscInt pt = 0; pt < npoints; pt++) p_jet[pt] += scale * h_jet[pt] * (points[pt * dim + v] -… in PetscDTPTrimmedEvalJet_Internal()
1354 p_jet = &p_i[l * npoints]; in PetscDTPTrimmedEvalJet_Internal()
1355 for (PetscInt pt = 0; pt < npoints; pt++) p_jet[pt] += scale * mult * h_jet[pt]; in PetscDTPTrimmedEvalJet_Internal()
1409 PetscErrorCode PetscDTPTrimmedEvalJet(PetscInt dim, PetscInt npoints, const PetscReal points[], Pet… in PetscDTPTrimmedEvalJet() argument
1412 …PetscCall(PetscDTPTrimmedEvalJet_Internal(dim, npoints, points, degree, formDegree, jetDegree, p)); in PetscDTPTrimmedEvalJet()
1552 static PetscErrorCode PetscDTGaussJacobiQuadrature_Newton_Internal(PetscInt npoints, PetscReal a, P… in PetscDTGaussJacobiQuadrature_Newton_Internal() argument
1564 a2 = PetscLGamma(a + npoints + 1); in PetscDTGaussJacobiQuadrature_Newton_Internal()
1565 a3 = PetscLGamma(b + npoints + 1); in PetscDTGaussJacobiQuadrature_Newton_Internal()
1566 a4 = PetscLGamma(a + b + npoints + 1); in PetscDTGaussJacobiQuadrature_Newton_Internal()
1567 a5 = PetscLGamma(npoints + 1); in PetscDTGaussJacobiQuadrature_Newton_Internal()
1578 for (k = 0; k < ia; k++) gf *= (npoints + 1. + k) / (npoints + b + 1. + k); in PetscDTGaussJacobiQuadrature_Newton_Internal()
1580 for (k = 0; k < ib; k++) gf *= (npoints + 1. + k) / (npoints + a + 1. + k); in PetscDTGaussJacobiQuadrature_Newton_Internal()
1588 for (k = 0; k < npoints; ++k) { in PetscDTGaussJacobiQuadrature_Newton_Internal()
1589 …PetscReal r = PetscCosReal(PETSC_PI * (1. - (4. * k + 3. + 2. * b) / (4. * npoints + 2. * (a + b +… in PetscDTGaussJacobiQuadrature_Newton_Internal()
1598 PetscCall(PetscDTComputeJacobi(a, b, npoints, r, &f)); in PetscDTGaussJacobiQuadrature_Newton_Internal()
1599 PetscCall(PetscDTComputeJacobiDerivative(a, b, npoints, r, 1, &fp)); in PetscDTGaussJacobiQuadrature_Newton_Internal()
1605 PetscCall(PetscDTComputeJacobiDerivative(a, b, npoints, x[k], 1, &dP)); in PetscDTGaussJacobiQuadrature_Newton_Internal()
1629 static PetscErrorCode PetscDTGaussJacobiQuadrature_GolubWelsch_Internal(PetscInt npoints, PetscReal… in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal() argument
1660 PetscCall(PetscMalloc2(npoints, &diag, npoints, &subdiag)); in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1661 PetscCall(PetscMalloc1(npoints * npoints, &V)); in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1662 PetscCall(PetscDTJacobiMatrix_Internal(npoints, a, b, diag, subdiag)); in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1663 for (i = 0; i < npoints - 1; i++) subdiag[i] = PetscSqrtReal(subdiag[i]); in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1664 PetscCall(PetscDTSymmetricTridiagonalEigensolve(npoints, diag, subdiag, x, V)); in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1665 for (i = 0; i < npoints; i++) w[i] = PetscSqr(PetscRealPart(V[i * npoints])) * mu0; in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1677 PetscCall(PetscSortedReal(npoints, x, &sorted)); in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1682 PetscCall(PetscMalloc2(npoints, &order, npoints, &tmp)); in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1683 for (i = 0; i < npoints; i++) order[i] = i; in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1684 PetscCall(PetscSortRealWithPermutation(npoints, x, order)); in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1685 PetscCall(PetscArraycpy(tmp, x, npoints)); in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1686 for (i = 0; i < npoints; i++) x[i] = tmp[order[i]]; in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1687 PetscCall(PetscArraycpy(tmp, w, npoints)); in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1688 for (i = 0; i < npoints; i++) w[i] = tmp[order[i]]; in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1695 static PetscErrorCode PetscDTGaussJacobiQuadrature_Internal(PetscInt npoints, PetscReal alpha, Pets… in PetscDTGaussJacobiQuadrature_Internal() argument
1698 …PetscCheck(npoints >= 1, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Number of points must be posi… in PetscDTGaussJacobiQuadrature_Internal()
1703 if (newton) PetscCall(PetscDTGaussJacobiQuadrature_Newton_Internal(npoints, alpha, beta, x, w)); in PetscDTGaussJacobiQuadrature_Internal()
1704 else PetscCall(PetscDTGaussJacobiQuadrature_GolubWelsch_Internal(npoints, alpha, beta, x, w)); in PetscDTGaussJacobiQuadrature_Internal()
1707 for (i = 0; i < (npoints + 1) / 2; i++) { in PetscDTGaussJacobiQuadrature_Internal()
1708 PetscInt j = npoints - 1 - i; in PetscDTGaussJacobiQuadrature_Internal()
1746 PetscErrorCode PetscDTGaussJacobiQuadrature(PetscInt npoints, PetscReal a, PetscReal b, PetscReal a… in PetscDTGaussJacobiQuadrature() argument
1751 …PetscCall(PetscDTGaussJacobiQuadrature_Internal(npoints, alpha, beta, x, w, PetscDTGaussQuadrature… in PetscDTGaussJacobiQuadrature()
1753 for (i = 0; i < npoints; i++) { in PetscDTGaussJacobiQuadrature()
1761 static PetscErrorCode PetscDTGaussLobattoJacobiQuadrature_Internal(PetscInt npoints, PetscReal alph… in PetscDTGaussLobattoJacobiQuadrature_Internal() argument
1766 …PetscCheck(npoints >= 2, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Number of points must be posi… in PetscDTGaussLobattoJacobiQuadrature_Internal()
1772 x[npoints - 1] = 1.; in PetscDTGaussLobattoJacobiQuadrature_Internal()
1773 …if (npoints > 2) PetscCall(PetscDTGaussJacobiQuadrature_Internal(npoints - 2, alpha + 1., beta + 1… in PetscDTGaussLobattoJacobiQuadrature_Internal()
1774 for (i = 1; i < npoints - 1; i++) w[i] /= (1. - x[i] * x[i]); in PetscDTGaussLobattoJacobiQuadrature_Internal()
1775 …PetscCall(PetscDTGaussLobattoJacobiEndweights_Internal(npoints, alpha, beta, &w[0], &w[npoints - 1… in PetscDTGaussLobattoJacobiQuadrature_Internal()
1803 PetscErrorCode PetscDTGaussLobattoJacobiQuadrature(PetscInt npoints, PetscReal a, PetscReal b, Pets… in PetscDTGaussLobattoJacobiQuadrature() argument
1808 …PetscCall(PetscDTGaussLobattoJacobiQuadrature_Internal(npoints, alpha, beta, x, w, PetscDTGaussQua… in PetscDTGaussLobattoJacobiQuadrature()
1810 for (i = 0; i < npoints; i++) { in PetscDTGaussLobattoJacobiQuadrature()
1839 PetscErrorCode PetscDTGaussQuadrature(PetscInt npoints, PetscReal a, PetscReal b, PetscReal x[], Pe… in PetscDTGaussQuadrature() argument
1844 …PetscCall(PetscDTGaussJacobiQuadrature_Internal(npoints, 0., 0., x, w, PetscDTGaussQuadratureNewto… in PetscDTGaussQuadrature()
1846 for (i = 0; i < npoints; i++) { in PetscDTGaussQuadrature()
1881 PetscErrorCode PetscDTGaussLobattoLegendreQuadrature(PetscInt npoints, PetscGaussLobattoLegendreCre… in PetscDTGaussLobattoLegendreQuadrature() argument
1886 …PetscCheck(npoints >= 2, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Must provide at least 2 grid … in PetscDTGaussLobattoLegendreQuadrature()
1888 PetscCall(PetscDTGaussLobattoJacobiQuadrature_Internal(npoints, 0., 0., x, w, newton)); in PetscDTGaussLobattoLegendreQuadrature()
1911 PetscErrorCode PetscDTGaussTensorQuadrature(PetscInt dim, PetscInt Nc, PetscInt npoints, PetscReal … in PetscDTGaussTensorQuadrature() argument
1914 …PetscInt totpoints = dim > 1 ? dim > 2 ? npoints * PetscSqr(npoints) : PetscSqr(npoints) : n… in PetscDTGaussTensorQuadrature()
1934 PetscCall(PetscMalloc1(npoints, &ww)); in PetscDTGaussTensorQuadrature()
1935 PetscCall(PetscDTGaussQuadrature(npoints, a, b, x, ww)); in PetscDTGaussTensorQuadrature()
1936 for (PetscInt i = 0; i < npoints; ++i) in PetscDTGaussTensorQuadrature()
1942 PetscCall(PetscMalloc2(npoints, &xw, npoints, &ww)); in PetscDTGaussTensorQuadrature()
1943 PetscCall(PetscDTGaussQuadrature(npoints, a, b, xw, ww)); in PetscDTGaussTensorQuadrature()
1944 for (PetscInt i = 0; i < npoints; ++i) { in PetscDTGaussTensorQuadrature()
1945 for (PetscInt j = 0; j < npoints; ++j) { in PetscDTGaussTensorQuadrature()
1946 x[(i * npoints + j) * dim + 0] = xw[i]; in PetscDTGaussTensorQuadrature()
1947 x[(i * npoints + j) * dim + 1] = xw[j]; in PetscDTGaussTensorQuadrature()
1948 for (PetscInt c = 0; c < Nc; ++c) w[(i * npoints + j) * Nc + c] = ww[i] * ww[j]; in PetscDTGaussTensorQuadrature()
1955 PetscCall(PetscMalloc2(npoints, &xw, npoints, &ww)); in PetscDTGaussTensorQuadrature()
1956 PetscCall(PetscDTGaussQuadrature(npoints, a, b, xw, ww)); in PetscDTGaussTensorQuadrature()
1957 for (PetscInt i = 0; i < npoints; ++i) { in PetscDTGaussTensorQuadrature()
1958 for (PetscInt j = 0; j < npoints; ++j) { in PetscDTGaussTensorQuadrature()
1959 for (PetscInt k = 0; k < npoints; ++k) { in PetscDTGaussTensorQuadrature()
1960 x[((i * npoints + j) * npoints + k) * dim + 0] = xw[i]; in PetscDTGaussTensorQuadrature()
1961 x[((i * npoints + j) * npoints + k) * dim + 1] = xw[j]; in PetscDTGaussTensorQuadrature()
1962 x[((i * npoints + j) * npoints + k) * dim + 2] = xw[k]; in PetscDTGaussTensorQuadrature()
1963 …for (PetscInt c = 0; c < Nc; ++c) w[((i * npoints + j) * npoints + k) * Nc + c] = ww[i] * ww[j] * … in PetscDTGaussTensorQuadrature()
1974 PetscCall(PetscQuadratureSetOrder(*q, 2 * npoints - 1)); in PetscDTGaussTensorQuadrature()
2002 PetscErrorCode PetscDTStroudConicalQuadrature(PetscInt dim, PetscInt Nc, PetscInt npoints, PetscRea… in PetscDTStroudConicalQuadrature() argument
2028 for (PetscInt i = 0; i < dim; ++i) totpoints *= npoints; in PetscDTStroudConicalQuadrature()
2031 PetscCall(PetscMalloc2(npoints, &p1, npoints, &w1)); in PetscDTStroudConicalQuadrature()
2033 for (PetscInt i = 0, totprev = 1, totrem = totpoints / npoints; i < dim; ++i) { in PetscDTStroudConicalQuadrature()
2037 PetscCall(PetscDTGaussJacobiQuadrature(npoints, -1., 1., i, 0.0, p1, w1)); in PetscDTStroudConicalQuadrature()
2039 for (PetscInt j = 0; j < npoints; j++) { in PetscDTStroudConicalQuadrature()
2047 totprev *= npoints; in PetscDTStroudConicalQuadrature()
2048 totrem /= npoints; in PetscDTStroudConicalQuadrature()
2053 PetscCall(PetscQuadratureSetOrder(*q, 2 * npoints - 1)); in PetscDTStroudConicalQuadrature()
2363 PetscInt K, k, npoints; in PetscDTTanhSinhTensorQuadrature() local
2389 npoints = 2 * K - 1; in PetscDTTanhSinhTensorQuadrature()
2390 PetscCall(PetscMalloc1(npoints * dim, &x)); in PetscDTTanhSinhTensorQuadrature()
2391 PetscCall(PetscMalloc1(npoints, &w)); in PetscDTTanhSinhTensorQuadrature()
2403 PetscCall(PetscQuadratureSetData(*q, dim, 1, npoints, x, w)); in PetscDTTanhSinhTensorQuadrature()