Commit 01bc2c35 authored by Daniel Wortmann's avatar Daniel Wortmann

Changed assignment of k1,k2,k3 arrays as there was a bug

parent 62accbea
......@@ -251,9 +251,9 @@ CONTAINS
ENDIF
IF (PRESENT(k1)) THEN
IF (.NOT.PRESENT(k2).OR..NOT.PRESENT(k3)) CALL juDFT_error("BUG: always read k1,k2,k3")
k1=d%eig_int(3+1:3+d%nmat,nrec)
k2=d%eig_int(3+d%nmat+1:3+2*d%nmat,nrec)
k3=d%eig_int(3+2*d%nmat+1:3+3*d%nmat,nrec)
k1=d%eig_int(3+1:3+size(k1),nrec)
k2=d%eig_int(3+d%nmat+1:3+d%nmat+size(k1),nrec)
k3=d%eig_int(3+2*d%nmat+1:3+2*d%nmat*size(k1),nrec)
ENDIF
IF (PRESENT(kveclo)) kveclo=d%eig_int(4+3*d%nmat:3+3*d%nmat+SIZE(kveclo),nrec)
......@@ -318,9 +318,9 @@ CONTAINS
IF (PRESENT(k1)) THEN
IF (.NOT.PRESENT(k2).OR..NOT.PRESENT(k3)) CALL juDFT_error("BUG: always write k1,k2,k3")
d%eig_int(3+1:3+d%nmat,nrec)=k1
d%eig_int(3+d%nmat+1:3+2*d%nmat,nrec)=k2
d%eig_int(3+2*d%nmat+1:3+3*d%nmat,nrec)=k3
d%eig_int(3+1:3+size(k1),nrec)=k1
d%eig_int(3+d%nmat+1:3+d%nmat+size(k1),nrec)=k2
d%eig_int(3+2*d%nmat+1:3+2*d%nmat+size(k1),nrec)=k3
ENDIF
IF (PRESENT(kveclo)) d%eig_int(4+3*d%nmat:3+3*d%nmat+SIZE(kveclo),nrec)=kveclo
......
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