def METHOD_NAME(self, residualFn=None):
    if residualFn is not None:
        return residualFn(self.var, self.matrix, self.RHSvector)
    else:
        return numerix.L2norm(self._calcResidualVector())