Lines Matching refs:q
60 PetscErrorCode PetscQuadratureCreate(MPI_Comm comm, PetscQuadrature *q) in PetscQuadratureCreate() argument
63 PetscAssertPointer(q, 2); in PetscQuadratureCreate()
65 …PetscCall(PetscHeaderCreate(*q, PETSCQUADRATURE_CLASSID, "PetscQuadrature", "Quadrature", "DT", co… in PetscQuadratureCreate()
66 (*q)->ct = DM_POLYTOPE_UNKNOWN; in PetscQuadratureCreate()
67 (*q)->dim = -1; in PetscQuadratureCreate()
68 (*q)->Nc = 1; in PetscQuadratureCreate()
69 (*q)->order = -1; in PetscQuadratureCreate()
70 (*q)->numPoints = 0; in PetscQuadratureCreate()
71 (*q)->points = NULL; in PetscQuadratureCreate()
72 (*q)->weights = NULL; in PetscQuadratureCreate()
91 PetscErrorCode PetscQuadratureDuplicate(PetscQuadrature q, PetscQuadrature *r) in PetscQuadratureDuplicate() argument
99 PetscAssertPointer(q, 1); in PetscQuadratureDuplicate()
100 PetscCall(PetscQuadratureCreate(PetscObjectComm((PetscObject)q), r)); in PetscQuadratureDuplicate()
101 PetscCall(PetscQuadratureGetCellType(q, &ct)); in PetscQuadratureDuplicate()
103 PetscCall(PetscQuadratureGetOrder(q, &order)); in PetscQuadratureDuplicate()
105 PetscCall(PetscQuadratureGetData(q, &dim, &Nc, &Nq, &points, &weights)); in PetscQuadratureDuplicate()
126 PetscErrorCode PetscQuadratureDestroy(PetscQuadrature *q) in PetscQuadratureDestroy() argument
129 if (!*q) PetscFunctionReturn(PETSC_SUCCESS); in PetscQuadratureDestroy()
130 PetscValidHeaderSpecific(*q, PETSCQUADRATURE_CLASSID, 1); in PetscQuadratureDestroy()
131 if (--((PetscObject)*q)->refct > 0) { in PetscQuadratureDestroy()
132 *q = NULL; in PetscQuadratureDestroy()
135 PetscCall(PetscFree((*q)->points)); in PetscQuadratureDestroy()
136 PetscCall(PetscFree((*q)->weights)); in PetscQuadratureDestroy()
137 PetscCall(PetscHeaderDestroy(q)); in PetscQuadratureDestroy()
156 PetscErrorCode PetscQuadratureGetCellType(PetscQuadrature q, DMPolytopeType *ct) in PetscQuadratureGetCellType() argument
159 PetscValidHeaderSpecific(q, PETSCQUADRATURE_CLASSID, 1); in PetscQuadratureGetCellType()
161 *ct = q->ct; in PetscQuadratureGetCellType()
178 PetscErrorCode PetscQuadratureSetCellType(PetscQuadrature q, DMPolytopeType ct) in PetscQuadratureSetCellType() argument
181 PetscValidHeaderSpecific(q, PETSCQUADRATURE_CLASSID, 1); in PetscQuadratureSetCellType()
182 q->ct = ct; in PetscQuadratureSetCellType()
201 PetscErrorCode PetscQuadratureGetOrder(PetscQuadrature q, PetscInt *order) in PetscQuadratureGetOrder() argument
204 PetscValidHeaderSpecific(q, PETSCQUADRATURE_CLASSID, 1); in PetscQuadratureGetOrder()
206 *order = q->order; in PetscQuadratureGetOrder()
223 PetscErrorCode PetscQuadratureSetOrder(PetscQuadrature q, PetscInt order) in PetscQuadratureSetOrder() argument
226 PetscValidHeaderSpecific(q, PETSCQUADRATURE_CLASSID, 1); in PetscQuadratureSetOrder()
227 q->order = order; in PetscQuadratureSetOrder()
249 PetscErrorCode PetscQuadratureGetNumComponents(PetscQuadrature q, PetscInt *Nc) in PetscQuadratureGetNumComponents() argument
252 PetscValidHeaderSpecific(q, PETSCQUADRATURE_CLASSID, 1); in PetscQuadratureGetNumComponents()
254 *Nc = q->Nc; in PetscQuadratureGetNumComponents()
274 PetscErrorCode PetscQuadratureSetNumComponents(PetscQuadrature q, PetscInt Nc) in PetscQuadratureSetNumComponents() argument
277 PetscValidHeaderSpecific(q, PETSCQUADRATURE_CLASSID, 1); in PetscQuadratureSetNumComponents()
278 q->Nc = Nc; in PetscQuadratureSetNumComponents()
307 PetscErrorCode PetscQuadratureGetData(PetscQuadrature q, PeOp PetscInt *dim, PeOp PetscInt *Nc, PeO… in PetscQuadratureGetData() argument
310 PetscValidHeaderSpecific(q, PETSCQUADRATURE_CLASSID, 1); in PetscQuadratureGetData()
313 *dim = q->dim; in PetscQuadratureGetData()
317 *Nc = q->Nc; in PetscQuadratureGetData()
321 *npoints = q->numPoints; in PetscQuadratureGetData()
325 *points = q->points; in PetscQuadratureGetData()
329 *weights = q->weights; in PetscQuadratureGetData()
495 PetscErrorCode PetscQuadraturePushForward(PetscQuadrature q, PetscInt imageDim, const PetscReal ori… in PetscQuadraturePushForward() argument
505 PetscValidHeaderSpecific(q, PETSCQUADRATURE_CLASSID, 1); in PetscQuadraturePushForward()
506 …PetscCheck(imageDim >= PetscAbsInt(formDegree), PetscObjectComm((PetscObject)q), PETSC_ERR_ARG_INC… in PetscQuadraturePushForward()
507 PetscCall(PetscQuadratureGetData(q, &dim, &Nc, &Npoints, &points, &weights)); in PetscQuadraturePushForward()
509 …PetscCheck(Nc % formSize == 0, PetscObjectComm((PetscObject)q), PETSC_ERR_ARG_INCOMP, "Number of c… in PetscQuadraturePushForward()
540 PetscCall(PetscQuadratureCreate(PetscObjectComm((PetscObject)q), Jinvstarq)); in PetscQuadraturePushForward()
566 PetscErrorCode PetscQuadratureSetData(PetscQuadrature q, PetscInt dim, PetscInt Nc, PetscInt npoint… in PetscQuadratureSetData() argument
569 PetscValidHeaderSpecific(q, PETSCQUADRATURE_CLASSID, 1); in PetscQuadratureSetData()
570 if (dim >= 0) q->dim = dim; in PetscQuadratureSetData()
571 if (Nc >= 0) q->Nc = Nc; in PetscQuadratureSetData()
572 if (npoints >= 0) q->numPoints = npoints; in PetscQuadratureSetData()
575 q->points = points; in PetscQuadratureSetData()
579 q->weights = weights; in PetscQuadratureSetData()
586 PetscInt q, d, c; in PetscQuadratureView_Ascii() local
595 for (q = 0; q < quad->numPoints; ++q) { in PetscQuadratureView_Ascii()
596 PetscCall(PetscViewerASCIIPrintf(v, "p%" PetscInt_FMT " (", q)); in PetscQuadratureView_Ascii()
600 PetscCall(PetscViewerASCIIPrintf(v, "%+g", (double)quad->points[q * quad->dim + d])); in PetscQuadratureView_Ascii()
603 if (quad->Nc > 1) PetscCall(PetscViewerASCIIPrintf(v, "w%" PetscInt_FMT " (", q)); in PetscQuadratureView_Ascii()
606 PetscCall(PetscViewerASCIIPrintf(v, "%+g", (double)quad->weights[q * quad->Nc + c])); in PetscQuadratureView_Ascii()
664 PetscErrorCode PetscQuadratureExpandComposite(PetscQuadrature q, PetscInt numSubelements, const Pet… in PetscQuadratureExpandComposite() argument
672 PetscValidHeaderSpecific(q, PETSCQUADRATURE_CLASSID, 1); in PetscQuadratureExpandComposite()
677 PetscCall(PetscQuadratureGetCellType(q, &ct)); in PetscQuadratureExpandComposite()
678 PetscCall(PetscQuadratureGetOrder(q, &order)); in PetscQuadratureExpandComposite()
679 PetscCall(PetscQuadratureGetData(q, &dim, &Nc, &npoints, &points, &weights)); in PetscQuadratureExpandComposite()
1911 …adrature(PetscInt dim, PetscInt Nc, PetscInt npoints, PetscReal a, PetscReal b, PetscQuadrature *q) in PetscDTGaussTensorQuadrature() argument
1972 PetscCall(PetscQuadratureCreate(PETSC_COMM_SELF, q)); in PetscDTGaussTensorQuadrature()
1973 PetscCall(PetscQuadratureSetCellType(*q, ct)); in PetscDTGaussTensorQuadrature()
1974 PetscCall(PetscQuadratureSetOrder(*q, 2 * npoints - 1)); in PetscDTGaussTensorQuadrature()
1975 PetscCall(PetscQuadratureSetData(*q, dim, Nc, totpoints, x, w)); in PetscDTGaussTensorQuadrature()
1976 PetscCall(PetscObjectChangeTypeName((PetscObject)*q, "GaussTensor")); in PetscDTGaussTensorQuadrature()
2002 …adrature(PetscInt dim, PetscInt Nc, PetscInt npoints, PetscReal a, PetscReal b, PetscQuadrature *q) in PetscDTStroudConicalQuadrature() argument
2051 PetscCall(PetscQuadratureCreate(PETSC_COMM_SELF, q)); in PetscDTStroudConicalQuadrature()
2052 PetscCall(PetscQuadratureSetCellType(*q, ct)); in PetscDTStroudConicalQuadrature()
2053 PetscCall(PetscQuadratureSetOrder(*q, 2 * npoints - 1)); in PetscDTStroudConicalQuadrature()
2054 PetscCall(PetscQuadratureSetData(*q, dim, Nc, totpoints, x, w)); in PetscDTStroudConicalQuadrature()
2055 PetscCall(PetscObjectChangeTypeName((PetscObject)*q, "StroudConical")); in PetscDTStroudConicalQuadrature()
2243 PetscQuadrature q; in PetscDTSimplexQuadrature() local
2326 PetscCall(PetscQuadratureCreate(PETSC_COMM_SELF, &q)); in PetscDTSimplexQuadrature()
2327 PetscCall(PetscQuadratureSetCellType(q, ct)); in PetscDTSimplexQuadrature()
2328 PetscCall(PetscQuadratureSetOrder(q, degree)); in PetscDTSimplexQuadrature()
2329 PetscCall(PetscQuadratureSetData(q, dim, 1, num_full_nodes, points, weights)); in PetscDTSimplexQuadrature()
2330 *quad = q; in PetscDTSimplexQuadrature()
2353 …anhSinhTensorQuadrature(PetscInt dim, PetscInt level, PetscReal a, PetscReal b, PetscQuadrature *q) in PetscDTTanhSinhTensorQuadrature() argument
2386 PetscCall(PetscQuadratureCreate(PETSC_COMM_SELF, q)); in PetscDTTanhSinhTensorQuadrature()
2387 PetscCall(PetscQuadratureSetCellType(*q, ct)); in PetscDTTanhSinhTensorQuadrature()
2388 PetscCall(PetscQuadratureSetOrder(*q, 2 * K + 1)); in PetscDTTanhSinhTensorQuadrature()
2403 PetscCall(PetscQuadratureSetData(*q, dim, 1, npoints, x, w)); in PetscDTTanhSinhTensorQuadrature()
2612 …ErrorCode PetscDTTensorQuadratureCreate(PetscQuadrature q1, PetscQuadrature q2, PetscQuadrature *q) in PetscDTTensorQuadratureCreate() argument
2624 PetscAssertPointer(q, 3); in PetscDTTensorQuadratureCreate()
2766 PetscCall(PetscQuadratureCreate(PETSC_COMM_SELF, q)); in PetscDTTensorQuadratureCreate()
2767 PetscCall(PetscQuadratureSetCellType(*q, ct)); in PetscDTTensorQuadratureCreate()
2768 PetscCall(PetscQuadratureSetOrder(*q, order)); in PetscDTTensorQuadratureCreate()
2778 PetscCall(PetscQuadratureSetData(*q, dim, Nc, Np, x, w)); in PetscDTTensorQuadratureCreate()
3379 PetscInt dim, qdim, d, Na, o, Nq, q, qp; in PetscQuadratureComputePermutations() local
3399 for (q = 0; q < Nq; ++q) { in PetscQuadratureComputePermutations()
3400 CoordinatesRefToReal(dim, dim, xi0, v0, J, &xq[q * dim], txq); in PetscQuadratureComputePermutations()
3407 …, PETSC_ERR_PLIB, "Transformed quad point %" PetscInt_FMT " does not match another quad point", q); in PetscQuadratureComputePermutations()
3408 idx[q] = qp; in PetscQuadratureComputePermutations()
3438 …e ct, PetscInt qorder, PetscDTSimplexQuadratureType qtype, PetscQuadrature *q, PetscQuadrature *fq) in PetscDTCreateQuadratureByCell() argument
3451 PetscCall(PetscDTGaussTensorQuadrature(dim, 1, quadPointsPerEdge, -1.0, 1.0, q)); in PetscDTCreateQuadratureByCell()
3456 PetscCall(PetscDTSimplexQuadrature(dim, 2 * qorder, qtype, q)); in PetscDTCreateQuadratureByCell()
3466 PetscCall(PetscDTTensorQuadratureCreate(q1, q2, q)); in PetscDTCreateQuadratureByCell()
3494 …reateDefaultQuadrature(DMPolytopeType ct, PetscInt qorder, PetscQuadrature *q, PetscQuadrature *fq) in PetscDTCreateDefaultQuadrature() argument
3497 PetscCall(PetscDTCreateQuadratureByCell(ct, qorder, PETSCDTSIMPLEXQUAD_DEFAULT, q, fq)); in PetscDTCreateDefaultQuadrature()