def _multi_digamma(self, a, p, name="<STR_LIT>"):
with self._name_scope(name):<EOL><INDENT>seq = self._multi_gamma_sequence(a, p)<EOL>return tf.reduce_sum(input_tensor=tf.math.digamma(seq), axis=[-<NUM_LIT:1>])<EOL><DEDENT>
Computes the multivariate digamma function; Psi_p(a).