def vv_full_head

bend-prot@0.1.0.0 · sim.bend:190

import Base
import bend-prot@0.1.0.0/geom.bend as G

190def vv_full_head(h: DynAtom, dt: F32) -> DynAtom:
191  match h:
192    case DA{serial, elem, pos, vel, +frc, +mass}:
193      acc = G.Vec3.scale(((1.0 / mass : F32) * acc_conv() : F32), frc)
194      DA{serial, elem, pos, G.Vec3.add(vel, G.Vec3.scale((0.5 * dt : F32), acc)), frc, mass}