Home
last modified time | relevance | path

Searched refs:A1 (Results 1 – 13 of 13) sorted by relevance

/petsc/src/mat/tests/
H A Dex195.c14 Mat A1, A2, A3, A4, A5, B, C, C1, nest; in main() local
30 PetscCall(MatCreate(comm, &A1)); in main()
32 PetscCall(MatSetSizes(A1, PETSC_DECIDE, PETSC_DECIDE, m * n, m * n)); in main()
33 PetscCall(MatSetFromOptions(A1)); in main()
34 PetscCall(MatSetUp(A1)); in main()
35 PetscCall(MatGetOwnershipRange(A1, &istart, &iend)); in main()
42 PetscCall(MatSetValues(A1, 1, &ii, 1, &J, &v, INSERT_VALUES)); in main()
46 PetscCall(MatSetValues(A1, 1, &ii, 1, &J, &v, INSERT_VALUES)); in main()
50 PetscCall(MatSetValues(A1, 1, &ii, 1, &J, &v, INSERT_VALUES)); in main()
54 PetscCall(MatSetValues(A1, 1, &ii, 1, &J, &v, INSERT_VALUES)); in main()
[all …]
H A Dex113.c10 Mat A, A1, A2, Mtmp, dstMat; in main() local
27 PetscCall(MatDuplicate(A, MAT_COPY_VALUES, &A1)); 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()
41 PetscCall(MatMatMult(A, A1, MAT_INITIAL_MATRIX, fill, &Mtmp)); in main()
49 PetscCall(MatDestroy(&A1)); in main()
/petsc/src/ts/tutorials/hybrid/
H A Dex1fwd.c71 PetscScalar tmp[2], A1[2][2], A2[2], denorm; in ShiftGradients() local
80 A1[0][0] = 110. * u[1] * (-actx->lambda1) / denorm + 1.; in ShiftGradients()
81 A1[1][0] = -110. * u[0] * (-actx->lambda1) / denorm; in ShiftGradients()
82 A1[0][1] = 110. * u[1] * 1. / denorm; in ShiftGradients()
83 A1[1][1] = -110. * u[0] * 1. / denorm + 1.; in ShiftGradients()
89 A1[0][0] = 110. * u[1] * (actx->lambda2) / denorm + 1; in ShiftGradients()
90 A1[1][0] = -110. * u[0] * (actx->lambda2) / denorm; in ShiftGradients()
91 A1[0][1] = -110. * u[1] * 1. / denorm; in ShiftGradients()
92 A1[1][1] = 110. * u[0] * 1. / denorm + 1.; in ShiftGradients()
101 tmp[0] = A1[0][0] * x[0] + A1[0][1] * x[1]; in ShiftGradients()
[all …]
H A Dex1fd.c49 PetscScalar tmp[2], A1[2][2], A2[2], denorm1, denorm2; in ShiftGradients() local
59 A1[0][0] = 110. * u[1] * (-actx->lambda1) / denorm1 + 1.; in ShiftGradients()
60 A1[0][1] = -110. * u[0] * (-actx->lambda1) / denorm1; in ShiftGradients()
61 A1[1][0] = 110. * u[1] * 1. / denorm1; in ShiftGradients()
62 A1[1][1] = -110. * u[0] * 1. / denorm1 + 1.; in ShiftGradients()
68 A1[0][0] = 110. * u[1] * (-actx->lambda1) / denorm2 + 1; in ShiftGradients()
69 A1[0][1] = -110. * u[0] * (-actx->lambda1) / denorm2; in ShiftGradients()
70 A1[1][0] = 110. * u[1] * 1. / denorm2; in ShiftGradients()
71 A1[1][1] = -110. * u[0] * 1. / denorm2 + 1.; in ShiftGradients()
81 tmp[0] = A1[0][0] * x[0] + A1[0][1] * x[1]; in ShiftGradients()
[all …]
H A Dex1adj.c49 PetscScalar tmp[2], A1[2][2], A2[2], denorm; in ShiftGradients() local
58 A1[0][0] = 110. * u[1] * (-actx->lambda1) / denorm + 1.; in ShiftGradients()
59 A1[0][1] = -110. * u[0] * (-actx->lambda1) / denorm; in ShiftGradients()
60 A1[1][0] = 110. * u[1] * 1. / denorm; in ShiftGradients()
61 A1[1][1] = -110. * u[0] * 1. / denorm + 1.; in ShiftGradients()
67 A1[0][0] = 110. * u[1] * (actx->lambda2) / denorm + 1; in ShiftGradients()
68 A1[0][1] = -110. * u[0] * (actx->lambda2) / denorm; in ShiftGradients()
69 A1[1][0] = -110. * u[1] * 1. / denorm; in ShiftGradients()
70 A1[1][1] = 110. * u[0] * 1. / denorm + 1.; in ShiftGradients()
80 tmp[0] = A1[0][0] * x[0] + A1[0][1] * x[1]; in ShiftGradients()
[all …]
/petsc/src/ts/impls/implicit/alpha/
H A Dalpha2.c26 Vec A0, Aa, A1; member
107 Vec X1 = X, V1 = th->V1, A1 = th->A1; in TSAlpha_StageVecs() local
118 PetscCall(VecWAXPY(A1, -1.0, X0, X1)); in TSAlpha_StageVecs()
119 PetscCall(VecAXPY(A1, -dt, V0)); in TSAlpha_StageVecs()
120 PetscCall(VecAXPBY(A1, -(1 - 2 * Beta) / (2 * Beta), 1 / (dt * dt * Beta), A0)); in TSAlpha_StageVecs()
122 PetscCall(VecWAXPY(V1, (1.0 - Gamma) / Gamma, A0, A1)); in TSAlpha_StageVecs()
131 PetscCall(VecWAXPY(Aa, -1.0, A0, A1)); in TSAlpha_StageVecs()
250 PetscCall(VecCopy(th->A1, th->A0)); in TSStep_Alpha()
422 PetscCall(VecDestroy(&th->A1)); in TSReset_Alpha()
456 PetscCall(VecDuplicate(ts->vec_sol, &th->A1)); in TSSetUp_Alpha()
/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
722 A1 = _mm256_loadu_pd(W + 4); in PetscKernel_A_gets_A_times_B_9()
728 S1 = _mm256_fmadd_pd(A1, B0, S1); in PetscKernel_A_gets_A_times_B_9()
731 S4 = _mm256_fmadd_pd(A1, B1, S4); in PetscKernel_A_gets_A_times_B_9()
734 S7 = _mm256_fmadd_pd(A1, B2, S7); in PetscKernel_A_gets_A_times_B_9()
770 A1 = _mm256_loadu_pd(W + 31); in PetscKernel_A_gets_A_times_B_9()
776 S1 = _mm256_fmadd_pd(A1, B6, S1); in PetscKernel_A_gets_A_times_B_9()
779 S4 = _mm256_fmadd_pd(A1, B7, S4); in PetscKernel_A_gets_A_times_B_9()
782 S7 = _mm256_fmadd_pd(A1, B8, S7); in PetscKernel_A_gets_A_times_B_9()
802 A1 = _mm256_loadu_pd(W + 49); in PetscKernel_A_gets_A_times_B_9()
[all …]
/petsc/src/ksp/ksp/impls/cg/pipecg2/
H A Dpipecg2.c555 Vec X, B, Z, P, W, Q, U, M, N, R, S, C, D, E, F, G[2], H[2], A1, B1; in KSPSolve_PIPECG2() local
586 A1 = ksp->work[17]; in KSPSolve_PIPECG2()
661 …PetscCall(VecMergedOpsShort_Private(X, R, Z, W, P, Q, C, D, G[0], H[0], G[1], H[1], S, A1, B1, E, … in KSPSolve_PIPECG2()
676 …PetscCall(VecMergedOps_Private(X, R, Z, W, P, Q, C, D, G[0], H[0], G[1], H[1], S, A1, B1, E, F, M,… in KSPSolve_PIPECG2()
/petsc/src/ts/tutorials/
H A Dex53.c691 PetscScalar A1 = 3.0 / (B * (1.0 + nu_u)); in mandel_drainage_pressure() local
700 u[0] = ((2.0 * P_0) / (a * A1)) * p; in mandel_drainage_pressure()
902 PetscScalar A1 = 3.0 / (B * (1.0 + nu_u)); in mandel_2d_p() local
912 u[0] = ((2.0 * F) / (a * A1)) * p; in mandel_2d_p()
1302 PetscScalar A1 = 3.0 / (B * (1.0 + nu_u));
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/share/petsc/matlab/
H A Dlaplacian.m54 % A1 = delsq(numgrid('S',32)); % input 'S' specifies square grid.
56 % norm(A1-A2,inf)
/petsc/src/ts/impls/multirate/
H A Dmprk.c69 …o, PetscInt s, const PetscReal Abase[], const PetscReal bbase[], PetscReal A1[], PetscReal b1[], P… in TSMPRKGenerateTableau2() argument
78 A1[(k * s + i) * ratio * s + k * s + j] = Abase[i * s + j]; in TSMPRKGenerateTableau2()
93 …o, PetscInt s, const PetscReal Abase[], const PetscReal bbase[], PetscReal A1[], PetscReal b1[], P… in TSMPRKGenerateTableau3() argument
102 … A1[((k * ratio + l) * s + i) * ratio * ratio * s + (k * ratio + l) * s + j] = Abase[i * s + j]; in TSMPRKGenerateTableau3()
/petsc/src/mat/impls/dense/seq/
H A Ddense.c1854 static PetscErrorCode MatEqual_SeqDense(Mat A1, Mat A2, PetscBool *flg) in MatEqual_SeqDense() argument
1856 Mat_SeqDense *mat1 = (Mat_SeqDense *)A1->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()
1870 PetscCall(MatDenseGetArrayRead(A1, &v1)); in MatEqual_SeqDense()
1872 for (i = 0; i < A1->cmap->n; i++) { in MatEqual_SeqDense()
1873 PetscCall(PetscArraycmp(v1, v2, A1->rmap->n, flg)); in MatEqual_SeqDense()
1878 PetscCall(MatDenseRestoreArrayRead(A1, &v1)); in MatEqual_SeqDense()
/petsc/doc/manual/
H A Dksp.md1853 KSPSetOperators(ksp,A1,A1);