Added ohdev() and ohdev_phase() numpy equivalents