\Mh dZddlZddlmZdgZededejddZdS) z5Bethe Hessian or deformed Laplacian matrix of graphs.N)not_implemented_forbethe_hessian_matrixdirected multigraphc Nddl}|t|}|Ztdtj|Dtdtj|Dz dz }tj||d}|j\}}|j|j |dd||d }|j|j ||d }|d zdz |z||zz |zS) uReturns the Bethe Hessian matrix of G. The Bethe Hessian is a family of matrices parametrized by r, defined as H(r) = (r^2 - 1) I - r A + D where A is the adjacency matrix, D is the diagonal matrix of node degrees, and I is the identify matrix. It is equal to the graph laplacian when the regularizer r = 1. The default choice of regularizer should be the ratio [2]_ .. math:: r_m = \left(\sum k_i \right)^{-1}\left(\sum k_i^2 \right) - 1 Parameters ---------- G : Graph A NetworkX graph r : float Regularizer parameter nodelist : list, optional The rows and columns are ordered according to the nodes in nodelist. If nodelist is None, then the ordering is produced by ``G.nodes()``. Returns ------- H : scipy.sparse.csr_array The Bethe Hessian matrix of `G`, with parameter `r`. Examples -------- >>> k = [3, 2, 2, 1, 0] >>> G = nx.havel_hakimi_graph(k) >>> H = nx.bethe_hessian_matrix(G) >>> H.toarray() array([[ 3.5625, -1.25 , -1.25 , -1.25 , 0. ], [-1.25 , 2.5625, -1.25 , 0. , 0. ], [-1.25 , -1.25 , 2.5625, 0. , 0. ], [-1.25 , 0. , 0. , 1.5625, 0. ], [ 0. , 0. , 0. , 0. , 0.5625]]) See Also -------- bethe_hessian_spectrum adjacency_matrix laplacian_matrix References ---------- .. [1] A. Saade, F. Krzakala and L. Zdeborová "Spectral Clustering of Graphs with the Bethe Hessian", Advances in Neural Information Processing Systems, 2014. .. [2] C. M. Le, E. Levina "Estimating the number of communities in networks by spectral methods" arXiv:1507.00827, 2015. rNc3&K|] \}}|dzV dS)N.0vds b/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/networkx/linalg/bethehessianmatrix.py z'bethe_hessian_matrix..Hs*..A1......c3 K|] \}}|V dS)Nr r s rrz'bethe_hessian_matrix..Hs&4P4P41aQ4P4P4P4P4P4Prcsr)nodelistformat)axis)rr ) scipylistsumnxdegreeto_scipy_sparse_arrayshapesparse csr_arrayspdiagseye) GrrspAnmDIs rrr st77y ..1... . .4P4P29Q<<4P4P4P1P1P PST T  XeDDDA 7DAq BI--aeeemmQ1U-SSTTA BIMM!QuM==>>A qD1H>AE !A %%r)NN)__doc__networkxrnetworkx.utilsr__all__ _dispatchablerr rrr0s;;...... ! "Z  \""C&C&C&#"! C&C&C&r