expand auto

Doc/DevRef/Multilayers.tex

@@ -718,6 +718,10 @@ In the case of a semi-infinite layer, the vertical component becomes zero and
 
 Last updated to reflect the actual code in May 2023.
 
+%===============================================================================
+\subsection{Call chain}
+%===============================================================================
+
 \footnotesize
 \parindent=0pt
 \parskip=1ex
@@ -758,6 +762,12 @@ compute up and down travelling amplitudes for each sample layer.
 
 Functions \T{fluxes} and \T{topLayerR} are implemented
 in files \SRC{Resample/Specular}{ComputeFluxScalar.cpp} and
-\SRC{Resample/Specular}{ComputeFluxMagnetic.cpp}.
+\SRC{Resample/Specular}{ComputeFluxMagnetic.cpp},
+where they share some local functions.
+
+%===============================================================================
+\subsection{Scalar fluxes}
+%===============================================================================
+
 
 \index{Multilayer|)}%

Resample/Specular/ComputeFluxMagnetic.cpp

@@ -183,7 +183,7 @@ Fluxes Compute::SpecularMagnetic::fluxes(const SliceStack& slices, const R3& k,
     ASSERT(slices.size() == kz.size());
 
     Fluxes result;
-    for (auto& coeff : ::computeFlux(slices, kz, slices.roughnessModel(), forward))
+    for (const MatrixFlux& coeff : ::computeFlux(slices, kz, slices.roughnessModel(), forward))
         result.emplace_back(new MatrixFlux(coeff));
 
     return result;