Lines Matching refs:i
42 PetscInt i; in Limit_Upwind() local
43 for (i = 0; i < info->m; i++) lmt[i] = 0; in Limit_Upwind()
47 PetscInt i; in Limit_LaxWendroff() local
48 for (i = 0; i < info->m; i++) lmt[i] = jR[i]; in Limit_LaxWendroff()
52 PetscInt i; in Limit_BeamWarming() local
53 for (i = 0; i < info->m; i++) lmt[i] = jL[i]; in Limit_BeamWarming()
57 PetscInt i; in Limit_Fromm() local
58 for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] + jR[i]); in Limit_Fromm()
62 PetscInt i; in Limit_Minmod() local
63 for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i], jR[i]); in Limit_Minmod()
67 PetscInt i; in Limit_Superbee() local
68 …for (i = 0; i < info->m; i++) lmt[i] = MaxMod2(MinMod2(jL[i], 2 * jR[i]), MinMod2(2 * jL[i], jR[i]… in Limit_Superbee()
72 PetscInt i; in Limit_MC() local
73 for (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i], 0.5 * (jL[i] + jR[i]), 2 * jR[i]); in Limit_MC()
77 PetscInt i; in Limit_VanLeer() local
78 …for (i = 0; i < info->m; i++) lmt[i] = (jL[i] * Abs(jR[i]) + Abs(jL[i]) * jR[i]) / (Abs(jL[i]) + A… in Limit_VanLeer()
82 PetscInt i; in Limit_VanAlbada() local
83 …for (i = 0; i < info->m; i++) lmt[i] = (jL[i] * Sqr(jR[i]) + Sqr(jL[i]) * jR[i]) / (Sqr(jL[i]) + S… in Limit_VanAlbada()
87 PetscInt i; in Limit_VanAlbadaTVD() local
88 …for (i = 0; i < info->m; i++) lmt[i] = (jL[i] * jR[i] < 0) ? 0 : (jL[i] * Sqr(jR[i]) + Sqr(jL[i]) … in Limit_VanAlbadaTVD()
92 PetscInt i; in Limit_Koren() local
93 …for (i = 0; i < info->m; i++) lmt[i] = ((jL[i] * Sqr(jR[i]) + 2 * Sqr(jL[i]) * jR[i]) / (2 * Sqr(j… in Limit_Koren()
97 PetscInt i; in Limit_KorenSym() local
98 …for (i = 0; i < info->m; i++) lmt[i] = (1.5 * (jL[i] * Sqr(jR[i]) + Sqr(jL[i]) * jR[i]) / (2 * Sqr… in Limit_KorenSym()
102 PetscInt i; in Limit_Koren3() local
103 for (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i], (jL[i] + 2 * jR[i]) / 3, 2 * jR[i]); in Limit_Koren3()
111 PetscInt i; in Limit_CadaTorrilhon2() local
112 for (i = 0; i < info->m; i++) lmt[i] = CadaTorrilhonPhiHatR_Eq13(jL[i], jR[i]); in Limit_CadaTorrilhon2()
118 PetscInt i; in Limit_CadaTorrilhon3R() local
119 for (i = 0; i < info->m; i++) { in Limit_CadaTorrilhon3R()
120 const PetscReal eta = (Sqr(jL[i]) + Sqr(jR[i])) / Sqr(r * hx); in Limit_CadaTorrilhon3R()
121 …i] = ((eta < 1 - eps) ? (jL[i] + 2 * jR[i]) / 3 : ((eta > 1 + eps) ? CadaTorrilhonPhiHatR_Eq13(jL[… in Limit_CadaTorrilhon3R()
145 PetscInt i; in Limit2_Upwind() local
146 for (i = 0; i < info->m; i++) lmt[i] = 0; in Limit2_Upwind()
150 PetscInt i; in Limit2_LaxWendroff() local
152 for (i = 0; i < info->m; i++) lmt[i] = jR[i] / info->hxs; in Limit2_LaxWendroff()
154 for (i = 0; i < info->m; i++) lmt[i] = jR[i] / (info->hxs / 2.0 + info->hxf / 2.0); in Limit2_LaxWendroff()
156 for (i = 0; i < info->m; i++) lmt[i] = jR[i] / info->hxf; in Limit2_LaxWendroff()
158 for (i = 0; i < info->m; i++) lmt[i] = jR[i] / info->hxf; in Limit2_LaxWendroff()
160 for (i = 0; i < info->m; i++) lmt[i] = jR[i] / (info->hxf / 2.0 + info->hxs / 2.0); in Limit2_LaxWendroff()
162 for (i = 0; i < info->m; i++) lmt[i] = jR[i] / info->hxs; in Limit2_LaxWendroff()
164 for (i = 0; i < info->m; i++) lmt[i] = jR[i] / info->hxs; in Limit2_LaxWendroff()
169 PetscInt i; in Limit2_BeamWarming() local
171 for (i = 0; i < info->m; i++) lmt[i] = jL[i] / info->hxs; in Limit2_BeamWarming()
173 for (i = 0; i < info->m; i++) lmt[i] = jL[i] / info->hxs; in Limit2_BeamWarming()
175 for (i = 0; i < info->m; i++) lmt[i] = jL[i] / (info->hxs / 2.0 + info->hxf / 2.0); in Limit2_BeamWarming()
177 for (i = 0; i < info->m; i++) lmt[i] = jL[i] / info->hxf; in Limit2_BeamWarming()
179 for (i = 0; i < info->m; i++) lmt[i] = jL[i] / info->hxf; in Limit2_BeamWarming()
181 for (i = 0; i < info->m; i++) lmt[i] = jL[i] / (info->hxf / 2.0 + info->hxs / 2.0); in Limit2_BeamWarming()
183 for (i = 0; i < info->m; i++) lmt[i] = jL[i] / info->hxs; in Limit2_BeamWarming()
188 PetscInt i; in Limit2_Fromm() local
190 for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] + jR[i]) / info->hxs; in Limit2_Fromm()
192 …for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] / info->hxs + jR[i] / (info->hxs / 2.0 + info-… in Limit2_Fromm()
194 …for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] / (info->hxs / 2.0 + info->hxf / 2.0) + jR[i] … in Limit2_Fromm()
196 for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] + jR[i]) / info->hxf; in Limit2_Fromm()
198 …for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] / info->hxf + jR[i] / (info->hxf / 2.0 + info-… in Limit2_Fromm()
200 …for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] / (info->hxf / 2.0 + info->hxs / 2.0) + jR[i] … in Limit2_Fromm()
202 for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] + jR[i]) / info->hxs; in Limit2_Fromm()
207 PetscInt i; in Limit2_Minmod() local
209 for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i], jR[i]) / info->hxs; in Limit2_Minmod()
211 …for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i] / info->hxs, jR[i] / (info->hxs / 2.0 + info-… in Limit2_Minmod()
213 …for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i] / (info->hxs / 2.0 + info->hxf / 2.0), jR[i] … in Limit2_Minmod()
215 for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i], jR[i]) / info->hxf; in Limit2_Minmod()
217 …for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i] / info->hxf, jR[i] / (info->hxf / 2.0 + info-… in Limit2_Minmod()
219 …for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i] / (info->hxf / 2.0 + info->hxs / 2.0), jR[i] … in Limit2_Minmod()
221 for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i], jR[i]) / info->hxs; in Limit2_Minmod()
226 PetscInt i; in Limit2_Superbee() local
228 …for (i = 0; i < info->m; i++) lmt[i] = MaxMod2(MinMod2(jL[i], 2 * jR[i]), MinMod2(2 * jL[i], jR[i]… in Limit2_Superbee()
230 …i = 0; i < info->m; i++) lmt[i] = MaxMod2(MinMod2(jL[i] / info->hxs, 2 * jR[i] / (info->hxs / 2.0 … in Limit2_Superbee()
232 …i = 0; i < info->m; i++) lmt[i] = MaxMod2(MinMod2(jL[i] / (info->hxs / 2.0 + info->hxf / 2.0), 2 *… in Limit2_Superbee()
234 …for (i = 0; i < info->m; i++) lmt[i] = MinMod2(MinMod2(jL[i], 2 * jR[i]), MinMod2(2 * jL[i], jR[i]… in Limit2_Superbee()
236 …i = 0; i < info->m; i++) lmt[i] = MinMod2(MinMod2(jL[i] / info->hxf, 2 * jR[i] / (info->hxf / 2.0 … in Limit2_Superbee()
238 …i = 0; i < info->m; i++) lmt[i] = MinMod2(MinMod2(jL[i] / (info->hxf / 2.0 + info->hxs / 2.0), 2 *… in Limit2_Superbee()
240 …for (i = 0; i < info->m; i++) lmt[i] = MinMod2(MinMod2(jL[i], 2 * jR[i]), MinMod2(2 * jL[i], jR[i]… in Limit2_Superbee()
245 PetscInt i; in Limit2_MC() local
247 …for (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i], 0.5 * (jL[i] + jR[i]), 2 * jR[i]) / info… in Limit2_MC()
249 …r (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxs, 0.5 * (jL[i] / info->hxs + jR[… in Limit2_MC()
251 …i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxs / 2.0 + info->hxf / 2.0), 0.5 * (… in Limit2_MC()
253 …for (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i], 0.5 * (jL[i] + jR[i]), 2 * jR[i]) / info… in Limit2_MC()
255 …r (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxf, 0.5 * (jL[i] / info->hxf + jR[… in Limit2_MC()
257 …i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxf / 2.0 + info->hxs / 2.0), 0.5 * (… in Limit2_MC()
259 …for (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i], 0.5 * (jL[i] + jR[i]), 2 * jR[i]) / info… in Limit2_MC()
264 PetscInt i; in Limit2_Koren3() local
266 …for (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i], (jL[i] + 2 * jR[i]) / 3, 2 * jR[i]) / in… in Limit2_Koren3()
268 … (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxs, (jL[i] / info->hxs + 2 * jR[i] … in Limit2_Koren3()
270 …i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxs / 2.0 + info->hxf / 2.0), (jL[i] … in Limit2_Koren3()
272 …for (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i], (jL[i] + 2 * jR[i]) / 3, 2 * jR[i]) / in… in Limit2_Koren3()
274 … (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxf, (jL[i] / info->hxf + 2 * jR[i] … in Limit2_Koren3()
276 …i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxf / 2.0 + info->hxs / 2.0), (jL[i] … in Limit2_Koren3()
278 …for (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i], (jL[i] + 2 * jR[i]) / 3, 2 * jR[i]) / in… in Limit2_Koren3()
285 PetscInt i; in Limit3_Upwind() local
286 for (i = 0; i < info->m; i++) lmt[i] = 0; in Limit3_Upwind()
290 PetscInt i; in Limit3_LaxWendroff() local
292 for (i = 0; i < info->m; i++) lmt[i] = jR[i] / info->hxs; in Limit3_LaxWendroff()
294 for (i = 0; i < info->m; i++) lmt[i] = jR[i] / (info->hxs / 2.0 + info->hxm / 2.0); in Limit3_LaxWendroff()
296 for (i = 0; i < info->m; i++) lmt[i] = jR[i] / info->hxm; in Limit3_LaxWendroff()
298 for (i = 0; i < info->m; i++) lmt[i] = jR[i] / (info->hxm / 2.0 + info->hxf / 2.0); in Limit3_LaxWendroff()
300 for (i = 0; i < info->m; i++) lmt[i] = jR[i] / info->hxf; in Limit3_LaxWendroff()
305 PetscInt i; in Limit3_BeamWarming() local
307 for (i = 0; i < info->m; i++) lmt[i] = jL[i] / info->hxs; in Limit3_BeamWarming()
309 for (i = 0; i < info->m; i++) lmt[i] = jL[i] / (info->hxs / 2.0 + info->hxf / 2.0); in Limit3_BeamWarming()
311 for (i = 0; i < info->m; i++) lmt[i] = jL[i] / info->hxm; in Limit3_BeamWarming()
313 for (i = 0; i < info->m; i++) lmt[i] = jL[i] / (info->hxm / 2.0 + info->hxf / 2.0); in Limit3_BeamWarming()
315 for (i = 0; i < info->m; i++) lmt[i] = jL[i] / info->hxf; in Limit3_BeamWarming()
320 PetscInt i; in Limit3_Fromm() local
322 for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] + jR[i]) / info->hxs; in Limit3_Fromm()
324 …for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] / info->hxs + jR[i] / (info->hxs / 2.0 + info-… in Limit3_Fromm()
326 …for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] / (info->hxs / 2.0 + info->hxm / 2.0) + jR[i] … in Limit3_Fromm()
328 …for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] / info->hxm + jR[i] / (info->hxs / 2.0 + info-… in Limit3_Fromm()
330 …for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] / (info->hxm / 2.0 + info->hxs / 2.0) + jR[i] … in Limit3_Fromm()
332 for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] + jR[i]) / info->hxm; in Limit3_Fromm()
334 …for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] / info->hxm + jR[i] / (info->hxm / 2.0 + info-… in Limit3_Fromm()
336 …for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] / (info->hxm / 2.0 + info->hxf / 2.0) + jR[i] … in Limit3_Fromm()
338 …for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] / info->hxf + jR[i] / (info->hxf / 2.0 + info-… in Limit3_Fromm()
340 …for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] / (info->hxf / 2.0 + info->hxm / 2.0) + jR[i] … in Limit3_Fromm()
342 for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] + jR[i]) / info->hxf; in Limit3_Fromm()
347 PetscInt i; in Limit3_Minmod() local
349 for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i], jR[i]) / info->hxs; in Limit3_Minmod()
351 …for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i] / info->hxs, jR[i] / (info->hxs / 2.0 + info-… in Limit3_Minmod()
353 …for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i] / (info->hxs / 2.0 + info->hxf / 2.0), jR[i] … in Limit3_Minmod()
355 …for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i] / info->hxm, jR[i] / (info->hxm / 2.0 + info-… in Limit3_Minmod()
357 …for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i] / (info->hxm / 2.0 + info->hxs / 2.0), jR[i] … in Limit3_Minmod()
359 for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i], jR[i]) / info->hxm; in Limit3_Minmod()
361 …for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i] / info->hxm, jR[i] / (info->hxm / 2.0 + info-… in Limit3_Minmod()
363 …for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i] / (info->hxm / 2.0 + info->hxf / 2.0), jR[i] … in Limit3_Minmod()
365 …for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i] / info->hxf, jR[i] / (info->hxf / 2.0 + info-… in Limit3_Minmod()
367 …for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i] / (info->hxf / 2.0 + info->hxm / 2.0), jR[i] … in Limit3_Minmod()
369 for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] + jR[i]) / info->hxf; in Limit3_Minmod()
374 PetscInt i; in Limit3_Superbee() local
376 …for (i = 0; i < info->m; i++) lmt[i] = MaxMod2(MinMod2(jL[i], 2 * jR[i]), MinMod2(2 * jL[i], jR[i]… in Limit3_Superbee()
378 …i = 0; i < info->m; i++) lmt[i] = MaxMod2(MinMod2(jL[i] / info->hxs, 2 * jR[i] / (info->hxs / 2.0 … in Limit3_Superbee()
380 …i = 0; i < info->m; i++) lmt[i] = MaxMod2(MinMod2(jL[i] / (info->hxs / 2.0 + info->hxm / 2.0), 2 *… in Limit3_Superbee()
382 …i = 0; i < info->m; i++) lmt[i] = MinMod2(MinMod2(jL[i] / info->hxm, 2 * jR[i] / (info->hxm / 2.0 … in Limit3_Superbee()
384 …i = 0; i < info->m; i++) lmt[i] = MinMod2(MinMod2(jL[i] / (info->hxm / 2.0 + info->hxs / 2.0), 2 *… in Limit3_Superbee()
386 …for (i = 0; i < info->m; i++) lmt[i] = MaxMod2(MinMod2(jL[i], 2 * jR[i]), MinMod2(2 * jL[i], jR[i]… in Limit3_Superbee()
388 …i = 0; i < info->m; i++) lmt[i] = MaxMod2(MinMod2(jL[i] / info->hxm, 2 * jR[i] / (info->hxm / 2.0 … in Limit3_Superbee()
390 …i = 0; i < info->m; i++) lmt[i] = MaxMod2(MinMod2(jL[i] / (info->hxm / 2.0 + info->hxf / 2.0), 2 *… in Limit3_Superbee()
392 …i = 0; i < info->m; i++) lmt[i] = MinMod2(MinMod2(jL[i] / info->hxf, 2 * jR[i] / (info->hxf / 2.0 … in Limit3_Superbee()
394 …i = 0; i < info->m; i++) lmt[i] = MinMod2(MinMod2(jL[i] / (info->hxf / 2.0 + info->hxm / 2.0), 2 *… in Limit3_Superbee()
396 …for (i = 0; i < info->m; i++) lmt[i] = MaxMod2(MinMod2(jL[i], 2 * jR[i]), MinMod2(2 * jL[i], jR[i]… in Limit3_Superbee()
401 PetscInt i; in Limit3_MC() local
403 …for (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i], 0.5 * (jL[i] + jR[i]), 2 * jR[i]) / info… in Limit3_MC()
405 …r (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxs, 0.5 * (jL[i] / info->hxs + jR[… in Limit3_MC()
407 …i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxs / 2.0 + info->hxm / 2.0), 0.5 * (… in Limit3_MC()
409 …r (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxm, 0.5 * (jL[i] / info->hxm + jR[… in Limit3_MC()
411 …i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxm / 2.0 + info->hxs / 2.0), 0.5 * (… in Limit3_MC()
413 …for (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i], 0.5 * (jL[i] + jR[i]), 2 * jR[i]) / info… in Limit3_MC()
415 …r (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxm, 0.5 * (jL[i] / info->hxm + jR[… in Limit3_MC()
417 …i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxm / 2.0 + info->hxf / 2.0), 0.5 * (… in Limit3_MC()
419 …r (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxf, 0.5 * (jL[i] / info->hxf + jR[… in Limit3_MC()
421 …i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxf / 2.0 + info->hxm / 2.0), 0.5 * (… in Limit3_MC()
423 …for (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i], 0.5 * (jL[i] + jR[i]), 2 * jR[i]) / info… in Limit3_MC()
428 PetscInt i; in Limit3_Koren3() local
430 …for (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i], (jL[i] + 2 * jR[i]) / 3, 2 * jR[i]) / in… in Limit3_Koren3()
432 … (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxs, (jL[i] / info->hxs + 2 * jR[i] … in Limit3_Koren3()
434 …i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxs / 2.0 + info->hxm / 2.0), (jL[i] … in Limit3_Koren3()
436 … (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxm, (jL[i] / info->hxm + 2 * jR[i] … in Limit3_Koren3()
438 …i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxm / 2.0 + info->hxs / 2.0), (jL[i] … in Limit3_Koren3()
440 …for (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i], (jL[i] + 2 * jR[i]) / 3, 2 * jR[i]) / in… in Limit3_Koren3()
442 … (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxm, (jL[i] / info->hxm + 2 * jR[i] … in Limit3_Koren3()
444 …i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxm / 2.0 + info->hxf / 2.0), (jL[i] … in Limit3_Koren3()
446 … (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxf, (jL[i] / info->hxf + 2 * jR[i] … in Limit3_Koren3()
448 …i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxf / 2.0 + info->hxm / 2.0), (jL[i] … in Limit3_Koren3()
450 …for (i = 0; i < info->m; i++) lmt[i] = MinMod3(2 * jL[i], (jL[i] + 2 * jR[i]) / 3, 2 * jR[i]) / in… in Limit3_Koren3()
526 PetscInt i, j, k, Mx, dof, xs, xm; in FVRHSFunction() local
546 for (i = xs - 2; i < 0; i++) { in FVRHSFunction()
547 for (j = 0; j < dof; j++) x[i * dof + j] = x[j]; in FVRHSFunction()
549 for (i = Mx; i < xs + xm + 2; i++) { in FVRHSFunction()
550 for (j = 0; j < dof; j++) x[i * dof + j] = x[(xs + xm - 1) * dof + j]; in FVRHSFunction()
554 for (i = xs - 1; i < xs + xm + 1; i++) { in FVRHSFunction()
558 …cs.characteristic)(ctx->physics.user, dof, &x[i * dof], ctx->R, ctx->Rinv, ctx->speeds, ctx->xmin … in FVRHSFunction()
565 jmpL = x[(i + 0) * dof + j] - x[(i - 1) * dof + j]; in FVRHSFunction()
566 jmpR = x[(i + 1) * dof + j] - x[(i + 0) * dof + j]; in FVRHSFunction()
580 slope[i * dof + j] = tmp; in FVRHSFunction()
584 for (i = xs; i < xs + xm + 1; i++) { in FVRHSFunction()
590 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hx / 2; in FVRHSFunction()
591 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hx / 2; in FVRHSFunction()
593 …ann)(ctx->physics.user, dof, uL, uR, ctx->flux, &maxspeed, ctx->xmin + hx * i, ctx->xmin, ctx->xma… in FVRHSFunction()
595 if (i > xs) { in FVRHSFunction()
596 for (j = 0; j < dof; j++) f[(i - 1) * dof + j] -= ctx->flux[j] / hx; in FVRHSFunction()
598 if (i < xs + xm) { in FVRHSFunction()
599 for (j = 0; j < dof; j++) f[i * dof + j] += ctx->flux[j] / hx; in FVRHSFunction()
623 PetscInt i, j, k, dof, xs, xm, Mx; in FVSample() local
631 for (i = xs; i < xs + xm; i++) { in FVSample()
632 const PetscReal h = (ctx->xmax - ctx->xmin) / Mx, xi = ctx->xmin + h / 2 + i * h; in FVSample()
635 for (k = 0; k < dof; k++) u[i * dof + k] = 0; in FVSample()
639 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample()
652 PetscInt imin, imax, Mx, i, j, xs, xm, dof; in SolutionStatsView() local
667 for (i = xs; i < xs + xm; i++) { in SolutionStatsView()
668 for (j = 0; j < dof; j++) tvsum += PetscAbsScalar(x[i * dof + j] - x[(i - 1) * dof + j]); in SolutionStatsView()