Thank you Mark for your response. I only need to double check a few things for which using generic C kernels should suffice. I need to run only 1 or 2 steps and not in parallel.<br>Thanks again for your help.<br>Sapna<br><br>
<div class="gmail_quote">On Wed, Sep 15, 2010 at 9:43 PM, Mark Abraham <span dir="ltr"><<a href="mailto:mark.abraham@anu.edu.au">mark.abraham@anu.edu.au</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div><div></div><div class="h5"><br><br>----- Original Message -----<br>From: sapna sarupria <<a href="mailto:sapna.sarupria@gmail.com" target="_blank">sapna.sarupria@gmail.com</a>><br>Date: Thursday, September 16, 2010 1:35<br>
Subject: [gmx-users] Wrtiing pair interaction energy on the fly<br>To: Discussion list for GROMACS users <<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>Dear all<br>
<font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>I understand that pair interactions in gromacs are calculated through assembly loops. Can some one guide me to the resources using which I can write out the interaction energy/force for each pair directly from gromacs during the simulation (without using the rerun option). The rerun option is especially difficult when you want pair interaction between all atoms of one's system since one has to define that many index groups. I most cases rerun is good, but in my case I have to test interaction between several pairs. <br>
<br></div></div>You are limited to 256 energy groups with or without using -rerun. If you need more, then perhaps you should consider some other simulation software that is better suited to this task. With GROMACS, you could use an environment variable to force the use of the generic C kernels, and inside that kernel use printf to write out the energy for each interaction. This will be horribly slow, and will be vastly simpler if not attempted in parallel.<br>
<font color="#888888"><br>Mark
</font><br>--<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br></blockquote></div><br><br clear="all"><br>-- <br>Sapna Sarupria<br>Post-doctoral Researcher<br>
Princeton University<br>New Jersey 08540<br>U.S.A.<br><br>Life isn't about finding yourself. Life is about creating yourself. -- George Bernard Shaw. <br>Dare to Dream<br>