43   Blasius    *blasius;  in main()  local
52   PetscCall(PetscCalloc1(1, &blasius));  in main()
53   blasius->Ma    = 2;   /* Mach number */  in main()
54   blasius->Pr    = 0.7; /* Prandtl number */  in main()
55   blasius->h_0   = 2.;  /* relative temperature at the wall */  in main()
56   blasius->N     = 10;  /* Number of Chebyshev terms */  in main()
57   blasius->gamma = 1.4; /* specific heat ratio */  in main()
60 …PetscCall(PetscOptionsReal("-mach", "Mach number at freestream", "", blasius->Ma, &blasius->Ma, NU…  in main()
61   PetscCall(PetscOptionsReal("-prandtl", "Prandtl number", "", blasius->Pr, &blasius->Pr, NULL));  in main()
62 …PetscCall(PetscOptionsReal("-h_0", "Relative enthalpy at wall", "", blasius->h_0, &blasius->h_0, N…  in main()
63 …PetscCall(PetscOptionsReal("-gamma", "Ratio of specific heats", "", blasius->gamma, &blasius->gamm…  in main()
64 …PetscCall(PetscOptionsInt("-N", "Number of Chebyshev terms for f", "", blasius->N, &blasius->N, NU…  in main()
67   blasius->dx_deta = 2 / L; /* this helps to map [-1,1] to [0,L] */  in main()
68   PetscCall(PetscMalloc2(blasius->N - 3, &blasius->x, blasius->N - 3, &weight));  in main()
69   PetscCall(PetscDTGaussQuadrature(blasius->N - 3, -1., 1., blasius->x, weight));  in main()
83   PetscCall(VecSetSizes(x, PETSC_DECIDE, 2 * blasius->N - 1));  in main()
90   PetscCall(SNESSetFunction(snes, r, FormFunction, blasius));  in main()
116   PetscCall(PetscFree2(blasius->x, weight));  in main()
117   PetscCall(PetscFree(blasius));  in main()
162   Blasius           *blasius = (Blasius *)ctx;  in FormFunction()  local
165   PetscInt           N  = blasius->N;  in FormFunction()
166   PetscReal          Ma = blasius->Ma, Pr = blasius->Pr;  in FormFunction()
181   ChebyshevEval(N, Tf, -1., blasius->dx_deta, f);  in FormFunction()
184   ChebyshevEval(N, Tf, 1., blasius->dx_deta, f);  in FormFunction()
187     ChebyshevEval(N, Tf, blasius->x[i], blasius->dx_deta, f);  in FormFunction()
189     ChebyshevEval(N - 1, Th, blasius->x[i], blasius->dx_deta, h);  in FormFunction()
190     r[N + 2 + i] = h[2] + Pr * f[0] * h[1] + Pr * (blasius->gamma - 1) * PetscSqr(Ma * f[2]);  in FormFunction()
192   ChebyshevEval(N - 1, Th, -1., blasius->dx_deta, h);  in FormFunction()
193   r[N] = h[0] - blasius->h_0; /* Left end boundary condition */  in FormFunction()
194   ChebyshevEval(N - 1, Th, 1., blasius->dx_deta, h);  in FormFunction()