52     const PetscScalar hx = i + 1 == info.xs + info.xm ? x[i] - x[i - 1] : x[i + 1] - x[i];  in FormRHSFunction()  local
54     f[i].u  = hx * (u[i].v);  in FormRHSFunction()
55 …f[i].v  = -hx * (PetscSqr(user->gammaTilde) * u[i].u + (PetscSqr(user->gamma) / user->xi) * (u[i].…  in FormRHSFunction()
56     f[i].th = -hx * (u[i].v + 1) * (u[i].th + (1 + user->epsilon) * PetscLogScalar(u[i].v + 1));  in FormRHSFunction()
88       const PetscScalar hx = x[i + 1] - x[i];  in FormIFunction()  local
89       f[i].u               = hx * udot[i].u;  in FormIFunction()
90       f[i].v               = hx * udot[i].v - PetscSqr(user->c) * (u[i + 1].u - u[i].u) / hx;  in FormIFunction()
91       f[i].th              = hx * udot[i].th;  in FormIFunction()
93       const PetscScalar hx = x[i] - x[i - 1];  in FormIFunction()  local
94       f[i].u               = hx * udot[i].u;  in FormIFunction()
95       f[i].v               = hx * udot[i].v - PetscSqr(user->c) * (u[i - 1].u - u[i].u) / hx;  in FormIFunction()
96       f[i].th              = hx * udot[i].th;  in FormIFunction()
98       const PetscScalar hx = x[i + 1] - x[i];  in FormIFunction()  local
99       f[i].u               = hx * udot[i].u;  in FormIFunction()
100 …f[i].v               = hx * udot[i].v - PetscSqr(user->c) * (u[i - 1].u - 2. * u[i].u + u[i + 1].u…  in FormIFunction()
101       f[i].th              = hx * udot[i].th;  in FormIFunction()
133       const PetscScalar hx  = x[i + 1] - x[i];  in FormIJacobian()  local
135       const PetscScalar dxx0 = PetscSqr(user->c) / hx, dxxR = -PetscSqr(user->c) / hx;  in FormIJacobian()
137         {{a * hx, 0, 0},        {0, 0, 0}   },  in FormIJacobian()
138         {{0, a * hx + dxx0, 0}, {0, dxxR, 0}},  in FormIJacobian()
139         {{0, 0, a * hx},        {0, 0, 0}   }  in FormIJacobian()
144       const PetscScalar hx  = x[i + 1] - x[i];  in FormIJacobian()  local
146       const PetscScalar dxxL = -PetscSqr(user->c) / hx, dxx0 = PetscSqr(user->c) / hx;  in FormIJacobian()
148         {{0, 0, 0},    {a * hx, 0, 0}       },  in FormIJacobian()
149         {{0, dxxL, 0}, {0, a * hx + dxx0, 0}},  in FormIJacobian()
150         {{0, 0, 0},    {0, 0, a * hx}       }  in FormIJacobian()
155       const PetscScalar hx  = x[i + 1] - x[i];  in FormIJacobian()  local
157 … PetscScalar dxxL = -PetscSqr(user->c) / hx, dxx0 = 2. * PetscSqr(user->c) / hx, dxxR = -PetscSqr(…  in FormIJacobian()
159         {{0, 0, 0},    {a * hx, 0, 0},        {0, 0, 0}   },  in FormIJacobian()
160         {{0, dxxL, 0}, {0, a * hx + dxx0, 0}, {0, dxxR, 0}},  in FormIJacobian()
161         {{0, 0, 0},    {0, 0, a * hx},        {0, 0, 0}   }  in FormIJacobian()
214   PetscReal         ftime, hx, dt;  in main()  local
255   hx = 20.0 / (PetscReal)(mx - 1);  in main()
256   dt = 0.4 * PetscSqr(hx) / PetscSqr(user.c); /* Diffusive stability limit */  in main()