| /petsc/src/binding/petsc4py/test/ |
| H A D | test_nsp.py | 16 u1 = PETSc.Vec().createSeq(3)
18 u1[0], u1[1], u1[2] = [1, 2, 0]
19 u1.normalize()
22 basis = [u1, u2]
|
| /petsc/src/tao/complementarity/tutorials/ |
| H A D | minsurf1.c | 430 PetscReal u1, u2, nf1, nf2, njac11, njac12, njac21, njac22; in MSA_BoundaryConditions() local
472 u1 = xt; 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()
484 u1 = u1 - (njac22 * nf1 - njac12 * nf2) / det; in MSA_BoundaryConditions()
488 boundary[i] = u1 * u1 - u2 * u2; in MSA_BoundaryConditions()
|
| /petsc/src/tao/unconstrained/tutorials/ |
| H A D | minsurf1.c | 478 PetscReal u1, u2, nf1, nf2, njac11, njac12, njac21, njac22; in MSA_BoundaryConditions() local
520 u1 = xt; 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()
532 u1 = u1 - (njac22 * nf1 - njac12 * nf2) / det; in MSA_BoundaryConditions()
536 boundary[i] = u1 * u1 - u2 * u2; in MSA_BoundaryConditions()
|
| H A D | minsurf2.c | 695 PetscReal u1, u2, nf1, nf2, njac11, njac12, njac21, njac22; in MSA_BoundaryConditions() local
742 u1 = xt; 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()
754 u1 = u1 - (njac22 * nf1 - njac12 * nf2) / det; in MSA_BoundaryConditions()
758 boundary[i] = u1 * u1 - u2 * u2; in MSA_BoundaryConditions()
|
| /petsc/src/snes/tutorials/ |
| H A D | ex58.c | 468 PetscScalar u1, u2, nf1, nf2, njac11, njac12, njac21, njac22; in FormBoundaryConditions() local
522 u1 = xt; 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()
534 u1 = u1 - (njac22 * nf1 - njac12 * nf2) / det; in FormBoundaryConditions()
538 boundary[i] = u1 * u1 - u2 * u2; in FormBoundaryConditions()
|
| H A D | ex4.c | 464 PetscScalar u1, u2, nf1, nf2, njac11, njac12, njac21, njac22; in FormBoundaryConditions_Enneper() local
512 u1 = xt; 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()
524 u1 = u1 - (njac22 * nf1 - njac12 * nf2) / det; in FormBoundaryConditions_Enneper()
528 boundary[i] = u1 * u1 - u2 * u2; in FormBoundaryConditions_Enneper()
|
| H A D | ex69.c | 198 PetscReal u1, u2, u3, u4, u5, u6; in SolKxSolution() local
508 u1 = -km * (_PC1 * t2 + _PC2 * t4 + _PC3 * t2 + _PC4 * t4 + t10 * AA * t13 + t10 * BB * t16); in SolKxSolution()
635 u1 *= PetscCosReal(km * z); /* x velocity */ in SolKxSolution()
636 sum1 += u1; in SolKxSolution()
769 PetscReal u1, u2, u3, u4, u5, u6, Z, x = pos[0], z = pos[1]; in SolCxSolution() local
2800 …u1 = (t55 + t90 + t120 + t155) / (0.4e1 * t158 * t19 * t2 + 0.8e1 * t158 * t19 * t1 * t13 + 0.4e1 … in SolCxSolution()
2944 u1 *= PetscCosReal(nz * PETSC_PI * z); /* x velocity */ in SolCxSolution()
2945 sum1 += u1; in SolCxSolution()
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaijfact10.c | 13 MatScalar *u, *d, *w, *wp, u0, u1, u2, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() local
67 u1 = u[1]; 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 D | sbaijfact9.c | 12 MatScalar *u, *d, *w, *wp, u0, u1, u2, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12; in MatCholeskyFactorNumeric_SeqSBAIJ_6() local
99 u1 = u[1]; 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 D | sbaijfact12.c | 13 MatScalar *u, *d, *w, *wp, u0, u1, u2, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering() local
63 u1 = u[1]; 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 D | sbaijfact11.c | 12 MatScalar *u, *d, *w, *wp, u0, u1, u2, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12; in MatCholeskyFactorNumeric_SeqSBAIJ_7() local
100 u1 = u[1]; 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 D | plate2f.F90 | 520 PetscReal yt, hx, hy, u1, u2, nf1, nf2
586 u1 = xt
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
601 u1 = u1 - (njac22*nf1 - njac12*nf2)/det
609 boundary_v(1 + i) = u1*u1 - u2*u2
|
| H A D | plate2.c | 611 PetscReal u1, u2, nf1, nf2, njac11, njac12, njac21, njac22; in MSA_BoundaryConditions() local
664 u1 = xt; 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()
676 u1 = u1 - (njac22 * nf1 - njac12 * nf2) / det; in MSA_BoundaryConditions()
680 boundary[i] = u1 * u1 - u2 * u2; in MSA_BoundaryConditions()
|
| /petsc/src/ts/tests/ |
| H A D | ex35.c | 122 PetscReal u1, u2, mag, zx, zy; in main() local
124 u1 = PetscRealPart(randVecNums[p * dim]); in main()
127 x[p * dim] = u1; in main()
130 mag = PetscSqrtReal(-2.0 * PetscLogReal(u1)); in main()
|
| H A D | ex30.c | 597 PetscReal u1, u2; in go() local
600 PetscCall(PetscRandomGetValueReal(rand, &u1)); in go()
601 } while (u1 == 0); in go()
604 …PetscReal mag = sigma * PetscSqrtReal(-2.0 * PetscLogReal(u1)); // is this the same scale gr… in go()
|
| /petsc/src/sys/classes/random/tutorials/ |
| H A D | ex2.c | 98 PetscScalar u1, u2; in stdNormalArray() local
103 PetscCall(PetscRandomGetValue(ran, &u1)); in stdNormalArray()
106 t = PetscSqrtReal(-2 * PetscLogReal(PetscRealPart(u1))); in stdNormalArray()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex43-solcx.h | 7 PetscReal Z, u1, u2, u3, u4, u5, u6, ZA, ZB; in evaluate_solCx() local
2093 …u1 = (t55 + t90 + t120 + t155) / (0.4e1 * t158 * t19 * t2 + 0.8e1 * t158 * t19 * t1 * t13 + 0.4e1 … in evaluate_solCx()
2233 u1 *= PetscCosReal(n * PETSC_PI * z); /* x velocity */ in evaluate_solCx()
2234 sum1 += u1; in evaluate_solCx()
|