Long-memory initialisation for RENs

[nic-barbara]

The long-memory initialiser for RENs will work as intended for contracting RENs, but won't for any of the other REN parametrisations because the H matrix also depends on the output variables. We should fix this.