Simplify and fix bem with pinned phi1 values

* Much easier to assume that there is only one pinned value and it is
  0. Since we only need this to fix the singularity of the Jacobian that
  is all we need.

* Included paranoid check for this.

* Probably can merge the two bem handlers now.