Double gamma function