'Variance of the inverse Wishart distribution
Only valid if the degrees of freedom are greater than the dimension of
the scale matrix plus three.
Parameters
%(_doc_default_callparams)s
Returns
var : float
The variance of the distribution'

def var(self, df, scale):

    (dim, df, scale) = self._process_parameters(df, scale)
    out = self._var(dim, df, scale)
    return (_squeeze_output(out) if (out is not None) else out)
