def lj_fat_pbc

bend-prot@0.1.0.0 · pbc.bend:52

import Base
import bend-prot@0.1.0.0/geom.bend as G
import bend-prot@0.1.0.0/force.bend as F

52def lj_fat_pbc(+pi: G.Vec3, +pj: G.Vec3, +eps: F32, +sig2: F32, +c2: F32, +box: G.Vec3) -> G.Vec3:
53  +dv = mic_vec(G.Vec3.sub(pi, pj), box)
54  +d2 = G.Vec3.norm2(dv)
55  F.lj_fgo(d2, dv, eps, sig2, F32.is_lt(d2, c2))