On 06/06/12 13:37, Paul Emsley wrote:
On 06/06/12 12:27, Nicholas Firth wrote:
Hi All,
I'm having some issues with getting some more 'interesting' ring systems to kekulize. I'm creating the molecule in RDKit by individually adding all the atoms and bonds from a graph. However I have found similar issue when trying to initialise a molecule from smiles.
Two examples of molecules which won't kekulize are:
n1nnnn1
c1cccp1
however, when adding hydrogens to these smiles, so
n1n[nH]nn1
c1ccc[pH]1
These will work perfectly fine. This is a problem as the graph structure I have does not include any hydrogens. For a fix I have tried to addHs to my molecule, however I'm having issue with these functions because I'm writing in C++ and want to return this value, so don't have any pointers as this results in a memory leak.



I guess I'm missing something. addHs() returns an ROMol - it is not clear to me why that gives you pointer issues/problems.


Argh.  Just missed the Cancel Button in time.

Yes, ROMol does indeed return a pointer. Why does that cause a problem? I guess we need to look at the details...?


Paul.


------------------------------------------------------------------------------
Live Security Virtual Conference
Exclusive live event will cover all the ways today's security and 
threat landscape has changed and how IT managers can respond. Discussions 
will include endpoint security, mobile security and the latest in malware 
threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
_______________________________________________
Rdkit-discuss mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/rdkit-discuss

Reply via email to