def sum_vels_ds
import bend-prot@0.1.0.0/sim.bend as Simimport Base
import bend-prot@0.1.0.0/geom.bend as G
81def sum_vels_ds(ds: List<&2, DynAtom>) -> G.Vec3:
82 match ds:
83 case Nil{}:
84 G.V3{0.0, 0.0, 0.0}
85 case h <> t:
86 match h:
87 case DA{serial, elem, pos, vel, frc, mass}:
88 G.Vec3.add(vel, sum_vels_ds(t))