Calculate and store the random effects on the B scale for each evalution.

They were being calculated anyway but in storage allocated in each call to
`updateμ!` so allocating only once is actually faster.