Multistep initialization set up
Vectorization of code is running successfully, but has big error
problems. I think these have to do with the hanging node constraints.
Constraints bug fixed. Seems to work properly!
Fixed a silly typo in the BDF2 code.