Lines Matching refs:info
40 void Limit_Upwind(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, PetscScalar *lmt) in Limit_Upwind() argument
43 for (i = 0; i < info->m; i++) lmt[i] = 0; in Limit_Upwind()
45 void Limit_LaxWendroff(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, PetscScalar *l… in Limit_LaxWendroff() argument
48 for (i = 0; i < info->m; i++) lmt[i] = jR[i]; in Limit_LaxWendroff()
50 void Limit_BeamWarming(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, PetscScalar *l… in Limit_BeamWarming() argument
53 for (i = 0; i < info->m; i++) lmt[i] = jL[i]; in Limit_BeamWarming()
55 void Limit_Fromm(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, PetscScalar *lmt) in Limit_Fromm() argument
58 for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] + jR[i]); in Limit_Fromm()
60 void Limit_Minmod(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, PetscScalar *lmt) in Limit_Minmod() argument
63 for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i], jR[i]); in Limit_Minmod()
65 void Limit_Superbee(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, PetscScalar *lmt) in Limit_Superbee() argument
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()
70 void Limit_MC(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, PetscScalar *lmt) in Limit_MC() argument
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()
75 void Limit_VanLeer(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, PetscScalar *lmt) in Limit_VanLeer() argument
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()
80 void Limit_VanAlbada(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, PetscScalar *lmt… in Limit_VanAlbada() argument
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()
85 void Limit_VanAlbadaTVD(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, PetscScalar *… in Limit_VanAlbadaTVD() argument
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()
90 void Limit_Koren(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, PetscScalar *lmt) /*… in Limit_Koren() argument
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()
95 void Limit_KorenSym(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, PetscScalar *lmt)… in Limit_KorenSym() argument
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()
100 void Limit_Koren3(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, PetscScalar *lmt) in Limit_Koren3() argument
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()
109 void Limit_CadaTorrilhon2(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, PetscScalar… in Limit_CadaTorrilhon2() argument
112 for (i = 0; i < info->m; i++) lmt[i] = CadaTorrilhonPhiHatR_Eq13(jL[i], jR[i]); in Limit_CadaTorrilhon2()
114 void Limit_CadaTorrilhon3R(PetscReal r, LimitInfo info, const PetscScalar *jL, const PetscScalar *j… in Limit_CadaTorrilhon3R() argument
117 const PetscReal eps = 1e-7, hx = info->hx; in Limit_CadaTorrilhon3R()
119 for (i = 0; i < info->m; i++) { in Limit_CadaTorrilhon3R()
124 void Limit_CadaTorrilhon3R0p1(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, PetscSc… in Limit_CadaTorrilhon3R0p1() argument
126 Limit_CadaTorrilhon3R(0.1, info, jL, jR, lmt); in Limit_CadaTorrilhon3R0p1()
128 void Limit_CadaTorrilhon3R1(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, PetscScal… in Limit_CadaTorrilhon3R1() argument
130 Limit_CadaTorrilhon3R(1, info, jL, jR, lmt); in Limit_CadaTorrilhon3R1()
132 void Limit_CadaTorrilhon3R10(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, PetscSca… in Limit_CadaTorrilhon3R10() argument
134 Limit_CadaTorrilhon3R(10, info, jL, jR, lmt); in Limit_CadaTorrilhon3R10()
136 void Limit_CadaTorrilhon3R100(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, PetscSc… in Limit_CadaTorrilhon3R100() argument
138 Limit_CadaTorrilhon3R(100, info, jL, jR, lmt); in Limit_CadaTorrilhon3R100()
143 void Limit2_Upwind(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, const PetscInt sf,… in Limit2_Upwind() argument
146 for (i = 0; i < info->m; i++) lmt[i] = 0; in Limit2_Upwind()
148 void Limit2_LaxWendroff(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, const PetscIn… in Limit2_LaxWendroff() argument
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()
167 void Limit2_BeamWarming(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, const PetscIn… in Limit2_BeamWarming() argument
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()
186 void Limit2_Fromm(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, const PetscInt sf, … in Limit2_Fromm() argument
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()
205 void Limit2_Minmod(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, const PetscInt sf,… in Limit2_Minmod() argument
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()
224 void Limit2_Superbee(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, const PetscInt s… in Limit2_Superbee() argument
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 …info->m; i++) lmt[i] = MaxMod2(MinMod2(jL[i] / info->hxs, 2 * jR[i] / (info->hxs / 2.0 + info->hxf… in Limit2_Superbee()
232 …info->m; i++) lmt[i] = MaxMod2(MinMod2(jL[i] / (info->hxs / 2.0 + info->hxf / 2.0), 2 * jR[i] / in… 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 …info->m; i++) lmt[i] = MinMod2(MinMod2(jL[i] / info->hxf, 2 * jR[i] / (info->hxf / 2.0 + info->hxs… in Limit2_Superbee()
238 …info->m; i++) lmt[i] = MinMod2(MinMod2(jL[i] / (info->hxf / 2.0 + info->hxs / 2.0), 2 * jR[i] / in… 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()
243 void Limit2_MC(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, const PetscInt sf, con… in Limit2_MC() argument
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 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxs, 0.5 * (jL[i] / info->hxs + jR[i] / (info->hx… in Limit2_MC()
251 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxs / 2.0 + info->hxf / 2.0), 0.5 * (jL[i] / (in… 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 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxf, 0.5 * (jL[i] / info->hxf + jR[i] / (info->hx… in Limit2_MC()
257 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxf / 2.0 + info->hxs / 2.0), 0.5 * (jL[i] / (in… 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()
262 void Limit2_Koren3(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, const PetscInt sf,… in Limit2_Koren3() argument
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 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxs, (jL[i] / info->hxs + 2 * jR[i] / (info->hxf … in Limit2_Koren3()
270 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxs / 2.0 + info->hxf / 2.0), (jL[i] / (info->hx… 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 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxf, (jL[i] / info->hxf + 2 * jR[i] / (info->hxf … in Limit2_Koren3()
276 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxf / 2.0 + info->hxs / 2.0), (jL[i] / (info->hx… 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()
283 void Limit3_Upwind(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, const PetscInt sm,… in Limit3_Upwind() argument
286 for (i = 0; i < info->m; i++) lmt[i] = 0; in Limit3_Upwind()
288 void Limit3_LaxWendroff(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, const PetscIn… in Limit3_LaxWendroff() argument
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()
303 void Limit3_BeamWarming(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, const PetscIn… in Limit3_BeamWarming() argument
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()
318 void Limit3_Fromm(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, const PetscInt sm, … in Limit3_Fromm() argument
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()
345 void Limit3_Minmod(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, const PetscInt sm,… in Limit3_Minmod() argument
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()
372 void Limit3_Superbee(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, const PetscInt s… in Limit3_Superbee() argument
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 …info->m; i++) lmt[i] = MaxMod2(MinMod2(jL[i] / info->hxs, 2 * jR[i] / (info->hxs / 2.0 + info->hxm… in Limit3_Superbee()
380 …info->m; i++) lmt[i] = MaxMod2(MinMod2(jL[i] / (info->hxs / 2.0 + info->hxm / 2.0), 2 * jR[i] / in… in Limit3_Superbee()
382 …info->m; i++) lmt[i] = MinMod2(MinMod2(jL[i] / info->hxm, 2 * jR[i] / (info->hxm / 2.0 + info->hxs… in Limit3_Superbee()
384 …info->m; i++) lmt[i] = MinMod2(MinMod2(jL[i] / (info->hxm / 2.0 + info->hxs / 2.0), 2 * jR[i] / in… 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 …info->m; i++) lmt[i] = MaxMod2(MinMod2(jL[i] / info->hxm, 2 * jR[i] / (info->hxm / 2.0 + info->hxf… in Limit3_Superbee()
390 …info->m; i++) lmt[i] = MaxMod2(MinMod2(jL[i] / (info->hxm / 2.0 + info->hxf / 2.0), 2 * jR[i] / in… in Limit3_Superbee()
392 …info->m; i++) lmt[i] = MinMod2(MinMod2(jL[i] / info->hxf, 2 * jR[i] / (info->hxf / 2.0 + info->hxm… in Limit3_Superbee()
394 …info->m; i++) lmt[i] = MinMod2(MinMod2(jL[i] / (info->hxf / 2.0 + info->hxm / 2.0), 2 * jR[i] / in… 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()
399 void Limit3_MC(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, const PetscInt sm, con… in Limit3_MC() argument
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 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxs, 0.5 * (jL[i] / info->hxs + jR[i] / (info->hx… in Limit3_MC()
407 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxs / 2.0 + info->hxm / 2.0), 0.5 * (jL[i] / (in… in Limit3_MC()
409 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxm, 0.5 * (jL[i] / info->hxm + jR[i] / (info->hx… in Limit3_MC()
411 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxm / 2.0 + info->hxs / 2.0), 0.5 * (jL[i] / (in… 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 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxm, 0.5 * (jL[i] / info->hxm + jR[i] / (info->hx… in Limit3_MC()
417 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxm / 2.0 + info->hxf / 2.0), 0.5 * (jL[i] / (in… in Limit3_MC()
419 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxf, 0.5 * (jL[i] / info->hxf + jR[i] / (info->hx… in Limit3_MC()
421 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxf / 2.0 + info->hxm / 2.0), 0.5 * (jL[i] / (in… 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()
426 void Limit3_Koren3(LimitInfo info, const PetscScalar *jL, const PetscScalar *jR, const PetscInt sm,… in Limit3_Koren3() argument
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 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxs, (jL[i] / info->hxs + 2 * jR[i] / (info->hxm … in Limit3_Koren3()
434 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxs / 2.0 + info->hxm / 2.0), (jL[i] / (info->hx… in Limit3_Koren3()
436 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxm, (jL[i] / info->hxm + 2 * jR[i] / (info->hxm … in Limit3_Koren3()
438 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxm / 2.0 + info->hxs / 2.0), (jL[i] / (info->hx… 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 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxm, (jL[i] / info->hxm + 2 * jR[i] / (info->hxm … in Limit3_Koren3()
444 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxm / 2.0 + info->hxf / 2.0), (jL[i] / (info->hx… in Limit3_Koren3()
446 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / info->hxf, (jL[i] / info->hxf + 2 * jR[i] / (info->hxf … in Limit3_Koren3()
448 …info->m; i++) lmt[i] = MinMod3(2 * jL[i] / (info->hxf / 2.0 + info->hxm / 2.0), (jL[i] / (info->hx… 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()
555 struct _LimitInfo info; in FVRHSFunction() local
573 info.m = dof; in FVRHSFunction()
574 info.hx = hx; in FVRHSFunction()
575 (*ctx->limit)(&info, cjmpL, cjmpR, ctx->cslope); in FVRHSFunction()