def coul_row

bend-prot@0.1.0.0 · force.bend:112

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

112def coul_row(
113  +pi: G.Vec3, +qi: F32, +ke: F32, +c2: F32,
114  ys: List<&2, P.Atom>, qs: List<&2, F32>
115) -> F32:
116  match ys qs:
117    case Nil{} Nil{}:
118      0.0
119    case Nil{} qh <> qt:
120      0.0
121    case h <> t Nil{}:
122      0.0
123    case h <> t qh <> qt:
124      (coul_row(pi, qi, ke, c2, t, qt) + coul_head(pi, qi, ke, c2, h, qh) : F32)