Assembly of non-symmetric system matrices does not work

[kmoks]

The subroutine ADDEM2 for assembling dense system matrices does not account for possibly non-symmetric element matrices when the model also has linear couplings (constraint equations). Off-diagonal terms that should end up on the diagonal are then either ignored or counted twice, instead of being summed.