def lj_frow_pbc

bend-prot@0.1.0.0 · pbc.bend:62

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

62def lj_frow_pbc(+pi: G.Vec3, +eps: F32, +sig2: F32, +c2: F32, +box: G.Vec3, ys: List<&2, P.Atom>) -> G.Vec3:
63  match ys:
64    case Nil{}:
65      G.V3{0.0, 0.0, 0.0}
66    case h <> t:
67      G.Vec3.add(lj_frow_pbc(pi, eps, sig2, c2, box, t), lj_fhead_pbc(pi, eps, sig2, c2, box, h))