Home
last modified time | relevance | path

Searched refs:p2 (Results 1 – 25 of 29) 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/vec/is/is/utils/
H A Disdiff.c124 PetscInt n1, n2, n3, p1, p2, *iout; in ISSum() local
146 p2 = 0; in ISSum()
150 n3 += n2 - p2; in ISSum()
153 while (p2 < n2 && i2[p2] < i1[p1]) { in ISSum()
155 p2++; in ISSum()
157 if (p2 == n2) { in ISSum()
162 if (i2[p2] == i1[p1]) { in ISSum()
165 p2++; in ISSum()
169 if (p2 == n2) { in ISSum()
174 while (p1 < n1 && i1[p1] < i2[p2]) { in ISSum()
[all …]
/petsc/share/petsc/matlab/
H A Dlaplacian.m340 p2 = mod(p-p1,u(1)*u(2)); variable
341 p3 = (p - p2 - p1)/(u(1)*u(2)) + 1;
342 p2 = p2/u(1) + 1; variable
344 V2 = sin(kron((1:u(2))'*(pi/(u(2)+1)),p2))*(2/(u(2)+1))^0.5;
346 V2 = cos(kron((0.5:1:u(2)-0.5)'*(pi/u(2)),p2-1))*(2/u(2))^0.5;
347 V2(:,p2==1) = 1/u(2)^0.5;
349 V2 = sin(kron((1:u(2))'*(pi/2/(u(2)+0.5)),2*p2 - 1))...
352 V2 = cos(kron((0.5:1:u(2)-0.5)'*(pi/2/(u(2)+0.5)),2*p2 - 1))...
357 V2(:,mod(p2,2)==1) = cos(a*(pi/u(2)*(p2(mod(p2,2)==1)-1)))...
359 pp = p2(mod(p2,2)==0);
[all …]
/petsc/src/snes/tutorials/
H A Dex74ftube.h1 double precision p1, p2, p3, p4, r1, r2, r3, r4 double precision e1, e2, e3, e4, v1, v2, v3, v4 dou…
6 …p2, p3, p4, r1, r2, r3, r4, e1, e2, e3, e4, v1, v2, v3, v4, ru1, ru2, ru3, ru4, et1, et2, et3, et4…
H A Dex48.c1091 …p2 = _mm_set1_pd(2), p05 = _mm_set1_pd(0.5), p42 = _mm_setr_pd(4, 2), p24 = _mm_shuffle_pd(p42, p4… local
1092 p4dv1p2du0 = _mm_add_pd(_mm_mul_pd(p4, dv1), _mm_mul_pd(p2, du0)), /* 4 dv1 + 2 du0 */
/petsc/src/mat/impls/baij/seq/
H A Dbaijfact13.c20 MatScalar p1, p2, p3, p4, m1, m2, m3, m4, m5, m6, m7, m8, m9, x1, x2, x3, x4; in MatILUFactorNumeric_SeqBAIJ_3_inplace() local
65 p2 = pc[1]; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
73 …if (p1 != 0.0 || p2 != 0.0 || p3 != 0.0 || p4 != 0.0 || p5 != 0.0 || p6 != 0.0 || p7 != 0.0 || p8 … in MatILUFactorNumeric_SeqBAIJ_3_inplace()
86 pc[1] = m2 = p2 * x1 + p5 * x2 + p8 * x3; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
90 pc[4] = m5 = p2 * x4 + p5 * x5 + p8 * x6; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
94 pc[7] = m8 = p2 * x7 + p5 * x8 + p8 * x9; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
280 MatScalar p1, p2, p3, p4, m1, m2, m3, m4, m5, m6, m7, m8, m9, x1, x2, x3, x4; in MatILUFactorNumeric_SeqBAIJ_3_NaturalOrdering_inplace() local
322 p2 = pc[1]; in MatILUFactorNumeric_SeqBAIJ_3_NaturalOrdering_inplace()
330 …if (p1 != 0.0 || p2 != 0.0 || p3 != 0.0 || p4 != 0.0 || p5 != 0.0 || p6 != 0.0 || p7 != 0.0 || p8 … in MatILUFactorNumeric_SeqBAIJ_3_NaturalOrdering_inplace()
343 pc[1] = m2 = p2 * x1 + p5 * x2 + p8 * x3; in MatILUFactorNumeric_SeqBAIJ_3_NaturalOrdering_inplace()
[all …]
H A Dbaijfact11.c20 MatScalar p1, p2, p3, p4, m1, m2, m3, m4, m5, m6, m7, m8, m9, x1, x2, x3, x4; in MatILUFactorNumeric_SeqBAIJ_4_inplace() local
76 p2 = pc[1]; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
91 …if (p1 != 0.0 || p2 != 0.0 || p3 != 0.0 || p4 != 0.0 || p5 != 0.0 || p6 != 0.0 || p7 != 0.0 || p8 … in MatILUFactorNumeric_SeqBAIJ_4_inplace()
111 pc[1] = m2 = p2 * x1 + p6 * x2 + p10 * x3 + p14 * x4; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
116 pc[5] = m6 = p2 * x5 + p6 * x6 + p10 * x7 + p14 * x8; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
121 pc[9] = m10 = p2 * x9 + p6 * x10 + p10 * x11 + p14 * x12; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
126 pc[13] = m14 = p2 * x13 + p6 * x14 + p10 * x15 + p14 * x16; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
351 MatScalar p1, p2, p3, p4, m1, m2, m3, m4, m5, m6, m7, m8, m9, x1, x2, x3, x4; in MatILUFactorNumeric_SeqBAIJ_4_NaturalOrdering_inplace() local
404 p2 = pc[1]; in MatILUFactorNumeric_SeqBAIJ_4_NaturalOrdering_inplace()
419 …if (p1 != 0.0 || p2 != 0.0 || p3 != 0.0 || p4 != 0.0 || p5 != 0.0 || p6 != 0.0 || p7 != 0.0 || p8 … in MatILUFactorNumeric_SeqBAIJ_4_NaturalOrdering_inplace()
[all …]
H A Dbaijfact7.c17 MatScalar p1, p2, p3, p4, m1, m2, m3, m4, m5, m6, m7, m8, m9, x1, x2, x3, x4; in MatILUFactorNumeric_SeqBAIJ_6_inplace() local
99 p2 = pc[1]; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
134 …if (p1 != 0.0 || p2 != 0.0 || p3 != 0.0 || p4 != 0.0 || p5 != 0.0 || p6 != 0.0 || p7 != 0.0 || p8 … in MatILUFactorNumeric_SeqBAIJ_6_inplace()
174 pc[1] = m2 = p2 * x1 + p8 * x2 + p14 * x3 + p20 * x4 + p26 * x5 + p32 * x6; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
181 pc[7] = m8 = p2 * x7 + p8 * x8 + p14 * x9 + p20 * x10 + p26 * x11 + p32 * x12; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
188 pc[13] = m14 = p2 * x13 + p8 * x14 + p14 * x15 + p20 * x16 + p26 * x17 + p32 * x18; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
195 pc[19] = m20 = p2 * x19 + p8 * x20 + p14 * x21 + p20 * x22 + p26 * x23 + p32 * x24; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
202 pc[25] = m26 = p2 * x25 + p8 * x26 + p14 * x27 + p20 * x28 + p26 * x29 + p32 * x30; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
209 pc[31] = m32 = p2 * x31 + p8 * x32 + p14 * x33 + p20 * x34 + p26 * x35 + p32 * x36; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
480 MatScalar p1, p2, p3, p4, p5, p6, p7, p8, p9, p10, p11, p12, p13, p14, p15; in MatILUFactorNumeric_SeqBAIJ_6_NaturalOrdering_inplace() local
[all …]
H A Dbaijfact9.c21 MatScalar p1, p2, p3, p4, m1, m2, m3, m4, m5, m6, m7, m8, m9, x1, x2, x3, x4; in MatILUFactorNumeric_SeqBAIJ_5_inplace() local
98 p2 = pc[1]; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
122 …if (p1 != 0.0 || p2 != 0.0 || p3 != 0.0 || p4 != 0.0 || p5 != 0.0 || p6 != 0.0 || p7 != 0.0 || p8 … in MatILUFactorNumeric_SeqBAIJ_5_inplace()
151 pc[1] = m2 = p2 * x1 + p7 * x2 + p12 * x3 + p17 * x4 + p22 * x5; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
157 pc[6] = m7 = p2 * x6 + p7 * x7 + p12 * x8 + p17 * x9 + p22 * x10; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
163 pc[11] = m12 = p2 * x11 + p7 * x12 + p12 * x13 + p17 * x14 + p22 * x15; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
169 pc[16] = m17 = p2 * x16 + p7 * x17 + p12 * x18 + p17 * x19 + p22 * x20; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
175 pc[21] = m22 = p2 * x21 + p7 * x22 + p12 * x23 + p17 * x24 + p22 * x25; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
424 MatScalar p1, p2, p3, p4, p5, p6, p7, p8, p9, p10, p11, p12, p13, p14, p15; in MatILUFactorNumeric_SeqBAIJ_5_NaturalOrdering_inplace() local
481 p2 = pc[1]; in MatILUFactorNumeric_SeqBAIJ_5_NaturalOrdering_inplace()
[all …]
H A Dbaijfact5.c17 MatScalar p1, p2, p3, p4, m1, m2, m3, m4, m5, m6, m7, m8, m9, x1, x2, x3, x4; in MatILUFactorNumeric_SeqBAIJ_7_inplace() local
116 p2 = pc[1]; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
164 …if (p1 != 0.0 || p2 != 0.0 || p3 != 0.0 || p4 != 0.0 || p5 != 0.0 || p6 != 0.0 || p7 != 0.0 || p8 … in MatILUFactorNumeric_SeqBAIJ_7_inplace()
217 pc[1] = m2 = p2 * x1 + p9 * x2 + p16 * x3 + p23 * x4 + p30 * x5 + p37 * x6 + p44 * x7; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
225 pc[8] = m9 = p2 * x8 + p9 * x9 + p16 * x10 + p23 * x11 + p30 * x12 + p37 * x13 + p44 * x14; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
233 … pc[15] = m16 = p2 * x15 + p9 * x16 + p16 * x17 + p23 * x18 + p30 * x19 + p37 * x20 + p44 * x21; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
241 … pc[22] = m23 = p2 * x22 + p9 * x23 + p16 * x24 + p23 * x25 + p30 * x26 + p37 * x27 + p44 * x28; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
249 … pc[29] = m30 = p2 * x29 + p9 * x30 + p16 * x31 + p23 * x32 + p30 * x33 + p37 * x34 + p44 * x35; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
257 … pc[36] = m37 = p2 * x36 + p9 * x37 + p16 * x38 + p23 * x39 + p30 * x40 + p37 * x41 + p44 * x42; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
265 … pc[43] = m44 = p2 * x43 + p9 * x44 + p16 * x45 + p23 * x46 + p30 * x47 + p37 * x48 + p44 * x49; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
[all …]
H A Dbaijfact.c234 MatScalar p1, p2, p3, p4; in MatILUFactorNumeric_SeqBAIJ_2_inplace() local
273 p2 = pc[1]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
276 …if (p1 != (PetscScalar)0.0 || p2 != (PetscScalar)0.0 || p3 != (PetscScalar)0.0 || p4 != (PetscScal… in MatILUFactorNumeric_SeqBAIJ_2_inplace()
284 pc[1] = m2 = p2 * x1 + p4 * x2; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
286 pc[3] = m4 = p2 * x3 + p4 * x4; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
345 MatScalar p1, p2, p3, p4, m1, m2, m3, m4, x1, x2, x3, x4; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace() local
380 p2 = pc[1]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
383 …if (p1 != (PetscScalar)0.0 || p2 != (PetscScalar)0.0 || p3 != (PetscScalar)0.0 || p4 != (PetscScal… in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
391 pc[1] = m2 = p2 * x1 + p4 * x2; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
393 pc[3] = m4 = p2 * x3 + p4 * x4; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
/petsc/src/dm/impls/swarm/tests/output/
H A Dex6_bsi_2d_multiple_2.out28 t 0.1: p2 error [0.011 0.39] sol [(0.827175 -0.561339) (17.735719 26.180494)] exact [(0.821001 -0.5…
52 t 0.1: p2 error [0.0018 0.057] sol [(0.821827 -0.569355) (18.009015 25.997706)] exact [(0.821001 -0…
76 t 0.1: p2 error [0.00067 0.024] sol [(0.820503 -0.571367) (18.076299 25.953113)] exact [(0.821001 -…
H A Dex6_bsi_2d_multiple_3.out28 t 0.01: p2 error [3.6e-05 0.0073] sol [(-0.951726 0.307065) (-9.716591 -30.092732)] exact [(-0.9516…
52 t 0.01: p2 error [3.6e-05 0.0073] sol [(-0.951726 0.307064) (-9.716664 -30.092823)] exact [(-0.9516…
76 t 0.01: p2 error [3.6e-05 0.0073] sol [(-0.951726 0.307064) (-9.716699 -30.092866)] exact [(-0.9516…
H A Dex6_im_2d_multiple_0.out28 t 0.1: p2 error [0.027 0.87] sol [(0.804764 -0.593029) (18.759203 25.457908)] exact [(0.821001 -0.5…
52 t 0.1: p2 error [0.008 0.26] sol [(0.816216 -0.577366) (18.266044 25.817982)] exact [(0.821001 -0.5…
76 t 0.1: p2 error [0.0031 0.1] sol [(0.819099 -0.573375) (18.140635 25.908236)] exact [(0.821001 -0.5…
H A Dex6_bsi_2d_multiple_1.out28 t 0.1: p2 error [0.014 0.42] sol [(0.823708 -0.557243) (17.844321 26.328908)] exact [(0.821001 -0.5…
52 t 0.1: p2 error [0.0029 0.088] sol [(0.820572 -0.568044) (18.047976 26.050124)] exact [(0.821001 -0…
76 t 0.1: p2 error [0.0011 0.042] sol [(0.819927 -0.570810) (18.093911 25.976577)] exact [(0.821001 -0…
H A Dex5_bsi_2d_multiple_3.out45 t 0.1: p2 error [1.8e-08 5.8e-07] sol [(0.821001 -0.570926) (18.054271 25.962343)] exact [(0.821001…
106 t 0.1: p2 error [2.3e-09 7.2e-08] sol [(0.821001 -0.570926) (18.054271 25.962344)] exact [(0.821001…
207 t 0.1: p2 error [2.8e-10 9e-09] sol [(0.821001 -0.570926) (18.054271 25.962344)] exact [(0.821001 -…
H A Dex5_im_2d_multiple_0.out45 t 0.001: p2 error [2.6e-10 8.3e-09] sol [(-0.827197 0.561912) (-17.769227 -26.158260)] exact [(-0.8…
106 t 0.001: p2 error [6.6e-11 2.1e-09] sol [(-0.827197 0.561912) (-17.769227 -26.158260)] exact [(-0.8…
207 t 0.001: p2 error [1.6e-11 5.2e-10] sol [(-0.827197 0.561912) (-17.769227 -26.158260)] exact [(-0.8…
H A Dex5_bsi_2d_multiple_2.out45 t 0.001: p2 error [2.6e-10 1.7e-08] sol [(-0.827197 0.561912) (-17.769227 -26.158260)] exact [(-0.8…
106 t 0.001: p2 error [6.6e-11 4.2e-09] sol [(-0.827197 0.561912) (-17.769227 -26.158260)] exact [(-0.8…
207 t 0.001: p2 error [1.6e-11 1e-09] sol [(-0.827197 0.561912) (-17.769227 -26.158260)] exact [(-0.827…
H A Dex5_bsi_2d_multiple_1.out45 t 0.001: p2 error [5e-06 0.00016] sol [(-0.827193 0.561909) (-17.769320 -26.158388)] exact [(-0.827…
106 t 0.001: p2 error [2.5e-06 7.9e-05] sol [(-0.827195 0.561911) (-17.769273 -26.158324)] exact [(-0.8…
207 t 0.001: p2 error [1.2e-06 4e-05] sol [(-0.827196 0.561912) (-17.769250 -26.158292)] exact [(-0.827…
/petsc/src/snes/utils/dm/
H A Ddminterpolatesnes.c680 const PetscScalar p2 = ref[2]; in HexMap_Private() local
682 …0] = x0 + f_1 * p0 + f_3 * p1 + f_4 * p2 + f_01 * p0 * p1 + f_12 * p1 * p2 + f_02 * p0 * p2 + f_01… in HexMap_Private()
683 …+ g_1 * p0 + g_3 * p1 + g_4 * p2 + g_01 * p0 * p1 + g_01 * p0 * p1 + g_12 * p1 * p2 + g_02 * p0 * … in HexMap_Private()
684 …+ h_1 * p0 + h_3 * p1 + h_4 * p2 + h_01 * p0 * p1 + h_01 * p0 * p1 + h_12 * p1 * p2 + h_02 * p0 * … in HexMap_Private()
/petsc/src/ts/tutorials/power_grid/
H A Dex8.c188 …p, PetscScalar x, PetscScalar y, PetscInt i, PetscInt j, PetscInt N, PetscScalar *p2, AppCtx *user) in adv2() argument
197 *p2 = fpos * (p[j][i] - p[j - 1][i]) / user->dy + fneg * (p[j + 1][i] - p[j][i]) / user->dy; in adv2()
199 …*p2 = fpos * (3 * p[j][i] - 4 * p[j - 1][i] + p[j - 2][i]) / (2 * user->dy) + fneg * (-p[j + 2][i]… in adv2()
201 …*p2 = fpos * (2 * p[j + 1][i] + 3 * p[j][i] - 6 * p[j - 1][i] + p[j - 2][i]) / (6 * user->dy) + fn… in adv2()
H A Dex7.c218 …(PetscScalar **p, PetscScalar x, PetscInt i, PetscInt j, PetscInt N, PetscScalar *p2, AppCtx *user) in adv2() argument
227 *p2 = fpos * (p[j][i] - p[j - 1][i]) / user->dy + fneg * (p[j + 1][i] - p[j][i]) / user->dy; in adv2()
229 …*p2 = fpos * (3 * p[j][i] - 4 * p[j - 1][i] + p[j - 2][i]) / (2 * user->dy) + fneg * (-p[j + 2][i]… in adv2()
231 …*p2 = fpos * (2 * p[j + 1][i] + 3 * p[j][i] - 6 * p[j - 1][i] + p[j - 2][i]) / (6 * user->dy) + fn… in adv2()
H A Dex6.c188 …(PetscScalar **p, PetscScalar x, PetscInt i, PetscInt j, PetscInt N, PetscScalar *p2, AppCtx *user) in adv2() argument
207 *p2 = f * (p[j + 1][i] - p[j - 1][i]) / (2 * user->dy); in adv2()
/petsc/src/mat/interface/
H A Dmatproduct.c1224 MatProductType p1, p2; in MatProductSymbolic_ABC_Basic() local
1238 p2 = MATPRODUCT_AtB; in MatProductSymbolic_ABC_Basic()
1246 p2 = MATPRODUCT_AB; in MatProductSymbolic_ABC_Basic()
1254 p2 = MATPRODUCT_AB; in MatProductSymbolic_ABC_Basic()
1278 PetscCall(MatProductSetType(mmabc->ABC, p2)); in MatProductSymbolic_ABC_Basic()
1283 …P, "Symbolic ProductType %s not supported with %s and %s", MatProductTypes[p2], ((PetscObject)A)->… in MatProductSymbolic_ABC_Basic()
/petsc/src/vec/is/sf/impls/basic/
H A Dsfpack.c58 Type *p = (Type *)packed, *p2; \
64 else if (opt) { /* has optimizations available */ p2 = p; \
71 PetscCall(PetscArraycpy(p2, u2 + (X * Y * k + X * j) * MBS, opt->dx[r] * MBS)); \
72 p2 += opt->dx[r] * MBS; \

12