law bond_forces_cons

bend-prot@0.1.0.0 · LAWS.bend:799

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/bonded.bend as B

799# Bond forces split head atom from tail atoms.
800law bond_forces_cons:
801  for h: P.Atom
802  for t: List<&2, P.Atom>
803  for bonds: List<&2, U32>
804  for k: F32
805  {B.bond_forces(h <> t, bonds, k) == B.bond_force_list(P.Atom.serial(h), bonds, h <> t, k) <> B.bond_forces(t, bonds, k) : List<&2, G.Vec3>}