def coul_at
import bend-prot@0.1.0.0/force.bend as Forceimport Base
import bend-prot@0.1.0.0/geom.bend as G
101def coul_at(+pi: G.Vec3, +pj: G.Vec3, +qq: F32, +ke: F32, +c2: F32) -> F32:
102 +d2 = G.Vec3.dist2(pi, pj)
103 coul_pair_go(d2, qq, ke, F32.is_lt(d2, c2))