Bonded.lj_row_elem_excl: def of bend-prot@0.1.0.0
Import: import bend-prot@0.1.0.0/bonded.bend as Bonded
Lines 224-233 of bend-prot@0.1.0.0/bonded.bend: https://jend.leeeet.dev/read.txt?path=bend-prot@0.1.0.0/bonded.bend&start_line=224&end_line=233

import Base
import bend-prot@0.1.0.0/geom.bend as G
import bend-prot@0.1.0.0/protein.bend as P
import bend-prot@0.1.0.0/force.bend as F

# Per-element rows/totals that also skip 1-2 and 1-3 pairs: the physically
# meaningful combination (uniform rows blow up on bonded neighbors).
def lj_row_elem_excl(+si: U32, +se: U32, +pi: G.Vec3, +ys: List<&2, P.Atom>, +c2: F32, +bonds: List<&2, U32>) -> F32:
  match ys:
    case Nil{}:
      0.0
    case h <> t:
      match h:
        case P.Atom{serial, +elem, +pj}:
          (Bool.pick(F32, excluded(si, serial, bonds), 0.0, F.lj_pair_at(pi, pj, F.lj_combine_eps(F.lj_eps(se), F.lj_eps(elem)), F.lj_combine_sig2(F.lj_sig(se), F.lj_sig(elem)), c2)) + lj_row_elem_excl(si, se, pi, t, c2, bonds) : F32)

Uses (13):
  https://jend.leeeet.dev/definition.txt?path=Base/base.bend&name=Bool.pick
  https://jend.leeeet.dev/definition.txt?path=Base/base.bend&name=F32
  https://jend.leeeet.dev/definition.txt?path=Base/base.bend&name=F32.add
  https://jend.leeeet.dev/definition.txt?path=Base/base.bend&name=List
  https://jend.leeeet.dev/definition.txt?path=Base/base.bend&name=U32
  https://jend.leeeet.dev/definition.txt?path=bend-prot@0.1.0.0/bonded.bend&name=excluded
  https://jend.leeeet.dev/definition.txt?path=bend-prot@0.1.0.0/force.bend&name=lj_combine_eps
  https://jend.leeeet.dev/definition.txt?path=bend-prot@0.1.0.0/force.bend&name=lj_combine_sig2
  https://jend.leeeet.dev/definition.txt?path=bend-prot@0.1.0.0/force.bend&name=lj_eps
  https://jend.leeeet.dev/definition.txt?path=bend-prot@0.1.0.0/force.bend&name=lj_pair_at
  https://jend.leeeet.dev/definition.txt?path=bend-prot@0.1.0.0/force.bend&name=lj_sig
  https://jend.leeeet.dev/definition.txt?path=bend-prot@0.1.0.0/geom.bend&name=Vec3
  https://jend.leeeet.dev/definition.txt?path=bend-prot@0.1.0.0/protein.bend&name=Atom