|Subject:||Re: [ESPResSo] Pressure calculation|
|Date:||Mon, 28 Jan 2008 17:29:46 +0100|
I have a problem with the pressure calculation in Espresso.
After Torstens email from last week about the bug in the pressure
calculation with tabulated bonded potentials I checked out the latest
release. This release gives idential results to my old runs if I do not
execute the analyze pressure command, so everything I observe seems to
be related to the pressure calculation:
As soon as I execute analyze pressure I get a segmentation fault with
the preceding error message:
add_bonded_virials: WARNING: Bond type 9277 , atom 8925 unhandled, Atom
It seems that at some point while looping through the atoms and the
bonds in add_bonded_virials a pointer goes astray. But I could not
identify why at this atom and at this bond. My system has 1323 molecules
with 9 atoms each. Each molecule has 8 bonds, 7 angles and 6 dihedral
angles. Some of the bonds and angles are tabulated. The error occurs
sort of unmotivated it seems.
The debugger has just helped to identify that the error occurs in pressure.h
356 p2 = local_particles[p1->bl.e[i++]];
Could it be that it has sth to do with the change in pressure.h where
parts from add_bonded_virials have been taken out and are now executed
in calc_bonded_force (I do not know who mad that change) ?
BTW: The error is not caused by Torsten's change to calculate the
pressure for tabulated bonds, if I take that out I still get the
Has anybody observed anything similar?
MPI for Polymer Research
55128 Mainz / Germany
ESPResSo mailing list
This email was Anti Virus checked by Astaro Security Gateway. http://www.astaro.com
|[Prev in Thread]||Current Thread||[Next in Thread]|