Home
last modified time | relevance | path

Searched +full:- +full:b (Results 1 – 25 of 1115) sorted by relevance

12345678910>>...45

/petsc/src/dm/dt/fe/tests/output/
H A Dex2_2.out2 B[0]: 1. 0. 0.
3 B[1]: -1. 0. 0.
4 B[2]: -1. 0. 0.
5 B[3]: -1. 0. 0.
6 B[4]: 0. 1. 0.
7 B[5]: 0. -1. 0.
8 B[6]: 0. -1. 0.
9 B[7]: 0. -1. 0.
10 B[8]: 0. 0. 1.
11 B[9]: 0. 0. -1.
[all …]
H A Dex2_1.out2 B[0]: 1. 0. 0.
3 B[1]: -1. 0. 0.
4 B[2]: -1. 0. 0.
5 B[3]: -1. 0. 0.
6 B[4]: 0. 1. 0.
7 B[5]: 0. -1. 0.
8 B[6]: 0. -1. 0.
9 B[7]: 0. -1. 0.
10 B[8]: 0. 0. 1.
11 B[9]: 0. 0. -1.
[all …]
/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 …]
/petsc/src/snes/tutorials/output/
H A Dex12_p4est_fas_q2_nonconformal_serial.out6 … resid norm 2.499198301910e+02 true resid norm 3.216249248292e+04 ||r(i)||/||b|| 1.000000000000e+00
7 … resid norm 8.614011391108e+01 true resid norm 1.952597363694e+04 ||r(i)||/||b|| 6.071038694314e-01
8 … resid norm 2.745949105361e+01 true resid norm 1.930071384658e+04 ||r(i)||/||b|| 6.001000655293e-01
9 … resid norm 2.559253117247e+01 true resid norm 1.639100195786e+04 ||r(i)||/||b|| 5.096309611751e-01
10 … resid norm 2.282154921037e+01 true resid norm 9.122321919069e+03 ||r(i)||/||b|| 2.836323062932e-01
11 … resid norm 9.015546483100e+00 true resid norm 2.894802732061e+03 ||r(i)||/||b|| 9.000554710113e-02
12 … resid norm 6.311361207223e+00 true resid norm 1.903443566068e+03 ||r(i)||/||b|| 5.918209128471e-02
13 … resid norm 3.316691807429e+00 true resid norm 7.510656749157e+02 ||r(i)||/||b|| 2.335222232277e-02
14 … resid norm 2.729920508905e+00 true resid norm 6.941012731568e+02 ||r(i)||/||b|| 2.158107844177e-02
15 … resid norm 1.834338478600e+00 true resid norm 2.992015148139e+02 ||r(i)||/||b|| 9.302808697827e-03
[all …]
H A Dex12_mg_newton_coarse_2.out4 Number of 0-cells per rank: 4225
5 Number of 1-cells per rank: 12416
6 Number of 2-cells per rank: 8192
12 Number of 0-cells per rank: 1224
13 Number of 1-cells per rank: 3543
14 Number of 2-cells per rank: 2320
19 Number of 0-cells per rank: 408
20 Number of 1-cells per rank: 1154
21 Number of 2-cells per rank: 747
26 Number of 0-cells per rank: 136
[all …]
/petsc/src/mat/impls/baij/seq/
H A Dbaij.h4 #include <../src/mat/impls/baij/seq/ftn-kernels/fsolvebaij.h>
14 MATSEQBAIJ format - Block compressed row storage. The i[] and j[]
43 PETSC_INTERN PetscErrorCode MatSeqBAIJSetPreallocation_SeqBAIJ(Mat B, PetscInt bs, PetscInt nz, con…
266 PetscKernel_A_gets_A_times_B_2: A = A * B with size bs=2
269 + A,B - square bs by bs arrays stored in column major order
270 - W - bs*bs work array
273 . A = A * B
276 static inline PetscErrorCode PetscKernel_A_gets_A_times_B_2(PetscScalar *A, const PetscScalar *B, P… in PetscKernel_A_gets_A_times_B_2() argument
279 A[0] = W[0] * B[0] + W[2] * B[1]; in PetscKernel_A_gets_A_times_B_2()
280 A[1] = W[1] * B[0] + W[3] * B[1]; in PetscKernel_A_gets_A_times_B_2()
[all …]
H A Daijbaij.c5 Mat B; in MatConvert_SeqBAIJ_SeqAIJ() local
6 Mat_SeqAIJ *b; in MatConvert_SeqBAIJ_SeqAIJ() local
8 Mat_SeqBAIJ *a = (Mat_SeqBAIJ *)A->data; in MatConvert_SeqBAIJ_SeqAIJ()
9 PetscInt bs = A->rmap->bs, *ai = a->i, *aj = a->j, n = A->rmap->N / bs, i, j, k; in MatConvert_SeqBAIJ_SeqAIJ()
11 MatScalar *aa = a->a; in MatConvert_SeqBAIJ_SeqAIJ()
15 B = *newmat; in MatConvert_SeqBAIJ_SeqAIJ()
16 for (i = 0; i < n; i++) maxlen = PetscMax(maxlen, ai[i + 1] - ai[i]); in MatConvert_SeqBAIJ_SeqAIJ()
20 maxlen = PetscMax(maxlen, ai[i + 1] - ai[i]); in MatConvert_SeqBAIJ_SeqAIJ()
21 for (j = 0; j < bs; j++) rowlengths[i * bs + j] = bs * (ai[i + 1] - ai[i]); in MatConvert_SeqBAIJ_SeqAIJ()
23 PetscCall(MatCreate(PetscObjectComm((PetscObject)A), &B)); in MatConvert_SeqBAIJ_SeqAIJ()
[all …]
/petsc/lib/petsc/bin/maint/abi-compliance-checker/doc/
H A DChangelog.html1 <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xht…
4 <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
10 margin-top: 1.0em;
11 background-color: #deeef7;
12 font-family: Helvetica, Arial, FreeSans, san-serif;
19 h1 { font-size: 3.8em; color: #211108; margin-bottom: 3px;margin-top:0px;padding-top:0px;}
20 h1 .small { font-size: 0.4em; }
21 h1 a { text-decoration: none }
22 h2 { font-size: 1.5em; color: #211108; }
23 h3 { text-align: center; color: #211108; }
[all …]
/petsc/src/ksp/ksp/utils/lmvm/
H A Dlmvmimpl.c39 … " title = {Representations of quasi-Newton matrices and their use in limited memory methods},"
41 " ISSN = {1436-4646},"
44 " number = {1-3},"
47 … " author = {Byrd, Richard H. and Nocedal, Jorge and Schnabel, Robert B.},"
50 " pages = {129-156}"
53 PETSC_INTERN PetscErrorCode MatReset_LMVM(Mat B, MatLMVMResetMode mode) in MatReset_LMVM() argument
55 Mat_LMVM *lmvm = (Mat_LMVM *)B->data; in MatReset_LMVM()
58 lmvm->k = 0; in MatReset_LMVM()
59 lmvm->prev_set = PETSC_FALSE; in MatReset_LMVM()
60 lmvm->shift = 0.0; in MatReset_LMVM()
[all …]
H A Dlmvmutils.c7 MatLMVMUpdate - Adds (X-Xprev) and (F-Fprev) updates to a `MATLMVM` matrix.
10 + B - A `MATLMVM` matrix
11 . X - Solution vector
12 - F - Function vector
24 …If the sizes of `B` have not been specified (using `MatSetSizes()` or `MatSetLayouts()`) before `M…
25 …called, the row size and layout of `B` will be set to match `F` and the column size and layout of …
30 PetscErrorCode MatLMVMUpdate(Mat B, Vec X, Vec F) in MatLMVMUpdate() argument
36 PetscValidHeaderSpecific(B, MAT_CLASSID, 1); in MatLMVMUpdate()
39 PetscCall(PetscObjectBaseTypeCompare((PetscObject)B, MATLMVM, &same)); in MatLMVMUpdate()
41 /* If B has specified layouts, this will check X and F are compatible; in MatLMVMUpdate()
[all …]
/petsc/src/sys/tests/linter/output/
H A DtestDuplicates.patch1 @@ -4,5 +4,5 @@
4 -PetscErrorCode testDuplicatesWithChanges(PetscInt *a, PetscScalar *b)
5 +static PetscErrorCode testDuplicatesWithChanges(PetscInt *a, PetscScalar *b)
8 @@ -9,5 +9,5 @@
11 - PetscAssertPointer(a, 1);
14 PetscAssertPointer(b, 5);
15 @@ -11,5 +11,5 @@
18 - PetscAssertPointer(b, 5);
19 + PetscAssertPointer(b, 2);
21 PetscAssertPointer(b, 7);
[all …]
/petsc/share/petsc/datafiles/meshes/
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
29 4 3 4 30 2f 4 b
31 4 2b 1 7 35 8 33
38 4 3 2f 39 9 a b
39 4 3 4 a 9 b 35
41 4 4 30 3a a b c
42 4 4 5 b a c 36
44 4 5 31 3b b c d
45 4 5 6 c b d 37
[all …]
/petsc/src/ksp/ksp/utils/lmvm/diagbrdn/
H A Ddiagbrdn.c4 static PetscErrorCode MatSolve_DiagBrdn(Mat B, Vec F, Vec dX) in MatSolve_DiagBrdn() argument
6 Mat_LMVM *lmvm = (Mat_LMVM *)B->data; in MatSolve_DiagBrdn()
9 PetscCall(MatSolve(lmvm->J0, F, dX)); in MatSolve_DiagBrdn()
13 static PetscErrorCode MatMult_DiagBrdn(Mat B, Vec X, Vec Z) in MatMult_DiagBrdn() argument
15 Mat_LMVM *lmvm = (Mat_LMVM *)B->data; in MatMult_DiagBrdn()
18 PetscCall(MatMult(lmvm->J0, X, Z)); in MatMult_DiagBrdn()
22 static PetscErrorCode MatUpdate_DiagBrdn(Mat B, Vec X, Vec F) in MatUpdate_DiagBrdn() argument
24 Mat_LMVM *lmvm = (Mat_LMVM *)B->data; in MatUpdate_DiagBrdn()
27 if (!lmvm->m) PetscFunctionReturn(PETSC_SUCCESS); in MatUpdate_DiagBrdn()
28 if (lmvm->prev_set) { in MatUpdate_DiagBrdn()
[all …]
/petsc/src/ksp/ksp/utils/lmvm/dfp/
H A Ddfp.c5 …B_{k+1} = (I - y_k (y_k^T s_k)^{-1} s_k^T) B_k (I - s_k (y_k^T s_k)^{-1} y_k^T) + y_k (y_k^T s_k)^…
9 a_k = (y_k^T s_k)^{-1} y_k^T x
10 g = B_k (x - a s_k) <--- recursion
11 B_{k+1}x = g + y_k(a - (y_k^T s_k)^{-1} s_k^T g)
13 PETSC_INTERN PetscErrorCode DFPKernel_Recursive(Mat B, MatLMVMMode mode, Vec X, Vec BX) in DFPKernel_Recursive() argument
24 PetscCall(MatLMVMGetRange(B, &oldest, &next)); in DFPKernel_Recursive()
25 PetscCall(MatLMVMGetUpdatedBasis(B, S_t, &S, NULL, NULL)); in DFPKernel_Recursive()
26 PetscCall(MatLMVMGetUpdatedBasis(B, Y_t, &Y, NULL, NULL)); in DFPKernel_Recursive()
27 PetscCall(MatLMVMGetUpdatedProducts(B, Y_t, S_t, LMBLOCK_DIAGONAL, &YtS)); in DFPKernel_Recursive()
31 PetscCall(PetscMalloc1(next - oldest, &alpha)); in DFPKernel_Recursive()
[all …]
/petsc/src/mat/impls/maij/
H A Dmaij.c5 MatMAIJGetAIJ - Get the `MATAIJ` matrix describing the blockwise action of the `MATMAIJ` matrix
10 . A - the `MATMAIJ` matrix
13 . B - the `MATAIJ` matrix
22 PetscErrorCode MatMAIJGetAIJ(Mat A, Mat *B) in MatMAIJGetAIJ() argument
30 Mat_MPIMAIJ *b = (Mat_MPIMAIJ *)A->data; in MatMAIJGetAIJ() local
32 *B = b->A; in MatMAIJGetAIJ()
34 Mat_SeqMAIJ *b = (Mat_SeqMAIJ *)A->data; in MatMAIJGetAIJ() local
36 *B = b->AIJ; in MatMAIJGetAIJ()
38 *B = A; in MatMAIJGetAIJ()
44 …MatMAIJRedimension - Get a new `MATMAIJ` matrix with the same action, but for a different block si…
[all …]
/petsc/src/mat/impls/aij/seq/bas/
H A Dbasfactor.c7 Mat_SeqAIJ *a = (Mat_SeqAIJ *)A->data; in MatICCFactorSymbolic_SeqAIJ_Bas()
8 Mat_SeqSBAIJ *b; in MatICCFactorSymbolic_SeqAIJ_Bas() local
10 PetscInt reallocs = 0, i, *ai = a->i, *aj = a->j, am = A->rmap->n, *ui; in MatICCFactorSymbolic_SeqAIJ_Bas()
14 PetscReal fill = info->fill, levels = info->levels; in MatICCFactorSymbolic_SeqAIJ_Bas()
19 …->rmap->n == A->cmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Must be square matrix, rows %" Pet… in MatICCFactorSymbolic_SeqAIJ_Bas()
30 for (i = 0; i < am; i++) ui[i + 1] = ui[i] + ai[i + 1] - adiag[i]; in MatICCFactorSymbolic_SeqAIJ_Bas()
34 aj = a->j + adiag[i]; in MatICCFactorSymbolic_SeqAIJ_Bas()
35 ncols = ui[i + 1] - ui[i]; in MatICCFactorSymbolic_SeqAIJ_Bas()
62 b = (Mat_SeqSBAIJ *)fact->data; in MatICCFactorSymbolic_SeqAIJ_Bas()
63 PetscCall(PetscMalloc1(ui[am], &b->a)); in MatICCFactorSymbolic_SeqAIJ_Bas()
[all …]
/petsc/src/ksp/ksp/utils/lmvm/bfgs/
H A Dbfgs.c7 …B_{k+1} = B_k + y_k (y_k^T s_k)^{-1} y_k^T - B_k s_k (s_k^T B_k s_k)^{-1} s_k^T B_k + y_k (y_k^T s…
11 B_{k+1} = B_0 + \sum_i B_i y_i (y_i^T s_i)^{-1} y_i^T - s_i (s_i^T B_i s_i)^{-1} s_i^T B_i
15 static PetscErrorCode BFGSKernel_Recursive_Inner(Mat B, MatLMVMMode mode, PetscInt oldest, PetscInt… in BFGSKernel_Recursive_Inner() argument
17 Mat_LMVM *lmvm = (Mat_LMVM *)B->data; in BFGSKernel_Recursive_Inner()
18 Mat_SymBrdn *lsb = (Mat_SymBrdn *)lmvm->ctx; in BFGSKernel_Recursive_Inner()
20 LMBasis BkS = lsb->basis[LMVMModeMap(SYMBROYDEN_BASIS_BKS, mode)]; in BFGSKernel_Recursive_Inner()
22 LMProducts StBkS = lsb->products[LMVMModeMap(SYMBROYDEN_PRODUCTS_STBKS, mode)]; in BFGSKernel_Recursive_Inner()
27 PetscCall(MatLMVMGetUpdatedBasis(B, Y_t, &Y, NULL, NULL)); in BFGSKernel_Recursive_Inner()
28 PetscCall(MatLMVMGetUpdatedProducts(B, LMBASIS_Y, LMBASIS_S, LMBLOCK_DIAGONAL, &YtS)); in BFGSKernel_Recursive_Inner()
29 PetscCall(MatLMVMGetWorkRow(B, &StBkX)); in BFGSKernel_Recursive_Inner()
[all …]
/petsc/src/ksp/ksp/utils/lmvm/brdn/
H A Dbadbrdn.c5 static PetscErrorCode MatSolve_LMVMBadBrdn_Recursive(Mat B, Vec F, Vec dX) in MatSolve_LMVMBadBrdn_Recursive() argument
8 PetscCall(BroydenKernel_Recursive(B, MATLMVM_MODE_DUAL, F, dX)); in MatSolve_LMVMBadBrdn_Recursive()
12 static PetscErrorCode MatSolve_LMVMBadBrdn_CompactDense(Mat B, Vec F, Vec dX) in MatSolve_LMVMBadBrdn_CompactDense() argument
15 PetscCall(BroydenKernel_CompactDense(B, MATLMVM_MODE_DUAL, F, dX)); in MatSolve_LMVMBadBrdn_CompactDense()
19 PETSC_UNUSED static PetscErrorCode MatSolve_LMVMBadBrdn_Dense(Mat B, Vec F, Vec dX) in MatSolve_LMVMBadBrdn_Dense() argument
22 PetscCall(BroydenKernel_Dense(B, MATLMVM_MODE_DUAL, F, dX)); in MatSolve_LMVMBadBrdn_Dense()
26 static PetscErrorCode MatSolveHermitianTranspose_LMVMBadBrdn_Recursive(Mat B, Vec F, Vec dX) in MatSolveHermitianTranspose_LMVMBadBrdn_Recursive() argument
29 PetscCall(BroydenKernelHermitianTranspose_Recursive(B, MATLMVM_MODE_DUAL, F, dX)); in MatSolveHermitianTranspose_LMVMBadBrdn_Recursive()
33 static PetscErrorCode MatSolveHermitianTranspose_LMVMBadBrdn_CompactDense(Mat B, Vec F, Vec dX) in MatSolveHermitianTranspose_LMVMBadBrdn_CompactDense() argument
36 PetscCall(BroydenKernelHermitianTranspose_CompactDense(B, MATLMVM_MODE_DUAL, F, dX)); in MatSolveHermitianTranspose_LMVMBadBrdn_CompactDense()
[all …]
H A Dbrdn.c6 static PetscErrorCode MatSolve_LMVMBrdn_Recursive(Mat B, Vec F, Vec dX) in MatSolve_LMVMBrdn_Recursive() argument
9 PetscCall(BadBroydenKernel_Recursive(B, MATLMVM_MODE_DUAL, F, dX)); in MatSolve_LMVMBrdn_Recursive()
13 static PetscErrorCode MatSolveHermitianTranspose_LMVMBrdn_Recursive(Mat B, Vec F, Vec dX) in MatSolveHermitianTranspose_LMVMBrdn_Recursive() argument
16 PetscCall(BadBroydenKernelHermitianTranspose_Recursive(B, MATLMVM_MODE_DUAL, F, dX)); in MatSolveHermitianTranspose_LMVMBrdn_Recursive()
20 static PetscErrorCode MatSolve_LMVMBrdn_CompactDense(Mat B, Vec F, Vec dX) in MatSolve_LMVMBrdn_CompactDense() argument
23 PetscCall(BadBroydenKernel_CompactDense(B, MATLMVM_MODE_DUAL, F, dX)); in MatSolve_LMVMBrdn_CompactDense()
27 static PetscErrorCode MatSolveHermitianTranspose_LMVMBrdn_CompactDense(Mat B, Vec F, Vec dX) in MatSolveHermitianTranspose_LMVMBrdn_CompactDense() argument
30 PetscCall(BadBroydenKernelHermitianTranspose_CompactDense(B, MATLMVM_MODE_DUAL, F, dX)); in MatSolveHermitianTranspose_LMVMBrdn_CompactDense()
36 B_k = B_{k-1} + (y_{k-1} - B_{k-1} s_{k-1}) (s_{k-1}^H s_{k-1})^{-1} s_{k-1}^H
40 alpha_k = (s_{k-1}^H s_{k-1})^{-1} s_{k-1}^H x
[all …]
/petsc/src/ksp/ksp/utils/lmvm/symbrdn/
H A Dsymbadbrdn.c3 static PetscErrorCode MatMult_LMVMSymBadBrdn_Recursive(Mat B, Vec X, Vec Y) in MatMult_LMVMSymBadBrdn_Recursive() argument
5 Mat_LMVM *lmvm = (Mat_LMVM *)B->data; in MatMult_LMVMSymBadBrdn_Recursive()
6 Mat_SymBrdn *lsb = (Mat_SymBrdn *)lmvm->ctx; in MatMult_LMVMSymBadBrdn_Recursive()
9 …PetscCheck(lsb->psi_scalar != PETSC_DETERMINE, PetscObjectComm((PetscObject)B), PETSC_ERR_ARG_WRON… in MatMult_LMVMSymBadBrdn_Recursive()
10 if (lsb->psi_scalar == 0.0) { in MatMult_LMVMSymBadBrdn_Recursive()
11 PetscCall(DFPKernel_Recursive(B, MATLMVM_MODE_PRIMAL, X, Y)); in MatMult_LMVMSymBadBrdn_Recursive()
12 } else if (lsb->psi_scalar == 1.0) { in MatMult_LMVMSymBadBrdn_Recursive()
13 PetscCall(BFGSKernel_Recursive(B, MATLMVM_MODE_PRIMAL, X, Y)); in MatMult_LMVMSymBadBrdn_Recursive()
15 PetscCall(SymBroydenKernel_Recursive(B, MATLMVM_MODE_PRIMAL, X, Y, PETSC_TRUE)); in MatMult_LMVMSymBadBrdn_Recursive()
20 static PetscErrorCode MatSolve_LMVMSymBadBrdn_Recursive(Mat B, Vec X, Vec Y) in MatSolve_LMVMSymBadBrdn_Recursive() argument
[all …]
H A Dsymbrdn.c8 PetscScalar numer = (1.0 - psi) * PetscRealPart(PetscConj(yts) * yts); in PhiFromPsi()
16 B_{k+1} = B_k + [ B_k s_k | y_k ] [-----+-----] [-----------]
22 B_{k+1} = B_0 + \sum_i [ B_i s_i | y_i ] [-----+-----] [-----------]
28 static PetscErrorCode SymBroydenKernel_Recursive_Inner(Mat B, MatLMVMMode mode, PetscInt oldest, Pe… in SymBroydenKernel_Recursive_Inner() argument
30 Mat_LMVM *lmvm = (Mat_LMVM *)B->data; in SymBroydenKernel_Recursive_Inner()
31 Mat_SymBrdn *lsb = (Mat_SymBrdn *)lmvm->ctx; in SymBroydenKernel_Recursive_Inner()
33 LMBasis BkS = lsb->basis[LMVMModeMap(SYMBROYDEN_BASIS_BKS, mode)]; in SymBroydenKernel_Recursive_Inner()
34 LMProducts M00 = lsb->products[LMVMModeMap(SYMBROYDEN_PRODUCTS_M00, mode)]; in SymBroydenKernel_Recursive_Inner()
35 LMProducts M01 = lsb->products[LMVMModeMap(SYMBROYDEN_PRODUCTS_M01, mode)]; in SymBroydenKernel_Recursive_Inner()
36 LMProducts M11 = lsb->products[LMVMModeMap(SYMBROYDEN_PRODUCTS_M11, mode)]; in SymBroydenKernel_Recursive_Inner()
[all …]
/petsc/src/mat/impls/sell/mpi/mpihip/
H A Dmpisellhip.hip.cxx5 static PetscErrorCode MatMPISELLSetPreallocation_MPISELLHIP(Mat B, PetscInt d_rlenmax, const PetscI… in MatMPISELLSetPreallocation_MPISELLHIP() argument
7 Mat_MPISELL *b = (Mat_MPISELL *)B->data; in MatMPISELLSetPreallocation_MPISELLHIP() local
10 PetscCall(PetscLayoutSetUp(B->rmap)); in MatMPISELLSetPreallocation_MPISELLHIP()
11 PetscCall(PetscLayoutSetUp(B->cmap)); in MatMPISELLSetPreallocation_MPISELLHIP()
13 if (!B->preallocated) { in MatMPISELLSetPreallocation_MPISELLHIP()
15 PetscCall(MatCreate(PETSC_COMM_SELF, &b->A)); in MatMPISELLSetPreallocation_MPISELLHIP()
16 PetscCall(MatBindToCPU(b->A, B->boundtocpu)); in MatMPISELLSetPreallocation_MPISELLHIP()
17 PetscCall(MatSetSizes(b->A, B->rmap->n, B->cmap->n, B->rmap->n, B->cmap->n)); in MatMPISELLSetPreallocation_MPISELLHIP()
18 PetscCall(MatSetType(b->A, MATSEQSELLHIP)); in MatMPISELLSetPreallocation_MPISELLHIP()
19 PetscCall(MatCreate(PETSC_COMM_SELF, &b->B)); in MatMPISELLSetPreallocation_MPISELLHIP()
[all …]
/petsc/src/mat/impls/sell/mpi/mpicuda/
H A Dmpisellcuda.cu5 static PetscErrorCode MatMPISELLSetPreallocation_MPISELLCUDA(Mat B, PetscInt d_rlenmax, const Petsc… in MatMPISELLSetPreallocation_MPISELLCUDA() argument
7 Mat_MPISELL *b = (Mat_MPISELL *)B->data; in MatMPISELLSetPreallocation_MPISELLCUDA() local
10 PetscCall(PetscLayoutSetUp(B->rmap)); in MatMPISELLSetPreallocation_MPISELLCUDA()
11 PetscCall(PetscLayoutSetUp(B->cmap)); in MatMPISELLSetPreallocation_MPISELLCUDA()
13 if (!B->preallocated) { in MatMPISELLSetPreallocation_MPISELLCUDA()
15 PetscCall(MatCreate(PETSC_COMM_SELF, &b->A)); in MatMPISELLSetPreallocation_MPISELLCUDA()
16 PetscCall(MatBindToCPU(b->A, B->boundtocpu)); in MatMPISELLSetPreallocation_MPISELLCUDA()
17 PetscCall(MatSetSizes(b->A, B->rmap->n, B->cmap->n, B->rmap->n, B->cmap->n)); in MatMPISELLSetPreallocation_MPISELLCUDA()
18 PetscCall(MatSetType(b->A, MATSEQSELLCUDA)); in MatMPISELLSetPreallocation_MPISELLCUDA()
19 PetscCall(MatCreate(PETSC_COMM_SELF, &b->B)); in MatMPISELLSetPreallocation_MPISELLCUDA()
[all …]
/petsc/src/mat/impls/sell/mpi/
H A Dmmsell.c11 Thus certain data-structures must be rebuilt.
18 Mat_MPISELL *sell = (Mat_MPISELL *)A->data; in MatDisAssemble_MPISELL()
19 Mat B = sell->B, Bnew; in MatDisAssemble_MPISELL() local
20 Mat_SeqSELL *Bsell = (Mat_SeqSELL *)B->data; in MatDisAssemble_MPISELL()
21 PetscInt i, j, totalslices, N = A->cmap->N, row; in MatDisAssemble_MPISELL()
25 /* free stuff related to matrix-vec multiply */ in MatDisAssemble_MPISELL()
26 PetscCall(VecDestroy(&sell->lvec)); in MatDisAssemble_MPISELL()
27 PetscCall(VecScatterDestroy(&sell->Mvctx)); in MatDisAssemble_MPISELL()
28 if (sell->colmap) { in MatDisAssemble_MPISELL()
30 PetscCall(PetscHMapIDestroy(&sell->colmap)); in MatDisAssemble_MPISELL()
[all …]

12345678910>>...45