Home
last modified time | relevance | path

Searched refs:u2 (Results 1 – 19 of 19) sorted by relevance

/petsc/src/binding/petsc4py/test/
H A Dtest_nsp.py17 u2 = PETSc.Vec().createSeq(3)
20 u2[0], u2[1], u2[2] = [2, -1, 0]
21 u2.normalize()
22 basis = [u1, u2]
/petsc/src/tao/complementarity/tutorials/
H A Dminsurf1.c430 PetscReal u1, u2, nf1, nf2, njac11, njac12, njac21, njac22; in MSA_BoundaryConditions() local
473 u2 = -yt; in MSA_BoundaryConditions()
475 nf1 = u1 + u1 * u2 * u2 - u1 * u1 * u1 / three - xt; in MSA_BoundaryConditions()
476 nf2 = -u2 - u1 * u1 * u2 + u2 * u2 * u2 / three - yt; in MSA_BoundaryConditions()
479 njac11 = one + u2 * u2 - u1 * u1; in MSA_BoundaryConditions()
480 njac12 = two * u1 * u2; in MSA_BoundaryConditions()
481 njac21 = -two * u1 * u2; in MSA_BoundaryConditions()
482 njac22 = -one - u1 * u1 + u2 * u2; in MSA_BoundaryConditions()
485 u2 = u2 - (njac11 * nf2 - njac21 * nf1) / det; in MSA_BoundaryConditions()
488 boundary[i] = u1 * u1 - u2 * u2; in MSA_BoundaryConditions()
/petsc/src/tao/unconstrained/tutorials/
H A Dminsurf1.c478 PetscReal u1, u2, nf1, nf2, njac11, njac12, njac21, njac22; in MSA_BoundaryConditions() local
521 u2 = -yt; in MSA_BoundaryConditions()
523 nf1 = u1 + u1 * u2 * u2 - u1 * u1 * u1 / three - xt; in MSA_BoundaryConditions()
524 nf2 = -u2 - u1 * u1 * u2 + u2 * u2 * u2 / three - yt; in MSA_BoundaryConditions()
527 njac11 = one + u2 * u2 - u1 * u1; in MSA_BoundaryConditions()
528 njac12 = two * u1 * u2; in MSA_BoundaryConditions()
529 njac21 = -two * u1 * u2; in MSA_BoundaryConditions()
530 njac22 = -one - u1 * u1 + u2 * u2; in MSA_BoundaryConditions()
533 u2 = u2 - (njac11 * nf2 - njac21 * nf1) / det; in MSA_BoundaryConditions()
536 boundary[i] = u1 * u1 - u2 * u2; in MSA_BoundaryConditions()
H A Dminsurf2.c695 PetscReal u1, u2, nf1, nf2, njac11, njac12, njac21, njac22; in MSA_BoundaryConditions() local
743 u2 = -yt; in MSA_BoundaryConditions()
745 nf1 = u1 + u1 * u2 * u2 - u1 * u1 * u1 / three - xt; in MSA_BoundaryConditions()
746 nf2 = -u2 - u1 * u1 * u2 + u2 * u2 * u2 / three - yt; in MSA_BoundaryConditions()
749 njac11 = one + u2 * u2 - u1 * u1; in MSA_BoundaryConditions()
750 njac12 = two * u1 * u2; in MSA_BoundaryConditions()
751 njac21 = -two * u1 * u2; in MSA_BoundaryConditions()
752 njac22 = -one - u1 * u1 + u2 * u2; in MSA_BoundaryConditions()
755 u2 = u2 - (njac11 * nf2 - njac21 * nf1) / det; in MSA_BoundaryConditions()
758 boundary[i] = u1 * u1 - u2 * u2; in MSA_BoundaryConditions()
/petsc/src/snes/tutorials/
H A Dex58.c468 PetscScalar u1, u2, nf1, nf2, njac11, njac12, njac21, njac22; in FormBoundaryConditions() local
523 u2 = -yt; in FormBoundaryConditions()
525 nf1 = u1 + u1 * u2 * u2 - u1 * u1 * u1 / three - xt; in FormBoundaryConditions()
526 nf2 = -u2 - u1 * u1 * u2 + u2 * u2 * u2 / three - yt; in FormBoundaryConditions()
529 njac11 = one + u2 * u2 - u1 * u1; in FormBoundaryConditions()
530 njac12 = two * u1 * u2; in FormBoundaryConditions()
531 njac21 = -two * u1 * u2; in FormBoundaryConditions()
532 njac22 = -one - u1 * u1 + u2 * u2; in FormBoundaryConditions()
535 u2 = u2 - (njac11 * nf2 - njac21 * nf1) / det; in FormBoundaryConditions()
538 boundary[i] = u1 * u1 - u2 * u2; in FormBoundaryConditions()
H A Dex4.c464 PetscScalar u1, u2, nf1, nf2, njac11, njac12, njac21, njac22; in FormBoundaryConditions_Enneper() local
513 u2 = -yt; in FormBoundaryConditions_Enneper()
515 nf1 = u1 + u1 * u2 * u2 - u1 * u1 * u1 / three - xt; in FormBoundaryConditions_Enneper()
516 nf2 = -u2 - u1 * u1 * u2 + u2 * u2 * u2 / three - yt; in FormBoundaryConditions_Enneper()
519 njac11 = one + u2 * u2 - u1 * u1; in FormBoundaryConditions_Enneper()
520 njac12 = two * u1 * u2; in FormBoundaryConditions_Enneper()
521 njac21 = -two * u1 * u2; in FormBoundaryConditions_Enneper()
522 njac22 = -one - u1 * u1 + u2 * u2; in FormBoundaryConditions_Enneper()
525 u2 = u2 - (njac11 * nf2 - njac21 * nf1) / det; in FormBoundaryConditions_Enneper()
528 boundary[i] = u1 * u1 - u2 * u2; in FormBoundaryConditions_Enneper()
H A Dex69.c198 PetscReal u1, u2, u3, u4, u5, u6; in SolKxSolution() local
518 …u2 = UU * _PC1 * t3 + UU * _PC2 * t6 - _PC1 * t6 * Rm + _PC2 * t3 * Rm - VV * _PC3 * t3 - VV * _P… in SolKxSolution()
637 u2 *= PetscSinReal(km * z); /* z velocity */ in SolKxSolution()
638 sum2 += u2; in SolKxSolution()
769 PetscReal u1, u2, u3, u4, u5, u6, Z, x = pos[0], z = pos[1]; in SolCxSolution() local
2841 …u2 = (-0.4e1 * t6 * t7 * t8 + 0.2e1 * t14 * t15 * t19 - 0.2e1 * t23 * t7 * t24 + 0.2e1 * t28 * t29… in SolCxSolution()
2938 …u5 = (-2 * Z * nz * PETSC_PI * u2 - u3 * 2 * nz * PETSC_PI) * PetscCosReal(nz * PETSC_PI * z); /* … in SolCxSolution()
2940 …u6 = (u3 * 2 * nz * PETSC_PI + 4 * Z * nz * PETSC_PI * u2) * PetscCosReal(nz * PETSC_PI * z); /* z… in SolCxSolution()
2946 u2 *= PetscSinReal(nz * PETSC_PI * z); /* z velocity */ in SolCxSolution()
2947 sum2 += u2; in SolCxSolution()
/petsc/src/mat/impls/sbaij/seq/
H A Dsbaijfact10.c13 MatScalar *u, *d, *w, *wp, u0, u1, u2, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() local
68 u2 = u[2]; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
140 m0 = uik[0] = -(d0 * u0 + d6 * u1 + d12 * u2 + d18 * u3 + d24 * u4 + d30 * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
141 m1 = uik[1] = -(d1 * u0 + d7 * u1 + d13 * u2 + d19 * u3 + d25 * u4 + d31 * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
142 m2 = uik[2] = -(d2 * u0 + d8 * u1 + d14 * u2 + d20 * u3 + d26 * u4 + d32 * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
143 m3 = uik[3] = -(d3 * u0 + d9 * u1 + d15 * u2 + d21 * u3 + d27 * u4 + d33 * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
144 m4 = uik[4] = -(d4 * u0 + d10 * u1 + d16 * u2 + d22 * u3 + d28 * u4 + d34 * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
145 m5 = uik[5] = -(d5 * u0 + d11 * u1 + d17 * u2 + d23 * u3 + d29 * u4 + d35 * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
183 dk[0] += m0 * u0 + m1 * u1 + m2 * u2 + m3 * u3 + m4 * u4 + m5 * u5; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
184 dk[1] += m6 * u0 + m7 * u1 + m8 * u2 + m9 * u3 + m10 * u4 + m11 * u5; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
[all …]
H A Dsbaijfact9.c12 MatScalar *u, *d, *w, *wp, u0, u1, u2, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12; in MatCholeskyFactorNumeric_SeqSBAIJ_6() local
100 u2 = u[2]; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
135 uik[0] = -(d[0] * u0 + d[6] * u1 + d[12] * u2 + d[18] * u3 + d[24] * u4 + d[30] * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
136 uik[1] = -(d[1] * u0 + d[7] * u1 + d[13] * u2 + d[19] * u3 + d[25] * u4 + d[31] * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
137 uik[2] = -(d[2] * u0 + d[8] * u1 + d[14] * u2 + d[20] * u3 + d[26] * u4 + d[32] * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
138 uik[3] = -(d[3] * u0 + d[9] * u1 + d[15] * u2 + d[21] * u3 + d[27] * u4 + d[33] * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
139 uik[4] = -(d[4] * u0 + d[10] * u1 + d[16] * u2 + d[22] * u3 + d[28] * u4 + d[34] * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
140 uik[5] = -(d[5] * u0 + d[11] * u1 + d[17] * u2 + d[23] * u3 + d[29] * u4 + d[35] * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
178 dk[0] += uik[0] * u0 + uik[1] * u1 + uik[2] * u2 + uik[3] * u3 + uik[4] * u4 + uik[5] * u5; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
179 dk[1] += uik[6] * u0 + uik[7] * u1 + uik[8] * u2 + uik[9] * u3 + uik[10] * u4 + uik[11] * u5; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
[all …]
H A Dsbaijfact12.c13 MatScalar *u, *d, *w, *wp, u0, u1, u2, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering() local
64 u2 = u[2]; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
112 …uik[0] = -(d[0] * u0 + d[7] * u1 + d[14] * u2 + d[21] * u3 + d[28] * u4 + d[35] * u5 + d[42] * u6); in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
113 …uik[1] = -(d[1] * u0 + d[8] * u1 + d[15] * u2 + d[22] * u3 + d[29] * u4 + d[36] * u5 + d[43] * u6); in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
114 …uik[2] = -(d[2] * u0 + d[9] * u1 + d[16] * u2 + d[23] * u3 + d[30] * u4 + d[37] * u5 + d[44] * u6); in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
115 …uik[3] = -(d[3] * u0 + d[10] * u1 + d[17] * u2 + d[24] * u3 + d[31] * u4 + d[38] * u5 + d[45] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
116 …uik[4] = -(d[4] * u0 + d[11] * u1 + d[18] * u2 + d[25] * u3 + d[32] * u4 + d[39] * u5 + d[46] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
117 …uik[5] = -(d[5] * u0 + d[12] * u1 + d[19] * u2 + d[26] * u3 + d[33] * u4 + d[40] * u5 + d[47] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
118 …uik[6] = -(d[6] * u0 + d[13] * u1 + d[20] * u2 + d[27] * u3 + d[34] * u4 + d[41] * u5 + d[48] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
169 …dk[0] += uik[0] * u0 + uik[1] * u1 + uik[2] * u2 + uik[3] * u3 + uik[4] * u4 + uik[5] * u5 + uik[6… in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
[all …]
H A Dsbaijfact11.c12 MatScalar *u, *d, *w, *wp, u0, u1, u2, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12; in MatCholeskyFactorNumeric_SeqSBAIJ_7() local
101 u2 = u[2]; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
149 …uik[0] = -(d[0] * u0 + d[7] * u1 + d[14] * u2 + d[21] * u3 + d[28] * u4 + d[35] * u5 + d[42] * u6); in MatCholeskyFactorNumeric_SeqSBAIJ_7()
150 …uik[1] = -(d[1] * u0 + d[8] * u1 + d[15] * u2 + d[22] * u3 + d[29] * u4 + d[36] * u5 + d[43] * u6); in MatCholeskyFactorNumeric_SeqSBAIJ_7()
151 …uik[2] = -(d[2] * u0 + d[9] * u1 + d[16] * u2 + d[23] * u3 + d[30] * u4 + d[37] * u5 + d[44] * u6); in MatCholeskyFactorNumeric_SeqSBAIJ_7()
152 …uik[3] = -(d[3] * u0 + d[10] * u1 + d[17] * u2 + d[24] * u3 + d[31] * u4 + d[38] * u5 + d[45] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7()
153 …uik[4] = -(d[4] * u0 + d[11] * u1 + d[18] * u2 + d[25] * u3 + d[32] * u4 + d[39] * u5 + d[46] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7()
154 …uik[5] = -(d[5] * u0 + d[12] * u1 + d[19] * u2 + d[26] * u3 + d[33] * u4 + d[40] * u5 + d[47] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7()
155 …uik[6] = -(d[6] * u0 + d[13] * u1 + d[20] * u2 + d[27] * u3 + d[34] * u4 + d[41] * u5 + d[48] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7()
206 …dk[0] += uik[0] * u0 + uik[1] * u1 + uik[2] * u2 + uik[3] * u3 + uik[4] * u4 + uik[5] * u5 + uik[6… in MatCholeskyFactorNumeric_SeqSBAIJ_7()
[all …]
/petsc/src/tao/bound/tutorials/
H A Dplate2f.F90520 PetscReal yt, hx, hy, u1, u2, nf1, nf2
587 u2 = -yt
592 nf1 = u1 + u1*u2*u2 - u1*u1*u1/three - xt
593 nf2 = -u2 - u1*u1*u2 + u2*u2*u2/three - yt
596 njac11 = one + u2*u2 - u1*u1
597 njac12 = two*u1*u2
598 njac21 = -two*u1*u2
599 njac22 = -one - u1*u1 + u2*u2
602 u2 = u2 - (njac11*nf2 - njac21*nf1)/det
609 boundary_v(1 + i) = u1*u1 - u2*u2
H A Dplate2.c611 PetscReal u1, u2, nf1, nf2, njac11, njac12, njac21, njac22; in MSA_BoundaryConditions() local
665 u2 = -yt; in MSA_BoundaryConditions()
667 nf1 = u1 + u1 * u2 * u2 - u1 * u1 * u1 / three - xt; in MSA_BoundaryConditions()
668 nf2 = -u2 - u1 * u1 * u2 + u2 * u2 * u2 / three - yt; in MSA_BoundaryConditions()
671 njac11 = one + u2 * u2 - u1 * u1; in MSA_BoundaryConditions()
672 njac12 = two * u1 * u2; in MSA_BoundaryConditions()
673 njac21 = -two * u1 * u2; in MSA_BoundaryConditions()
674 njac22 = -one - u1 * u1 + u2 * u2; in MSA_BoundaryConditions()
677 u2 = u2 - (njac11 * nf2 - njac21 * nf1) / det; in MSA_BoundaryConditions()
680 boundary[i] = u1 * u1 - u2 * u2; in MSA_BoundaryConditions()
/petsc/src/ts/tests/
H A Dex35.c122 PetscReal u1, u2, mag, zx, zy; in main() local
125 u2 = PetscRealPart(randVecNums[p * dim + 1]); in main()
128 x[p * dim + 1] = u2; in main()
132 zx = mag * PetscCosReal(2 * PETSC_PI * u2) + 0; in main()
133 zy = mag * PetscSinReal(2 * PETSC_PI * u2) + 0; in main()
H A Dex30.c597 PetscReal u1, u2; in go() local
602 PetscCall(PetscRandomGetValueReal(rand, &u2)); in go()
605 xx_t[grid][tid][pp] = mag * PetscCosReal(2.0 * PETSC_PI * u2); in go()
606 yy_t[grid][tid][pp] = mag * PetscSinReal(2.0 * PETSC_PI * u2); in go()
/petsc/src/sys/classes/random/tutorials/
H A Dex2.c98 PetscScalar u1, u2; in stdNormalArray() local
104 PetscCall(PetscRandomGetValue(ran, &u2)); in stdNormalArray()
107 eps[i] = t * PetscCosReal(2 * PETSC_PI * PetscRealPart(u2)); in stdNormalArray()
108 eps[i + 1] = t * PetscSinReal(2 * PETSC_PI * PetscRealPart(u2)); in stdNormalArray()
/petsc/src/vec/is/sf/impls/basic/
H A Dsfpack.c57 const Type *u = (const Type *)unpacked, *u2; \
66 u2 = u + opt->start[r] * MBS; \
71 PetscCall(PetscArraycpy(p2, u2 + (X * Y * k + X * j) * MBS, opt->dx[r] * MBS)); \
98 Type *u = (Type *)unpacked, *u2; \
109 u2 = u + opt->start[r] * MBS; \
114 PetscCall(PetscArraycpy(u2 + (X * Y * k + X * j) * MBS, p, opt->dx[r] * MBS)); \
139 Type *u = (Type *)unpacked, *u2; \
152 u2 = u + opt->start[r] * MBS; \
157 … for (i = 0; i < opt->dx[r] * MBS; i++) OpApply(Op, u2[(X * Y * k + X * j) * MBS + i], p[i]); \
/petsc/src/ksp/ksp/tutorials/
H A Dex43-solcx.h7 PetscReal Z, u1, u2, u3, u4, u5, u6, ZA, ZB; in evaluate_solCx() local
2133 …u2 = (-0.4e1 * t6 * t7 * t8 + 0.2e1 * t14 * t15 * t19 - 0.2e1 * t23 * t7 * t24 + 0.2e1 * t28 * t29… in evaluate_solCx()
2227 …u5 = (PetscReal)(-2 * Z * n * PETSC_PI * u2 - u3 * 2 * n * PETSC_PI) * PetscCosReal(n * PETSC_PI *… in evaluate_solCx()
2229 …u6 = (PetscReal)(u3 * 2 * n * PETSC_PI + 4 * Z * n * PETSC_PI * u2) * PetscCosReal(n * PETSC_PI * … in evaluate_solCx()
2235 u2 *= PetscSinReal(n * PETSC_PI * z); /* z velocity */ in evaluate_solCx()
2236 sum2 += u2; in evaluate_solCx()
/petsc/src/ts/utils/dmplexlandau/
H A Dplexland.c60 double u2 = 0, c02 = (double)*c2_0_arr, xx; in gamma_m1_f() local
64 for (PetscInt i = 0; i < dim; ++i) u2 += x[i] * x[i]; in gamma_m1_f()
66 xx = u2 / c02; in gamma_m1_f()