def sasa_of_atom

bend-prot@0.1.0.0 · sasa.bend:55

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

55def sasa_of_atom(
56  h: P.Atom, +offs: List<&2, G.Vec3>, +rad: F32, +ys: List<&2, P.Atom>,
57  +rc2: F32, +area_pt: F32
58) -> F32:
59  match h:
60    case P.Atom{+serial, elem, +pos}:
61      (area_pt * F32.from_nat(sasa_cover(pos, offs, rad, ys, rc2, serial)) : F32)