'Derivative of the variance function v\'(mu)
May be undefined at zero.'

def deriv(self, mu):

    der = (self.power * (np.fabs(mu) ** (self.power - 1)))
    ii = np.flatnonzero((mu < 0))
    der[ii] *= (-1)
    return der
