Hi, 
   I noticed this while going through the source code:

file: fe_interface.C
line: 169

it reads:

return FEXYZ<1>::n_shape_functions(t, o);

should be:

return FEXYZ<3>::n_shape_functions(t, o)

I hope this corrects it.

Manav