| /petsc/src/dm/dt/tests/output/ |
| H A D | ex1_1.out | 2 degree 1 at 0: B= 0 D= 1 D2= 0 3 degree 2 at 0: B= -0.5 D= 0 D2= 3 4 degree 3 at 0: B= 0 D= -1.5 D2= 0 5 degree 4 at 0: B= 0.375 D= 0 D2= -7.5 6 degree 5 at 0: B= 0 D= 1.875 D2= 0 7 degree 1 at 0.2: B= 0.2 D= 1 D2= 0 8 degree 2 at 0.2: B= -0.44 D= 0.6 D2= 3 9 degree 3 at 0.2: B= -0.28 D= -1.2 D2= 3 10 degree 4 at 0.2: B= 0.232 D= -1.36 D2= -5.4 11 degree 5 at 0.2: B= 0.3075 D= 0.888 D2= -9.24 [all …]
|
| /petsc/src/dm/dt/space/impls/ptrimmed/tests/output/ |
| H A D | ex1_1.out | 4 Polynomial space of degree 1 12 Polynomial space of degree 1 16 Polynomial space of degree 2 24 Polynomial space of degree 2 28 Polynomial space of degree 3 36 Polynomial space of degree 3 40 Polynomial space of degree 4 48 Polynomial space of degree 4 52 Polynomial space of degree 0 60 Polynomial space of degree 0 [all …]
|
| /petsc/src/dm/dt/space/impls/ptrimmed/tests/ |
| H A D | ex1.c | 5 static PetscErrorCode test(PetscInt dim, PetscInt formDegree, PetscInt degree, PetscInt nCopies) in test() argument 24 PetscCall(PetscSpaceSetDegree(sp, degree, PETSC_DETERMINE)); in test() 29 …PetscCall(PetscDTPTrimmedSize(dim, formDegree == 0 ? degree : degree + 1, PetscAbsInt(formDegree),… in test() 34 maxDexp = (PetscAbsInt(formDegree) == dim || formDegree == 0) ? degree : degree + 1; in test() 36 …PetscCheck(degree == d, comm, PETSC_ERR_PLIB, "Space degree mismatch, %" PetscInt_FMT " != %" Pets… in test() 62 for (PetscInt degree = 0; degree <= 4; degree++) { in main() local 63 if (formDegree == 0 && degree == 0) continue; in main() 64 …ies = 1; nCopies <= PetscMax(2, dim); nCopies++) PetscCall(test(dim, formDegree, degree, nCopies)); in main()
|
| /petsc/src/dm/dt/dualspace/impls/lagrange/tests/output/ |
| H A D | ex1_3_simplex_discontinuous_full.out | 1 New space: tetrahedron, order 0, trimmed 0, tensor 0, continuous 0, form degree -2 10 New space: tetrahedron, order 0, trimmed 0, tensor 0, continuous 0, form degree -1 19 New space: tetrahedron, order 0, trimmed 0, tensor 0, continuous 0, form degree 0 26 New space: tetrahedron, order 0, trimmed 0, tensor 0, continuous 0, form degree 1 35 New space: tetrahedron, order 0, trimmed 0, tensor 0, continuous 0, form degree 2 44 New space: tetrahedron, order 0, trimmed 0, tensor 0, continuous 0, form degree 3 51 New space: tetrahedron, order 1, trimmed 0, tensor 0, continuous 0, form degree -2 69 New space: tetrahedron, order 1, trimmed 0, tensor 0, continuous 0, form degree -1 87 New space: tetrahedron, order 1, trimmed 0, tensor 0, continuous 0, form degree 0 97 New space: tetrahedron, order 1, trimmed 0, tensor 0, continuous 0, form degree 1 [all …]
|
| H A D | ex1_3_simplex_continuous_full.out | 1 New space: tetrahedron, order 0, trimmed 0, tensor 0, continuous 1, form degree -2 10 New space: tetrahedron, order 0, trimmed 0, tensor 0, continuous 1, form degree -1 19 New space: tetrahedron, order 0, trimmed 0, tensor 0, continuous 1, form degree 0 26 New space: tetrahedron, order 0, trimmed 0, tensor 0, continuous 1, form degree 1 35 New space: tetrahedron, order 0, trimmed 0, tensor 0, continuous 1, form degree 2 44 New space: tetrahedron, order 0, trimmed 0, tensor 0, continuous 1, form degree 3 51 New space: tetrahedron, order 1, trimmed 0, tensor 0, continuous 1, form degree -2 68 New space: tetrahedron, order 1, trimmed 0, tensor 0, continuous 1, form degree -1 85 New space: tetrahedron, order 1, trimmed 0, tensor 0, continuous 1, form degree 0 94 New space: tetrahedron, order 1, trimmed 0, tensor 0, continuous 1, form degree 1 [all …]
|
| H A D | ex1_3_wedge_discontinuous_full.out | 1 New space: triangular_prism, order 0, trimmed 0, tensor 2, continuous 0, form degree -2 10 New space: triangular_prism, order 0, trimmed 0, tensor 2, continuous 0, form degree -1 19 New space: triangular_prism, order 0, trimmed 0, tensor 2, continuous 0, form degree 0 26 New space: triangular_prism, order 0, trimmed 0, tensor 2, continuous 0, form degree 1 35 New space: triangular_prism, order 0, trimmed 0, tensor 2, continuous 0, form degree 2 44 New space: triangular_prism, order 0, trimmed 0, tensor 2, continuous 0, form degree 3 51 New space: triangular_prism, order 1, trimmed 0, tensor 2, continuous 0, form degree -2 75 New space: triangular_prism, order 1, trimmed 0, tensor 2, continuous 0, form degree -1 99 New space: triangular_prism, order 1, trimmed 0, tensor 2, continuous 0, form degree 0 111 New space: triangular_prism, order 1, trimmed 0, tensor 2, continuous 0, form degree 1 [all …]
|
| H A D | ex1_3_tensor_discontinuous_full.out | 1 New space: hexahedron, order 0, trimmed 0, tensor 1, continuous 0, form degree -2 10 New space: hexahedron, order 0, trimmed 0, tensor 1, continuous 0, form degree -1 19 New space: hexahedron, order 0, trimmed 0, tensor 1, continuous 0, form degree 0 26 New space: hexahedron, order 0, trimmed 0, tensor 1, continuous 0, form degree 1 35 New space: hexahedron, order 0, trimmed 0, tensor 1, continuous 0, form degree 2 44 New space: hexahedron, order 0, trimmed 0, tensor 1, continuous 0, form degree 3 51 New space: hexahedron, order 1, trimmed 0, tensor 1, continuous 0, form degree -2 81 New space: hexahedron, order 1, trimmed 0, tensor 1, continuous 0, form degree -1 111 New space: hexahedron, order 1, trimmed 0, tensor 1, continuous 0, form degree 0 125 New space: hexahedron, order 1, trimmed 0, tensor 1, continuous 0, form degree 1 [all …]
|
| H A D | ex1_3_wedge_continuous_trimmed.out | 1 New space: triangular_prism, order 1, trimmed 1, tensor 2, continuous 1, form degree -2 11 New space: triangular_prism, order 1, trimmed 1, tensor 2, continuous 1, form degree -1 25 New space: triangular_prism, order 1, trimmed 1, tensor 2, continuous 1, form degree 0 36 New space: triangular_prism, order 1, trimmed 1, tensor 2, continuous 1, form degree 1 50 New space: triangular_prism, order 1, trimmed 1, tensor 2, continuous 1, form degree 2 60 New space: triangular_prism, order 1, trimmed 1, tensor 2, continuous 1, form degree 3 67 New space: triangular_prism, order 2, trimmed 1, tensor 2, continuous 1, form degree -2 100 New space: triangular_prism, order 2, trimmed 1, tensor 2, continuous 1, form degree -1 144 New space: triangular_prism, order 2, trimmed 1, tensor 2, continuous 1, form degree 0 167 New space: triangular_prism, order 2, trimmed 1, tensor 2, continuous 1, form degree 1 [all …]
|
| H A D | ex1_3_tensor_continuous_full.out | 1 New space: hexahedron, order 0, trimmed 0, tensor 1, continuous 1, form degree -2 10 New space: hexahedron, order 0, trimmed 0, tensor 1, continuous 1, form degree -1 19 New space: hexahedron, order 0, trimmed 0, tensor 1, continuous 1, form degree 0 26 New space: hexahedron, order 0, trimmed 0, tensor 1, continuous 1, form degree 1 35 New space: hexahedron, order 0, trimmed 0, tensor 1, continuous 1, form degree 2 44 New space: hexahedron, order 0, trimmed 0, tensor 1, continuous 1, form degree 3 51 New space: hexahedron, order 1, trimmed 0, tensor 1, continuous 1, form degree -2 80 New space: hexahedron, order 1, trimmed 0, tensor 1, continuous 1, form degree -1 109 New space: hexahedron, order 1, trimmed 0, tensor 1, continuous 1, form degree 0 122 New space: hexahedron, order 1, trimmed 0, tensor 1, continuous 1, form degree 1 [all …]
|
| H A D | ex1_3_wedge_discontinuous_trimmed.out | 1 New space: triangular_prism, order 1, trimmed 1, tensor 2, continuous 0, form degree -2 12 New space: triangular_prism, order 1, trimmed 1, tensor 2, continuous 0, form degree -1 27 New space: triangular_prism, order 1, trimmed 1, tensor 2, continuous 0, form degree 0 39 New space: triangular_prism, order 1, trimmed 1, tensor 2, continuous 0, form degree 1 54 New space: triangular_prism, order 1, trimmed 1, tensor 2, continuous 0, form degree 2 65 New space: triangular_prism, order 1, trimmed 1, tensor 2, continuous 0, form degree 3 72 New space: triangular_prism, order 2, trimmed 1, tensor 2, continuous 0, form degree -2 103 New space: triangular_prism, order 2, trimmed 1, tensor 2, continuous 0, form degree -1 145 New space: triangular_prism, order 2, trimmed 1, tensor 2, continuous 0, form degree 0 169 New space: triangular_prism, order 2, trimmed 1, tensor 2, continuous 0, form degree 1 [all …]
|
| H A D | ex1_3_simplex_discontinuous_trimmed.out | 1 New space: tetrahedron, order 1, trimmed 1, tensor 0, continuous 0, form degree -2 11 New space: tetrahedron, order 1, trimmed 1, tensor 0, continuous 0, form degree -1 23 New space: tetrahedron, order 1, trimmed 1, tensor 0, continuous 0, form degree 0 33 New space: tetrahedron, order 1, trimmed 1, tensor 0, continuous 0, form degree 1 45 New space: tetrahedron, order 1, trimmed 1, tensor 0, continuous 0, form degree 2 55 New space: tetrahedron, order 1, trimmed 1, tensor 0, continuous 0, form degree 3 62 New space: tetrahedron, order 2, trimmed 1, tensor 0, continuous 0, form degree -2 83 New space: tetrahedron, order 2, trimmed 1, tensor 0, continuous 0, form degree -1 109 New space: tetrahedron, order 2, trimmed 1, tensor 0, continuous 0, form degree 0 125 New space: tetrahedron, order 2, trimmed 1, tensor 0, continuous 0, form degree 1 [all …]
|
| H A D | ex1_3_simplex_continuous_trimmed.out | 1 New space: tetrahedron, order 1, trimmed 1, tensor 0, continuous 1, form degree -2 10 New space: tetrahedron, order 1, trimmed 1, tensor 0, continuous 1, form degree -1 21 New space: tetrahedron, order 1, trimmed 1, tensor 0, continuous 1, form degree 0 30 New space: tetrahedron, order 1, trimmed 1, tensor 0, continuous 1, form degree 1 41 New space: tetrahedron, order 1, trimmed 1, tensor 0, continuous 1, form degree 2 50 New space: tetrahedron, order 1, trimmed 1, tensor 0, continuous 1, form degree 3 57 New space: tetrahedron, order 2, trimmed 1, tensor 0, continuous 1, form degree -2 80 New space: tetrahedron, order 2, trimmed 1, tensor 0, continuous 1, form degree -1 105 New space: tetrahedron, order 2, trimmed 1, tensor 0, continuous 1, form degree 0 120 New space: tetrahedron, order 2, trimmed 1, tensor 0, continuous 1, form degree 1 [all …]
|
| H A D | ex1_3_wedge_continuous_full.out | 1 New space: triangular_prism, order 0, trimmed 0, tensor 2, continuous 1, form degree -2 10 New space: triangular_prism, order 0, trimmed 0, tensor 2, continuous 1, form degree -1 19 New space: triangular_prism, order 0, trimmed 0, tensor 2, continuous 1, form degree 0 26 New space: triangular_prism, order 0, trimmed 0, tensor 2, continuous 1, form degree 1 35 New space: triangular_prism, order 0, trimmed 0, tensor 2, continuous 1, form degree 2 44 New space: triangular_prism, order 0, trimmed 0, tensor 2, continuous 1, form degree 3 51 New space: triangular_prism, order 1, trimmed 0, tensor 2, continuous 1, form degree -2 74 New space: triangular_prism, order 1, trimmed 0, tensor 2, continuous 1, form degree -1 97 New space: triangular_prism, order 1, trimmed 0, tensor 2, continuous 1, form degree 0 108 New space: triangular_prism, order 1, trimmed 0, tensor 2, continuous 1, form degree 1 [all …]
|
| H A D | ex1_3_tensor_continuous_trimmed.out | 1 New space: hexahedron, order 1, trimmed 1, tensor 1, continuous 1, form degree -2 12 New space: hexahedron, order 1, trimmed 1, tensor 1, continuous 1, form degree -1 29 New space: hexahedron, order 1, trimmed 1, tensor 1, continuous 1, form degree 0 42 New space: hexahedron, order 1, trimmed 1, tensor 1, continuous 1, form degree 1 59 New space: hexahedron, order 1, trimmed 1, tensor 1, continuous 1, form degree 2 70 New space: hexahedron, order 1, trimmed 1, tensor 1, continuous 1, form degree 3 77 New space: hexahedron, order 2, trimmed 1, tensor 1, continuous 1, form degree -2 121 New space: hexahedron, order 2, trimmed 1, tensor 1, continuous 1, form degree -1 183 New space: hexahedron, order 2, trimmed 1, tensor 1, continuous 1, form degree 0 218 New space: hexahedron, order 2, trimmed 1, tensor 1, continuous 1, form degree 1 [all …]
|
| H A D | ex1_3_tensor_discontinuous_trimmed.out | 1 New space: hexahedron, order 1, trimmed 1, tensor 1, continuous 0, form degree -2 13 New space: hexahedron, order 1, trimmed 1, tensor 1, continuous 0, form degree -1 31 New space: hexahedron, order 1, trimmed 1, tensor 1, continuous 0, form degree 0 45 New space: hexahedron, order 1, trimmed 1, tensor 1, continuous 0, form degree 1 63 New space: hexahedron, order 1, trimmed 1, tensor 1, continuous 0, form degree 2 75 New space: hexahedron, order 1, trimmed 1, tensor 1, continuous 0, form degree 3 82 New space: hexahedron, order 2, trimmed 1, tensor 1, continuous 0, form degree -2 124 New space: hexahedron, order 2, trimmed 1, tensor 1, continuous 0, form degree -1 184 New space: hexahedron, order 2, trimmed 1, tensor 1, continuous 0, form degree 0 217 New space: hexahedron, order 2, trimmed 1, tensor 1, continuous 0, form degree 1 [all …]
|
| H A D | ex1_1_continuous_full.out | 1 New space: segment, order 0, trimmed 0, tensor 0, continuous 1, form degree 0 16 New space: segment, order 0, trimmed 0, tensor 0, continuous 1, form degree 1 31 New space: segment, order 1, trimmed 0, tensor 0, continuous 1, form degree 0 38 New space: segment, order 1, trimmed 0, tensor 0, continuous 1, form degree 1 56 New space: segment, order 2, trimmed 0, tensor 0, continuous 1, form degree 0 75 New space: segment, order 2, trimmed 0, tensor 0, continuous 1, form degree 1 96 New space: segment, order 3, trimmed 0, tensor 0, continuous 1, form degree 0 118 New space: segment, order 3, trimmed 0, tensor 0, continuous 1, form degree 1
|
| H A D | ex1_1_discontinuous_full.out | 1 New space: segment, order 0, trimmed 0, tensor 0, continuous 0, form degree 0 16 New space: segment, order 0, trimmed 0, tensor 0, continuous 0, form degree 1 31 New space: segment, order 1, trimmed 0, tensor 0, continuous 0, form degree 0 49 New space: segment, order 1, trimmed 0, tensor 0, continuous 0, form degree 1 67 New space: segment, order 2, trimmed 0, tensor 0, continuous 0, form degree 0 88 New space: segment, order 2, trimmed 0, tensor 0, continuous 0, form degree 1 109 New space: segment, order 3, trimmed 0, tensor 0, continuous 0, form degree 0 133 New space: segment, order 3, trimmed 0, tensor 0, continuous 0, form degree 1
|
| H A D | ex1_1_continuous_trimmed.out | 1 New space: segment, order 1, trimmed 1, tensor 0, continuous 1, form degree 0 8 New space: segment, order 1, trimmed 1, tensor 0, continuous 1, form degree 1 23 New space: segment, order 2, trimmed 1, tensor 0, continuous 1, form degree 0 42 New space: segment, order 2, trimmed 1, tensor 0, continuous 1, form degree 1 60 New space: segment, order 3, trimmed 1, tensor 0, continuous 1, form degree 0 82 New space: segment, order 3, trimmed 1, tensor 0, continuous 1, form degree 1
|
| H A D | ex1_1_discontinuous_trimmed.out | 1 New space: segment, order 1, trimmed 1, tensor 0, continuous 0, form degree 0 19 New space: segment, order 1, trimmed 1, tensor 0, continuous 0, form degree 1 34 New space: segment, order 2, trimmed 1, tensor 0, continuous 0, form degree 0 55 New space: segment, order 2, trimmed 1, tensor 0, continuous 0, form degree 1 73 New space: segment, order 3, trimmed 1, tensor 0, continuous 0, form degree 0 97 New space: segment, order 3, trimmed 1, tensor 0, continuous 0, form degree 1
|
| /petsc/src/dm/dt/tests/ |
| H A D | ex15.c | 7 static PetscErrorCode testQuadrature(PetscInt dim, PetscInt degree, PetscDTSimplexQuadratureType ty… in testQuadrature() argument 12 PetscInt p_degree = (degree + 1) / 2; in testQuadrature() 13 PetscInt p_degree_min = degree - p_degree; in testQuadrature() 19 PetscCall(PetscDTSimplexQuadrature(dim, degree, type, &quad)); in testQuadrature() 33 …ror in <P_PKD(%" PetscInt_FMT "),P_PKD%" PetscInt_FMT "d)> = %g\n", dim, degree, PetscDTSimplexQua… in testQuadrature() 34 …error in <P_PKD(%" PetscInt_FMT "),P_PKD(%" PetscInt_FMT ")> = %g", dim, degree, PetscDTSimplexQua… in testQuadrature()
|
| /petsc/src/dm/dt/space/impls/ptrimmed/ |
| H A D | spaceptrimmed.c | 21 tdegree = f == 0 ? sp->degree : sp->degree + 1; in PetscSpacePTrimmedView_Ascii() 22 …Int_FMT "/\\%" PetscInt_FMT ")\n", PetscAbsInt(f), f < 0 ? "*" : "", sp->degree, tdegree, PetscAbs… in PetscSpacePTrimmedView_Ascii() 100 if (sp->degree == PETSC_DEFAULT) sp->degree = 0; in PetscSpaceSetUp_Ptrimmed() 101 …etscCheck(sp->degree >= 0, PetscObjectComm((PetscObject)sp), PETSC_ERR_ARG_OUTOFRANGE, "Invalid ne… in PetscSpaceSetUp_Ptrimmed() 102 …ree = (pt->formDegree == 0 || PetscAbsInt(pt->formDegree) == sp->Nv) ? sp->degree : sp->degree + 1; in PetscSpaceSetUp_Ptrimmed() 123 PetscCall(PetscDTPTrimmedSize(sp->Nv, f == 0 ? sp->degree : sp->degree + 1, pt->formDegree, dim)); in PetscSpaceGetDimension_Ptrimmed() 138 PetscInt jet, degree, Nf, Ncopies, Njet; in PetscSpaceEvaluate_Ptrimmed() local 159 degree = f == 0 ? sp->degree : sp->degree + 1; in PetscSpaceEvaluate_Ptrimmed() 164 PetscCall(PetscDTPTrimmedSize(dim, degree, f, &Nb)); in PetscSpaceEvaluate_Ptrimmed() 166 PetscCall(PetscDTPTrimmedEvalJet(dim, npoints, points, degree, f, jet, eval)); in PetscSpaceEvaluate_Ptrimmed() [all …]
|
| /petsc/src/dm/impls/plex/tests/ |
| H A D | ex75.c | 8 PetscInt dim = 1, d, cStart, cEnd, c, q, degree = 2, coordSize, offset; in main() local 39 PetscCall(PetscSectionSetDof(coordSec, c, dim * (degree + 1))); in main() 40 PetscCall(PetscSectionSetFieldDof(coordSec, c, 0, dim * (degree + 1))); in main() 52 for (q = 0; q < degree + 1; ++q) { in main() 54 …oords[offset + dim * q + d] = 100. + (PetscScalar)c + (1.0 / (PetscScalar)degree) * (PetscScalar)q; in main()
|
| /petsc/src/vec/is/sf/tutorials/ |
| H A D | ex1.c | 388 const PetscInt *degree; in main() local 389 PetscCall(PetscSFComputeDegreeBegin(sf, °ree)); in main() 390 PetscCall(PetscSFComputeDegreeEnd(sf, °ree)); in main() 392 PetscCall(PetscIntView(nrootsalloc, degree, PETSC_VIEWER_STDOUT_WORLD)); in main() 413 const PetscInt *degree; in main() local 415 PetscCall(PetscSFComputeDegreeBegin(sf, °ree)); in main() 416 PetscCall(PetscSFComputeDegreeEnd(sf, °ree)); in main() 417 for (i = 0, inedges = 0; i < nrootsalloc; i++) inedges += degree[i]; in main() 429 const PetscInt *degree; in main() local 431 PetscCall(PetscSFComputeDegreeBegin(sf, °ree)); in main() [all …]
|
| /petsc/src/dm/dt/space/impls/poly/ |
| H A D | spacepoly.c | 19 …ce of degree %" PetscInt_FMT "\n", poly->tensor ? "Tensor polynomial" : "Polynomial", sp->degree)); in PetscSpacePolynomialView_Ascii() 63 PetscInt degree = sp->degree; in PetscSpaceSetUp_Polynomial() local 82 PetscCall(PetscSpaceSetDegree(subsp, degree, PETSC_DETERMINE)); in PetscSpaceSetUp_Polynomial() 98 …scCheck(sp->degree >= 0, PetscObjectComm((PetscObject)sp), PETSC_ERR_ARG_OUTOFRANGE, "Negative deg… in PetscSpaceSetUp_Polynomial() 99 sp->maxDegree = sp->degree; in PetscSpaceSetUp_Polynomial() 106 PetscInt deg = sp->degree; in PetscSpaceGetDimension_Polynomial() 150 PetscCall(PetscDTBinomialInt(dim + sp->degree, dim, &Nb)); in PetscSpaceEvaluate_Polynomial() 163 if (sp->degree == 1) { in PetscSpaceEvaluate_Polynomial() 166 PetscCall(PetscDTPKDEvalJet(dim, npoints, points, sp->degree, jet, pScalar)); in PetscSpaceEvaluate_Polynomial()
|
| /petsc/src/dm/dt/space/interface/ |
| H A D | space.c | 208 …(maximally included) polynomial degree", "PetscSpaceSetDegree", sp->degree, &sp->degree, NULL, 0)); in PetscSpaceSetFromOptions() 294 s->degree = 0; in PetscSpaceCreate() 349 if (minDegree) *minDegree = sp->degree; in PetscSpaceGetDegree() 366 PetscErrorCode PetscSpaceSetDegree(PetscSpace sp, PetscInt degree, PetscInt maxDegree) in PetscSpaceSetDegree() argument 370 sp->degree = degree; in PetscSpaceSetDegree()
|