law coul_forces_cons

bend-prot@0.1.0.0 · LAWS.bend:689

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

689# Coulomb forces over cons cells split head from tail.
690law coul_forces_cons:
691  for h: P.Atom
692  for t: List<&2, P.Atom>
693  for qh: F32
694  for qt: List<&2, F32>
695  for ys: List<&2, P.Atom>
696  for qs_e: List<&2, F32>
697  for ke: F32
698  for c2: F32
699  {F.coul_forces(h <> t, qh <> qt, ys, qs_e, ke, c2) == F.coul_fself_atom(h, qh, ys, qs_e, ke, c2) <> F.coul_forces(t, qt, ys, qs_e, ke, c2) : List<&2, G.Vec3>}