@@ -363,23 +363,23 @@ <h3>Order of guessing<a class="headerlink" href="#order-of-guessing" title="Link
363363< span class ="g g-Whitespace "> </ span > < span class ="mi "> 165</ span > < span class ="c1 "> # check if the topology already has the attribute to partially guess it</ span >
364364< span class ="g g-Whitespace "> </ span > < span class ="mi "> 166</ span > < span class ="k "> if</ span > < span class ="nb "> hasattr</ span > < span class ="p "> (</ span > < span class ="bp "> self</ span > < span class ="o "> .</ span > < span class ="n "> _universe</ span > < span class ="o "> .</ span > < span class ="n "> atoms</ span > < span class ="p "> ,</ span > < span class ="n "> attr_to_guess</ span > < span class ="p "> )</ span > < span class ="ow "> and</ span > < span class ="ow "> not</ span > < span class ="n "> force_guess</ span > < span class ="p "> :</ span >
365365
366- < span class ="nn "> File ~/micromamba/envs/mda-user-guide/lib/python3.12/site-packages/MDAnalysis/guesser/default_guesser.py:438 ,</ span > in < span class ="ni "> DefaultGuesser.guess_bonds</ span > < span class ="nt "> (self, atoms, coords)</ span >
367- < span class ="g g-Whitespace "> </ span > < span class ="mi "> 436 </ span > < span class ="c1 "> # check that all types have a defined vdw</ span >
368- < span class ="g g-Whitespace "> </ span > < span class ="mi "> 437 </ span > < span class ="k "> if</ span > < span class ="ow "> not</ span > < span class ="nb "> all</ span > < span class ="p "> (</ span > < span class ="n "> val</ span > < span class ="ow "> in</ span > < span class ="n "> vdwradii</ span > < span class ="k "> for</ span > < span class ="n "> val</ span > < span class ="ow "> in</ span > < span class ="nb "> set</ span > < span class ="p "> (</ span > < span class ="n "> atomtypes</ span > < span class ="p "> )):</ span >
369- < span class ="ne "> --> </ span > < span class ="mi "> 438 </ span > < span class ="k "> raise</ span > < span class ="ne "> ValueError</ span > < span class ="p "> (</ span >
370- < span class ="g g-Whitespace "> </ span > < span class ="mi "> 439 </ span > < span class ="p "> (</ span >
371- < span class ="g g-Whitespace "> </ span > < span class ="mi "> 440 </ span > < span class ="s2 "> "vdw radii for types: "</ span >
372- < span class ="g g-Whitespace "> </ span > < span class ="mi "> 441 </ span > < span class ="o "> +</ span > < span class ="s2 "> ", "</ span > < span class ="o "> .</ span > < span class ="n "> join</ span > < span class ="p "> (</ span >
373- < span class ="g g-Whitespace "> </ span > < span class ="mi "> 442 </ span > < span class ="p "> [</ span > < span class ="n "> t</ span > < span class ="k "> for</ span > < span class ="n "> t</ span > < span class ="ow "> in</ span > < span class ="nb "> set</ span > < span class ="p "> (</ span > < span class ="n "> atomtypes</ span > < span class ="p "> )</ span > < span class ="k "> if</ span > < span class ="n "> t</ span > < span class ="ow "> not</ span > < span class ="ow "> in</ span > < span class ="n "> vdwradii</ span > < span class ="p "> ]</ span >
374- < span class ="g g-Whitespace "> </ span > < span class ="mi "> 443 </ span > < span class ="p "> )</ span >
375- < span class ="g g-Whitespace "> </ span > < span class ="mi "> 444 </ span > < span class ="o "> +</ span > < span class ="s2 "> ". These can be defined manually using the"</ span >
376- < span class ="g g-Whitespace "> </ span > < span class ="mi "> 445 </ span > < span class ="o "> +</ span > < span class ="sa "> f</ span > < span class ="s2 "> " keyword 'vdwradii'"</ span >
377- < span class ="g g-Whitespace "> </ span > < span class ="mi "> 446 </ span > < span class ="p "> )</ span >
378- < span class ="g g-Whitespace "> </ span > < span class ="mi "> 447 </ span > < span class ="p "> )</ span >
379- < span class ="g g-Whitespace "> </ span > < span class ="mi "> 449 </ span > < span class ="n "> lower_bound</ span > < span class ="o "> =</ span > < span class ="bp "> self</ span > < span class ="o "> .</ span > < span class ="n "> _kwargs</ span > < span class ="o "> .</ span > < span class ="n "> get</ span > < span class ="p "> (</ span > < span class ="s2 "> "lower_bound"</ span > < span class ="p "> ,</ span > < span class ="mf "> 0.1</ span > < span class ="p "> )</ span >
380- < span class ="g g-Whitespace "> </ span > < span class ="mi "> 451 </ span > < span class ="n "> box</ span > < span class ="o "> =</ span > < span class ="bp "> self</ span > < span class ="o "> .</ span > < span class ="n "> _kwargs</ span > < span class ="o "> .</ span > < span class ="n "> get</ span > < span class ="p "> (</ span > < span class ="s2 "> "box"</ span > < span class ="p "> ,</ span > < span class ="kc "> None</ span > < span class ="p "> )</ span >
381-
382- < span class ="ne "> ValueError</ span > : vdw radii for types: 29, 30, 10 , 1, 52 , 7, 57, 21, 26, 23, 5, 20, 3, 51, 55, 70, 72, 24, 56, 81, 9, 25, 54, 2, 22, 73, 6, 50, 31, 32 . These can be defined manually using the keyword 'vdwradii'
366+ < span class ="nn "> File ~/micromamba/envs/mda-user-guide/lib/python3.12/site-packages/MDAnalysis/guesser/default_guesser.py:441 ,</ span > in < span class ="ni "> DefaultGuesser.guess_bonds</ span > < span class ="nt "> (self, atoms, coords)</ span >
367+ < span class ="g g-Whitespace "> </ span > < span class ="mi "> 439 </ span > < span class ="c1 "> # check that all types have a defined vdw</ span >
368+ < span class ="g g-Whitespace "> </ span > < span class ="mi "> 440 </ span > < span class ="k "> if</ span > < span class ="ow "> not</ span > < span class ="nb "> all</ span > < span class ="p "> (</ span > < span class ="n "> val</ span > < span class ="ow "> in</ span > < span class ="n "> vdwradii</ span > < span class ="k "> for</ span > < span class ="n "> val</ span > < span class ="ow "> in</ span > < span class ="nb "> set</ span > < span class ="p "> (</ span > < span class ="n "> atomtypes</ span > < span class ="p "> )):</ span >
369+ < span class ="ne "> --> </ span > < span class ="mi "> 441 </ span > < span class ="k "> raise</ span > < span class ="ne "> ValueError</ span > < span class ="p "> (</ span >
370+ < span class ="g g-Whitespace "> </ span > < span class ="mi "> 442 </ span > < span class ="p "> (</ span >
371+ < span class ="g g-Whitespace "> </ span > < span class ="mi "> 443 </ span > < span class ="s2 "> "vdw radii for types: "</ span >
372+ < span class ="g g-Whitespace "> </ span > < span class ="mi "> 444 </ span > < span class ="o "> +</ span > < span class ="s2 "> ", "</ span > < span class ="o "> .</ span > < span class ="n "> join</ span > < span class ="p "> (</ span >
373+ < span class ="g g-Whitespace "> </ span > < span class ="mi "> 445 </ span > < span class ="p "> [</ span > < span class ="n "> t</ span > < span class ="k "> for</ span > < span class ="n "> t</ span > < span class ="ow "> in</ span > < span class ="nb "> set</ span > < span class ="p "> (</ span > < span class ="n "> atomtypes</ span > < span class ="p "> )</ span > < span class ="k "> if</ span > < span class ="n "> t</ span > < span class ="ow "> not</ span > < span class ="ow "> in</ span > < span class ="n "> vdwradii</ span > < span class ="p "> ]</ span >
374+ < span class ="g g-Whitespace "> </ span > < span class ="mi "> 446 </ span > < span class ="p "> )</ span >
375+ < span class ="g g-Whitespace "> </ span > < span class ="mi "> 447 </ span > < span class ="o "> +</ span > < span class ="s2 "> ". These can be defined manually using the"</ span >
376+ < span class ="g g-Whitespace "> </ span > < span class ="mi "> 448 </ span > < span class ="o "> +</ span > < span class ="sa "> f</ span > < span class ="s2 "> " keyword 'vdwradii'"</ span >
377+ < span class ="g g-Whitespace "> </ span > < span class ="mi "> 449 </ span > < span class ="p "> )</ span >
378+ < span class ="g g-Whitespace "> </ span > < span class ="mi "> 450 </ span > < span class ="p "> )</ span >
379+ < span class ="g g-Whitespace "> </ span > < span class ="mi "> 452 </ span > < span class ="n "> lower_bound</ span > < span class ="o "> =</ span > < span class ="bp "> self</ span > < span class ="o "> .</ span > < span class ="n "> _kwargs</ span > < span class ="o "> .</ span > < span class ="n "> get</ span > < span class ="p "> (</ span > < span class ="s2 "> "lower_bound"</ span > < span class ="p "> ,</ span > < span class ="mf "> 0.1</ span > < span class ="p "> )</ span >
380+ < span class ="g g-Whitespace "> </ span > < span class ="mi "> 454 </ span > < span class ="n "> box</ span > < span class ="o "> =</ span > < span class ="bp "> self</ span > < span class ="o "> .</ span > < span class ="n "> _kwargs</ span > < span class ="o "> .</ span > < span class ="n "> get</ span > < span class ="p "> (</ span > < span class ="s2 "> "box"</ span > < span class ="p "> ,</ span > < span class ="kc "> None</ span > < span class ="p "> )</ span >
381+
382+ < span class ="ne "> ValueError</ span > : vdw radii for types: 5, 20, 72 , 1, 73 , 7, 81, 29, 3, 30, 54, 22, 70, 31, 9, 25, 51, 56, 52, 57, 6, 10, 23, 32, 26, 24, 21, 55, 50, 2 . These can be defined manually using the keyword 'vdwradii'
383383</ pre > </ div >
384384</ div >
385385< p > However, the snippet below will re-guess the types, and now bond-guessing can work as the elements have vdW radii defined:</ p >
0 commit comments