[ghemical-devel] Duplicate Coordinates

Donald Ephraim Curtis donald-curtis at uiowa.edu
Fri Dec 15 18:56:52 EST 2006


What about just eliminating the divide by zero?  Anyone have any idea
where this is located in eng1_mm_tripos52.cpp?  I guess I'll look
through the code at some point this weekend.  I'm sure at some point i
can figure it out.  Thanks for the heads up Geoff.

-Donald

(Fri, Dec 15, 2006 at 06:38:01PM -0500) Geoffrey Hutchison <geoff at geoffhutchison.net>:
> 
> On Dec 15, 2006, at 4:41 PM, Donald Ephraim Curtis wrote:
> 
> >Before I add the duplicate atom coordinate "check" i wanted to make  
> >sure
> >we couldn't fix this problem in the libghemical tools.
> 
> I think this should ideally be fixed in the underlying libghemical.  
> Basically, there's a divide by zero and that kills everything. You  
> can see the same effect if you import a SMILES string (which will  
> have 0.0 for all coordinates).
> 
> So I think libghemical needs a quick check for coordinates which are  
> too close. (For example, sort by z-axis, and check for distances  
> under 0.4Ang / 40pm.)
> 
> Cheers,
> -Geoff


More information about the ghemical-devel mailing list