| /petsc/include/petsc/private/kernels/ |
| H A D | petscaxpy.h | 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); \ 100 #define PetscKernelAXPY3(U, a1, a2, a3, p1, p2, p3, n) \ argument 104 *U++ += a1 * *p1++ + a2 * *p2++ + a3 * *p3++; \ 106 *U++ += a1 * *p1++ + a2 * *p2++ + a3 * *p3++; \ 108 *U++ += a1 * *p1++ + a2 * *p2++ + a3 * *p3++; \ 114 U[0] += a1 * p1[0] + a2 * p2[0] + a3 * p3[0]; \ 115 U[1] += a1 * p1[1] + a2 * p2[1] + a3 * p3[1]; \ [all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex121.c | 23 PetscScalar *a, *a2, *a3; in main() local 97 PetscCall(VecGetArray(z2, &a3)); in main() 101 a3[i] = 0.0; in main() 106 a3[i] += a[xpInd] * a2[diffInd]; in main() 111 PetscCall(VecRestoreArray(z2, &a3)); in main()
|
| /petsc/src/snes/tutorials/ |
| H A D | ex74ftube.h | 1 …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/vec/vec/impls/seq/ftn-kernels/ |
| H A D | fmaxpy.F90 | 11 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))
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijfact9.c | 886 __m256d a0, a1, a2, a3, a4, a5, w0, w1, w2, w3, s0, s1, s2, v0, v1, v2, v3; in MatSolve_SeqBAIJ_9_NaturalOrdering() local 918 a3 = _mm256_loadu_pd(&v[9]); in MatSolve_SeqBAIJ_9_NaturalOrdering() 919 s0 = _mm256_fnmadd_pd(a3, w1, s0); in MatSolve_SeqBAIJ_9_NaturalOrdering() 934 a3 = _mm256_loadu_pd(&v[27]); in MatSolve_SeqBAIJ_9_NaturalOrdering() 935 s0 = _mm256_fnmadd_pd(a3, w3, s0); in MatSolve_SeqBAIJ_9_NaturalOrdering() 950 a3 = _mm256_loadu_pd(&v[45]); in MatSolve_SeqBAIJ_9_NaturalOrdering() 951 s0 = _mm256_fnmadd_pd(a3, w1, s0); in MatSolve_SeqBAIJ_9_NaturalOrdering() 966 a3 = _mm256_loadu_pd(&v[63]); in MatSolve_SeqBAIJ_9_NaturalOrdering() 967 s0 = _mm256_fnmadd_pd(a3, w3, s0); in MatSolve_SeqBAIJ_9_NaturalOrdering() 1010 a3 = _mm256_loadu_pd(&v[9]); in MatSolve_SeqBAIJ_9_NaturalOrdering() [all …]
|
| H A D | baij2.c | 679 __m256d a0, a1, a2, a3, a4, a5; in MatMult_SeqBAIJ_9_AVX2() local 743 a3 = _mm256_loadu_pd(&v[j * 81 + 18]); in MatMult_SeqBAIJ_9_AVX2() 744 z0 = _mm256_fmadd_pd(a3, w2, z0); in MatMult_SeqBAIJ_9_AVX2() 761 a3 = _mm256_loadu_pd(&v[j * 81 + 36]); in MatMult_SeqBAIJ_9_AVX2() 762 z0 = _mm256_fmadd_pd(a3, w0, z0); in MatMult_SeqBAIJ_9_AVX2() 788 a3 = _mm256_loadu_pd(&v[j * 81 + 63]); in MatMult_SeqBAIJ_9_AVX2() 789 z0 = _mm256_fmadd_pd(a3, w3, z0); in MatMult_SeqBAIJ_9_AVX2() 1333 __m256d a0, a1, a2, a3, a4, a5; in MatMult_SeqBAIJ_12_AVX2() local 1373 a3 = _mm256_loadu_pd(v + 12); in MatMult_SeqBAIJ_12_AVX2() 1374 z0 = _mm256_fmadd_pd(a3, w1, z0); in MatMult_SeqBAIJ_12_AVX2() [all …]
|
| /petsc/src/ts/impls/rosw/ |
| H A D | rosw.c | 961 PetscErrorCode TSRosWRegisterRos4(TSRosWType name, PetscReal gamma, PetscReal a2, PetscReal a3, Pet… in TSRosWRegisterRos4() argument 972 …if (a3 == (PetscReal)PETSC_DEFAULT || a3 == (PetscReal)PETSC_DETERMINE) a3 = (one / five - a2 / fo… in TSRosWRegisterRos4() 973 …a4 = a3; … in TSRosWRegisterRos4() 986 rhs[1] = one / three - a3 * a3 * b3; in TSRosWRegisterRos4() 987 rhs[2] = one / four - a3 * a3 * a3 * b3; in TSRosWRegisterRos4() 994 beta43 = (p56 - a2 * p43) / (b4 * a3 * a3 * (a3 - a2)); /* 7.21 */ in TSRosWRegisterRos4() 1000 M[1][0] = a4 * a4 * beta32beta2p - a3 * a3 * beta4jbetajp; in TSRosWRegisterRos4() 1003 M[2][0] = b4 * beta43 * a3 * a3 - p43; in TSRosWRegisterRos4() 1020 a32 = p42 / (b3 * a3 * beta2p + b4 * a4 * beta2p); in TSRosWRegisterRos4() 1031 A[2][0] = a3 - a32; in TSRosWRegisterRos4() [all …]
|
| /petsc/src/ts/adapt/impls/dsp/ |
| H A D | adaptdsp.c | 125 PetscReal a3 = dsp->Alpha[1]; in TSAdaptChoose_DSP() local 141 rho0 *= PetscPowReal(rho2, -a3); 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/binding/petsc4py/test/ |
| H A D | test_vec.py | 331 a3 = v3.getDict()['__array__'] 332 self.assertTrue(a is a3) 349 a3 = v3.getDict()['__array__'] 350 self.assertTrue(a is a3)
|
| /petsc/share/petsc/matlab/ |
| H A D | laplacian.m | 268 a3 = pi/2/(u(3)+1); variable 271 a3 = pi/2/u(3); variable 275 a3 = pi/4/(u(3)+0.5); variable 278 a3 = pi/u(3); variable 281 lambda3 = 4*sin(a3*N).^2;
|
| /petsc/include/ |
| H A D | petscmacros.h | 1203 #define PETSC_SELECT_16TH(a1, a2, a3, a4, a5, a6, a7, a8, a9, a10, a11, a12, a13, a14, a15, a16, ..… argument
|
| /petsc/src/dm/dt/interface/ |
| H A D | dt.c | 1563 PetscReal a2, a3, a4, a5; in PetscDTGaussJacobiQuadrature_Newton_Internal() local 1565 a3 = PetscLGamma(b + npoints + 1); in PetscDTGaussJacobiQuadrature_Newton_Internal() 1568 gf = PetscExpReal(a2 + a3 - (a4 + a5)); in PetscDTGaussJacobiQuadrature_Newton_Internal()
|
| /petsc/src/vec/vec/impls/seq/kokkos/ |
| H A D | veckok.kokkos.cxx | 762 …KokkosView yv, PetscScalar a0, PetscScalar a1, PetscScalar a2, PetscScalar a3, PetscScalar a4, Pet… in MAXPYFunctor() 768 a[3] = a3; in MAXPYFunctor()
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | testcase3D.cas | 7982 (wj-bsl-earsm/a3 1.8) 9033 1 a0 a3 2c 1 9 9115 a0 96 97 a3 1 0 9120 a0 a3 a2 9f 9 0 9173 2c 2b a2 a3 9 0 9184 2c a3 97 93 1 0
|