def coul_frow

bend-prot@0.1.0.0 · force.bend:326

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

326def coul_frow(+pi: G.Vec3, +qi: F32, +ke: F32, +c2: F32, ys: List<&2, P.Atom>, qs: List<&2, F32>) -> G.Vec3:
327  match ys qs:
328    case Nil{} Nil{}:
329      G.V3{0.0, 0.0, 0.0}
330    case Nil{} qh <> qt:
331      G.V3{0.0, 0.0, 0.0}
332    case h <> t Nil{}:
333      G.V3{0.0, 0.0, 0.0}
334    case h <> t qh <> qt:
335      G.Vec3.add(coul_frow(pi, qi, ke, c2, t, qt), coul_fhead(pi, qi, ke, c2, h, qh))