HyperbolicModule: Bugfix: only compute boundary d_ij on the upper triangular portion

[tamiko]

I just noticed this. We can avoid computing the lower triangular part of d_ij as it gets symmetrized (and overwritten) immediately afterwards.