method='Nelder-Mead', tol=1.e-12) F = sps.loggamma(X+Y*1j) M = np.abs(F) # n.b. this is phase of -log(gamma(z)), not log(gamma(z)) P = np.arctan2(-F.imag,-F.real)...
South Houses Caltech 1931a.png...
