24   PetscInt             N       = blasius->n_cheb;  in CompressibleBlasiusResidual()  local
33   Th = Tf + N;  in CompressibleBlasiusResidual()
37   ChebyshevEval(N, Tf, -1., blasius->eta_max, f);  in CompressibleBlasiusResidual()
42   ChebyshevEval(N, Tf, 1., blasius->eta_max, f);  in CompressibleBlasiusResidual()
45   for (int i = 0; i < N - 3; i++) {  in CompressibleBlasiusResidual()
46     ChebyshevEval(N, Tf, blasius->X[i], blasius->eta_max, f);  in CompressibleBlasiusResidual()
47     ChebyshevEval(N - 1, Th, blasius->X[i], blasius->eta_max, h);  in CompressibleBlasiusResidual()
59 …r[N + 2 + i] = (mu_rho_tilde[0] * h[2] + mu_rho_tilde[1] * h[1]) + Pr * f[0] * h[1] + Pr * (gamma …  in CompressibleBlasiusResidual()
63   ChebyshevEval(N - 1, Th, -1., blasius->eta_max, h);  in CompressibleBlasiusResidual()
64   r[N] = h[0] - blasius->T_wall / S_infty.Y.temperature;  in CompressibleBlasiusResidual()
67   ChebyshevEval(N - 1, Th, 1., blasius->eta_max, h);  in CompressibleBlasiusResidual()
68   r[N + 1] = h[0] - 1.;  in CompressibleBlasiusResidual()
80   PetscInt            N = blasius->n_cheb;  in ComputeChebyshevCoefficients()  local
86   PetscCall(PetscMalloc2(N - 3, &blasius->X, N - 3, &w));  in ComputeChebyshevCoefficients()
87   PetscCall(PetscDTGaussQuadrature(N - 3, -1., 1., blasius->X, w));  in ComputeChebyshevCoefficients()
92   PetscCall(VecSetSizes(sol, PETSC_DECIDE, 2 * N - 1));  in ComputeChebyshevCoefficients()
95   PetscCall(VecSetValue(sol, N, 1., INSERT_VALUES));  in ComputeChebyshevCoefficients()
106   for (int i = 0; i < N; i++) blasius->Tf_cheb[i] = cheb_coefs[i];  in ComputeChebyshevCoefficients()
107   for (int i = 0; i < N - 1; i++) blasius->Th_cheb[i] = cheb_coefs[i + N];  in ComputeChebyshevCoefficients()
166 static PetscErrorCode ModifyMesh(MPI_Comm comm, DM dm, PetscInt dim, PetscReal growth, PetscInt N, …  in ModifyMesh()  argument
195     PetscReal dy1 = refine_height * (growth - 1) / (pow(growth, N) - 1);  in ModifyMesh()
198     PetscReal logdy = (log(domain_max[1]) - log(refine_height)) / (faces[1] - N);  in ModifyMesh()
206       if (y_box_index <= N) {  in ModifyMesh()
210         PetscInt j   = y_box_index - N;  in ModifyMesh()
265   PetscInt   N                                    = 20;           // Number of Chebyshev terms  in NS_BLASIUS()  local
285   PetscCall(PetscOptionsInt("-n_chebyshev", "Number of Chebyshev terms", NULL, N, &N, NULL));  in NS_BLASIUS()
286 …scCheck(3 <= N && N <= BLASIUS_MAX_N_CHEBYSHEV, comm, PETSC_ERR_ARG_OUTOFRANGE, "-n_chebyshev %" P…  in NS_BLASIUS()
335   blasius_ctx->n_cheb   = N;  in NS_BLASIUS()