Add invdigamma function and tests