LAWS.coul_forces_cons: law of bend-prot@0.1.0.0
Import: import bend-prot@0.1.0.0/LAWS.bend as LAWS
Lines 689-699 of bend-prot@0.1.0.0/LAWS.bend: https://jend.leeeet.dev/read.txt?path=bend-prot@0.1.0.0/LAWS.bend&start_line=689&end_line=699

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

# Coulomb forces over cons cells split head from tail.
law coul_forces_cons:
  for h: P.Atom
  for t: List<&2, P.Atom>
  for qh: F32
  for qt: List<&2, F32>
  for ys: List<&2, P.Atom>
  for qs_e: List<&2, F32>
  for ke: F32
  for c2: F32
  {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>}

Uses (6):
  https://jend.leeeet.dev/definition.txt?path=Base/base.bend&name=F32
  https://jend.leeeet.dev/definition.txt?path=Base/base.bend&name=List
  https://jend.leeeet.dev/definition.txt?path=bend-prot@0.1.0.0/force.bend&name=coul_forces
  https://jend.leeeet.dev/definition.txt?path=bend-prot@0.1.0.0/force.bend&name=coul_fself_atom
  https://jend.leeeet.dev/definition.txt?path=bend-prot@0.1.0.0/geom.bend&name=Vec3
  https://jend.leeeet.dev/definition.txt?path=bend-prot@0.1.0.0/protein.bend&name=Atom