<div dir="ltr"><div dir="ltr">On Wed, Sep 4, 2019 at 4:29 PM Aritz Leonardo <aritz.leonardo@ehu.eus> wrote:</div><div dir="ltr"><br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
  

    
  
  <div bgcolor="#FFFFFF">If I print on the screen for a check , the values of evc and
      evc(igk_k) corresponding to the first 4 planewaves, (1st band )<br>
    
    <p>   WRITE(*,*) evc(1,1), evc(igk_k(1,ik),1)<br></p></div></blockquote><div>no, no!  evc(i,n) = i-th component of the n-th band; the i-th component corrsponds to (k+G)(i) = k(ik)+G(igk_k(i,ik))   where ik is the index of k-points. G is the array of G-vectors<br></div><div><br></div><div>Paolo<br></div></div><br>-- <br><div dir="ltr" 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>