136   for (CeedInt n = 0; n < nmodes; n++) {  in CalcSpectrum()  local
137     const CeedScalar dkappa = n == 0 ? kappa[0] : kappa[n] - kappa[n - 1];  in CalcSpectrum()
138     qn[n]                   = Calc_qn(kappa[n], dkappa, keta, kcut, ke, 1.0);  in CalcSpectrum()
139     Ektot += qn[n];  in CalcSpectrum()
143   for (CeedInt n = 0; n < nmodes; n++) qn[n] /= Ektot;  in CalcSpectrum()  local
167   CeedPragmaSIMD for (CeedInt n = 0; n < nmodes; n++) {  in StgShur14Calc()  local
168     xhat[0] = (X[0] - stg_ctx->u0 * t) * Max(2 * kappa[0] / kappa[n], 0.1);  in StgShur14Calc()
170     for (CeedInt i = 0; i < 3; i++) xdotd += d[i * nmodes + n] * xhat[i];  in StgShur14Calc()
171     const CeedScalar cos_kxdp = cos(kappa[n] * xdotd + phi[n]);  in StgShur14Calc()
172     vp[0] += sqrt(qn[n]) * sigma[0 * nmodes + n] * cos_kxdp;  in StgShur14Calc()
173     vp[1] += sqrt(qn[n]) * sigma[1 * nmodes + n] * cos_kxdp;  in StgShur14Calc()
174     vp[2] += sqrt(qn[n]) * sigma[2 * nmodes + n] * cos_kxdp;  in StgShur14Calc()
212   CeedPragmaSIMD for (CeedInt n = 0; n < nmodes; n++) {  in StgShur14Calc_PrecompEktot()  local
213     xhat[0] = (X[0] - stg_ctx->u0 * t) * Max(2 * kappa[0] / kappa[n], 0.1);  in StgShur14Calc_PrecompEktot()
215     for (CeedInt i = 0; i < 3; i++) xdotd += d[i * nmodes + n] * xhat[i];  in StgShur14Calc_PrecompEktot()
216     const CeedScalar cos_kxdp = cos(kappa[n] * xdotd + phi[n]);  in StgShur14Calc_PrecompEktot()
217     const CeedScalar dkappa   = n == 0 ? kappa[0] : kappa[n] - kappa[n - 1];  in StgShur14Calc_PrecompEktot()
218     const CeedScalar qn       = Calc_qn(kappa[n], dkappa, keta, kcut, ke, Ektot);  in StgShur14Calc_PrecompEktot()
219     vp[0] += sqrt(qn) * sigma[0 * nmodes + n] * cos_kxdp;  in StgShur14Calc_PrecompEktot()
220     vp[1] += sqrt(qn) * sigma[1 * nmodes + n] * cos_kxdp;  in StgShur14Calc_PrecompEktot()
221     vp[2] += sqrt(qn) * sigma[2 * nmodes + n] * cos_kxdp;  in StgShur14Calc_PrecompEktot()
270     CeedPragmaSIMD for (CeedInt n = 0; n < nmodes; n++) {  in StgShur14Preprocess()  local
271       const CeedScalar dkappa = n == 0 ? kappa[0] : kappa[n] - kappa[n - 1];  in StgShur14Preprocess()
272       Ek_tot += Calc_qn(kappa[n], dkappa, keta, kcut, ke, 1.0);  in StgShur14Preprocess()