<div dir="ltr"><div><div><div>Dear Stefano,<br><br></div>thank you so much. This is clear to me now.<br><br></div>Best,<br></div>Fan<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Nov 22, 2015 at 7:59 PM, stefano de gironcoli <span dir="ltr"><<a href="mailto:degironc@sissa.it" target="_blank">degironc@sissa.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
  
    
  
  <div text="#000000" bgcolor="#FFFFFF">
    <div>Dear Fan,<span class=""><br>
      <br>
      On 22/11/2015 23:59, Fan Zheng wrote:<br>
    </span></div><span class="">
    <blockquote type="cite">
      <div dir="ltr">
        <div>
          <div>
            <div>
              <div>
                <div>Dear all,<br>
                  <br>
                </div>
                I have a question about the electron phonon coupling
                calculation in the subroutine "elphel" (in the file
                elphonon.f90).<br>
                <br>
              </div>
              <div>This subroutine computes <n,k|dV|m,k+q>.
                |n,k> is the wavefunction for band "m" with
                wavevector "k".<br>
                <br>
              </div>
              <div>First, the program computes <n,k|dV.  Then, it is
                multiplied by |m,k+q><br>
              </div>
              <div><br>
              </div>
              <div>If q=0, it only needs to read to evc as:<br>
                ...<br>
                        IF (lgamma) THEN<br>
                           CALL get_buffer(evc, lrwfc, iuwfc, ikk)<br>
                        ELSE<br>
                           CALL get_buffer (evc, lrwfc, iuwfc, ikk)<br>
                           CALL get_buffer (evq, lrwfc, iuwfc, ikq)<br>
                        ENDIF<br>
                ...<br>
              </div>
              My question is in the q=0 case, array "evq" is used to
              multiply the computed <n,k|dV instead of evc. (line
              401-406 in version 5.2.1)<br>
              <br>
            </div>
            Can anyone explain  why "evc" array is not used in q=0 case
            for these line (401-406)? <br>
            <br>
          </div>
        </div>
      </div>
    </blockquote>
    <br></span>
    because in allocate_phq there is the following bit<br>
     if (lgamma) then<br>
         !<br>
         !  q=0  : evq and igkq are pointers to evc and igk<br>
         !<br>
         evq  => evc<br>
         igkq => igk<br>
      else<br>
         !<br>
         !  q!=0 : evq, igkq are allocated and calculated at point k+q<br>
         !<br>
         allocate (evq ( npwx*npol , nbnd))<br>
         allocate (igkq ( npwx))<br>
      endif<br>
    <br>
    therefore in the lgamma case evc and evq (and igkq and igk) are the
    same.<br>
    <br>
    best <br>
     stefano<br>
    <br>
    <br>
    <blockquote type="cite">
      <div dir="ltr">
        <div>Thanks,<br>
        </div>
        Fan<br>
      </div>
      <br>
      <fieldset></fieldset>
      <br>
      <pre>_______________________________________________
Pw_forum mailing list
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a></pre>
    </blockquote>
    <br>
  </div>

<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br></div>