def coul_head

bend-prot@0.1.0.0 · force.bend:105

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

105def coul_head(
106  +pi: G.Vec3, +qi: F32, +ke: F32, +c2: F32, h: P.Atom, qj: F32
107) -> F32:
108  match h:
109    case P.Atom{serial, elem, +pj}:
110      coul_at(pi, pj, (qi * qj : F32), ke, c2)