implement flu_d() and flu_vd(), closes gh-1