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

Reply via email to