Contact dynamics, contact Poisson bracket, and symplectic integrator -- Even Arnold nodds

By introducing an integration factor to the differential one-form of contact dynamics, equations of motion are derived variationally, and contact Poisson bracket and contact Lagrangian are formulated. Discrete symplectic integrator, named hybrid leap-frog method, is found to be a numerical solver of equations of motion for contact dissipative systems.