* allow setting a fixed value anywhere in the domain
* do not use boundaryType but a separate function to check for fixed cell
* set 1.0 on the diagonal in the system matrix and 0.0 for the off-diagonal entries
  to improve numerical accuracy