synctools_v0.py 文件源码

def get_couling_derivate_matrix(self, h, twist_number, s):
        if type(twist_number) == InPhase:
            G = networkx.complete_graph(self.n)
            d = np.zeros((self.n, self.n))
            for ir in range(self.n):
                ir_neigh = G.neighbors(ir)
                d[ir, ir_neigh] = 1
                d[ir, :] = d[ir, :] / np.sum(d[ir, :])
            return d
        else:
            raise Exception('Topology not compatible with state')