def lj_total_elem_excl
import bend-prot@0.1.0.0/bonded.bend as Bondedimport Base
import bend-prot@0.1.0.0/protein.bend as P
235def lj_total_elem_excl(xs: List<&2, P.Atom>, +ys: List<&2, P.Atom>, +c2: F32, +bonds: List<&2, U32>) -> F32:
236 match xs:
237 case Nil{}:
238 0.0
239 case h <> t:
240 match h:
241 case P.Atom{serial, elem, pos}:
242 (lj_total_elem_excl(t, ys, c2, bonds) + lj_row_elem_excl(serial, elem, pos, ys, c2, bonds) : F32)