law coul_total_cons

bend-prot@0.1.0.0 · LAWS.bend:167

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

167# A Coulomb total splits head charge from tail charges.
168law coul_total_cons:
169  for h: P.Atom
170  for t: List<&2, P.Atom>
171  for qh: F32
172  for qt: List<&2, F32>
173  for ys: List<&2, P.Atom>
174  for qs_e: List<&2, F32>
175  for ke: F32
176  for c2: F32
177  {F.coul_total(h <> t, qh <> qt, ys, qs_e, ke, c2) == (F.coul_total(t, qt, ys, qs_e, ke, c2) + F.coul_self(h, qh, ys, qs_e, ke, c2) : F32) : F32}