Hello all, I found either bug in GetSubstructMatches() or I don't understand its behaviour. First case: a=Chem.MolFromSmiles("CC(C)=O") b=Chem.MolFromSmarts("CC(C)=[S,O]") >>> a.GetSubstructMatches(b) ((0, 1, 2, 3),)
This is exactly what I expected but when I add some hydrogens to SMARTS, function behaviour is weird to me: >>> a=Chem.MolFromSmiles("[H]CC(C)=O") >>> b=Chem.MolFromSmarts("[H]CC(C)=[S,O]") >>> a.GetSubstructMatches(b) () Is it bug or feature? Regards, Rafal ps. This what I trying to do is to set some property to specific hydrogens in molecule. In python I have something like this: nmrHshift={'HCCO':1.2, 'HC(C)O:3.7, 'HOCC':4.8 } Therefore I use this strange notation with some explicite added Hs. Is there any better way to address this problem? ------------------------------------------------------------------------------ What NetFlow Analyzer can do for you? Monitors network bandwidth and traffic patterns at an interface-level. Reveals which users, apps, and protocols are consuming the most bandwidth. Provides multi-vendor support for NetFlow, J-Flow, sFlow and other flows. Make informed decisions using capacity planning reports. https://ad.doubleclick.net/ddm/clk/305295220;132659582;e _______________________________________________ Rdkit-discuss mailing list Rdkit-discuss@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/rdkit-discuss