32   PetscInt   N;  member
56   blasius->N     = 10;  /* Number of Chebyshev terms */  in main()
64 …PetscCall(PetscOptionsInt("-N", "Number of Chebyshev terms for f", "", blasius->N, &blasius->N, NU…  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()
129 static void ChebyshevEval(PetscInt N, const PetscScalar *Tf, PetscReal x, PetscReal dx_deta, PetscS…  in ChebyshevEval()  argument
138   for (int i = 3; i < N; i++) {  in ChebyshevEval()
165   PetscInt           N  = blasius->N;  in FormFunction()  local
177   Th = Tf + N;  in FormFunction()
181   ChebyshevEval(N, Tf, -1., blasius->dx_deta, f);  in FormFunction()
184   ChebyshevEval(N, Tf, 1., blasius->dx_deta, f);  in FormFunction()
186   for (int i = 0; i < N - 3; i++) {  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()
195   r[N + 1] = h[0] - 1; /* Left end boundary condition */  in FormFunction()