[gmx-users] gromacs 4.6 segfault

2013-01-15 Thread sebastian
Hey GROMACS users, using mdrun (version 4.6-beta3) on a GPU node (1 nvidia K10 with cuda drivers and runtime 4.2 + 2 times intel 6 core E5 with hyper threading and SSE4.1) I get allways after a few or few 100 ns the following segfault: line 15: 28957 Segmentation fault mdrun -deffnm

[gmx-users] do_dssp

2013-01-15 Thread James Starlight
Dear Gromacs Users! I'd like to know if it's possible to obtain some values from do_dssp except the graph diagrams of the ss dynamics ? For example I wounder to know values of how many residues were in helix, sheet or coil conformation during my simulation ? By the way during analysis of the

[gmx-users] suffice of position restraints

2013-01-15 Thread Ahmet yıldırım
Dear users, I usually apply the positions restraint of 100 ps on system. Does it produce a problem to apply the positions restraint for long time (100 or 200 ps)? cant it cause strain on structure? Thanks in advance -- Ahmet Yıldırım -- gmx-users mailing listgmx-users@gromacs.org

Re: [gmx-users] Re: Turn on pdb output after every step

2013-01-15 Thread Bogdan Costescu
On Mon, Jan 14, 2013 at 9:23 PM, santhosh mesan...@gmail.com wrote: I want to output pdb during the run. Before every integration step, I want to output the current configuration of the system. Before an integration step is after the previous integration step, which is where the normal

Re: [gmx-users] do_dssp

2013-01-15 Thread Justin Lemkul
On 1/15/13 3:19 AM, James Starlight wrote: Dear Gromacs Users! I'd like to know if it's possible to obtain some values from do_dssp except the graph diagrams of the ss dynamics ? For example I wounder to know values of how many residues were in helix, sheet or coil conformation during my

Re: [gmx-users] suffice of position restraints

2013-01-15 Thread Justin Lemkul
On 1/15/13 3:28 AM, Ahmet yıldırım wrote: Dear users, I usually apply the positions restraint of 100 ps on system. Does it produce a problem to apply the positions restraint for long time (100 or 200 ps)? cant it cause strain on structure? Unlikely. The amount of equilibration necessary

Re: [gmx-users] MD Simulation of Calcium with protein

2013-01-15 Thread Justin Lemkul
On 1/15/13 5:49 AM, Devika N T wrote: HI I would like to know the protocol to be followed for performing MD simulation for calcium with protein (Calmodulin) Can I follow the same protocol which is followed for a protein? Probably, as long as the force field you choose can deal with

[gmx-users] 60% slowdown with GPU / verlet and sd integrator

2013-01-15 Thread Floris Buelens
Hi, I'm seeing MD simulation running a lot slower with the sd integrator than with md - ca. 10 vs. 30 ns/day for my 47000 atom system. I found no documented indication that this should be the case. Timings and logs pasted in below - wall time seems to be accumulating up in Update and Rest,

Re: [gmx-users] Re: gromacs 4.6 segfault

2013-01-15 Thread Justin Lemkul
On 1/15/13 7:06 AM, Dr. Vitaly Chaban wrote: using mdrun (version 4.6-beta3) on a GPU node (1 nvidia K10 with cuda drivers and runtime 4.2 + 2 times intel 6 core E5 with hyper threading and SSE4.1) I get allways after a few or few 100 ns the following segfault: line 15: 28957 Segmentation

Re: [gmx-users] do_dssp

2013-01-15 Thread James Starlight
Justin, I want to obtain timescale on X as well as number of residues on Y on the xmp graph . By default that graph has not legend so it's hard to analyse it . By the way have you forced with some problems with dppc ? Simetimes I've obtain error like There are 286 residues in your selected

Re: [gmx-users] Re: gromacs 4.6 segfault

2013-01-15 Thread sebastian
On 01/15/2013 01:09 PM, Justin Lemkul wrote: On 1/15/13 7:06 AM, Dr. Vitaly Chaban wrote: using mdrun (version 4.6-beta3) on a GPU node (1 nvidia K10 with cuda drivers and runtime 4.2 + 2 times intel 6 core E5 with hyper threading and SSE4.1) I get allways after a few or few 100 ns the

Re: [gmx-users] Re: gromacs 4.6 segfault

2013-01-15 Thread Dr. Vitaly Chaban
On Tue, Jan 15, 2013 at 1:09 PM, Justin Lemkul jalem...@vt.edu wrote: On 1/15/13 7:06 AM, Dr. Vitaly Chaban wrote: using mdrun (version 4.6-beta3) on a GPU node (1 nvidia K10 with cuda drivers and runtime 4.2 + 2 times intel 6 core E5 with hyper threading and SSE4.1) I get allways after a

Re: [gmx-users] do_dssp

2013-01-15 Thread Justin Lemkul
On 1/15/13 7:13 AM, James Starlight wrote: Justin, I want to obtain timescale on X as well as number of residues on Y on the xmp graph . By default that graph has not legend so it's hard to analyse it . The legends are printed in the .xvg header. I've never had a problem plotting

Re: [gmx-users] Re: gromacs 4.6 segfault

2013-01-15 Thread Justin Lemkul
On 1/15/13 7:16 AM, Dr. Vitaly Chaban wrote: On Tue, Jan 15, 2013 at 1:09 PM, Justin Lemkul jalem...@vt.edu wrote: On 1/15/13 7:06 AM, Dr. Vitaly Chaban wrote: using mdrun (version 4.6-beta3) on a GPU node (1 nvidia K10 with cuda drivers and runtime 4.2 + 2 times intel 6 core E5 with

Re: [gmx-users] Re: gromacs 4.6 segfault

2013-01-15 Thread Dr. Vitaly Chaban
On Tue, Jan 15, 2013 at 1:09 PM, Justin Lemkul jalem...@vt.edu wrote: On 1/15/13 7:06 AM, Dr. Vitaly Chaban wrote: using mdrun (version 4.6-beta3) on a GPU node (1 nvidia K10 with cuda drivers and runtime 4.2 + 2 times intel 6 core E5 with hyper threading and SSE4.1) I get allways after a

[gmx-users] Re: MD Simulation of Calcium with protein

2013-01-15 Thread Thomas Schlesier
Am 15.01.2013 12:52, schrieb gmx-users-requ...@gromacs.org: On 1/15/13 5:49 AM, Devika N T wrote: HI I would like to know the protocol to be followed for performing MD simulation for calcium with protein (Calmodulin) Can I follow the same protocol which is followed for a protein? Probably,

Re: [gmx-users] Re: MD Simulation of Calcium with protein

2013-01-15 Thread Justin Lemkul
On 1/15/13 8:12 AM, Thomas Schlesier wrote: Am 15.01.2013 12:52, schrieb gmx-users-requ...@gromacs.org: On 1/15/13 5:49 AM, Devika N T wrote: HI I would like to know the protocol to be followed for performing MD simulation for calcium with protein (Calmodulin) Can I follow the same

[gmx-users] Re: MD Simulation of Calcium with protein

2013-01-15 Thread Thomas Schlesier
Am 15.01.2013 14:17, schrieb gmx-users-requ...@gromacs.org: On 1/15/13 8:12 AM, Thomas Schlesier wrote: Am 15.01.2013 12:52, schriebgmx-users-requ...@gromacs.org: On 1/15/13 5:49 AM, Devika N T wrote: HI I would like to know the protocol to be followed for performing MD simulation for

Re: [gmx-users] do_dssp

2013-01-15 Thread James Starlight
Justin, thanks both options works perfect. James 2013/1/15 Justin Lemkul jalem...@vt.edu: On 1/15/13 7:13 AM, James Starlight wrote: Justin, I want to obtain timescale on X as well as number of residues on Y on the xmp graph . By default that graph has not legend so it's hard to

[gmx-users] Re: Umbrella sampling and free energy estimation

2013-01-15 Thread James Starlight
Restating of my question about ussage of essential eigenvectors as the pulling coordinates I'd like to introduce some example. The case study is the open to close transition in some protein. For instance by means of essential dynamics sampling technique (in targeted-contraction mode) I've

[gmx-users] Re: Turn on pdb output after every step

2013-01-15 Thread santhosh
@Bogdan : It differs from the normal gromacs setup as this PDB write is expected to happen during the run. (While the simulation is still running). I have posted this in gmx-developers post as I felt its more apt after Justins advice. -- View this message in context:

[gmx-users] Gromacs-4.6-beta3 compile warnings intel-suite 2011 and 2013

2013-01-15 Thread Richard Broadbent
Dear All, I've just installed 4.6-beta3 on my ubuntu linux (Intel Xeon [sandy bridge]) box using both intel-suite/64/2011.10/319, and intel-suite/64/2013.0/079 with mkl Using either compiler I received several hundred warnings of type #120, #167, and #556 (see bellow for examples). I

Re: [gmx-users] gromacs 4.6 segfault

2013-01-15 Thread Szilárd Páll
Hi, Please try to run the very same simulation that crashes on CPUs only (-nb cpu). If that also crashes, please file a bug report because we surely have a software bug. If that does not crash, we will have to dig a bit deeper into what could be causing the issue. Log files from the two runs

Re: [gmx-users] 60% slowdown with GPU / verlet and sd integrator

2013-01-15 Thread Szilárd Páll
Hi Floris, Great feedback, this needs to be looked into. Could you please file a bug report, preferably with a tpr (and/or all inputs) as well as log files. Thanks, -- Szilárd On Tue, Jan 15, 2013 at 3:50 AM, Floris Buelens floris_buel...@yahoo.comwrote: Hi, I'm seeing MD simulation

Re: [gmx-users] Gromacs-4.6-beta3 compile warnings intel-suite 2011 and 2013

2013-01-15 Thread Roland Schulz
Hi, could you check if you get these warnings also with the latest version from git? We have changed quite a bit sense then. git clone https://gerrit.gromacs.org/p/gromacs git checkout release-4-6 Roland On Tue, Jan 15, 2013 at 9:34 AM, Richard Broadbent richard.broadben...@imperial.ac.uk

[gmx-users] How do forces from PME get calculated?

2013-01-15 Thread Schlesier, Thomas
Maybe a somewhat dumb question: How do forces from PME get calculated? In the manual (and also some textbooks) i have only found expressions for the potential for PME (or ewald-summation), but no information how the forces get calculated. One idea i have would be dF_i = - [ V(r_i(t)) -

[gmx-users] 1,4-interactions and pairs...again

2013-01-15 Thread Robson da Silva
Hi gmx users, I'm writing a MD code and, comparing with gromacs, I have some doubts. When I put gen-pairs=yes, fudges=0.5 and no [pairs] or [pairtypes] directives, the program (gmx) would generate 1,4 pairs with the parameters described in the [atomtypes] directive, right? But, at the end of the

AW: [gmx-users] How do forces from PME get calculated?

2013-01-15 Thread Florian Dommert
Hi Thomas, There are different possibilities to derive the forces from the PME method. Usually with SPME, the so called analytical differentiation scheme is applied, where the gradient of the reciprocal sum is directly calculated. This is computationally very efficient, because it requires only