def psicomputations(variance, lengthscale, Z, variational_posterior):
        # here are the "statistics" for psi0, psi1 and psi2
        # Produced intermediate results:
        # _psi1                NxM
        mu = variational_posterior.mean
        S = variational_posterior.variance
        gamma = variational_posterior.binary_prob

        psi0 = np.empty(mu.shape[0])
        psi0[:] = variance
        psi1 = _psi1computations(variance, lengthscale, Z, mu, S, gamma)
        psi2 = _psi2computations(variance, lengthscale, Z, mu, S, gamma)
        return psi0, psi1, psi2
Z - MxQ
dataset/ETHPy150Open SheffieldML/GPy/GPy/kern/src/psi_comp/ssrbf_psi_comp.py/psicomputations