def coul_fself_atom
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
337def coul_fself_atom(h: P.Atom, +qh: F32, +ys: List<&2, P.Atom>, +qs_e: List<&2, F32>, +ke: F32, +c2: F32) -> G.Vec3:
338 match h:
339 case P.Atom{serial, elem, +pi}:
340 coul_frow(pi, qh, ke, c2, ys, qs_e)