def angle_total
import bend-prot@0.1.0.0/bonded.bend as Bondedimport Base
import bend-prot@0.1.0.0/protein.bend as P
117def angle_total(+xs: List<&2, P.Atom>, triples: List<&1, U32 & U32 & U32>, +k: F32, +eq: F32) -> F32:
118 match triples:
119 case Nil{}:
120 0.0
121 case (sa, sb, sc) <> t:
122 (angle_triple(xs, sa, sb, sc, k, eq) + angle_total(xs, t, k, eq) : F32)