<div dir="ltr"><div>Hi users,</div><div><br></div><div>In section 7.1 of the developer's manual, the following example code is shown:</div><div><br></div><div>

        
                
                <span></span>
        
        
                <div class="gmail-page" title="Page 17">
                        <div class="gmail-section">
                                <div class="gmail-layoutArea">
                                        <div class="gmail-column">
                                                <pre><span style="font-size:12pt;font-family:"CMTT12"">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, &
</span></pre>
                                                <pre><span style="font-size:12pt;font-family:"CMTT12"">               npwx, swfc, npwx, (0.d0, 0.d0), ps, m)
   CALL mp_sum(ps,intra_bgrp_group)
</span></pre>
                                        </div>
                                </div>
                        </div>
                </div>
        
</div><div>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?</div><div><br></div><div>Best regards,</div><div>Andrew<br></div></div>