def lj_forces

bend-prot@0.1.0.0 · force.bend:198

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

198def lj_forces(
199  xs: List<&2, P.Atom>, +ys: List<&2, P.Atom>, +eps: F32, +sig2: F32, +c2: F32
200) -> List<&2, G.Vec3>:
201  match xs:
202    case Nil{}:
203      Nil{}
204    case h <> t:
205      lj_fself_atom(h, ys, eps, sig2, c2) <> lj_forces(t, ys, eps, sig2, c2)