~/jendhub

def angle_e

bend-prot@0.1.0.0 · bonded.bend:89

import bend-prot@0.1.0.0/bonded.bend as Bonded
import Base
import bend-prot@0.1.0.0/geom.bend as G
import bend-prot@0.1.0.0/force.bend as F

89# Harmonic angle energy k*(theta - eq)^2 about the middle serial.
90def angle_e(+pa: G.Vec3, +pb: G.Vec3, +pc: G.Vec3, +k: F32, +eq: F32) -> F32:
91  (k * F.lj_sq((G.Vec3.angle(G.Vec3.sub(pa, pb), G.Vec3.sub(pc, pb)) - eq : F32)) : F32)

Uses (7)

Used by (1)