def maxforce2_of

bend-prot@0.1.0.0 · sim.bend:242

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

242def maxforce2_of(ds: List<&2, DynAtom>) -> F32:
243  match ds:
244    case Nil{}:
245      0.0
246    case h <> t:
247      match h:
248        case DA{serial, elem, pos, vel, frc, mass}:
249          F32.max(G.Vec3.norm2(frc), maxforce2_of(t))