Lines Matching refs:nw1
2574 integer :: r2mnunit, nw1, nw2, nkp, nn local
2581 do nw1 = 1, num_wann
2585 if (nw1 .eq. nw2) delta = 1.0_dp
2590 (2.0_dp*delta - real(m_matrix(nw1, nw2, nn, nkp) + &
2591 conjg(m_matrix(nw2, nw1, nn, nkp)), kind=dp))
2596 write (r2mnunit, '(2i6,f20.12)') nw1, nw2, r2ave_mn
3100 integer :: nn, nw1, nw2, nw3 local
3104 loop_nw1: do nw1 = 1, num_wann
3105 loop_nw2: do nw2 = nw1 + 1, num_wann
3109 a11 = a11 + (m_w(nw1, nw1, nn) - m_w(nw2, nw2, nn))**2
3110 a12 = a12 + m_w(nw1, nw2, nn)*(m_w(nw1, nw1, nn) - m_w(nw2, nw2, nn))
3111 a22 = a22 + m_w(nw1, nw2, nn)**2
3131 rtmp1 = m_w(nw3, nw1, nn)*cc + m_w(nw3, nw2, nn)*ss
3132 rtmp2 = -m_w(nw3, nw1, nn)*ss + m_w(nw3, nw2, nn)*cc
3133 m_w(nw3, nw1, nn) = rtmp1
3138 rtmp1 = cc*m_w(nw1, nw3, nn) + ss*m_w(nw2, nw3, nn)
3139 rtmp2 = -ss*m_w(nw1, nw3, nn) + cc*m_w(nw2, nw3, nn)
3140 m_w(nw1, nw3, nn) = rtmp1
3146 rtmp1 = ur_rot(nw3, nw1)*cc + ur_rot(nw3, nw2)*ss
3147 rtmp2 = -ur_rot(nw3, nw1)*ss + ur_rot(nw3, nw2)*cc
3148 ur_rot(nw3, nw1) = rtmp1