21   const BlasiusContext blasius = (BlasiusContext)ctx;  in CompressibleBlasiusResidual()  local
24   PetscInt             N       = blasius->n_cheb;  in CompressibleBlasiusResidual()
25   State                S_infty = blasius->S_infty;  in CompressibleBlasiusResidual()
29 …PetscScalar Ma = Mach(&blasius->newtonian_ctx, S_infty.Y.temperature, U_infty), Pr = Prandtl(&blas…  in CompressibleBlasiusResidual()
30               gamma = HeatCapacityRatio(&blasius->newtonian_ctx);  in CompressibleBlasiusResidual()
37   ChebyshevEval(N, Tf, -1., blasius->eta_max, f);  in CompressibleBlasiusResidual()
42   ChebyshevEval(N, Tf, 1., blasius->eta_max, f);  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()
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()
76 PetscErrorCode ComputeChebyshevCoefficients(BlasiusContext blasius) {  in ComputeChebyshevCoefficients()  argument
80   PetscInt            N = blasius->n_cheb;  in ComputeChebyshevCoefficients()
86   PetscCall(PetscMalloc2(N - 3, &blasius->X, N - 3, &w));  in ComputeChebyshevCoefficients()
87   PetscCall(PetscDTGaussQuadrature(N - 3, -1., 1., blasius->X, w));  in ComputeChebyshevCoefficients()
97   PetscCall(SNESSetFunction(snes, res, CompressibleBlasiusResidual, blasius));  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()
110   PetscCall(PetscFree2(blasius->X, w));  in ComputeChebyshevCoefficients()