def coul_k
import bend-prot@0.1.0.0/force.bend as Forceimport Base
295# --- Coulomb forces: F_i = ke*qi*qj/d2^1.5 * (pi - pj), repulsive for
296# like charges along +d. Self pairs (d2 == 0) and pairs beyond c2 give zero.
297def coul_k(+qq: F32, +ke: F32, +d2: F32) -> F32:
298 ((ke * qq : F32) / ((d2 * F32.sqrt(d2) : F32)) : F32)