Post by jonathan
I previously supposed that the number of lipids selected was constant to
simplify the problem.
The distance criteria gives a different number of lipids for each frame.
To end up with? a constant number of atoms in the output trajectory,
each lipid selected in a frame is written in a separate output frame.
I then end up with an output trajectory of a single lipid, with a number
of frames equal to the number of lipids selected over the whole
You seem to want a dynamic selection for each frame. Then, to split that selection into individual molecules, and then concentate a pseudo-trajectory from the set of all molecules found in that frame. Then concatenate a further pseudo-trajectory from the trajectory created from each frame.
In the fullness of time you'll be able to do the first part in one hit, but for the moment I think you need to run g_select -on, split the trajectory into frames (so you don't have to read the whole file for the frame of interest), and then apply trjconv (mode 2) with the matching index group for each frame. Then somehow split each resulting frame into molecules (e.g. a further trjconv mode 2), and re-concatenate. All of that is best done with a shell/Perl/python script calling the existing GROMACS tools. The result will be only slightly slower than a C solution, and much faster to develop.
Post by jonathan
As my analysis does not depend on time, loosing the correspondence
between the input trajectory frame and the output trajectory frame is
It becomes obvious that using the index from g_select is not appropriate
for this on large trajectories ...
I'm trying to avoid modifying the trjconv or g_select source to get this
working, so any suggestions are welcomed !
If there's no other way, I'll switch to gmx-dev list.
gmx-users mailing list??? gmx-users at gromacs.org
Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-request at gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
-------------- next part --------------
An HTML attachment was scrubbed...