'Returns the logarithm of the determinant of a
factor-structured matrix.'

def logdet(self):

    logdet = np.sum(np.log(self.diag))
    logdet += np.sum(np.log(self.scales))
    logdet += np.sum(np.log((1 + (1 / self.scales))))
    return logdet
