Home
last modified time | relevance | path

Searched refs:A2 (Results 1 – 24 of 24) sorted by relevance

/petsc/src/mat/tests/
H A Dex23.c11 Mat A, B, A2, B2, T; in main() local
211 Mat A2; in main() local
232 PetscCall(MatCreate(PETSC_COMM_WORLD, &A2)); in main()
233 PetscCall(MatSetType(A2, MATIS)); in main()
239 PetscCall(MatLoad(A2, view)); in main()
240 PetscCall(CheckMat(A, A2, PETSC_TRUE, "Load")); in main()
241 PetscCall(MatView(A2, PETSC_VIEWER_STDOUT_WORLD)); in main()
247 PetscCall(MatLoad(A2, view)); in main()
248 PetscCall(MatView(A2, PETSC_VIEWER_STDOUT_WORLD)); in main()
254 PetscCall(MatLoad(A2, view)); in main()
[all …]
H A Dex5k.kokkos.cxx48 Mat A, A2; in main() local
77 PetscCall(MatCreate(PETSC_COMM_WORLD, &A2)); in main()
78 PetscCall(MatSetType(A2, MATAIJ)); in main()
79 PetscCall(MatLoad(A2, fd)); in main()
80 PetscCall(MatCreateVecs(A2, &x2, &y2)); in main()
83 PetscCall(MatGetSize(A2, &M, &N)); in main()
84 PetscCall(MatGetInfo(A2, MAT_GLOBAL_SUM, &info)); in main()
88 PetscCall(MatDuplicate(A2, MAT_COPY_VALUES, &A)); in main()
99 PetscCall(MatMult(A2, x2, y2)); in main()
127 PetscCall(MatDestroy(&A2)); in main()
H A Dex113.c10 Mat A, A1, A2, Mtmp, dstMat; in main() local
28 PetscCall(MatDuplicate(A, MAT_COPY_VALUES, &A2)); in main()
31 PetscCall(MatMatMult(A1, A2, MAT_INITIAL_MATRIX, fill, &Mtmp)); in main()
36 PetscCall(MatMatMult(A1, A2, MAT_INITIAL_MATRIX, fill, &Mtmp)); in main()
50 PetscCall(MatDestroy(&A2)); in main()
H A Dex6k.kokkos.cxx51 …Mat A2, P2, C2; /* Shadow matrices (of MATAIJ) of A,P,C for initialization and validati… in main() local
97 PetscCall(MatCreate(PETSC_COMM_WORLD, &A2)); in main()
98 PetscCall(MatSetType(A2, MATAIJ)); in main()
99 PetscCall(MatLoad(A2, fdA)); in main()
102 PetscCall(MatGetSize(A2, &M, &N)); in main()
103 PetscCall(MatGetInfo(A2, MAT_GLOBAL_SUM, &info)); in main()
107 PetscCall(MatDuplicate(A2, MAT_COPY_VALUES, &A)); in main()
126 P2 = A2; in main()
131 PetscCall(MatProductCreate(A2, P2, NULL, &C2)); in main()
189 PetscCall(MatDestroy(&A2)); in main()
H A Dex222.c55 Mat A2; in main() local
59 PetscCall(MatConvert(A, etype, MAT_INITIAL_MATRIX, &A2)); in main()
60 PetscCall(MatAXPY(A2, -1.0, Ae, DIFFERENT_NONZERO_PATTERN)); in main()
61 PetscCall(MatNorm(A2, NORM_FROBENIUS, &err)); in main()
63 PetscCall(MatDestroy(&A2)); in main()
H A Dex303.c9 Mat A2, B2; // diag, offdiag of mat2 in main() local
34 …PetscCall(MatConvert(A, mat_type, MAT_INITIAL_MATRIX, &A2)); // Copy A, B to A2, B2 b… in main()
45 PetscCall(MatCreateMPIAIJWithSeqAIJ(PETSC_COMM_WORLD, M, N, A2, B2, garray_h, &mat2)); in main()
52 PetscCall(MatConvert(A, mat_type, MAT_INITIAL_MATRIX, &A2)); in main()
70 PetscCall(MatCreateMPIAIJWithSeqAIJ(PETSC_COMM_WORLD, M, N, A2, B2, NULL, &mat2)); in main()
H A Dbench_spmv.c128 Mat A2 = NULL; in TimedSpMV() local
147 PetscCall(MatDestroy(&A2)); in TimedSpMV()
148 PetscCall(MatDuplicate(A, MAT_COPY_VALUES, &A2)); in TimedSpMV()
149 PetscCall(MatSetType(A2, petscmatformat)); in TimedSpMV()
150 …PetscCall(MatSetFromOptions(A2)); // This allows to change parameters such as slice height in SpMV… in TimedSpMV()
151 } else A2 = A; in TimedSpMV()
155 PetscCall(MatMult(A2, b, u)); in TimedSpMV()
163 if (repetitions > 0 && use_gpu) PetscCall(MatDestroy(&A2)); in TimedSpMV()
169 Mat A2 = NULL; in WarmUpDevice() local
185 PetscCall(MatDuplicate(A, MAT_COPY_VALUES, &A2)); in WarmUpDevice()
[all …]
H A Dex195.c14 Mat A1, A2, A3, A4, A5, B, C, C1, nest; in main() local
63 PetscCall(MatDuplicate(A1, MAT_COPY_VALUES, &A2)); in main()
83 Mat mata[] = {A1, A2, A3, A4}; in main()
152 PetscCall(MatDestroy(&A2)); in main()
H A Dex221.c137 Mat A2; in main() local
145 PetscCall(MatCreateSubMatrix(A, r, c, MAT_INITIAL_MATRIX, &A2)); in main()
149 A = A2; in main()
330 Mat A2, S2; in main() local
332 PetscCall(MatDuplicate(A, MAT_COPY_VALUES, &A2)); in main()
336 A = A2; in main()
H A Dex70.c260 Mat A2; in main() local
262 PetscCall(MatDuplicate(A, MAT_COPY_VALUES, &A2)); in main()
264 PetscCall(MatMultEqual(A, A2, 10, &flg)); in main()
270 PetscCall(MatComputeOperator(A2, MATDENSE, &A2E)); in main()
278 PetscCall(MatDestroy(&A2)); in main()
/petsc/src/ksp/ksp/tests/
H A Dex49.c8 Mat A, A2; in main() local
75 A2 = NULL; in main()
76 if (convert) PetscCall(MatConvert(A, MATAIJ, MAT_INITIAL_MATRIX, &A2)); in main()
99 PetscCall(KSPSetOperators(ksp, A2 ? A2 : A, A)); in main()
136 PetscCall(MatDestroy(&A2)); in main()
/petsc/src/ksp/ksp/tutorials/
H A Dex6f.F9021 Mat A, A2
114 PetscCallA(solve1(ksp, A, x, b, u, i, nsteps, A2, ierr))
152 Mat A2
167 PetscCallA(MatConvert(A, MATSAME, MAT_INITIAL_MATRIX, A2, ierr))
195 PetscCallA(KSPSetOperators(ksp, A, A2, ierr))
201 if (count == nsteps) PetscCallA(MatDestroy(A2, ierr))
/petsc/src/ts/tutorials/hybrid/
H A Dex1adj.c49 PetscScalar tmp[2], A1[2][2], A2[2], denorm; in ShiftGradients() local
63 A2[0] = 110. * u[1] * (-u[0]) / denorm; in ShiftGradients()
64 A2[1] = -110. * u[0] * (-u[0]) / denorm; in ShiftGradients()
72 A2[0] = 0; in ShiftGradients()
73 A2[1] = 0; in ShiftGradients()
82 y[0] = y[0] + A2[0] * x[0] + A2[1] * x[1]; in ShiftGradients()
92 y[0] = y[0] + A2[0] * x[0] + A2[1] * x[1]; in ShiftGradients()
H A Dex1fd.c49 PetscScalar tmp[2], A1[2][2], A2[2], denorm1, denorm2; in ShiftGradients() local
64 A2[0] = 110. * u[1] * (-u[0]) / denorm2; in ShiftGradients()
65 A2[1] = -110. * u[0] * (-u[0]) / denorm2; in ShiftGradients()
73 A2[0] = 0; in ShiftGradients()
74 A2[1] = 0; in ShiftGradients()
83 y[0] = y[0] + A2[0] * x[0] + A2[1] * x[1]; in ShiftGradients()
93 y[0] = y[0] + A2[0] * x[0] + A2[1] * x[1]; in ShiftGradients()
H A Dex1fwd.c71 PetscScalar tmp[2], A1[2][2], A2[2], denorm; in ShiftGradients() local
85 A2[0] = 110. * u[1] * (-u[0]) / denorm; in ShiftGradients()
86 A2[1] = -110. * u[0] * (-u[0]) / denorm; in ShiftGradients()
94 A2[0] = 0; in ShiftGradients()
95 A2[1] = 0; in ShiftGradients()
117 x[0] = tmp[0] + A2[0]; in ShiftGradients()
118 x[1] = tmp[1] + A2[1]; in ShiftGradients()
/petsc/src/tao/leastsquares/tutorials/
H A DtomographyGenerateData.m22 [A2, b2, xGT2] = PetscBinaryRead('tomographyData_A_b_xGT');
23 difference(full(A2), full(L));
78 [A2, b2, xGT2] = PetscBinaryRead('cs1Data_A_b_xGT');
/petsc/src/tao/matrix/
H A Dadamat.c91 Mat A2; in MatDuplicate_ADA() local
96 PetscCall(MatDuplicate(ctx->A, op, &A2)); in MatDuplicate_ADA()
103 PetscCall(MatCreateADA(A2, D1b, D2b, M)); in MatDuplicate_ADA()
106 PetscCall(PetscObjectDereference((PetscObject)A2)); in MatDuplicate_ADA()
/petsc/src/ts/impls/multirate/
H A Dmprk.c69 …l Abase[], const PetscReal bbase[], PetscReal A1[], PetscReal b1[], PetscReal A2[], PetscReal b2[]) in TSMPRKGenerateTableau2() argument
79 A2[(k * s + i) * ratio * s + k * s + j] = Abase[i * s + j] / ratio; in TSMPRKGenerateTableau2()
84 for (j = 0; j < s; j++) A2[(k * s + i) * ratio * s + l * s + j] = bbase[j] / ratio; in TSMPRKGenerateTableau2()
93 …base[], const PetscReal bbase[], PetscReal A1[], PetscReal b1[], PetscReal A2[], PetscReal b2[], P… in TSMPRKGenerateTableau3() argument
103 …A2[((k * ratio + l) * s + i) * ratio * ratio * s + (k * ratio + l) * s + j] = Abase[i * s + j] / r… in TSMPRKGenerateTableau3()
111 …A2[((k * ratio + m) * s + i) * ratio * ratio * s + (l * ratio + n) * s + j] = bbase[j] / ratio / r… in TSMPRKGenerateTableau3()
/petsc/src/mat/impls/baij/seq/
H A Dbaij.h707 …__m256d S0, S1, S2, S3, S4, S5, S6, S7, S8, B0, B1, B2, B6, B7, B8, A0, A1, A2, A3, A4, A5, A6, A… in PetscKernel_A_gets_A_times_B_9() local
723 A2 = _mm256_loadu_pd(W + 8); in PetscKernel_A_gets_A_times_B_9()
729 S2 = _mm256_fmadd_pd(A2, B0, S2); in PetscKernel_A_gets_A_times_B_9()
732 S5 = _mm256_fmadd_pd(A2, B1, S5); in PetscKernel_A_gets_A_times_B_9()
735 S8 = _mm256_fmadd_pd(A2, B2, S8); in PetscKernel_A_gets_A_times_B_9()
771 A2 = _mm256_loadu_pd(W + 35); in PetscKernel_A_gets_A_times_B_9()
777 S2 = _mm256_fmadd_pd(A2, B6, S2); in PetscKernel_A_gets_A_times_B_9()
780 S5 = _mm256_fmadd_pd(A2, B7, S5); in PetscKernel_A_gets_A_times_B_9()
783 S8 = _mm256_fmadd_pd(A2, B8, S8); in PetscKernel_A_gets_A_times_B_9()
803 A2 = _mm256_loadu_pd(W + 53); in PetscKernel_A_gets_A_times_B_9()
[all …]
/petsc/src/mat/impls/is/
H A Dmatis.c471 Mat A2, A3; in MatMPIXAIJComputeLocalToGlobalMapping_Private() local
528 PetscCall(MatCreate(comm, &A2)); in MatMPIXAIJComputeLocalToGlobalMapping_Private()
529 PetscCall(MatSetType(A2, MATMPIAIJ)); in MatMPIXAIJComputeLocalToGlobalMapping_Private()
530 PetscCall(MatSetSizes(A2, A->rmap->n, A->cmap->n, A->rmap->N, A->cmap->N)); in MatMPIXAIJComputeLocalToGlobalMapping_Private()
531 PetscCall(MatMPIAIJSetPreallocation(A2, 0, dnz, 0, onz)); in MatMPIXAIJComputeLocalToGlobalMapping_Private()
532 PetscCall(MatSetOption(A2, MAT_NO_OFF_PROC_ENTRIES, PETSC_TRUE)); in MatMPIXAIJComputeLocalToGlobalMapping_Private()
533 for (i = 0, j = 0; i < A2->rmap->n; i++) { in MatMPIXAIJComputeLocalToGlobalMapping_Private()
534 PetscInt row = i + A2->rmap->rstart, s = dnz[i] + onz[i]; in MatMPIXAIJComputeLocalToGlobalMapping_Private()
539 PetscCall(MatSetValues(A2, 1, &row, s, idxs, vals, INSERT_VALUES)); in MatMPIXAIJComputeLocalToGlobalMapping_Private()
546 PetscCall(MatAssemblyBegin(A2, MAT_FINAL_ASSEMBLY)); in MatMPIXAIJComputeLocalToGlobalMapping_Private()
[all …]
/petsc/share/petsc/matlab/
H A Dlaplacian.m55 % [~,~,A2] = laplacian([30,30]);
56 % norm(A1-A2,inf)
/petsc/
H A D.mailmap35 Jed Brown <jed@jedbrown.org> <jed@59A2.org>
/petsc/src/ts/tutorials/
H A Dex53.c1303 PetscScalar A2 = (alpha * (1.0 - 2.0*nu)) / (1.0 - nu);
1312 y1 = PetscTanReal(a1) - PetscRealPart(A1/A2)*a1;
1313 y2 = PetscTanReal(a2) - PetscRealPart(A1/A2)*a2;
1315 ym = PetscTanReal(am) - PetscRealPart(A1/A2)*am;
1342 PetscScalar sigma_zz = -1.0*(F/aL) - ((2.0*F)/aL) * (A2/A1) * A_x + ((2.0*F)/aL) * B_x;
/petsc/src/mat/impls/dense/seq/
H A Ddense.c1854 static PetscErrorCode MatEqual_SeqDense(Mat A1, Mat A2, PetscBool *flg) in MatEqual_SeqDense() argument
1857 Mat_SeqDense *mat2 = (Mat_SeqDense *)A2->data; in MatEqual_SeqDense()
1862 if (A1->rmap->n != A2->rmap->n) { in MatEqual_SeqDense()
1866 if (A1->cmap->n != A2->cmap->n) { in MatEqual_SeqDense()
1871 PetscCall(MatDenseGetArrayRead(A2, &v2)); in MatEqual_SeqDense()
1879 PetscCall(MatDenseRestoreArrayRead(A2, &v2)); in MatEqual_SeqDense()