| /petsc/src/ts/tutorials/ |
| H A D | ex43.c | 8 PetscReal u0, v0; /* initial conditions */ member
12 static void Exact(PetscReal t, PetscReal omega, PetscReal xi, PetscReal u0, PetscReal v0, PetscReal… in Exact() argument
18 PetscReal C1 = (v0 + a * u0) / w; in Exact()
19 PetscReal C2 = u0; in Exact()
24 PetscReal C1 = (w * u0 + xi * u0 + v0) / (2 * w); in Exact()
25 PetscReal C2 = (w * u0 - xi * u0 - v0) / (2 * w); in Exact()
30 PetscReal C1 = v0 + a * u0; in Exact()
31 PetscReal C2 = u0; in Exact()
46 Exact(t, user->Omega, user->Xi, user->u0, user->v0, &u, &v); in Solution()
169 …PetscCall(PetscOptionsReal("-initial_u", "Initial displacement", __FILE__, user.u0, &user.u0, NULL… in main()
[all …]
|
| H A D | ex11.h | 588 static PetscScalar s, c0, p0, r0, u0, w0, x0, x2, ri, cx, sgn0, wsp0, gasc1, gasc2, gasc3, gasc4; in riemannsolver() local
619 u0 = cvmgm_(uxl, uxr, &sgn0); in riemannsolver()
635 wsp0 = u0 + sgn0 * c0; in riemannsolver()
644 ri = u0 - sgn0 * 2. * gasc4 * c0; in riemannsolver()
656 *uxm = cvmgp_(uxm, &u0, &d__1); in riemannsolver()
|
| H A D | ex11.c | 1625 PetscReal p0, u0, wcp[3], wc[3]; in initLinearWave() local
1636 u0 = 0.; in initLinearWave()
1638 ueq[1] = u0; in initLinearWave()
|
| /petsc/src/ts/tests/ |
| H A D | ex81.c | 7 PetscReal u0; /* initial condition */ member
11 static void Exact(PetscReal t, PetscReal v0, PetscReal u0, PetscScalar *ut) in Exact() argument
13 if (ut) *ut = u0 + v0 * t; in Exact()
71 …PetscCall(PetscOptionsReal("-initial_u", "Initial displacement", __FILE__, user.u0, &user.u0, NULL… in main()
93 u[0] = user.u0; in main()
101 Exact(t_final, user.v0, user.u0, &u_exact); in main()
|
| H A D | ex80.c | 7 PetscReal u0, v0; /* initial conditions */ member
11 static void Exact(PetscReal t, PetscReal a0, PetscReal u0, PetscReal v0, PetscScalar *ut, PetscScal… in Exact() argument
13 if (ut) *ut = u0 + v0 * t + 0.5 * a0 * t * t; in Exact()
72 …PetscCall(PetscOptionsReal("-initial_u", "Initial displacement", __FILE__, user.u0, &user.u0, NULL… in main()
97 u[0] = user.u0; in main()
108 Exact(t_final, user.a0, user.u0, user.v0, &u_exact, &v_exact); in main()
|
| /petsc/src/dm/tests/ |
| H A D | ex30.c | 12 PetscReal alpha = 1, K = 1, rho0 = 1, u0 = 0, sigma = 0.2; in main() local
33 PetscCall(PetscOptionsReal("-u0", "Reference velocity", "", u0, &u0, NULL)); in main()
94 v[0] = -0.5 * u0; in main()
98 v[4] = 0.5 * u0; in main()
101 v[7] = -0.5 * u0; in main()
105 v[11] = 0.5 * u0; in main()
|
| /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
66 u0 = u[0]; 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
98 u0 = u[0]; 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
62 u0 = u[0]; 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
99 u0 = u[0]; 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/ts/utils/ |
| H A D | tsconvest.c | 53 Vec u, u0; in PetscConvEstGetConvRateTS_Temporal_Private() local
61 PetscCall(TSGetSolution(ts, &u0)); in PetscConvEstGetConvRateTS_Temporal_Private()
62 PetscCall(PetscObjectReference((PetscObject)u0)); in PetscConvEstGetConvRateTS_Temporal_Private()
109 PetscCall(TSSetSolution(ts, u0)); in PetscConvEstGetConvRateTS_Temporal_Private()
110 PetscCall(PetscConvEstComputeInitialGuess(ce, 0, NULL, u0)); in PetscConvEstGetConvRateTS_Temporal_Private()
111 PetscCall(VecDestroy(&u0)); in PetscConvEstGetConvRateTS_Temporal_Private()
|