def bond_forces

bend-prot@0.1.0.0 · bonded.bend:308

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

308def bond_forces(+xs: List<&2, P.Atom>, +bonds: List<&2, U32>, +k: F32) -> List<&2, G.Vec3>:
309  match xs:
310    case Nil{}:
311      Nil{}
312    case +h <> t:
313      bond_force_list(P.Atom.serial(h), bonds, xs, k) <> bond_forces(t, bonds, k)