def lj_head

bend-prot@0.1.0.0 · force.bend:49

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

49def lj_head(+pi: G.Vec3, +eps: F32, +sig2: F32, +c2: F32, h: P.Atom) -> F32:
50  match h:
51    case P.Atom{serial, elem, +pj}:
52      lj_pair_at(pi, pj, eps, sig2, c2)