18   const BlasiusContext  blasius = (BlasiusContext)ctx;  in CompressibleBlasiusResidual()  local
21   PetscInt              N       = blasius->n_cheb;  in CompressibleBlasiusResidual()
22   State                 S_infty = blasius->S_infty;  in CompressibleBlasiusResidual()
24   NewtonianIGProperties gas     = blasius->newt_ctx.gas;  in CompressibleBlasiusResidual()
34   ChebyshevEval(N, Tf, -1., blasius->eta_max, f);  in CompressibleBlasiusResidual()
39   ChebyshevEval(N, Tf, 1., blasius->eta_max, f);  in CompressibleBlasiusResidual()
43     ChebyshevEval(N, Tf, blasius->X[i], blasius->eta_max, f);  in CompressibleBlasiusResidual()
44     ChebyshevEval(N - 1, Th, blasius->X[i], blasius->eta_max, h);  in CompressibleBlasiusResidual()
60   ChebyshevEval(N - 1, Th, -1., blasius->eta_max, h);  in CompressibleBlasiusResidual()
61   r[N] = h[0] - blasius->T_wall / S_infty.Y.temperature;  in CompressibleBlasiusResidual()
64   ChebyshevEval(N - 1, Th, 1., blasius->eta_max, h);  in CompressibleBlasiusResidual()
73 PetscErrorCode ComputeChebyshevCoefficients(BlasiusContext blasius) {  in ComputeChebyshevCoefficients()  argument
77   PetscInt            N = blasius->n_cheb;  in ComputeChebyshevCoefficients()
83   PetscCall(PetscMalloc2(N - 3, &blasius->X, N - 3, &w));  in ComputeChebyshevCoefficients()
84   PetscCall(PetscDTGaussQuadrature(N - 3, -1., 1., blasius->X, w));  in ComputeChebyshevCoefficients()
94   PetscCall(SNESSetFunction(snes, res, CompressibleBlasiusResidual, blasius));  in ComputeChebyshevCoefficients()
103   for (int i = 0; i < N; i++) blasius->Tf_cheb[i] = cheb_coefs[i];  in ComputeChebyshevCoefficients()
104   for (int i = 0; i < N - 1; i++) blasius->Th_cheb[i] = cheb_coefs[i + N];  in ComputeChebyshevCoefficients()
107   PetscCall(PetscFree2(blasius->X, w));  in ComputeChebyshevCoefficients()