<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">Dear Paolo<br>Thank you for reply. I didn't say that when I use 1 node each have 32 cores , I don't have any problem with run of pw1.in and pw2.in wiith nbnd=1024. But once I increase  the number of proccessors to 32 node each have 32 core this error happen.<br><br>Best regards<br><br>Ali Kazempour<br>
<br>
<br>
<br>
Fritz-Haber-Institut              fax   : ++49-30-8413 4701<br>
der Max-Planck-Gesellschaft<br>
Faradayweg 4-6                    e-mail: kazempou@fhi-berlin.mpg.de<br>
D-14 195 Berlin-Dahlem / German<br><br>--- On <b>Thu, 6/10/10, Paolo Giannozzi <i><giannozz@democritos.it></i></b> wrote:<br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><br>From: Paolo Giannozzi <giannozz@democritos.it><br>Subject: Re: [Pw_forum] error in gamma<br>To: "PWSCF Forum" <pw_forum@pwscf.org><br>Date: Thursday, June 10, 2010, 4:44 AM<br><br><div class="plainMail"><a ymailto="mailto:kazempoor@ph.iut.ac.ir" href="/mc/compose?to=kazempoor@ph.iut.ac.ir">kazempoor@ph.iut.ac.ir</a> wrote:<br><br>>      from  pzpotrf  : error #        47<br>>       problems computing cholesky decomposition<br>> <br>> This is very strange because when I set nbnd=512<br> > the pw1.in and pw2.in run without any problem .<br><br>it is not strange. The diagonalization algorithm in PWscf<br>is expected to work for a
 number of bands that is much<br>smaller than the dimension of the basis set. It usually<br>works great if you calculate occupied bands and a few<br>unoccupied bands. The more bands you add, the more you<br>are likely to run into numerical trouble. The error<br>here occurs in a Scalapack routine (pzpotrf). Finding<br>out why and how is a highly nontrivial task<br><br>Paolo<br>-- <br>Paolo Giannozzi, Democritos and University of Udine, Italy<br>_______________________________________________<br>Pw_forum mailing list<br><a ymailto="mailto:Pw_forum@pwscf.org" href="/mc/compose?to=Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br><a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br></div></blockquote></td></tr></table><br>