Now we store the H(-R) for positive directions
adn H(R) for negative directions and use the reverse
phases.

This means we get faster indexing when creating the
matrices for SE calculations. I have checked that
the results are the same.