Commit bd71079a authored by Matthias Redies's avatar Matthias Redies

ensure matmul has same type

parent 8f635e0f
...@@ -47,7 +47,7 @@ CONTAINS ...@@ -47,7 +47,7 @@ CONTAINS
COMPLEX, PARAMETER :: img = (0.0, 1.0) COMPLEX, PARAMETER :: img = (0.0, 1.0)
! - arrays - ! - arrays -
INTEGER :: rrot(3, 3), invrrot(3, 3) REAL :: rrot(3, 3), invrrot(3, 3)
INTEGER :: g(3), g1(3) INTEGER :: g(3), g1(3)
REAL :: tau1(3), rtaual(3), rkpt(3), rkpthlp(3), trans(3) REAL :: tau1(3), rtaual(3), rkpt(3), rkpthlp(3), trans(3)
COMPLEX :: cmthlp(2*atoms%lmaxd + 1) COMPLEX :: cmthlp(2*atoms%lmaxd + 1)
...@@ -56,21 +56,21 @@ CONTAINS ...@@ -56,21 +56,21 @@ CONTAINS
tpiimg = -tpi_const*img tpiimg = -tpi_const*img
if (l_real) THEN if (l_real) THEN
rrot = transpose(sym%mrot(:, :, sym%invtab(iop))) rrot = transpose(1.0 * sym%mrot(:, :, sym%invtab(iop)))
invrrot = transpose(sym%mrot(:, :, iop)) invrrot = transpose(1.0 * sym%mrot(:, :, iop))
trans = sym%tau(:, iop) trans = sym%tau(:, iop)
else else
IF (iop <= sym%nop) THEN IF (iop <= sym%nop) THEN
trs = .false. trs = .false.
rrot = transpose(sym%mrot(:, :, sym%invtab(iop))) rrot = transpose(1.0 * sym%mrot(:, :, sym%invtab(iop)))
invrrot = transpose(sym%mrot(:, :, iop)) invrrot = transpose(1.0 * sym%mrot(:, :, iop))
trans = sym%tau(:, iop) trans = sym%tau(:, iop)
ELSE ELSE
trs = .true. trs = .true.
iiop = iop - sym%nop iiop = iop - sym%nop
rrot = -transpose(sym%mrot(:, :, sym%invtab(iiop))) rrot = -transpose(1.0 * sym%mrot(:, :, sym%invtab(iiop)))
invrrot = -transpose(sym%mrot(:, :, iiop)) invrrot = -transpose(1.0 * sym%mrot(:, :, iiop))
trans = sym%tau(:, iiop) trans = sym%tau(:, iiop)
END IF END IF
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