Generative programming methods for parallel partial differential field equation solvers