Commit 06ceb878 authored by Gregor Michalicek's avatar Gregor Michalicek

Corrected wrong cotra3 to matmul translation in dos/evaldos.f90

parent 0b4b09dc
...@@ -475,18 +475,18 @@ ...@@ -475,18 +475,18 @@
OPEN (18,FILE='bands'//spin12(jspin)) OPEN (18,FILE='bands'//spin12(jspin))
ntb = minval(nevk(:)) ntb = minval(nevk(:))
kx(1) = 0.0 kx(1) = 0.0
vkr(:,1)=matmul(cell%bmat,vk(:,1)) vkr(:,1)=matmul(vk(:,1),cell%bmat)
DO k = 2, kpts%nkpt DO k = 2, kpts%nkpt
vkr(:,k)=matmul(cell%bmat,vk(:,k)) vkr(:,k)=matmul(vk(:,k),cell%bmat)
dk = ( vkr(1,k) - vkr(1,k-1) )**2 + ( vkr(2,k) - vkr(2,k-1) )**2 +& dk = (vkr(1,k)-vkr(1,k-1))**2 + (vkr(2,k)-vkr(2,k-1) )**2 + &
( vkr(3,k) - vkr(3,k-1) )**2 (vkr(3,k)-vkr(3,k-1))**2
kx(k) = kx(k-1) + sqrt(dk) kx(k) = kx(k-1) + sqrt(dk)
ENDDO ENDDO
DO i = 1, ntb DO i = 1, ntb
DO k = 1, kpts%nkpt DO k = 1, kpts%nkpt
write(18,'(2f15.9)') kx(k),ev(i,k) write(18,'(2f15.9)') kx(k),ev(i,k)
ENDDO ENDDO
ENDDO ENDDO
CLOSE (18) CLOSE (18)
ENDIF ENDIF
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment