def bond_force_list

bend-prot@0.1.0.0 · bonded.bend:299

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

299def bond_force_list(+s: U32, bs: List<&2, U32>, +xs: List<&2, P.Atom>, +k: F32) -> G.Vec3:
300  match bs:
301    case Nil{}:
302      G.V3{0.0, 0.0, 0.0}
303    case x <> Nil{}:
304      G.V3{0.0, 0.0, 0.0}
305    case x <> y <> t:
306      G.Vec3.add(bond_force_bond(xs, s, x, y, k), bond_force_list(s, t, xs, k))