Derived jacobians

This should be ready to pull now; assembly runs ten times faster and the results are within ~1e-10 of the numeric jacobians.  You'll need the libMesh master or the upcoming release to avoid the euler_solver.C bug though.