AbsoluteDualityGap.py 文件源码

def _solveHyperplaneProjection(self, points):
        m, n = self.A.shape
        errors = np.zeros(m)
        for i in range(m):
            if i in self.ban_constraints:
                errors[i] = 9999999
            else:
                ai = self.A[i] / np.linalg.norm(self.A[i].T, self.normalize_c)
                bi = self.b[i] / np.linalg.norm(self.A[i].T, self.normalize_c)
                errors[i] = np.sum([ai * pt - bi for pt in points])
        minInd = np.argmin(errors)
        self.c = self.A[minInd] / np.linalg.norm(self.A[minInd].T,
                                                 self.normalize_c)
        self.c = self.c.tolist()[0]
        self.error = errors[minInd]
        self.dual = np.zeros(m)
        self.dual[minInd] = 1 / np.linalg.norm(self.A[minInd].T,
                                               self.normalize_c)
        self._solved = True
        return errors[minInd]