body b1, b2
vector anchor1, anchor2, direction
scalar jointMin, jointMax

let dir = relative b1.rot direction in
let del = (b2.pos + relative b2.rot anchor2) - (b1.pos + relative b1.rot anchor1) in

{ del dot dir
  del cross dir }

constraint { jointMin -inf } { jointMax inf }