131   for (CeedInt n = 0; n < nmodes; n++) {  in CalcSpectrum()  local
132     const CeedScalar dkappa = n == 0 ? kappa[0] : kappa[n] - kappa[n - 1];  in CalcSpectrum()
133     qn[n]                   = Calc_qn(kappa[n], dkappa, keta, kcut, ke, 1.0);  in CalcSpectrum()
134     Ektot += qn[n];  in CalcSpectrum()
138   for (CeedInt n = 0; n < nmodes; n++) qn[n] /= Ektot;  in CalcSpectrum()  local
162   CeedPragmaSIMD for (CeedInt n = 0; n < nmodes; n++) {  in StgShur14Calc()  local
163     xhat[0] = (X[0] - stg_ctx->u0 * t) * Max(2 * kappa[0] / kappa[n], 0.1);  in StgShur14Calc()
165     for (CeedInt i = 0; i < 3; i++) xdotd += d[i * nmodes + n] * xhat[i];  in StgShur14Calc()
166     const CeedScalar cos_kxdp = cos(kappa[n] * xdotd + phi[n]);  in StgShur14Calc()
167     vp[0] += sqrt(qn[n]) * sigma[0 * nmodes + n] * cos_kxdp;  in StgShur14Calc()
168     vp[1] += sqrt(qn[n]) * sigma[1 * nmodes + n] * cos_kxdp;  in StgShur14Calc()
169     vp[2] += sqrt(qn[n]) * sigma[2 * nmodes + n] * cos_kxdp;  in StgShur14Calc()
207   CeedPragmaSIMD for (CeedInt n = 0; n < nmodes; n++) {  in StgShur14Calc_PrecompEktot()  local
208     xhat[0] = (X[0] - stg_ctx->u0 * t) * Max(2 * kappa[0] / kappa[n], 0.1);  in StgShur14Calc_PrecompEktot()
210     for (CeedInt i = 0; i < 3; i++) xdotd += d[i * nmodes + n] * xhat[i];  in StgShur14Calc_PrecompEktot()
211     const CeedScalar cos_kxdp = cos(kappa[n] * xdotd + phi[n]);  in StgShur14Calc_PrecompEktot()
212     const CeedScalar dkappa   = n == 0 ? kappa[0] : kappa[n] - kappa[n - 1];  in StgShur14Calc_PrecompEktot()
213     const CeedScalar qn       = Calc_qn(kappa[n], dkappa, keta, kcut, ke, Ektot);  in StgShur14Calc_PrecompEktot()
214     vp[0] += sqrt(qn) * sigma[0 * nmodes + n] * cos_kxdp;  in StgShur14Calc_PrecompEktot()
215     vp[1] += sqrt(qn) * sigma[1 * nmodes + n] * cos_kxdp;  in StgShur14Calc_PrecompEktot()
216     vp[2] += sqrt(qn) * sigma[2 * nmodes + n] * cos_kxdp;  in StgShur14Calc_PrecompEktot()
272     CeedPragmaSIMD for (CeedInt n = 0; n < nmodes; n++) {  in StgShur14Preprocess()  local
273       const CeedScalar dkappa = n == 0 ? kappa[0] : kappa[n] - kappa[n - 1];  in StgShur14Preprocess()
274       Ek_tot += Calc_qn(kappa[n], dkappa, keta, kcut, ke, 1.0);  in StgShur14Preprocess()