The relevant formula can be found in Computer Physics Communications 277 (2022) 108384.
Computer Physics Communications 277 (2022) 108384