PDE assembly in serial and 2D seems to be doing what it's supposed to when boundary elements are not involved.
branches/ripleygmg_from_3668/ripley/src/Brick.cpp