[QE-users] npw vs npwx

[Andrew Xu]

Hi users,

In section 7.1 of the developer's manual, the following example code is
shown:

COMPLEX, ALLOCATABLE :: ps(:,:), wfc(:,:), swfc(:,:)
   ALLOCATE (ps(m,m), wfc(npwx,m),swfc(npwx,m))
   CALL zgemm (’c’, ’n’, m, m, npw, (1.d0, 0.d0), wfc, &

               npwx, swfc, npwx, (0.d0, 0.d0), ps, m)
   CALL mp_sum(ps,intra_bgrp_group)

In the zgemm statement, I'm confused why npw and npwx is used, instead of
only npw or only npwx. Is this saying the number of columns of op(wfc)
(which is the number of rows of swfc) is npwx, but we only use the first
npw columns of op(wfc) (rows of swfc) because the entries in the last npwx
- npw columns of op(wfc) (rows of swfc) are all zero?

Best regards,
Andrew
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20201111/3b161492/attachment.html>