content_image/355696.jpg content_image/355697.jpg content_image/355698.jpg content_image/355699.jpg content_image/355702.jpg content_image/355701.jpg content_image/355700.jpg content_image/355703.jpg content_image/355704.jpg
Liouville operators as follows: 

<img src='content_image/355696.jpg'>

<img src='content_image/355697.jpg'>

<img src='content_image/355698.jpg'>

where 

Defining 

<img src='content_image/355699.jpg'>

<img src='content_image/355702.jpg'>

<img src='content_image/355701.jpg'>

<img src='content_image/355700.jpg'>

where α = A, B, C , one possible reversible measure-preserving integration algorithm for the BK thermostat is then 

<img src='content_image/355703.jpg'>

Using the so-called direct translation technique [36] we can expand the above symmetric break-up of the Liouville operator into a pseudo-code form, ready to be implemented on the computer: 

<img src='content_image/355704.jpg'>

