Home
last modified time | relevance | path

Searched refs:D (Results 1 – 25 of 417) sorted by relevance

12345678910>>...17

/petsc/share/petsc/datafiles/meshes/
H A Dnozzle.igs7 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
16 190 1 0 1 1 0 0 0D 10
[all …]
/petsc/src/dm/dt/space/impls/wxy/
H A Dspacewxy.c58 …cSpace sp, PetscInt npoints, const PetscReal points[], PetscReal B[], PetscReal D[], PetscReal H[]) in PetscSpaceEvaluate_WXY() argument
68 PetscCall(PetscSpaceEvaluate(sp, npoints, points, B, D, H)); in PetscSpaceEvaluate_WXY()
108 if (D) { in PetscSpaceEvaluate_WXY()
119 D[p * p_inc + 0 * b_inc + 0 * c_inc + 0] = 0.; in PetscSpaceEvaluate_WXY()
120 D[p * p_inc + 0 * b_inc + 0 * c_inc + 1] = 2. * z; in PetscSpaceEvaluate_WXY()
121 D[p * p_inc + 0 * b_inc + 0 * c_inc + 2] = 2. * y; in PetscSpaceEvaluate_WXY()
122 D[p * p_inc + 0 * b_inc + 1 * c_inc + 0] = 0.; in PetscSpaceEvaluate_WXY()
123 D[p * p_inc + 0 * b_inc + 1 * c_inc + 1] = 0.; in PetscSpaceEvaluate_WXY()
124 D[p * p_inc + 0 * b_inc + 1 * c_inc + 2] = 0.; in PetscSpaceEvaluate_WXY()
125 D[p * p_inc + 0 * b_inc + 2 * c_inc + 0] = 0.; in PetscSpaceEvaluate_WXY()
[all …]
/petsc/src/dm/dt/tests/output/
H A Dex1_1.out2 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 …]
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
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
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
30 IS Object: 1D, Simplex, Order 4, Point 0 Symmetry 0 1 MPI process
36 Vec Object: 1D, Simplex, Order 4, Point 0 Symmetry 0 1 MPI process
41 IS Object: 1D, Simplex, Order 5, Point 0 Symmetry -1 1 MPI process
48 Vec Object: 1D, Simplex, Order 5, Point 0 Symmetry -1 1 MPI process
[all …]
/petsc/src/mat/impls/aij/seq/
H A Dmatmatmatmult.c17 …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()
36 PetscCheck(!D->product->data, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Not yet coded"); in MatMatMatMultSymbolic_SeqAIJ_SeqAIJ_SeqAIJ()
39 D->product->data = matmatmatmult; in MatMatMatMultSymbolic_SeqAIJ_SeqAIJ_SeqAIJ()
40 D->product->destroy = MatProductCtxDestroy_SeqAIJ_MatMatMatMult; in MatMatMatMultSymbolic_SeqAIJ_SeqAIJ_SeqAIJ()
[all …]
/petsc/src/mat/tests/
H A Dex216.c7 Mat A, B, C, D, E; in main() local
17 PetscCall(MatCreateDense(PETSC_COMM_WORLD, PETSC_DECIDE, PETSC_DECIDE, M, M, NULL, &D)); in main()
18 PetscCall(MatSetRandom(D, NULL)); in main()
19 PetscCall(MatTranspose(D, MAT_INITIAL_MATRIX, &A)); in main()
20 PetscCall(MatAXPY(D, 1.0, A, SAME_NONZERO_PATTERN)); in main()
22 PetscCall(MatSetBlockSize(D, bs)); in main()
23 PetscCall(MatDuplicate(D, MAT_COPY_VALUES, &E)); in main()
24 PetscCall(MatViewFromOptions(D, NULL, "-input_dense")); in main()
26 PetscCall(MatConvert(D, MATAIJ, MAT_INITIAL_MATRIX, &A)); in main()
27 PetscCall(MatConvert(D, MATBAIJ, MAT_INITIAL_MATRIX, &B)); in main()
[all …]
H A Dex109.c7 Mat A, B, C, D, AT; in main() local
88 PetscCall(MatProductCreate(A, B, NULL, &D)); in main()
89 PetscCall(MatProductSetType(D, MATPRODUCT_AB)); in main()
90 PetscCall(MatProductSetAlgorithm(D, "default")); in main()
91 PetscCall(MatProductSetFill(D, fill)); in main()
92 PetscCall(MatProductSetFromOptions(D)); in main()
93 PetscCall(MatProductSymbolic(D)); in main()
94 for (i = 0; i < 2; i++) PetscCall(MatProductNumeric(D)); in main()
95 PetscCall(MatEqual(C, D, &equal)); in main()
97 PetscCall(MatDestroy(&D)); in main()
[all …]
H A Dex66.c32 Mat A, B, C, D; in main() local
203 PetscCall(MatDuplicate(B, MAT_COPY_VALUES, &D)); in main()
204 PetscCall(MatSetOption(D, MAT_SYMMETRIC, symm)); in main()
205 PetscCall(MatMultEqual(B, D, 10, &flg)); in main()
208 PetscCall(MatMultTransposeEqual(B, D, 10, &flg)); in main()
211 PetscCall(MatDestroy(&D)); in main()
242 PetscCall(MatMatMult(B, C, nt ? MAT_REUSE_MATRIX : MAT_INITIAL_MATRIX, PETSC_DETERMINE, &D)); in main()
243 PetscCall(MatViewFromOptions(D, NULL, "-bc_view")); in main()
244 PetscCall(PetscObjectBaseTypeCompareAny((PetscObject)D, &flg, MATSEQDENSE, MATMPIDENSE, "")); in main()
247 PetscCall(MatCreateVecs(D, NULL, &v)); in main()
[all …]
H A Dex104.c11 Mat A, B, C, D; in main() local
103 PetscCall(MatDuplicate(C, MAT_COPY_VALUES, &D)); in main()
105 PetscCall(MatDestroy(&D)); in main()
114 PetscCall(MatTransposeMatMult(A, A, MAT_INITIAL_MATRIX, fill, &D)); /* D = A^T*A */ in main()
115 PetscCall(MatTransposeMatMult(A, A, MAT_REUSE_MATRIX, fill, &D)); in main()
116 PetscCall(MatTransposeMatMultEqual(A, A, D, 10, &equal)); in main()
120 PetscCall(MatDuplicate(D, MAT_COPY_VALUES, &C)); in main()
124 PetscCall(MatSetRandom(D, NULL)); in main()
125 PetscCall(MatMatMult(A, D, MAT_INITIAL_MATRIX, fill, &E)); in main()
126 PetscCall(MatMatMult(A, D, MAT_REUSE_MATRIX, fill, &E)); in main()
[all …]
H A Dex237.c31 Mat A, C, D, E; in main() local
212 … PetscCall(MatCreateDense(PETSC_COMM_WORLD, bs[j] * m, PETSC_DECIDE, bs[j] * M, N[k], NULL, &D)); in main()
216 PetscCall(MatConvert(D, MATDENSECUDA, MAT_INPLACE_MATRIX, &D)); in main()
244 PetscCall(MatProductCreateWithMat(A, C, NULL, D)); in main()
245 PetscCall(MatProductSetType(D, MATPRODUCT_AB)); in main()
246 PetscCall(MatProductSetFromOptions(D)); in main()
247 PetscCall(MatProductSymbolic(D)); in main()
252 PetscCall(MatProductNumeric(D)); in main()
255 for (t = 0; t < trial; ++t) PetscCall(MatProductNumeric(D)); in main()
264 PetscCall(MatDenseGetLocalMatrix(D, &Dt)); in main()
[all …]
/petsc/include/petsc/private/cpp/
H A Dregister_finalize.hpp74 template <typename D>
76 inline PetscErrorCode RegisterFinalizeable<D>::do_finalize_(Self &&self, Args &&...args) noexcept in do_finalize_()
88 template <typename D>
90 inline PetscErrorCode RegisterFinalizeable<D>::do_register_finalize_(Self &&self, Args &&...args) n… in do_register_finalize_()
106 template <typename D>
108 inline constexpr PetscErrorCode RegisterFinalizeable<D>::finalize_(Args &&...) noexcept in finalize_()
113 template <typename D>
115 inline constexpr PetscErrorCode RegisterFinalizeable<D>::register_finalize_(Args &&...) noexcept in register_finalize_()
126 template <typename D>
127 inline bool RegisterFinalizeable<D>::registered() const noexcept in registered()
[all …]
H A Darray.hpp18 template <class D, class...>
20 using type = D;
26 template <class D, class... T>
27 using array_return_type = std::array<typename return_type_helper<D, T...>::type, sizeof...(T)>;
37 template <class D = void, class... T>
38 …r detail::array_return_type<D, T...> make_array(T &&...t) noexcept(std::is_nothrow_constructible<d… in make_array()
/petsc/include/petsc/private/
H A Dveccupmimpl.h256 template <device::cupm::DeviceType T, typename D>
258 class Vec_CUPMBase<T, D>::VectorArray : public device::cupm::impl::RestoreableArray<T, MT, MA> {
273 template <device::cupm::DeviceType T, typename D>
275 inline Vec_CUPMBase<T, D>::VectorArray<MT, MA>::VectorArray(PetscDeviceContext dctx, Vec v) noexcep… in VectorArray()
278 …PetscCallAbort(PETSC_COMM_SELF, Vec_CUPMBase<T, D>::template GetArray<MT, MA, true>(v, &this->ptr_…
282 template <device::cupm::DeviceType T, typename D>
284 inline Vec_CUPMBase<T, D>::VectorArray<MT, MA>::~VectorArray() noexcept in ~VectorArray()
287 …PetscCallAbort(PETSC_COMM_SELF, Vec_CUPMBase<T, D>::template RestoreArray<MT, MA>(v_, &this->ptr_,… in ~VectorArray()
295 template <device::cupm::DeviceType T, typename D>
296 inline PetscErrorCode Vec_CUPMBase<T, D>::ResetAllocatedDevicePtr_(PetscDeviceContext dctx, Vec v, … in ResetAllocatedDevicePtr_()
[all …]
/petsc/src/ksp/pc/impls/sor/
H A Dsor.tex11 Consider the matrix problem $ A x = b$, where $A = L + U + D$.
38 x^{1/2} = \omega (L + D)^{-1} b
41 x = (1 - \omega) x^{1/2} + \omega (U + D)^{-1}(b - L x^{1/2}) = x^{1/2} + \omega (U+D)^{-1}(b - A x…
86 x = \omega (L + D)^{-1}b
97 x = (L + D/\omega)^{-1}b
100 …econditioner obtained by apply the two step process $ (L + D/\omega)^{-1} $ and then $ (U + D/\ome…
103 (L + D/\omega)^{-1} A (U + D/\omega)^{-1} y = (L + D/\omega)^{-1} b.
105 Then after this system is solved, $ x = (U + D/\omega)^{-1} y$. If an initial guess that is nonzero…
106 initial guess for $ y$ must be computed via $ y = (U + D/\omega) x$.
108 …(L + D/\omega)^{-1} A (U + D/\omega)^{-1} & = & (L + D/\omega)^{-1} (L + D + U) (U + D/\omega)^{-…
[all …]
/petsc/src/dm/impls/stag/tests/output/
H A Dex42_dim-3.out8 1 DOF per vertex (0D)
9 1 DOF per edge (1D)
10 1 DOF per face (2D)
11 1 DOF per element (3D)
92 1 DOF per vertex (0D)
93 0 DOF per edge (1D)
94 0 DOF per face (2D)
95 0 DOF per element (3D)
249 0 DOF per vertex (0D)
250 1 DOF per edge (1D)
[all …]
H A Dex16_2.out8 2 DOF per vertex (0D)
9 1 DOF per face (1D)
10 1 DOF per element (2D)
22 2 DOF per vertex (0D)
23 1 DOF per face (1D)
24 1 DOF per element (2D)
36 2 DOF per vertex (0D)
37 1 DOF per face (1D)
38 1 DOF per element (2D)
H A Dex42_dim-2.out8 1 DOF per vertex (0D)
9 1 DOF per face (1D)
10 1 DOF per element (2D)
43 1 DOF per vertex (0D)
44 0 DOF per face (1D)
45 0 DOF per element (2D)
83 0 DOF per vertex (0D)
84 1 DOF per face (1D)
85 0 DOF per element (2D)
112 0 DOF per vertex (0D)
[all …]
/petsc/src/sys/objects/device/impls/
H A Dimpldevicebase.hpp93 template <typename D>
94 inline PetscErrorCode DeviceBase<D>::getDevice(PetscDevice device, PetscInt id) noexcept in getDevice()
106 template <typename D>
107 inline PetscErrorCode DeviceBase<D>::configureDevice(PetscDevice device) noexcept in configureDevice()
114 template <typename D>
115 inline PetscErrorCode DeviceBase<D>::viewDevice(PetscDevice device, PetscViewer viewer) noexcept in viewDevice()
122 template <typename D>
123 inline PetscErrorCode DeviceBase<D>::getAttribute(PetscDevice device, PetscDeviceAttribute attr, vo… in getAttribute()
130 template <typename D>
132 inline PetscErrorCode DeviceBase<D>::PetscOptionDevice(F &&OptionsFunction, PetscOptionItems PetscO… in PetscOptionDevice()
[all …]
/petsc/src/dm/impls/stag/tutorials/output/
H A Dex1_periodic_seq.out7 1 DOF per vertex (0D)
8 1 DOF per element (1D)
18 1 DOF per vertex (0D)
19 1 DOF per element (1D)
29 1 DOF per vertex (0D)
30 0 DOF per element (1D)
40 1 DOF per vertex (0D)
41 0 DOF per element (1D)
51 0 DOF per vertex (0D)
52 1 DOF per element (1D)
/petsc/src/mat/impls/aij/mpi/
H A Dmpimatmatmatmult.c49 …tscErrorCode MatMatMatMultSymbolic_MPIAIJ_MPIAIJ_MPIAIJ(Mat A, Mat B, Mat C, PetscReal fill, Mat D) in MatMatMatMultSymbolic_MPIAIJ_MPIAIJ_MPIAIJ() argument
56 MatCheckProduct(D, 5); in MatMatMatMultSymbolic_MPIAIJ_MPIAIJ_MPIAIJ()
57 …PetscCheck(!D->product->data, PetscObjectComm((PetscObject)D), PETSC_ERR_PLIB, "Product data not e… in MatMatMatMultSymbolic_MPIAIJ_MPIAIJ_MPIAIJ()
58 product = D->product; in MatMatMatMultSymbolic_MPIAIJ_MPIAIJ_MPIAIJ()
65 PetscCall(MatMatMultSymbolic_MPIAIJ_MPIAIJ(A, BC, fill, D)); in MatMatMatMultSymbolic_MPIAIJ_MPIAIJ_MPIAIJ()
69 PetscCall(MatMatMultSymbolic_MPIAIJ_MPIAIJ_nonscalable(A, BC, fill, D)); in MatMatMatMultSymbolic_MPIAIJ_MPIAIJ_MPIAIJ()
74 D->ops->matmatmultnumeric = MatMatMatMultNumeric_MPIAIJ_MPIAIJ_MPIAIJ; in MatMatMatMultSymbolic_MPIAIJ_MPIAIJ_MPIAIJ()
78 PetscErrorCode MatMatMatMultNumeric_MPIAIJ_MPIAIJ_MPIAIJ(Mat A, Mat B, Mat C, Mat D) in MatMatMatMultNumeric_MPIAIJ_MPIAIJ_MPIAIJ() argument
84 MatCheckProduct(D, 4); in MatMatMatMultNumeric_MPIAIJ_MPIAIJ_MPIAIJ()
85 …PetscCheck(D->product->data, PetscObjectComm((PetscObject)D), PETSC_ERR_PLIB, "Product data empty"… in MatMatMatMultNumeric_MPIAIJ_MPIAIJ_MPIAIJ()
[all …]
/petsc/src/dm/tutorials/output/
H A Dex19_1.out1 1D -- DM_BOUNDARY_NONE
6 1D -- DM_BOUNDARY_GHOSTED
11 1D -- DM_BOUNDARY_PERIODIC
16 2D -- DM_BOUNDARY_NONE
21 2D -- DM_BOUNDARY_GHOSTED
26 2D -- DM_BOUNDARY_PERIODIC
31 3D -- DM_BOUNDARY_NONE
36 3D -- DM_BOUNDARY_GHOSTED
41 3D -- DM_BOUNDARY_PERIODIC
/petsc/src/binding/petsc4py/test/
H A Dtest_mat_py.py240 self.D = mat.createVecLeft()
243 self.D.destroy()
247 self.D.scale(a)
250 self.D.shift(a)
253 self.D.zeroEntries()
256 y.pointwiseMult(x, self.D)
262 dctx.D = self.D.duplicate()
264 self.D.copy(dctx.D)
269 self.D.copy(vd)
275 self.D.axpy(1, vd)
[all …]
/petsc/src/mat/impls/transpose/
H A Dhtransm.c6 Mat B, D; member
20 PetscCall(PetscObjectQuery((PetscObject)data->D, "MatProductCtx_HT", (PetscObject *)&container)); in MatProductCtxDestroy_HT()
22 PetscCall(PetscObjectCompose((PetscObject)data->D, "MatProductCtx_HT", NULL)); in MatProductCtxDestroy_HT()
27 static PetscErrorCode MatProductNumeric_HT(Mat D) in MatProductNumeric_HT() argument
35 MatCheckProduct(D, 1); in MatProductNumeric_HT()
36 …PetscCheck(D->product->data, PetscObjectComm((PetscObject)D), PETSC_ERR_PLIB, "Product data empty"… in MatProductNumeric_HT()
37 product = D->product; in MatProductNumeric_HT()
38 PetscCall(PetscObjectQuery((PetscObject)D, "MatProductCtx_HT", (PetscObject *)&container)); in MatProductNumeric_HT()
39 …PetscCheck(container, PetscObjectComm((PetscObject)D), PETSC_ERR_PLIB, "MatProductCtx_HT missing"); in MatProductNumeric_HT()
49 PetscCall((*data->numeric)(D)); in MatProductNumeric_HT()
[all …]
H A Dtransm.c317 Mat D; member
328 …PetscCall(PetscObjectQuery((PetscObject)data->D, "MatProductCtx_Transpose", (PetscObject *)&contai… in MatProductCtxDestroy_Transpose()
330 PetscCall(PetscObjectCompose((PetscObject)data->D, "MatProductCtx_Transpose", NULL)); in MatProductCtxDestroy_Transpose()
335 static PetscErrorCode MatProductNumeric_Transpose(Mat D) in MatProductNumeric_Transpose() argument
342 MatCheckProduct(D, 1); in MatProductNumeric_Transpose()
343 …PetscCheck(D->product->data, PetscObjectComm((PetscObject)D), PETSC_ERR_PLIB, "Product data empty"… in MatProductNumeric_Transpose()
344 product = D->product; in MatProductNumeric_Transpose()
345 PetscCall(PetscObjectQuery((PetscObject)D, "MatProductCtx_Transpose", (PetscObject *)&container)); in MatProductNumeric_Transpose()
346 …PetscCheck(container, PetscObjectComm((PetscObject)D), PETSC_ERR_PLIB, "MatProductCtx_Transpose mi… in MatProductNumeric_Transpose()
350 PetscCall((*data->numeric)(D)); in MatProductNumeric_Transpose()
[all …]
/petsc/src/snes/tutorials/
H A Dex46.c13 PetscReal D; /* The diffusion coefficient */ member
42 user.D = 1.0; in main()
43 PetscCall(PetscOptionsReal("-D", "The diffusion coefficient D", __FILE__, user.D, &user.D, NULL)); in main()
125 PetscReal D, K, hx, hy, hxdhy, hydhx; in FormFunctionLocal() local
129 D = user->D; in FormFunctionLocal()
150 …f[j][i] = D * (uxx + uyy) - (K * funcA(x[j][i], user) * PetscSqrtScalar(ux * ux + uy * uy) + funcU… in FormFunctionLocal()
166 PetscScalar D, K, A, v[5], hx, hy, hxdhy, hydhx, ux, uy; in FormJacobianLocal() local
171 D = user->D; in FormJacobianLocal()
205 v[0] = -D * hxdhy; in FormJacobianLocal()
208 v[1] = -D * hydhx; in FormJacobianLocal()
[all …]

12345678910>>...17