law bond_total_cons

bend-prot@0.1.0.0 · LAWS.bend:526

import Base
import bend-prot@0.1.0.0/protein.bend as P
import bend-prot@0.1.0.0/bonded.bend as B

526# A bond total splits head bond from tail bonds.
527law bond_total_cons:
528  for xs: List<&2, P.Atom>
529  for s1: U32
530  for s2: U32
531  for t: List<&1, U32 & U32>
532  for k: F32
533  {B.bond_total(xs, (s1, s2) <> t, k) == (B.bond_cov_pair(xs, s1, s2, k) + B.bond_total(xs, t, k) : F32) : F32}