Bugfix/petsc linear operator setup

Fixed an issue with preallocation in petsc matrix in case that the blocking of the domain and range space are different.