def lj_pair_at
import bend-prot@0.1.0.0/force.bend as Forceimport Base
import bend-prot@0.1.0.0/geom.bend as G
45def lj_pair_at(+pi: G.Vec3, +pj: G.Vec3, +eps: F32, +sig2: F32, +c2: F32) -> F32:
46 +d2 = G.Vec3.dist2(pi, pj)
47 lj_pair_go(d2, eps, sig2, F32.is_lt(d2, c2))