Fix FEMFunction projection bug reported by Manav Bhatia.

Discussed in [this](http://sourceforge.net/p/libmesh/mailman/message/31428231/) thread. To summarize, need to get quadrature rule info after reinit elements. This only broke on 3D faces because QGauss init's automatically for 1D rules. Example code provided by Manav Bhatia now runs successfully. Also tested in GRINS to be double sure.