layers.py 文件源码

def solve(self):

        DI = np.diag_indices_from(self.SII)

        D = self.SII[DI]
        #for i in xrange(1, self.SII.shape[0]): self.SII[i:,i-1] = self.SII[i-1,i:]

        self.SII[DI] += np.mean(D)
        _, self.W, _ = sposv(self.SII, self.SIO, overwrite_a=True, overwrite_b=False) #, lower=1)

        #self.SII[DI] = D