<div dir="ltr"><div>Paolo Umari has fixed the problem. Please try the attached patch (that includes the one I sent a few days ago)<br><br></div>Paolo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Oct 9, 2016 at 10:10 PM, Paolo Giannozzi <span dir="ltr"><<a href="mailto:p.giannozzi@gmail.com" target="_blank">p.giannozzi@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">On Sun, Oct 9, 2016 at 3:41 PM, Евгений Васильевич Тихонов <span dir="ltr"><<a href="mailto:e.tikhonov@physics.msu.ru" target="_blank">e.tikhonov@physics.msu.ru</a>></span> wrote:<br><br><div class="gmail_extra"><div class="gmail_quote"><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex" class="gmail_quote"><div dir="ltr"><div><div><div>[...] please check if GW examples are still working as intended.<br></div></div></div></div></blockquote><div><br></div><div>unfortunately they aren't working at all. The crash you get is easily fixed (see patch below) but there is another problem<br><br>Paolo<br><br>Index: GWW/head/lanczos_k.f90<br>==============================<wbr>==============================<wbr>=======<br>--- GWW/head/lanczos_k.f90 (revision 13098)<br>+++ GWW/head/lanczos_k.f90 (working copy)<br>@@ -19,7 +19,7 @@<br> USE gvect<br> USE constants, ONLY : e2, pi, tpi, fpi<br> USE cell_base, ONLY: at, alat, tpiba, omega, tpiba2<br>- USE wvfct, ONLY : g2kin, npwx, npw, nbnd<br>+ USE wvfct, ONLY : g2kin, npwx, nbnd<br> USE wavefunctions_module, ONLY : evc, psic<br> USE mp, ONLY : mp_sum, mp_barrier, mp_bcast<br> USE mp_world, ONLY : mpime, world_comm<br>@@ -26,7 +26,7 @@<br> USE gvecs, ONLY : nls, nlsm, doublegrid<br> USE g_psi_mod, ONLY : h_diag, s_diag<br> USE uspp, ONLY : vkb, nkb, okvan<br>- USE klist, ONLY : xk,igk_k<br>+ USE klist, ONLY : xk,igk_k, ngk<br> USE noncollin_module, ONLY : noncolin, npol<br> USE qpoint, ONLY : igkq<br><br>@@ -51,7 +51,7 @@<br> REAL(kind=DP), ALLOCATABLE :: c(:)<br><br><br>- INTEGER :: is,ig,ii,jj,it,ipol<br>+ INTEGER :: npw, is,ig,ii,jj,it,ipol<br> INTEGER :: iunlan<br> COMPLEX(kind=DP) :: csca<br><br>@@ -61,6 +61,8 @@<br> allocate(c(nstates))<br> allocate(spsi(npwx,nstates))<br><br>+ npw = ngk(ik)<br>+ allocate(igkq(1:npw))<br> igkq(1:npw) = igk_k(1:npw,ik)<br><br> t_out(:,:,:)=(0.d0,0.d0)<br>@@ -295,6 +297,7 @@<br><br> enddo<br><br>+ deallocate(igkq)<br> deallocate(psi_1,psi_2,psi_3)<br> deallocate(u_0,u_1)<br> deallocate(alpha,beta,gamma,n_<wbr>1)<span class="HOEnZb"><font color="#888888"><br><br><br><br><br></font></span></div></div><span class="HOEnZb"><font color="#888888"><br clear="all"><br>-- <br><div class="m_-3887732428851829641gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone <a href="tel:%2B39-0432-558216" value="+390432558216" target="_blank">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222" target="_blank">+39-0432-558222</a><br><br></div></div></div></div></div>
</font></span></div></div>
</blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</div>