0 is an acceptable error so the inequalities needed changed for the assertion in output.cc
Added cases to the setup so that simulation_type==2 does not have Interface style bcs. This makes sense because it is only to test moving domain problems with given boundary conditions.
Added placeholder values to the error returned by structure so that the error is not always zero. This prevents python from throwing a log(0/0) error. Also, added some comments for the second partials matrix in building the RHS.