Replaces the per-time-period R loop with a single C++ call that loops internally, reducing R-to-C++ transition overhead.
rZ_nrm_batch_cpp(Z, EZ, rho, s2, Y)List with updated Z and E_nrm (residuals)
3D array of latent values (n x n x T)
3D array of expected values (n x n x T)
Dyadic correlation parameter
Dyadic variance
3D array of observed values (n x n x T)