@@ -97,11 +97,7 @@ <h1>Source code for pymatgen.analysis.diffusion.aimd.clustering</h1><div class="
9797< span class ="sd "> provide the index for each point, and ss in the final sum squared</ span >
9898< span class ="sd "> distances.</ span >
9999< span class ="sd "> """</ span >
100- < span class ="n "> centroids</ span > < span class ="o "> =</ span > < span class ="p "> (</ span >
101- < span class ="n "> np</ span > < span class ="o "> .</ span > < span class ="n "> array</ span > < span class ="p "> (</ span > < span class ="n "> random</ span > < span class ="o "> .</ span > < span class ="n "> sample</ span > < span class ="p "> (</ span > < span class ="nb "> list</ span > < span class ="p "> (</ span > < span class ="n "> points</ span > < span class ="p "> ),</ span > < span class ="n "> k</ span > < span class ="p "> ))</ span >
102- < span class ="k "> if</ span > < span class ="n "> initial_centroids</ span > < span class ="ow "> is</ span > < span class ="kc "> None</ span >
103- < span class ="k "> else</ span > < span class ="n "> initial_centroids</ span >
104- < span class ="p "> )</ span >
100+ < span class ="n "> centroids</ span > < span class ="o "> =</ span > < span class ="n "> np</ span > < span class ="o "> .</ span > < span class ="n "> array</ span > < span class ="p "> (</ span > < span class ="n "> random</ span > < span class ="o "> .</ span > < span class ="n "> sample</ span > < span class ="p "> (</ span > < span class ="nb "> list</ span > < span class ="p "> (</ span > < span class ="n "> points</ span > < span class ="p "> ),</ span > < span class ="n "> k</ span > < span class ="p "> ))</ span > < span class ="k "> if</ span > < span class ="n "> initial_centroids</ span > < span class ="ow "> is</ span > < span class ="kc "> None</ span > < span class ="k "> else</ span > < span class ="n "> initial_centroids</ span >
105101
106102 < span class ="c1 "> # Initialize book keeping vars.</ span >
107103 < span class ="n "> iterations</ span > < span class ="o "> =</ span > < span class ="mi "> 0</ span >
@@ -243,9 +239,7 @@ <h1>Source code for pymatgen.analysis.diffusion.aimd.clustering</h1><div class="
243239 < span class ="k "> if</ span > < span class ="nb "> len</ span > < span class ="p "> (</ span > < span class ="n "> ind</ span > < span class ="p "> )</ span > < span class ="o "> ></ span > < span class ="mi "> 0</ span > < span class ="p "> :</ span >
244240 < span class ="n "> c</ span > < span class ="o "> =</ span > < span class ="n "> np</ span > < span class ="o "> .</ span > < span class ="n "> zeros</ span > < span class ="p "> (</ span > < span class ="n "> n</ span > < span class ="p "> )</ span >
245241 < span class ="k "> for</ span > < span class ="n "> j</ span > < span class ="ow "> in</ span > < span class ="n "> ind</ span > < span class ="p "> :</ span >
246- < span class ="n "> dist</ span > < span class ="p "> ,</ span > < span class ="n "> image</ span > < span class ="o "> =</ span > < span class ="bp "> self</ span > < span class ="o "> .</ span > < span class ="n "> lattice</ span > < span class ="o "> .</ span > < span class ="n "> get_distance_and_image</ span > < span class ="p "> (</ span >
247- < span class ="n "> centroids</ span > < span class ="p "> [</ span > < span class ="n "> i</ span > < span class ="p "> ],</ span > < span class ="n "> points</ span > < span class ="p "> [</ span > < span class ="n "> j</ span > < span class ="p "> ]</ span >
248- < span class ="p "> )</ span >
242+ < span class ="n "> dist</ span > < span class ="p "> ,</ span > < span class ="n "> image</ span > < span class ="o "> =</ span > < span class ="bp "> self</ span > < span class ="o "> .</ span > < span class ="n "> lattice</ span > < span class ="o "> .</ span > < span class ="n "> get_distance_and_image</ span > < span class ="p "> (</ span > < span class ="n "> centroids</ span > < span class ="p "> [</ span > < span class ="n "> i</ span > < span class ="p "> ],</ span > < span class ="n "> points</ span > < span class ="p "> [</ span > < span class ="n "> j</ span > < span class ="p "> ])</ span >
249243 < span class ="n "> c</ span > < span class ="o "> +=</ span > < span class ="n "> points</ span > < span class ="p "> [</ span > < span class ="n "> j</ span > < span class ="p "> ]</ span > < span class ="o "> +</ span > < span class ="n "> image</ span >
250244 < span class ="n "> c</ span > < span class ="o "> /=</ span > < span class ="nb "> len</ span > < span class ="p "> (</ span > < span class ="n "> ind</ span > < span class ="p "> )</ span >
251245 < span class ="n "> c</ span > < span class ="o "> =</ span > < span class ="n "> np</ span > < span class ="o "> .</ span > < span class ="n "> mod</ span > < span class ="p "> (</ span > < span class ="n "> c</ span > < span class ="p "> ,</ span > < span class ="mi "> 1</ span > < span class ="p "> )</ span >
@@ -272,10 +266,7 @@ <h1>Source code for pymatgen.analysis.diffusion.aimd.clustering</h1><div class="
272266 < span class ="k "> return</ span > < span class ="kc "> True</ span >
273267 < span class ="k "> if</ span > < span class ="n "> old_centroids</ span > < span class ="ow "> is</ span > < span class ="kc "> None</ span > < span class ="p "> :</ span >
274268 < span class ="k "> return</ span > < span class ="kc "> False</ span >
275- < span class ="k "> return</ span > < span class ="nb "> all</ span > < span class ="p "> (</ span >
276- < span class ="n "> np</ span > < span class ="o "> .</ span > < span class ="n "> allclose</ span > < span class ="p "> (</ span > < span class ="n "> pbc_diff</ span > < span class ="p "> (</ span > < span class ="n "> c1</ span > < span class ="p "> ,</ span > < span class ="n "> c2</ span > < span class ="p "> ),</ span > < span class ="p "> [</ span > < span class ="mi "> 0</ span > < span class ="p "> ,</ span > < span class ="mi "> 0</ span > < span class ="p "> ,</ span > < span class ="mi "> 0</ span > < span class ="p "> ])</ span >
277- < span class ="k "> for</ span > < span class ="n "> c1</ span > < span class ="p "> ,</ span > < span class ="n "> c2</ span > < span class ="ow "> in</ span > < span class ="nb "> zip</ span > < span class ="p "> (</ span > < span class ="n "> old_centroids</ span > < span class ="p "> ,</ span > < span class ="n "> centroids</ span > < span class ="p "> )</ span >
278- < span class ="p "> )</ span > </ div >
269+ < span class ="k "> return</ span > < span class ="nb "> all</ span > < span class ="p "> (</ span > < span class ="n "> np</ span > < span class ="o "> .</ span > < span class ="n "> allclose</ span > < span class ="p "> (</ span > < span class ="n "> pbc_diff</ span > < span class ="p "> (</ span > < span class ="n "> c1</ span > < span class ="p "> ,</ span > < span class ="n "> c2</ span > < span class ="p "> ),</ span > < span class ="p "> [</ span > < span class ="mi "> 0</ span > < span class ="p "> ,</ span > < span class ="mi "> 0</ span > < span class ="p "> ,</ span > < span class ="mi "> 0</ span > < span class ="p "> ])</ span > < span class ="k "> for</ span > < span class ="n "> c1</ span > < span class ="p "> ,</ span > < span class ="n "> c2</ span > < span class ="ow "> in</ span > < span class ="nb "> zip</ span > < span class ="p "> (</ span > < span class ="n "> old_centroids</ span > < span class ="p "> ,</ span > < span class ="n "> centroids</ span > < span class ="p "> ))</ span > </ div >
279270</ div >
280271
281272
0 commit comments