Structure preserving discretizations of the Liouville equation and their numerical tests