<div dir="ltr">On Sun, Oct 9, 2016 at 3:41 PM, Евгений Васильевич Тихонов <span dir="ltr"><<a target="_blank" href="mailto:e.tikhonov@physics.msu.ru">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>===================================================================<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_1)<br><br><br><br><br></div></div><br clear="all"><br>-- <br><div class="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></div>