def coul_forces

bend-prot@0.1.0.0 · force.bend:342

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

342def coul_forces(xs: List<&2, P.Atom>, +qs_o: List<&2, F32>, +ys: List<&2, P.Atom>, +qs_e: List<&2, F32>, +ke: F32, +c2: F32) -> List<&2, G.Vec3>:
343  match xs qs_o:
344    case Nil{} Nil{}:
345      Nil{}
346    case Nil{} qh <> qt:
347      Nil{}
348    case h <> t Nil{}:
349      Nil{}
350    case h <> t qh <> qt:
351      coul_fself_atom(h, qh, ys, qs_e, ke, c2) <> coul_forces(t, qt, ys, qs_e, ke, c2)