def lj_frow

bend-prot@0.1.0.0 · force.bend:182

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

182def lj_frow(
183  +pi: G.Vec3, +eps: F32, +sig2: F32, +c2: F32, ys: List<&2, P.Atom>
184) -> G.Vec3:
185  match ys:
186    case Nil{}:
187      G.V3{0.0, 0.0, 0.0}
188    case h <> t:
189      G.Vec3.add(lj_frow(pi, eps, sig2, c2, t), lj_fhead(pi, eps, sig2, c2, h))