Bugfix: mb = nb when neither are specified to setup_dmat

This was need to make scalapack solves work on non-square grids.