* Set off-diagonal matrix blocks to zero if internal Dirichlet is set * use correct problem() function in md framework