Home
last modified time | relevance | path

Searched +full:- +full:d (Results 1 – 25 of 1044) sorted by relevance

12345678910>>...42

/petsc/src/dm/dt/tests/output/
H A Dex1_1.out1 User-provided points
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
[all …]
H A Dex4_0.out1 IS Object: 1D, Simplex, Order 3, Point 0 Symmetry -1 1 MPI process
6 Vec Object: 1D, Simplex, Order 3, Point 0 Symmetry -1 1 MPI process
9 -0.447214
10 IS Object: 1D, Simplex, Order 3, Point 0 Symmetry 0 1 MPI process
15 Vec Object: 1D, Simplex, Order 3, Point 0 Symmetry 0 1 MPI process
17 -0.447214
19 IS Object: 1D, Simplex, Order 4, Point 0 Symmetry -1 1 MPI process
25 Vec Object: 1D, Simplex, Order 4, Point 0 Symmetry -1 1 MPI process
29 -0.654654
30 IS Object: 1D, Simplex, Order 4, Point 0 Symmetry 0 1 MPI process
[all …]
/petsc/share/petsc/datafiles/meshes/
H A Dgmsh-3d-binary.msh239 …���?�������?A������?+����������?�������?������?,����������?�������?�������?-�����������������?�…
43 �����?P��<���?C����������?�f�x���?�f�x���?D����������?�MUUU��?�������?E����������?�������?�MUUU��?…
48 �����?P��<���?��������d���^������?^������?��������e���A������?������?��������f���������?������?…
59 …���+������(������\���������������,������(������\���*������������-������(������]����…
61 …���B������(���&���]���T������������C������(���*���Y���+������������D������(���X���Y���*�…
62 ���g������������Z������2���d��� ���g������������[������2��� ���d���j������������\������…
65 …d������2������b���m������������e������2���������h������������f������2������d���h����…
H A Dnozzle.igs5 1.000,6,1HM,1,1.000,15H20201228.095057,1.0e-008,0.00,6HNoname,6HNoname, G 3
7 314 1 0 1 0 0 0 000010001D 1
8 314 1 0 1 0 0 0 0D 2
9 116 2 0 1 0 0 0 000010001D 3
10 116 1 0 1 0 0 0 0D 4
11 123 3 0 1 0 0 0 000010201D 5
12 123 1 0 1 0 0 0 0D 6
13 123 4 0 1 0 0 0 000010201D 7
14 123 1 0 1 0 0 0 0D 8
15 190 5 0 1 0 0 0 000010001D 9
[all …]
H A Dsquare.cas13 2.5000000000e-01 -5.0000000000e-01
14 3.7500000000e-01 -5.0000000000e-01
15 3.1538547488e-01 -3.9957278878e-01
16 3.7500000000e-01 5.0000000000e-01
17 4.1801295515e-01 4.1820081129e-01
18 5.0000000000e-01 5.0000000000e-01
19 1.2500000000e-01 -5.0000000000e-01
20 1.8999095390e-01 -3.8996719959e-01
21 1.5256269259e-01 -6.0843155578e-02
22 2.7596690063e-01 1.4377567833e-03
[all …]
H A Dwedge_cylinder.cas1 (0 "fluent24.2.0 build-id: 10192")
26 4 2b 1 2d 2c 1 8
27 4 1 2 2e 2d 2 9
32 4 5 6 32 31 6 d
33 4 1 2d 37 7 8 9
44 4 5 31 3b b c d
45 4 5 6 c b d 37
47 4 6 32 3c c d e
51 4 35 7 d 3f f 3a
52 4 b c 3c 3b d 14
[all …]
/petsc/src/tao/leastsquares/tutorials/matlab/more_wild_probs/
H A Ddfovec.m5 % Benchmarking Derivative-Free Optimization Algorithms
7 % SIAM J. Optimization, Vol. 20 (1), pp.172-191, 2009.
45 v = [4.0d0,2.0d0,1.0d0,5.0d-1,2.5d-1,1.67d-1,1.25d-1,1.0d-1,8.33d-2,...
46 7.14d-2,6.25d-2];
47 y1 = [1.4d-1,1.8d-1,2.2d-1,2.5d-1,2.9d-1,3.2d-1,3.5d-1,3.9d-1,3.7d-1,...
48 5.8d-1,7.3d-1,9.6d-1,1.34d0,2.1d0,4.39d0];
49 y2 = [1.957d-1,1.947d-1,1.735d-1,1.6d-1,8.44d-2,6.27d-2,4.56d-2,3.42d-2,...
50 3.23d-2,2.35d-2,2.46d-2];
53 y4 = [8.44d-1,9.08d-1,9.32d-1,9.36d-1,9.25d-1,9.08d-1,8.81d-1,8.5d-1,...
54 8.18d-1,7.84d-1,7.51d-1,7.18d-1,6.85d-1,6.58d-1,6.28d-1,6.03d-1,...
[all …]
/petsc/src/dm/interface/
H A Ddmperiodicity.c6 DMGetPeriodicity - Get the description of mesh periodicity
11 . dm - The `DM` object
14 + maxCell - Over distances greater than this, we can assume a point has crossed over to another she…
15 . Lstart - If we assume the mesh is a torus, this is the start of each coordinate, or `NULL` for 0…
16 - L - If we assume the mesh is a torus, this is the length of each coordinate, otherwise it i…
26 if (maxCell) *maxCell = dm->maxCell; in DMGetPeriodicity()
27 if (Lstart) *Lstart = dm->Lstart; in DMGetPeriodicity()
28 if (L) *L = dm->L; in DMGetPeriodicity()
33 DMSetPeriodicity - Set the description of mesh periodicity
38 + dm - The `DM` object
[all …]
H A Ddmgeommodel.c19 for (PetscInt d = 0; d < dE; ++d) norm += PetscSqr(PetscRealPart(mcoords[d])); in DMSnapToGeomModelBall() local
21 for (PetscInt d = 0; d < dE; ++d) gcoords[d] = mcoords[d] / norm; in DMSnapToGeomModelBall() local
23 for (PetscInt d = 0; d < dE; ++d) gcoords[d] = mcoords[d]; in DMSnapToGeomModelBall() local
39 for (PetscInt d = 0; d < dE - 1; ++d) norm += PetscSqr(PetscRealPart(mcoords[d])); in DMSnapToGeomModelCylinder() local
41 for (PetscInt d = 0; d < dE - 1; ++d) gcoords[d] = mcoords[d] * gmax[0] / norm; in DMSnapToGeomModelCylinder() local
42 gcoords[dE - 1] = mcoords[dE - 1]; in DMSnapToGeomModelCylinder()
44 for (PetscInt d = 0; d < dE; ++d) gcoords[d] = mcoords[d]; in DMSnapToGeomModelCylinder() local
50 DMGeomModelRegisterAll - Registers all of the geometry model methods in the `DM` package.
73 DMGeomModelRegister - Adds a geometry model to `DM`
78 + sname - name of a new user-defined geometry model
[all …]
/petsc/src/dm/dt/space/impls/wxy/
H A Dspacewxy.c14 PetscCall(PetscViewerASCIIPrintf(v, "WXY space of degree %" PetscInt_FMT "\n", sp->degree)); in PetscSpacePolynomialView_Ascii()
32 PetscSpace_WXY *wxy = (PetscSpace_WXY *)sp->data; in PetscSpaceDestroy_WXY()
41 PetscSpace_WXY *wxy = (PetscSpace_WXY *)sp->data; in PetscSpaceSetUp_WXY()
44 if (wxy->setupcalled) PetscFunctionReturn(PETSC_SUCCESS); in PetscSpaceSetUp_WXY()
45 …PetscCheck(sp->degree >= 0, PetscObjectComm((PetscObject)sp), PETSC_ERR_ARG_OUTOFRANGE, "Negative … in PetscSpaceSetUp_WXY()
46 sp->maxDegree = sp->degree; in PetscSpaceSetUp_WXY()
47 wxy->setupcalled = PETSC_TRUE; in PetscSpaceSetUp_WXY()
58 …cSpace sp, PetscInt npoints, const PetscReal points[], PetscReal B[], PetscReal D[], PetscReal H[]) in PetscSpaceEvaluate_WXY() argument
60 PetscSpace_WXY *wxy = (PetscSpace_WXY *)sp->data; in PetscSpaceEvaluate_WXY()
61 PetscInt dim = sp->Nv; in PetscSpaceEvaluate_WXY()
[all …]
/petsc/src/dm/impls/swarm/
H A Ddata_ex.c6 This is a doubly-connected edge list which enumerates the
12 proc ---------> proc
13 0 <-------- 1
33 |--------|-----------------|--|
44 Phase D: Send data
63 +Phase D
87 DMSwarmDataEx d; in DMSwarmDataExCreate() local
90 PetscCall(PetscNew(&d)); in DMSwarmDataExCreate()
91 PetscCallMPI(MPI_Comm_dup(comm, &d->comm)); in DMSwarmDataExCreate()
92 PetscCallMPI(MPI_Comm_rank(d->comm, &d->rank)); in DMSwarmDataExCreate()
[all …]
/petsc/src/mat/impls/sbaij/seq/
H A Dsbaijfact12.c9 Mat_SeqSBAIJ *a = (Mat_SeqSBAIJ *)A->data, *b = (Mat_SeqSBAIJ *)C->data; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
10 PetscInt i, j, mbs = a->mbs, *bi = b->i, *bj = b->j; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
12 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
13 MatScalar *u, *d, *w, *wp, u0, u1, u2, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering() local
17 PetscReal shift = info->shiftamount; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
22 allowzeropivot = PetscNot(A->erroriffailure); in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
29 ai = a->i; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
30 aj = a->j; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
31 aa = a->a; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
35 /*initialize k-th row with elements nonzero in row k of A */ in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
[all …]
H A Dsbaijfact11.c7 Mat_SeqSBAIJ *a = (Mat_SeqSBAIJ *)A->data, *b = (Mat_SeqSBAIJ *)C->data; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
8 IS perm = b->row; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
9 const PetscInt *ai, *aj, *perm_ptr, *bi = b->i, *bj = b->j; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
10 PetscInt *a2anew, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili, i, j, mbs = a->mbs; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
11 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
12 MatScalar *u, *d, *w, *wp, u0, u1, u2, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12; in MatCholeskyFactorNumeric_SeqSBAIJ_7() local
16 PetscReal shift = info->shiftamount; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
21 allowzeropivot = PetscNot(A->erroriffailure); in MatCholeskyFactorNumeric_SeqSBAIJ_7()
31 if (!a->permute) { in MatCholeskyFactorNumeric_SeqSBAIJ_7()
32 ai = a->i; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
[all …]
H A Dsbaijfact9.c7 Mat_SeqSBAIJ *a = (Mat_SeqSBAIJ *)A->data, *b = (Mat_SeqSBAIJ *)C->data; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
8 IS perm = b->row; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
9 const PetscInt *ai, *aj, *perm_ptr, mbs = a->mbs, *bi = b->i, *bj = b->j; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
11 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
12 MatScalar *u, *d, *w, *wp, u0, u1, u2, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12; in MatCholeskyFactorNumeric_SeqSBAIJ_6() local
15 PetscReal shift = info->shiftamount; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
20 allowzeropivot = PetscNot(A->erroriffailure); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
30 if (!a->permute) { in MatCholeskyFactorNumeric_SeqSBAIJ_6()
31 ai = a->i; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
32 aj = a->j; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
[all …]
H A Dsbaijfact8.c9 Mat_SeqSBAIJ *a = (Mat_SeqSBAIJ *)A->data, *b = (Mat_SeqSBAIJ *)C->data; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
10 PetscInt i, j, mbs = a->mbs, *bi = b->i, *bj = b->j; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
12 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
13 MatScalar *u, *d, *rtmp, *rtmp_ptr, work[25]; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering() local
14 PetscReal shift = info->shiftamount; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
19 allowzeropivot = PetscNot(A->erroriffailure); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
26 ai = a->i; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
27 aj = a->j; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
28 aa = a->a; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
32 /*initialize k-th row with elements nonzero in row k of A */ in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
[all …]
H A Dsbaijfact7.c7 Mat_SeqSBAIJ *a = (Mat_SeqSBAIJ *)A->data, *b = (Mat_SeqSBAIJ *)C->data; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
8 IS perm = b->row; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
9 const PetscInt *ai, *aj, *perm_ptr, mbs = a->mbs, *bi = b->i, *bj = b->j; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
11 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
12 MatScalar *u, *d, *rtmp, *rtmp_ptr, work[25]; in MatCholeskyFactorNumeric_SeqSBAIJ_5() local
13 PetscReal shift = info->shiftamount; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
18 allowzeropivot = PetscNot(A->erroriffailure); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
28 if (!a->permute) { in MatCholeskyFactorNumeric_SeqSBAIJ_5()
29 ai = a->i; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
30 aj = a->j; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
[all …]
/petsc/src/dm/impls/stag/
H A Dstag2d.c1 /* Functions specific to the 2-dimensional implementation of DMStag */
5 …DMStagCreate2d - Create an object to manage data living on the elements, faces, and vertices of a …
10 + comm - MPI communicator
11 . bndx - x boundary type, `DM_BOUNDARY_NONE`, `DM_BOUNDARY_PERIODIC`, or
13 . bndy - y boundary type, `DM_BOUNDARY_NONE`, `DM_BOUNDARY_PERIODIC`, or `DM_BOUNDARY_GHOST…
14 . M - global number of elements in x direction
15 . N - global number of elements in y direction
16 . m - number of ranks in the x direction (may be `PETSC_DECIDE`)
17 . n - number of ranks in the y direction (may be `PETSC_DECIDE`)
18 . dof0 - number of degrees of freedom per vertex/0-cell
[all …]
H A Dstag.c2 Implementation of DMStag, defining dimension-independent functions in the
3 DM API. stag1d.c, stag2d.c, and stag3d.c may include dimension-specific
5 DMStag-specific API functions, as well as internal functions.
12 PetscInt f0, f1, f2, f3, dof0, dof1, dof2, dof3, n_entries, k, d, cnt, n_fields, dim; in DMCreateFieldDecomposition_Stag() local
38 for (d = 0; d < dof0; ++d) { in DMCreateFieldDecomposition_Stag()
39 stencil0[dof0 * k + d].i = 0; in DMCreateFieldDecomposition_Stag()
40 stencil0[dof0 * k + d].j = 0; in DMCreateFieldDecomposition_Stag()
41 stencil0[dof0 * k + d].j = 0; in DMCreateFieldDecomposition_Stag()
45 for (d = 0; d < dof1; ++d) { in DMCreateFieldDecomposition_Stag()
46 stencil1[dof1 * k + d].i = 0; in DMCreateFieldDecomposition_Stag()
[all …]
H A Dstag3d.c1 /* Functions specific to the 3-dimensional implementation of DMStag */
5 …DMStagCreate3d - Create an object to manage data living on the elements, faces, edges, and vertice…
10 + comm - MPI communicator
11 . bndx - x boundary type, `DM_BOUNDARY_NONE`, `DM_BOUNDARY_PERIODIC`, or `DM_BOUNDARY_GHOST…
12 . bndy - y boundary type, `DM_BOUNDARY_NONE`, `DM_BOUNDARY_PERIODIC`, or `DM_BOUNDARY_GHOST…
13 . bndz - z boundary type, `DM_BOUNDARY_NONE`, `DM_BOUNDARY_PERIODIC`, or `DM_BOUNDARY_GHOST…
14 . M - global number of elements in x direction
15 . N - global number of elements in y direction
16 . P - global number of elements in z direction
17 . m - number of ranks in the x direction (may be `PETSC_DECIDE`)
[all …]
/petsc/src/dm/impls/plex/tests/
H A Dex8.c29 PetscCall(DMViewFromOptions(*dm, NULL, "-dm_view")); in ReadMesh()
39 options->runType = RUN_REFERENCE; in ProcessOptions()
40 options->transform = PETSC_FALSE; in ProcessOptions()
43 run = options->runType; in ProcessOptions()
44 …PetscCall(PetscOptionsEList("-run_type", "The run type", "ex8.c", runTypes, 3, runTypes[options->r… in ProcessOptions()
45 options->runType = (RunType)run; in ProcessOptions()
46 …PetscCall(PetscOptionsBool("-transform", "Use random transforms", "ex8.c", options->transform, &op… in ProcessOptions()
48 if (options->runType == RUN_FILE) { in ProcessOptions()
52 PetscCall(ReadMesh(PETSC_COMM_WORLD, options, &options->dm)); in ProcessOptions()
53 PetscCall(DMGetDimension(options->dm, &dim)); in ProcessOptions()
[all …]
/petsc/src/dm/impls/plex/tutorials/
H A Dex20.c24 …PetscCall(PetscOptionsFList("-dm_type", "DM implementation on which to define field", "ex20.c", DM… in main()
27 PetscCall(PetscOptionsGetBool(NULL, NULL, "-ex20_print", &print, NULL)); in main()
28 PetscCall(PetscOptionsGetInt(NULL, NULL, "-dm_plex_dim", &dim, NULL)); in main()
31 PetscCall(PetscOptionsGetIntArray(NULL, NULL, "-dm_plex_box_faces", faces, &n, NULL)); in main()
34 PetscCall(PetscOptionsGetRealArray(NULL, NULL, "-dm_plex_box_lower", lower, &n, NULL)); in main()
37 PetscCall(PetscOptionsGetRealArray(NULL, NULL, "-dm_plex_box_upper", upper, &n, NULL)); in main()
59 for (PetscInt d = 0; d < dim; d++) { in main() local
66 PetscCall(DMCreateGlobalVector(dm, &Vecxyz[d])); in main()
67 …PetscCall(PetscObjectSetName((PetscObject)Vecxyz[d], d == 0 ? "x_coordinate" : (d == 1 ? "y_coordi… in main()
81 xArray[sOff] = coordArray[cOff + d]; in main()
[all …]
/petsc/src/ts/tutorials/advection-diffusion-reaction/
H A Dshashi.F908 ! ------------------------------------------------------------------------
10 ! FormFunction - Evaluates nonlinear function, F(x).
13 ! snes - the SNES context
14 ! x - input vector
15 ! dummy - optional user-defined context (not used here)
18 ! f - function vector
31 ! - For default PETSc vectors, VecGetArray() returns a pointer to
33 ! - You MUST call VecRestoreArray() when you no longer need access to
35 ! - Note that the Fortran interface to VecGetArray() differs from the
46 ! ---------------------------------------------------------------------
[all …]
/petsc/src/dm/impls/stag/tests/
H A Dex7.c1 static char help[] = "Test DMStag 3d periodic and ghosted boundary conditions\n\n";
10 …PetscInt startx, starty, startz, nx, ny, nz, i, j, k, d, is, js, ks, dof0, dof1, dof2, dof3,… in main() local
44 for (d = 0; d < dofTotal; ++d) { in main()
45 …][j][i][d] != 1.0) PetscCall(PetscPrintf(PETSC_COMM_SELF, "[%d] Unexpected value %g (expecting %g)… in main()
46 a2[k][j][i][d] = 0.0; in main()
47 for (ks = -stencilWidth; ks <= stencilWidth; ++ks) { in main()
48 for (js = -stencilWidth; js <= stencilWidth; ++js) { in main()
49 …for (is = -stencilWidth; is <= stencilWidth; ++is) a2[k][j][i][d] += a1[k + ks][j + js][i + is][d]; in main()
62 /* For the all-periodic case, all values are the same . Otherwise, just check the local version */ in main()
67 for (d = 0; d < 3; ++d) expected *= (2 * stencilWidth + 1); in main()
[all …]
/petsc/src/mat/impls/aij/seq/
H A Dmatmatmatmult.c2 Defines matrix-matrix-matrix product routines for SeqAIJ matrices
3 D = A * B * C
12 PetscCall(MatDestroy(&matmatmatmult->BC)); in MatProductCtxDestroy_SeqAIJ_MatMatMatMult()
17 …tscErrorCode MatMatMatMultSymbolic_SeqAIJ_SeqAIJ_SeqAIJ(Mat A, Mat B, Mat C, PetscReal fill, Mat D) in MatMatMatMultSymbolic_SeqAIJ_SeqAIJ_SeqAIJ() argument
24 MatCheckProduct(D, 5); in MatMatMatMultSymbolic_SeqAIJ_SeqAIJ_SeqAIJ()
25 …PetscCheck(!D->product->data, PetscObjectComm((PetscObject)D), PETSC_ERR_PLIB, "Product data not e… in MatMatMatMultSymbolic_SeqAIJ_SeqAIJ_SeqAIJ()
29 PetscCall(PetscStrallocpy(D->product->alg, &alg)); in MatMatMatMultSymbolic_SeqAIJ_SeqAIJ_SeqAIJ()
30 PetscCall(MatProductSetAlgorithm(D, "sorted")); /* set alg for D = A*BC */ in MatMatMatMultSymbolic_SeqAIJ_SeqAIJ_SeqAIJ()
31 PetscCall(MatMatMultSymbolic_SeqAIJ_SeqAIJ(A, BC, fill, D)); in MatMatMatMultSymbolic_SeqAIJ_SeqAIJ_SeqAIJ()
32 PetscCall(MatProductSetAlgorithm(D, alg)); /* resume original algorithm */ in MatMatMatMultSymbolic_SeqAIJ_SeqAIJ_SeqAIJ()
[all …]
/petsc/src/dm/dt/tests/
H A Dex8.c9 PetscInt d, n, maxdim = 4; in main() local
16 for (d = 0; d <= maxdim; d++) { in main()
17 for (n = 0; n <= d + 2; n++) { in main()
20 PetscCall(PetscDTBinomialInt(d + n, d, &Nk)); in main()
24 PetscCall(PetscDTIndexToBary(d + 1, n, k, btup)); in main()
25 for (j = 0, sum = 0; j < d + 1; j++) { in main()
26 …C_ERR_PLIB, "PetscDTIndexToBary, d = %" PetscInt_FMT ", n = %" PetscInt_FMT ", k = %" PetscInt_FMT… in main()
29 …SC_ERR_PLIB, "PetscDTIndexToBary, d = %" PetscInt_FMT ", n = %" PetscInt_FMT ", k = %" PetscInt_FM… in main()
30 PetscCall(PetscDTBaryToIndex(d + 1, n, btup, &kchk)); in main()
31 …ETSC_ERR_PLIB, "PetscDTBaryToIndex, d = %" PetscInt_FMT ", n = %" PetscInt_FMT ", k = %" PetscInt_… in main()
[all …]

12345678910>>...42