def lj_total
import bend-prot@0.1.0.0/force.bend as Forceimport Base
import bend-prot@0.1.0.0/protein.bend as P
70def lj_total(
71 xs: List<&2, P.Atom>, +ys: List<&2, P.Atom>, +eps: F32, +sig2: F32, +c2: F32
72) -> F32:
73 match xs:
74 case Nil{}:
75 0.0
76 case h <> t:
77 (lj_total(t, ys, eps, sig2, c2) + lj_self_atom(h, ys, eps, sig2, c2) : F32)