<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 27, 2014 at 2:19 PM, Mendez Giraldez, Raul <span dir="ltr"><<a href="mailto:rmendez@email.unc.edu" target="_blank">rmendez@email.unc.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt">Dear Gromacs developers,
<div><br>
</div>
<div>I am trying to compute PCA on a <font face="Segoe UI, Helvetica, Arial, sans-serif">46464x46464 matrix (15488 atoms). g_covar produces segmentation fault. We thought it was a lack of physical memory, but when our computer cluster system administrator,
Shubin Liu, looked into the system logs, found that the maximum RAM used by the program was 16 Gb, while I had booked up to 800 Gb of RAM. Shubin realized that the actual problem arises at the eigensolver call within covar routine. So it looks like eigensolver
is unable to handle big matrices, but not because of lack of physical memory (is it based on LAPACK libraries, still 32-bit coded?).</font></div></div></div></blockquote><div><br></div><div>I've no idea if there are any limits, but you can choose the LAPACK implementation with which to link GROMACS (see install guide), so you should find out what you've been using (e.g. inspect CMakeCache.txt in your build tree), and explore alternatives.</div><div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt"><div><font face="Segoe UI, Helvetica, Arial, sans-serif"> I am using Gromacs </font><span style="font-size:10pt">4.6.3.</span></div>
<div><font face="Segoe UI, Helvetica, Arial, sans-serif"><br>
</font></div>
<div><font face="Segoe UI, Helvetica, Arial, sans-serif">Do you know an alternative to that? The reason I want to do that is to have some ion atoms in addition to my protein and see how they move along the low frequency modes (large variance PCs).</font></div></div></div></blockquote><div><br></div><div>How many of your protein atoms do you really need for that analysis? I presume you're already excluding the hydrogen atoms, since they're not doing much at low frequencies. Even if you think you really do need them, I suggest you start with something small that works.</div><div><br></div><div>Mark</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt">
<div><font face="Segoe UI, Helvetica, Arial, sans-serif"><br>
</font></div>
<div><font face="Segoe UI, Helvetica, Arial, sans-serif">Thank you so much,</font></div>
<div><font face="Segoe UI, Helvetica, Arial, sans-serif"><br>
</font></div>
<div><font face="Segoe UI, Helvetica, Arial, sans-serif">Raul</font></div>
<div><br>
<div style="font-family:Tahoma;font-size:13px">Raul Mendez Giraldez, PhD
<div>Dokholyan Lab</div>
<div>Dept. Biophysics & Biochemistry </div>
<div>Genetic Medicine</div>
<div>University of North Carolina</div>
<div>120 Mason Farm Road</div>
<div>Chapel Hill, NC 27599</div>
<div>US</div>
</div>
</div>
</div>
</div>
<br>--<br>
Gromacs Developers mailing list<br>
<br>
* Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List</a> before posting!<br>
<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>
<br>
* For (un)subscribe requests visit<br>
<a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-developers" target="_blank">https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-developers</a> or send a mail to <a href="mailto:gmx-developers-request@gromacs.org">gmx-developers-request@gromacs.org</a>.<br></blockquote></div><br></div></div>