def bond_total
import bend-prot@0.1.0.0/bonded.bend as Bondedimport Base
import bend-prot@0.1.0.0/protein.bend as P
82def bond_total(+xs: List<&2, P.Atom>, bonds: List<&1, U32 & U32>, +k: F32) -> F32:
83 match bonds:
84 case Nil{}:
85 0.0
86 case (s1, s2) <> t:
87 (bond_cov_pair(xs, s1, s2, k) + bond_total(xs, t, k) : F32)