def coul_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
321def coul_fhead(+pi: G.Vec3, +qi: F32, +ke: F32, +c2: F32, h: P.Atom, qj: F32) -> G.Vec3:
322 match h:
323 case P.Atom{serial, elem, +pj}:
324 coul_fat(pi, pj, (qi * qj : F32), ke, c2)