def ke_head
import bend-prot@0.1.0.0/sim.bend as Simimport Base
import bend-prot@0.1.0.0/geom.bend as G
99def ke_head(h: DynAtom) -> F32:
100 match h:
101 case DA{serial, elem, pos, vel, frc, mass}:
102 ((0.5 * mass : F32) * G.Vec3.norm2(vel) : F32)