[rand.eng.philox] Make the round states explicit.

This clarifies which state is the final result, and avoids the use of
the vaguely defined variable $X'$. It changes the index variable $q$
to be 1-based.

source/numerics.tex

@@ -3106,27 +3106,24 @@
 \begin{codeblock}
 @$i$@ = @$i$@ + 1
 if (@$i$@ == @$n$@) {
-  @$Y$@ = Philox(@$K$@, @$X$@) // \seebelow
-  @$Z$@ = @$Z$@ + 1
+  @$Y$@ = Philox(@$K$@, @$X$@)        // \seebelow
+  @$Z$@ = @$Z$@ + 1                // this updates $X$
   @$i$@ = 0
 }
 \end{codeblock}
 
 \pnum
 The \tcode{Philox} function maps the length-$n/2$ sequence $K$ and
-the length-$n$ sequence $X$ into a length-$n$ output sequence $Y$.
+the length-$n$ sequence $X$ into a length-$n$ output sequence.
 Philox applies an $r$-round substitution-permutation network to the values in $X$.
-A single round of the generation algorithm performs the following steps:
+That is, there are intermediate values $X^{(0)}, X^{(1)}, \dotsc, X^{(r)}$,
+where $X^{(0)} \cedef X$, and for each round $q$ (with $q = 1, \dotsc, r$),
+$X^{(q)}$ is computed from $X^{(q - 1)}$ as follows. The output sequence is $X^{(r)}$.
 \begin{itemize}
 \item
-The output sequence $X'$ of the previous round
-($X$ in case of the first round)
-is permuted to obtain the intermediate state $V$:
-\begin{codeblock}
-@$V_j = X'_{f_n(j)}$@
-\end{codeblock}
-where $j = 0, \dotsc, n - 1$ and
-$f_n(j)$ is defined in \tref{rand.eng.philox.f}.
+An intermediate state $V^{(q)}$ is obtained by permuting the previous output,
+$V^{(q)}_j \cedef X^{(q - 1)}_{f_n(j)}$,
+where $j = 0, \dotsc, n - 1$, and $f_n(j)$ is defined in \tref{rand.eng.philox.f}.
 
 \begin{floattable}{Values for the word permutation $\bm{f}_{\bm{n}}\bm{(j)}$}{rand.eng.philox.f}
 {l|l|l|l|l|l|l|l|l|l|l|l|l|l|l|l|l|l}
@@ -3144,12 +3141,13 @@
 \end{note}
 
 \item
-The following computations are applied to the elements of the $V$ sequence:
-\begin{codeblock}
-@$X_{2k + 0} = \mullo(V_{2k + 1}, M_{k}, w)$@
-@$X_{2k + 1} = \mulhi(V_{2k + 1}, M_{k}, w) \xor \mathit{key}^q_k \xor V_{2k}$@
-\end{codeblock}
-where:
+The next output $X^{(q)}$ is computed from the elements of the $V^{(q)}$ as follows.
+For $k = 0, \dotsc, n/2 - 1$,
+\begin{itemize}
+\item $X^{(q)}_{2k + 0} = \mullo(V^{(q)}_{2k + 1}, M_{k}, w)$, and
+\item $X^{(q)}_{2k + 1} = \mulhi(V^{(q)}_{2k + 1}, M_{k}, w) \xor \mathit{key}^q_k \xor V^{(q)}_{2k}$,
+\end{itemize}
+where
   \begin{itemize}
   \item
   $\mullo(\tcode{a}, \tcode{b}, \tcode{w})$ is
@@ -3161,15 +3159,9 @@
   the high half of the modular multiplication of \tcode{a} and \tcode{b}:
   $(\left\lfloor (\tcode{a} \cdot \tcode{b}) / 2^w \right\rfloor)$,
 
-  \item
-  $k = 0, \dotsc, n/2 - 1$ is the index in the sequences,
-
-  \item
-  $q = 0, \dotsc, r - 1$ is the index of the round,
-
   \item
   $\mathit{key}^q_k$ is the $k^\text{th}$ round key for round $q$,
-  $\mathit{key}^q_k \cedef (K_k + q \cdot C_k) \mod 2^w$,
+  $\mathit{key}^q_k \cedef (K_k + (q - 1) \cdot C_k) \mod 2^w$,
 
   \item
   $K_k$ are the elements of the key sequence $K$,
@@ -3182,10 +3174,6 @@
   \end{itemize}
 \end{itemize}
 
-\pnum
-After $r$ applications of the single-round function,
-\tcode{Philox} returns the sequence $Y = X'$.
-
 \indexlibraryglobal{philox_engine}%
 \indexlibrarymember{result_type}{philox_engine}%
 \begin{codeblock}