Sorry guys,<br><br>this email was not for you. :-)<br><br>cheers,<br><br>G.<br><br><div class="gmail_quote">On Thu, Mar 11, 2010 at 8:12 AM, Gianluca Giovannetti <span dir="ltr"><<a href="mailto:gianluca.giovannetti@gmail.com">gianluca.giovannetti@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Ciao Massimo,<br><br>dopo aver lottato a lungo con questa storia dell` interfaccia ho scritto a dei giapponesi che usano anche loro wannier90 per tirare fuori la matrice ridotta. ho dubbi che ci sia qualcosa di sbagliato nell`interpretazione delle uscite di Wannier.<br>

ti faccio sapere.<br><br>ciao<br><font color="#888888"><br>G.</font><div><div></div><div class="h5"><br><br><div class="gmail_quote">On Wed, Mar 10, 2010 at 2:39 PM, Gianluca Giovannetti <span dir="ltr"><<a href="mailto:gianluca.giovannetti@gmail.com" target="_blank">gianluca.giovannetti@gmail.com</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Dear All,<br><br>in the file plot.F90 of Wannier i see a piece:<br><br><br>    do loop_kpt=1,total_pts<br>

       ham_kprm=cmplx_0<br>       if (index(bands_plot_mode,'s-k').ne.0) then<br>          do loop_rpt=1,nrpts<br>

             rdotk=twopi*dot_product(plot_kpoint(:,loop_kpt),irvec(:,loop_rpt))<br>             fac=exp(cmplx_i*rdotk)/real(ndegen(loop_rpt),dp)<br>             ham_kprm=ham_kprm+fac*ham_r(:,:,loop_rpt)<br>          end do<br>



       elseif (index(bands_plot_mode,'cut').ne.0) then<br>          do loop_rpt=1,nrpts_cut<br>             rdotk=twopi*dot_product(plot_kpoint(:,loop_kpt),irvec_cut(:,loop_rpt))<br>!!$[aam] check divide by ndegen?<br>



             fac=exp(cmplx_i*rdotk)<br>             ham_kprm=ham_kprm+fac*ham_r_cut(:,:,loop_rpt)<br>          end do<br>       else<br><br><br>what it is "ndegen" variable?<br><br>Once i get the H(n,m,R) in the file _hr.dat (in which are printed  irvec(:,loop_rpt) and ham_r(:,:,loop_rpt)  as coded  in hamiltonian.F90)  i usually work  postprocessing it and  i can get the correct Fourier Trasform (i can check the bands structure for example) only if i use ndegen as above (the values of this variable can be printed out with "iprint =3").<br>


<br>this variable ndegen is updated in hamiltonian.F90  at the subroutine hamiltonian_wigner_seitz(count_pts)  but i have problem to understand what it is.<br>could you give me some more details?<br><br>thank you.<br><font color="#888888"><br>


Gianluca<br><br><br><br>
</font></blockquote></div><br>
</div></div></blockquote></div><br>