law coul_trunc_charges

bend-prot@0.1.0.0 · LAWS.bend:264

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

264# Coulomb with no charges answers zero, however many atoms remain.
265law coul_trunc_charges:
266  for h: P.Atom
267  for t: List<&2, P.Atom>
268  for ys: List<&2, P.Atom>
269  for qs: List<&2, F32>
270  for ke: F32
271  for c2: F32
272  {F.coul_total(h <> t, Nil{}, ys, qs, ke, c2) == 0.0 : F32}