law angle_total_cons

bend-prot@0.1.0.0 · LAWS.bend:541

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

541# An angle total splits head triple from tail triples.
542law angle_total_cons:
543  for xs: List<&2, P.Atom>
544  for sa: U32
545  for sb: U32
546  for sc: U32
547  for t: List<&1, U32 & U32 & U32>
548  for k: F32
549  for eq: F32
550  {B.angle_total(xs, (sa, sb, sc) <> t, k, eq) == (B.angle_triple(xs, sa, sb, sc, k, eq) + B.angle_total(xs, t, k, eq) : F32) : F32}