law coul_row_cons

bend-prot@0.1.0.0 · LAWS.bend:155

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

155# A Coulomb row splits head charge from tail charges.
156law coul_row_cons:
157  for pi: G.Vec3
158  for qi: F32
159  for ke: F32
160  for c2: F32
161  for h: P.Atom
162  for t: List<&2, P.Atom>
163  for qh: F32
164  for qt: List<&2, F32>
165  {F.coul_row(pi, qi, ke, c2, h <> t, qh <> qt) == (F.coul_row(pi, qi, ke, c2, t, qt) + F.coul_head(pi, qi, ke, c2, h, qh) : F32) : F32}