def lj_fhead
import bend-prot@0.1.0.0/force.bend as Forceimport Base
import bend-prot@0.1.0.0/geom.bend as G
import bend-prot@0.1.0.0/protein.bend as P
177def lj_fhead(+pi: G.Vec3, +eps: F32, +sig2: F32, +c2: F32, h: P.Atom) -> G.Vec3:
178 match h:
179 case P.Atom{serial, elem, +pj}:
180 lj_fat(pi, pj, eps, sig2, c2)