Home
last modified time | relevance | path

Searched refs:a2 (Results 1 – 25 of 31) sorted by relevance

12

/petsc/include/petsc/private/kernels/
H A Dpetscaxpy.h38 #define PetscKernelAXPY2(U, a1, a2, p1, p2, n) \ argument
40 fortranmaxpy2_(U, &a1, &a2, p1, p2, &n); \
42 #define PetscKernelAXPY3(U, a1, a2, a3, p1, p2, p3, n) \ argument
44 fortranmaxpy3_(U, &a1, &a2, &a3, p1, p2, p3, &n); \
46 #define PetscKernelAXPY4(U, a1, a2, a3, a4, p1, p2, p3, p4, n) \ argument
48 fortranmaxpy4_(U, &a1, &a2, &a3, &a4, p1, p2, p3, p4, &n); \
76 #define PetscKernelAXPY2(U, a1, a2, p1, p2, n) \ argument
80 *U++ += a1 * *p1++ + a2 * *p2++; \
82 *U++ += a1 * *p1++ + a2 * *p2++; \
84 *U++ += a1 * *p1++ + a2 * *p2++; \
[all …]
/petsc/src/mat/impls/diagonal/tests/
H A Dex1.c8 Vec a, a2, b, b2, c, c2, A_diag, A_inv_diag; in main() local
24 PetscCall(VecDuplicate(a, &a2)); in main()
25 PetscCall(VecCopy(a, a2)); in main()
30 PetscCall(MatCreateDiagonal(a2, &A)); in main()
32 PetscCall(VecDestroy(&a2)); in main()
35 PetscCall(VecDuplicate(a, &a2)); in main()
85 PetscCall(MatGetDiagonal(A_dup, a2)); in main()
86 PetscCall(VecAXPY(a2, -1.0, a)); in main()
87 PetscCall(VecNorm(a2, NORM_INFINITY, &vec_norm)); in main()
93 PetscCall(MatGetDiagonal(A, a2)); in main()
[all …]
/petsc/src/dm/impls/stag/tests/
H A Dex7.c9 PetscScalar *a, ****a1, ****a2, expected; in main() local
40 PetscCall(DMStagVecGetArray(dm, vecLocal2, &a2)); in main()
46 a2[k][j][i][d] = 0.0; in main()
49 …for (is = -stencilWidth; is <= stencilWidth; ++is) a2[k][j][i][d] += a1[k + ks][j + js][i + is][d]; in main()
57 PetscCall(DMStagVecRestoreArray(dm, vecLocal2, &a2)); in main()
73 PetscCall(DMStagVecGetArrayRead(dm, vecLocal2, &a2)); in main()
91 if (a2[k][j][i][d] != expected) { in main()
92 …pected value %g (expecting %g)\n", rank, i, j, k, d, (double)PetscRealPart(a2[k][j][i][d]), (doubl… in main()
100 if (a2[k][j][i][d] != expected) { in main()
101 …pected value %g (expecting %g)\n", rank, i, j, k, d, (double)PetscRealPart(a2[k][j][i][d]), (doubl… in main()
[all …]
H A Dex10.c9 PetscScalar *a, ***a1, ***a2, expected; in main() local
39 PetscCall(DMStagVecGetArray(dm, vecLocal2, &a2)); in main()
44 a2[j][i][d] = 0.0; in main()
46 for (is = -stencilWidth; is <= stencilWidth; ++is) a2[j][i][d] += a1[j + js][i + is][d]; in main()
52 PetscCall(DMStagVecRestoreArray(dm, vecLocal2, &a2)); in main()
68 PetscCall(DMStagVecGetArrayRead(dm, vecLocal2, &a2)); in main()
83 if (a2[j][i][d] != expected) { in main()
84 …nexpected value %g (expecting %g)\n", rank, i, j, d, (double)PetscRealPart(a2[j][i][d]), (double)P… in main()
92 if (a2[j][i][d] != expected) { in main()
93 …nexpected value %g (expecting %g)\n", rank, i, j, d, (double)PetscRealPart(a2[j][i][d]), (double)P… in main()
[all …]
H A Dex9.c9 PetscScalar *a, ****a1, ****a2, expected, sum; in main() local
40 PetscCall(DMStagVecGetArray(dm, vecLocal2, &a2)); in main()
46 a2[k][j][i][d] = 0.0; in main()
47 for (ks = -stencilWidth; ks <= stencilWidth; ++ks) a2[k][j][i][d] += a1[k + ks][j][i][d]; in main()
48 for (js = -stencilWidth; js <= stencilWidth; ++js) a2[k][j][i][d] += a1[k][j + js][i][d]; in main()
49 for (is = -stencilWidth; is <= stencilWidth; ++is) a2[k][j][i][d] += a1[k][j][i + is][d]; in main()
50 a2[k][j][i][d] -= 2.0 * a1[k][j][i][d]; in main()
56 PetscCall(DMStagVecRestoreArray(dm, vecLocal2, &a2)); in main()
H A Dex12.c9 PetscScalar *a, ***a1, ***a2, expected, sum; in main() local
39 PetscCall(DMStagVecGetArray(dm, vecLocal2, &a2)); in main()
44 a2[j][i][d] = 0.0; in main()
45 for (js = -stencilWidth; js <= stencilWidth; ++js) a2[j][i][d] += a1[j + js][i][d]; in main()
46 for (is = -stencilWidth; is <= stencilWidth; ++is) a2[j][i][d] += a1[j][i + is][d]; in main()
47 a2[j][i][d] -= a1[j][i][d]; in main()
52 PetscCall(DMStagVecRestoreArray(dm, vecLocal2, &a2)); in main()
/petsc/src/vec/is/is/utils/
H A Discomp.c37 PetscInt sz1, sz2, *a1, *a2; in ISEqual() local
72 PetscCall(PetscMalloc1(sz2, &a2)); in ISEqual()
75 PetscCall(PetscArraycpy(a2, ptr2, sz2)); in ISEqual()
78 PetscCall(PetscIntSortSemiOrdered(sz2, a2)); in ISEqual()
79 PetscCall(PetscArraycmp(a1, a2, sz1, &flag)); in ISEqual()
85 PetscCall(PetscFree(a2)); in ISEqual()
/petsc/src/vec/vec/impls/seq/ftn-kernels/
H A Dfmaxpy.F9011 PetscScalar, intent(in) :: a0, a1, a2, a3
22 x(1:n) = x(1:n) + (a0*y0(1:n) + a1*y1(1:n) + a2*y2(1:n) + a3*y3(1:n))
28 PetscScalar, intent(in) :: a0, a1, a2
38 x(1:n) = x(1:n) + (a0*y0(1:n) + a1*y1(1:n) + a2*y2(1:n))
/petsc/src/snes/tutorials/
H A Dex74ftube.h1 …ble precision ru1, ru2, ru3, ru4, et1, et2, et3, et4 double precision csubv, gamma, a1, a2, a3, a4,
6 … e3, e4, v1, v2, v3, v4, ru1, ru2, ru3, ru4, et1, et2, et3, et4, csubv, gamma, a1, a2, a3, a4, sspd
/petsc/src/sys/classes/draw/impls/image/
H A Ddrawimage.h87 int x, y, s, a2 = w * w, b2 = h * h; in PetscImageDrawEllipse() local
88 for (x = 0, y = h, s = 2 * b2 + a2 * (1 - 2 * h); b2 * x <= a2 * y; x++) { in PetscImageDrawEllipse()
92 s += 4 * a2 * (1 - y); in PetscImageDrawEllipse()
97 for (x = w, y = 0, s = 2 * a2 + b2 * (1 - 2 * w); a2 * y <= b2 * x; y++) { in PetscImageDrawEllipse()
104 s += a2 * ((4 * y) + 6); in PetscImageDrawEllipse()
/petsc/src/mat/tests/
H A Dex121.c23 PetscScalar *a, *a2, *a3; in main() local
96 PetscCall(VecGetArray(w, &a2)); in main()
106 a3[i] += a[xpInd] * a2[diffInd]; in main()
110 PetscCall(VecRestoreArray(w, &a2)); in main()
H A Dex90.c50 PetscScalar a1[] = {1, 2, 4, 1, 2}, a2[] = {2, 4, 1, 2, 1}; in main() local
63 …comm, 2, 2, PETSC_DETERMINE, PETSC_DETERMINE, rank ? i2 : i1, rank ? j2 : j1, rank ? a2 : a1, &A)); in main()
/petsc/src/mat/impls/baij/seq/
H A Dbaijfact9.c886 __m256d a0, a1, a2, a3, a4, a5, w0, w1, w2, w3, s0, s1, s2, v0, v1, v2, v3; in MatSolve_SeqBAIJ_9_NaturalOrdering() local
914 a2 = _mm256_loadu_pd(&v[8]); in MatSolve_SeqBAIJ_9_NaturalOrdering()
915 s2 = _mm256_fnmadd_pd(a2, w0, s2); in MatSolve_SeqBAIJ_9_NaturalOrdering()
930 a2 = _mm256_loadu_pd(&v[26]); in MatSolve_SeqBAIJ_9_NaturalOrdering()
931 s2 = _mm256_fnmadd_pd(a2, w2, s2); in MatSolve_SeqBAIJ_9_NaturalOrdering()
946 a2 = _mm256_loadu_pd(&v[44]); in MatSolve_SeqBAIJ_9_NaturalOrdering()
947 s2 = _mm256_fnmadd_pd(a2, w0, s2); in MatSolve_SeqBAIJ_9_NaturalOrdering()
962 a2 = _mm256_loadu_pd(&v[62]); in MatSolve_SeqBAIJ_9_NaturalOrdering()
963 s2 = _mm256_fnmadd_pd(a2, w2, s2); in MatSolve_SeqBAIJ_9_NaturalOrdering()
978 a2 = _mm256_maskload_pd(v + 80, _mm256_set_epi64x(0LL, 0LL, 0LL, 1LL << 63)); in MatSolve_SeqBAIJ_9_NaturalOrdering()
[all …]
H A Dbaij2.c679 __m256d a0, a1, a2, a3, a4, a5; in MatMult_SeqBAIJ_9_AVX2() local
729 a2 = _mm256_loadu_pd(&v[j * 81 + 8]); in MatMult_SeqBAIJ_9_AVX2()
730 z2 = _mm256_fmadd_pd(a2, w0, z2); in MatMult_SeqBAIJ_9_AVX2()
738 a2 = _mm256_loadu_pd(&v[j * 81 + 17]); in MatMult_SeqBAIJ_9_AVX2()
739 z2 = _mm256_fmadd_pd(a2, w1, z2); in MatMult_SeqBAIJ_9_AVX2()
756 a2 = _mm256_loadu_pd(&v[j * 81 + 35]); in MatMult_SeqBAIJ_9_AVX2()
757 z2 = _mm256_fmadd_pd(a2, w3, z2); in MatMult_SeqBAIJ_9_AVX2()
774 a2 = _mm256_loadu_pd(&v[j * 81 + 53]); in MatMult_SeqBAIJ_9_AVX2()
775 z2 = _mm256_fmadd_pd(a2, w1, z2); in MatMult_SeqBAIJ_9_AVX2()
783 a2 = _mm256_loadu_pd(&v[j * 81 + 62]); in MatMult_SeqBAIJ_9_AVX2()
[all …]
/petsc/src/ts/impls/rosw/
H A Drosw.c961 PetscErrorCode TSRosWRegisterRos4(TSRosWType name, PetscReal gamma, PetscReal a2, PetscReal a3, Pet… in TSRosWRegisterRos4() argument
972 …FAULT || a3 == (PetscReal)PETSC_DETERMINE) a3 = (one / five - a2 / four) / (one / four - a2 / thre… in TSRosWRegisterRos4()
980 M[1][1] = a2 * a2; in TSRosWRegisterRos4()
983 M[2][1] = a2 * a2 * a2; in TSRosWRegisterRos4()
994 beta43 = (p56 - a2 * p43) / (b4 * a3 * a3 * (a3 - a2)); /* 7.21 */ in TSRosWRegisterRos4()
1001 M[1][1] = a2 * a2 * beta4jbetajp; in TSRosWRegisterRos4()
1002 M[1][2] = -a2 * a2 * beta32beta2p; in TSRosWRegisterRos4()
1004 M[2][1] = -b4 * beta43 * a2 * a2; in TSRosWRegisterRos4()
1008 rhs[2] = -a2 * a2 * p32; in TSRosWRegisterRos4()
1027 A[1][0] = a2; in TSRosWRegisterRos4()
[all …]
/petsc/src/dm/dt/tests/
H A Dex1.c207 PetscReal a1, b1, a2, b2; in main() local
212 a2 = 2.2; in main()
217 a2 = 2.; in main()
222 PetscCall(CheckJacobiQuadrature(i, a2, b2, PetscDTGaussJacobiQuadrature, 2 * i - 1)); in main()
226 PetscCall(CheckJacobiQuadrature(i, a2, b2, PetscDTGaussLobattoJacobiQuadrature, 2 * i - 3)); in main()
/petsc/src/ts/adapt/impls/dsp/
H A Dadaptdsp.c124 PetscReal a2 = dsp->Alpha[0]; in TSAdaptChoose_DSP() local
140 rho0 *= PetscPowReal(rho1, -a2); in TSAdaptChoose_DSP()
187 double a2 = (double)dsp->Alpha[0], a3 = (double)dsp->Alpha[1]; in TSAdaptView_DSP() local
189 …ewerASCIIPrintf(viewer, "filter parameters kBeta=[%g,%g,%g] Alpha=[%g,%g]\n", b1, b2, b3, a2, a3)); in TSAdaptView_DSP()
/petsc/src/snes/impls/ls/
H A Dls.c90 PetscReal a1, a2; in SNESNEWTONLSCheckResidual_Private() local
108 PetscCall(VecNorm(W2, NORM_2, &a2)); in SNESNEWTONLSCheckResidual_Private()
109 …o(snes, "||J^T(F-Ax)||/||F-AX|| %14.12e near zero implies inconsistent rhs\n", (double)(a2 / a1))); in SNESNEWTONLSCheckResidual_Private()
/petsc/src/ksp/ksp/tutorials/
H A Dex54f.F9036 PetscReal::thk, a1, a2
128 a2 = (dd(2, 1)*shp(1, kj) + dd(2, 2)*shp(2, kj))*xsj
134 ss(i1, j1) = ss(i1, j1) + a1*shp(1, ki) + a2*shp(2, ki)
/petsc/src/binding/petsc4py/test/
H A Dtest_vec.py329 a2 = v2.getDict()['__array__']
330 self.assertTrue(a is a2)
347 a2 = v2.getDict()['__array__']
348 self.assertTrue(a is a2)
/petsc/src/ts/tutorials/
H A Dex53.c1267 PetscReal a1, a2, am;
1309 a2 = a1 + PETSC_PI/2;
1313 y2 = PetscTanReal(a2) - PetscRealPart(A1/A2)*a2;
1314 am = (a1 + a2)/2.0;
1320 a2 = am;
1324 am = a2;
1515 PetscReal a1, a2, am; in mandelZeros() local
1540 a2 = a1 + PETSC_PI / 2; in mandelZeros()
1544 y2 = PetscTanReal(a2) - PetscRealPart((1.0 - nu) / (nu_u - nu)) * a2; in mandelZeros()
1545 am = (a1 + a2) / 2.0; in mandelZeros()
[all …]
/petsc/share/petsc/matlab/
H A Dlaplacian.m248 a2 = pi/2/(u(2)+1); variable
251 a2 = pi/2/u(2); variable
255 a2 = pi/4/(u(2)+0.5); variable
258 a2 = pi/u(2); variable
261 lambda2 = 4*sin(a2*N).^2;
/petsc/doc/tutorials/
H A Din_person.md42 [Slides](https://59a2.org/files/20130618-PRACE.pdf)
63 [Slides](https://59a2.org/files/20120220-TACC.pdf).
/petsc/include/
H A Dpetscmacros.h1203 #define PETSC_SELECT_16TH(a1, a2, a3, a4, a5, a6, a7, a8, a9, a10, a11, a12, a13, a14, a15, a16, ..… argument
/petsc/lib/petsc/bin/maint/petsclinter/petsclinter/classes/docs/
H A D_doc_section_base.py900 i for i, (a1, a2) in enumerate(itertools.zip_longest(text, fixed)) if a1 != a2

12