Home
last modified time | relevance | path

Searched refs:a1 (Results 1 – 25 of 26) sorted by relevance

12

/petsc/include/petsc/private/kernels/
H A Dpetscaxpy.h32 #define PetscKernelAXPY(U, a1, p1, n) \ argument
36 PetscCallBLAS("BLASaxpy", BLASaxpy_(&nn, &a1, p1, &one, U, &one)); \
38 #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++; \
[all …]
/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))
44 PetscScalar, intent(in) :: a0, a1
53 x(1:n) = x(1:n) + (a0*y0(1:n) + a1*y1(1:n))
/petsc/src/dm/impls/stag/tests/
H A Dex9.c9 PetscScalar *a, ****a1, ****a2, expected, sum; in main() local
39 PetscCall(DMStagVecGetArrayRead(dm, vecLocal1, &a1)); in main()
45 …if (a1[k][j][i][d] != 1.0) PetscCall(PetscPrintf(PETSC_COMM_SELF, "[%d] Unexpected value %g (expec… 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()
55 PetscCall(DMStagVecRestoreArrayRead(dm, vecLocal1, &a1)); in main()
H A Dex12.c9 PetscScalar *a, ***a1, ***a2, expected, sum; in main() local
38 PetscCall(DMStagVecGetArrayRead(dm, vecLocal1, &a1)); in main()
43 …if (a1[j][i][d] != 1.0) PetscCall(PetscPrintf(PETSC_COMM_SELF, "[%d] Unexpected value %g (expectin… 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()
51 PetscCall(DMStagVecRestoreArrayRead(dm, vecLocal1, &a1)); in main()
H A Dex10.c9 PetscScalar *a, ***a1, ***a2, expected; in main() local
38 PetscCall(DMStagVecGetArrayRead(dm, vecLocal1, &a1)); in main()
43 …if (a1[j][i][d] != 1.0) PetscCall(PetscPrintf(PETSC_COMM_SELF, "[%d] Unexpected value %g (expectin… in main()
46 for (is = -stencilWidth; is <= stencilWidth; ++is) a2[j][i][d] += a1[j + js][i + is][d]; in main()
51 PetscCall(DMStagVecRestoreArrayRead(dm, vecLocal1, &a1)); in main()
H A Dex7.c9 PetscScalar *a, ****a1, ****a2, expected; in main() local
39 PetscCall(DMStagVecGetArrayRead(dm, vecLocal1, &a1)); in main()
45 …if (a1[k][j][i][d] != 1.0) PetscCall(PetscPrintf(PETSC_COMM_SELF, "[%d] Unexpected value %g (expec… in main()
49 …for (is = -stencilWidth; is <= stencilWidth; ++is) a2[k][j][i][d] += a1[k + ks][j + js][i + is][d]; in main()
56 PetscCall(DMStagVecRestoreArrayRead(dm, vecLocal1, &a1)); in main()
/petsc/src/vec/is/is/utils/
H A Discomp.c37 PetscInt sz1, sz2, *a1, *a2; in ISEqual() local
71 PetscCall(PetscMalloc1(sz1, &a1)); in ISEqual()
74 PetscCall(PetscArraycpy(a1, ptr1, sz1)); in ISEqual()
77 PetscCall(PetscIntSortSemiOrdered(sz1, a1)); in ISEqual()
79 PetscCall(PetscArraycmp(a1, a2, sz1, &flag)); in ISEqual()
84 PetscCall(PetscFree(a1)); in ISEqual()
/petsc/src/snes/impls/ls/
H A Dls.c48 PetscReal a1; in SNESNEWTONLSCheckLocalMin_Private() local
62 PetscCall(VecNorm(W, NORM_2, &a1)); in SNESNEWTONLSCheckLocalMin_Private()
63 …cInfo(snes, "|| J^T F|| %14.12e near zero implies found a local minimum\n", (double)(a1 / fnorm))); in SNESNEWTONLSCheckLocalMin_Private()
64 if (a1 / fnorm < 1.e-4) *ismin = PETSC_TRUE; in SNESNEWTONLSCheckLocalMin_Private()
76 a1 = PetscAbsScalar(result) / (fnorm * wnorm); in SNESNEWTONLSCheckLocalMin_Private()
77 …F^T J random)/(|| F ||*||random|| %14.12e near zero implies found a local minimum\n", (double)a1)); in SNESNEWTONLSCheckLocalMin_Private()
78 if (a1 < 1.e-4) *ismin = PETSC_TRUE; in SNESNEWTONLSCheckLocalMin_Private()
90 PetscReal a1, a2; in SNESNEWTONLSCheckResidual_Private() local
107 PetscCall(VecNorm(W1, NORM_2, &a1)); in SNESNEWTONLSCheckResidual_Private()
109 …if (a1 != 0.0) PetscCall(PetscInfo(snes, "||J^T(F-Ax)||/||F-AX|| %14.12e near zero implies inconsi… in SNESNEWTONLSCheckResidual_Private()
/petsc/src/ts/tutorials/
H A Dex24.c166 PetscScalar a, a0, a1; in FormIFunction() local
169 a1 = 1.; in FormIFunction()
171 f[i + 1] += 200. * a * a1; in FormIFunction()
213 PetscScalar v[2][2], a, a0, a1, a00, a01, a10, a11; in FormIJacobian() local
220 a1 = 1.; in FormIJacobian()
224 v[0][1] = 200. * (a * a01 + a1 * a0); in FormIJacobian()
225 v[1][0] = 200. * (a * a10 + a0 * a1); in FormIJacobian()
226 v[1][1] = 200. * (a * a11 + a1 * a1); in FormIJacobian()
H A Dex53.c1267 PetscReal a1, a2, am;
1308 a1 = ((PetscReal) i - 1.0) * PETSC_PI * PETSC_PI / 4.0 + EPS;
1309 a2 = a1 + PETSC_PI/2;
1312 y1 = PetscTanReal(a1) - PetscRealPart(A1/A2)*a1;
1314 am = (a1 + a2)/2.0;
1318 a1 = am;
1515 PetscReal a1, a2, am; in mandelZeros() local
1539 a1 = ((PetscReal)i - 1.0) * PETSC_PI * PETSC_PI / 4.0 + EPS; in mandelZeros()
1540 a2 = a1 + PETSC_PI / 2; in mandelZeros()
1541 am = a1; in mandelZeros()
[all …]
/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/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
912 a1 = _mm256_loadu_pd(&v[4]); in MatSolve_SeqBAIJ_9_NaturalOrdering()
913 s1 = _mm256_fnmadd_pd(a1, w0, s1); in MatSolve_SeqBAIJ_9_NaturalOrdering()
928 a1 = _mm256_loadu_pd(&v[22]); in MatSolve_SeqBAIJ_9_NaturalOrdering()
929 s1 = _mm256_fnmadd_pd(a1, w2, s1); in MatSolve_SeqBAIJ_9_NaturalOrdering()
944 a1 = _mm256_loadu_pd(&v[40]); in MatSolve_SeqBAIJ_9_NaturalOrdering()
945 s1 = _mm256_fnmadd_pd(a1, w0, s1); in MatSolve_SeqBAIJ_9_NaturalOrdering()
960 a1 = _mm256_loadu_pd(&v[58]); in MatSolve_SeqBAIJ_9_NaturalOrdering()
961 s1 = _mm256_fnmadd_pd(a1, w2, s1); in MatSolve_SeqBAIJ_9_NaturalOrdering()
976 a1 = _mm256_loadu_pd(&v[76]); in MatSolve_SeqBAIJ_9_NaturalOrdering()
[all …]
H A Dbaij2.c679 __m256d a0, a1, a2, a3, a4, a5; in MatMult_SeqBAIJ_9_AVX2() local
727 a1 = _mm256_loadu_pd(&v[j * 81 + 4]); in MatMult_SeqBAIJ_9_AVX2()
728 z1 = _mm256_fmadd_pd(a1, w0, z1); in MatMult_SeqBAIJ_9_AVX2()
736 a1 = _mm256_loadu_pd(&v[j * 81 + 13]); in MatMult_SeqBAIJ_9_AVX2()
737 z1 = _mm256_fmadd_pd(a1, w1, z1); in MatMult_SeqBAIJ_9_AVX2()
754 a1 = _mm256_loadu_pd(&v[j * 81 + 31]); in MatMult_SeqBAIJ_9_AVX2()
755 z1 = _mm256_fmadd_pd(a1, w3, z1); in MatMult_SeqBAIJ_9_AVX2()
772 a1 = _mm256_loadu_pd(&v[j * 81 + 49]); in MatMult_SeqBAIJ_9_AVX2()
773 z1 = _mm256_fmadd_pd(a1, w1, z1); in MatMult_SeqBAIJ_9_AVX2()
781 a1 = _mm256_loadu_pd(&v[j * 81 + 58]); in MatMult_SeqBAIJ_9_AVX2()
[all …]
/petsc/src/mat/tests/
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/dm/dt/tests/
H A Dex1.c207 PetscReal a1, b1, a2, b2; in main() local
210 a1 = -0.6; in main()
215 a1 = 0.; in main()
221 PetscCall(CheckJacobiQuadrature(i, a1, b1, PetscDTGaussJacobiQuadrature, 2 * i - 1)); in main()
225 PetscCall(CheckJacobiQuadrature(i, a1, b1, PetscDTGaussLobattoJacobiQuadrature, 2 * i - 3)); in main()
/petsc/src/snes/impls/vi/
H A Dvi.c154 PetscReal a1; in SNESVICheckLocalMin_Private() local
163 PetscCall(VecNorm(W, NORM_2, &a1)); in SNESVICheckLocalMin_Private()
164 …(PetscInfo(snes, "|| J^T F|| %g near zero implies found a local minimum\n", (double)(a1 / fnorm))); in SNESVICheckLocalMin_Private()
165 if (a1 / fnorm < 1.e-4) *ismin = PETSC_TRUE; in SNESVICheckLocalMin_Private()
177 a1 = PetscAbsScalar(result) / (fnorm * wnorm); in SNESVICheckLocalMin_Private()
178 …s, "(F^T J random)/(|| F ||*||random|| %g near zero implies found a local minimum\n", (double)a1)); in SNESVICheckLocalMin_Private()
179 if (a1 < 1.e-4) *ismin = PETSC_TRUE; in SNESVICheckLocalMin_Private()
/petsc/src/ksp/ksp/tutorials/
H A Dex54f.F9036 PetscReal::thk, a1, a2
127 a1 = (dd(1, 1)*shp(1, kj) + dd(1, 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.py327 a1 = v1.getDict()['__array__']
328 self.assertTrue(a is a1)
345 a1 = v1.getDict()['__array__']
346 self.assertTrue(a is a1)
/petsc/share/petsc/matlab/
H A Dlaplacian.m230 a1 = pi/2/(u(1)+1); variable
233 a1 = pi/2/u(1); variable
237 a1 = pi/4/(u(1)+0.5); variable
240 a1 = pi/u(1); variable
244 lambda1 = 4*sin(a1*N).^2;
/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
/petsc/src/mat/impls/sbaij/mpi/
H A Dmpisbaij.c1321 PetscInt sz1, sz2, *a1, *a2, i, j, k, nmatch; in ISEqual_private() local
1333 PetscCall(PetscMalloc1(sz1, &a1)); in ISEqual_private()
1335 PetscCall(PetscArraycpy(a1, ptr1, sz1)); in ISEqual_private()
1337 PetscCall(PetscSortInt(sz1, a1)); in ISEqual_private()
1344 if (a1[i] == a2[j]) { in ISEqual_private()
1353 PetscCall(PetscFree(a1)); in ISEqual_private()
/petsc/src/vec/is/section/interface/
H A Dsection.c3996 TYPE *a0 = (TYPE *)origArray, *a1; \
3998 PetscCall(PetscMalloc1(size, &a1)); \
4003 PetscCall(PetscMemcpy(&a1[o1], &a0[o0], n * unitsize)); \
4005 *newArray = (void *)a1; \
/petsc/src/dm/dt/interface/
H A Ddt.c1556 PetscReal a1, a6, gf; in PetscDTGaussJacobiQuadrature_Newton_Internal() local
1560 a1 = PetscPowReal(2.0, a + b + 1); in PetscDTGaussJacobiQuadrature_Newton_Internal()
1585 a6 = a1 * gf; in PetscDTGaussJacobiQuadrature_Newton_Internal()
/petsc/src/vec/vec/impls/seq/kokkos/
H A Dveckok.kokkos.cxx762 …MAXPYFunctor(PetscScalarKokkosView yv, PetscScalar a0, PetscScalar a1, PetscScalar a2, PetscScalar… in MAXPYFunctor()
766 a[1] = a1; in MAXPYFunctor()

12