libvb/Constellation: with the Newton algorithm so fast, separate imaging from solving