Compare two sample by coordinate

Initial code for Liouville integrators