<div dir="ltr"><div>Are you referring to code like the following?</div><div>        CALL zgemm( 'N', 'N', npw, nbnd_occ, nkb, &</div>             (0.d0,0.5d0), vkb(1,1), npwx, ps2(1,1,1), nkb, (1.d0,0.0d0), &<br><div>             ppsi, npwx )</div><div>It's old-style fortran: ps2 is not passed as an array but as a pointer. ps2(1,1,1) points to the first element of the array. In this particular case, it is the same as ps2 with no indices</div><div><br></div><div>Paolo<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jul 23, 2020 at 10:49 PM Lun Yue <<a href="mailto:xifntx@gmail.com">xifntx@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi QE developers,<br>
<br>
I have a question that I hope you have the time for.<br>
<br>
The momentum matrix elements "write_p_avg.f90" includes the commutator <br>
[r, V_nonlocal], which is implemented in "compute_ppsi.f90".<br>
In this module, the call to "zgemm" performs the sum between p*psi and <br>
the nonlocal part included in the array "ps2". I was wondering why only <br>
a single element of "ps2" is added, even though it is huge array <br>
calculated in the previous step (it is used in zgemm and deallocated <br>
right afterwards). Is this a bug, or am I missing something?<br>
<br>
Also, is this part of the code well tested?<br>
<br>
<br>
Best regards,<br>
<br>
Lun Yue<br>
<br>
Louisiana State University<br>
<br>
_______________________________________________<br>
developers mailing list<br>
<a href="mailto:developers@lists.quantum-espresso.org" target="_blank">developers@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/developers" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/developers</a><br>
</blockquote></div><br clear="all"><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>