Lines Matching refs:quad
584 static PetscErrorCode PetscQuadratureView_Ascii(PetscQuadrature quad, PetscViewer v) in PetscQuadratureView_Ascii() argument
590 if (quad->Nc > 1) in PetscQuadratureView_Ascii()
591 …" PetscInt_FMT " components\n", DMPolytopeTypes[quad->ct], quad->order, quad->numPoints, quad->dim… in PetscQuadratureView_Ascii()
592 … " points (dim %" PetscInt_FMT ")\n", DMPolytopeTypes[quad->ct], quad->order, quad->numPoints, qua… in PetscQuadratureView_Ascii()
595 for (q = 0; q < quad->numPoints; ++q) { in PetscQuadratureView_Ascii()
598 for (d = 0; d < quad->dim; ++d) { 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()
604 for (c = 0; c < quad->Nc; ++c) { in PetscQuadratureView_Ascii()
606 PetscCall(PetscViewerASCIIPrintf(v, "%+g", (double)quad->weights[q * quad->Nc + c])); in PetscQuadratureView_Ascii()
608 if (quad->Nc > 1) PetscCall(PetscViewerASCIIPrintf(v, ")")); in PetscQuadratureView_Ascii()
628 PetscErrorCode PetscQuadratureView(PetscQuadrature quad, PetscViewer viewer) in PetscQuadratureView() argument
633 PetscValidHeader(quad, 1); in PetscQuadratureView()
635 if (!viewer) PetscCall(PetscViewerASCIIGetStdout(PetscObjectComm((PetscObject)quad), &viewer)); in PetscQuadratureView()
638 if (isascii) PetscCall(PetscQuadratureView_Ascii(quad, viewer)); in PetscQuadratureView()
2140 …Quadrature(PetscInt dim, PetscInt degree, PetscDTSimplexQuadratureType type, PetscQuadrature *quad) in PetscDTSimplexQuadrature() argument
2150 PetscCall(PetscDTStroudConicalQuadrature(dim, 1, points_per_dim, -1, 1, quad)); in PetscDTSimplexQuadrature()
2225 PetscCall(PetscDTSimplexQuadrature(dim, degree, PETSCDTSIMPLEXQUAD_CONIC, quad)); in PetscDTSimplexQuadrature()
2330 *quad = q; in PetscDTSimplexQuadrature()
3375 PetscErrorCode PetscQuadratureComputePermutations(PetscQuadrature quad, PeOp PetscInt *Np, IS *perm… in PetscQuadratureComputePermutations() argument
3382 PetscCall(PetscQuadratureGetData(quad, &qdim, NULL, &Nq, &xq, &wq)); in PetscQuadratureComputePermutations()
3383 PetscCall(PetscQuadratureGetCellType(quad, &ct)); in PetscQuadratureComputePermutations()