Fix incorrect diagnostics output after Prasad's changes to element matrix setup
Fix build failure after latest merge
Only use current element number of rows when calculating Jacobian with finite differencing. This was causing build failures in combination with the half-merge of Prasad's element matrix setup changes. That was fixed but this is still a good idea.