def kinetic_energy(self, pos, mom, cache={}): return 0.5 * mom.dot(la.cho_solve( (self.mass_matrix_chol, True), mom))