'Parameters
dim : int
Dimension of the scale matrix
%(_doc_default_callparams)s
Notes
As this function does no argument checking, it should not be
called directly; use \'var\' instead.'

def _var(self, dim, df, scale):

    var = (scale ** 2)
    diag = scale.diagonal()
    var += np.outer(diag, diag)
    var *= df
    return var
