Used in rolling re-runs: pads old_trg with one extra row and column of
NA, then fills in the gaps from new_trg (which already contains the
newly observed diagonal).
Usage
diag_writer(old_trg, new_trg)
Arguments
old_trg
A (n) x (n) matrix - the old triangle.
new_trg
A (n+1) x (n+1) matrix - the new triangle including the
freshly observed diagonal.