Lattice relaxation does not work with domain decomposition

The code is not currently working, and fails inside par_vec_init. This needs to be fixed an a test needs to be added.