[Open Babel] Atom charge lost in 2.3.x SDF-to-SMILES conversion

2012-05-16 Thread Craig James
Here's another problem in 2.3.x that's not in 2.2.x. I realize these are bogus molecules, but it illustrates a serious problem. $ echo "[Fe+2][Fe-2].[Fe+2].[Fe-2]" | babel -i smi -o can [Fe-2][Fe+2].[Fe-2].[Fe+2] Everything looks good -- charges are maintained. Now go through the SDF forma

Re: [Open Babel] Periodic boundaries and atom typing large systems

2012-05-16 Thread Geoff Hutchison
> To make the OB force fields work with large amounts of atoms, > the non-bonded pairs should be detected dynamically. Well, I'm not sure if that's really a performance win. Probably these need to be *updated* every 10 steps or something. > force field calculations. I noticed that periodic bound

Re: [Open Babel] Periodic boundaries and atom typing large systems

2012-05-16 Thread Tim Vandermeersch
Hi, The problem is that all vdw and electrostatic pairs are precomputed in the current implementation. Setting cut-off values doesn't change this. To make the OB force fields work with large amounts of atoms, the non-bonded pairs should be detected dynamically. Tim On Wed, May 16, 2012 at 6:17 P

Re: [Open Babel] Periodic boundaries and atom typing large systems

2012-05-16 Thread Tom Daff
> Isn't there some way to turn down the VdW via a distance cutoff or so? (EnableCutOff and SetVdWCutoff) I've never used these but hopefully they will prevent the pairwise generation of all of the VdW terms. Curiously, this doesn't work. Enabling the cutoffs and setting to very small values still

Re: [Open Babel] Periodic boundaries and atom typing large systems

2012-05-16 Thread Geoff Hutchison
> Isn't there some way to turn down the VdW via a distance cutoff or so? > (EnableCutOff and SetVdWCutoff) I've never used these but hopefully they will > prevent the pairwise generation of all of the VdW terms. Yes, exactly. -Geoff --

Re: [Open Babel] Best approach for removing a molecule from a complex

2012-05-16 Thread inviso
Oops, meant to reply to the list. That indeed was the problem. Building a list of atoms to remove and then removing then in a separate loop worked like a charm. Thanks, Greg On Wed, May 16, 2012 at 7:16 AM, Noel O'Boyle wrote: > >> Any thoughts on why pybel.formula is not updating when I use

Re: [Open Babel] Best approach for removing a molecule from a complex

2012-05-16 Thread Noel O'Boyle
> Any thoughts on why pybel.formula is not updating when I use the > underlying OBMol.DeleteAtom? > This is a FAQ. Don't delete atoms from whether a loop using an Open Babel iterator as it invalidates the iterator. Instead, just create a list of OBAtoms, and in a separate loop delete each OBAtom i

[Open Babel] Compile OpenBabel-2.3.1 with VS2010 errors

2012-05-16 Thread matchcoder
I am a new to OpenBabel. I try to compile OpenBabel-2.3.1 with VS2010, But I get some errors,can somebody help me, thanks in advance. Error 1: alias.cpp D:\MinGW\home\matchcoder\openbabel-2.3.1\src\alias.cpp(75) : error C2001: newlin e in constant D:\MinGW\home\matchcoder\openbabel-2.3.1\src\alias

Re: [Open Babel] FP2/FP3 parser

2012-05-16 Thread Noel O'Boyle
I don't have any C code but the output is a series of hexadecimal numbers listed in reverse order (i.e the last number is for the first 32 bits, the hexadecimal number itself is not backwards or anything). If you read all of those numbers into a list, and then pass it to a function like pybel's _fi

Re: [Open Babel] Periodic boundaries and atom typing large systems

2012-05-16 Thread Noel O'Boyle
Isn't there some way to turn down the VdW via a distance cutoff or so? (EnableCutOff and SetVdWCutoff) I've never used these but hopefully they will prevent the pairwise generation of all of the VdW terms. - Noel On 15 May 2012 22:40, Tom Daff wrote: > I work with molecular crystals and was ver