def coul_self

bend-prot@0.1.0.0 · force.bend:126

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

126def coul_self(
127  h: P.Atom, +qh: F32, +ys: List<&2, P.Atom>, +qs_e: List<&2, F32>,
128  +ke: F32, +c2: F32
129) -> F32:
130  match h:
131    case P.Atom{serial, elem, +pi}:
132      coul_row(pi, qh, ke, c2, ys, qs_e)