Commit 65a10945 authored by Matthias Redies's avatar Matthias Redies

fix a bunch of warnings

parent cce6330b
......@@ -29,7 +29,7 @@ CONTAINS
DO j = 1, i
dg = gpt(:, j) - gpt(:, i)
g = gptnorm(dg, cell%bmat)
IF (g == 0) THEN
IF (abs(g) < 1e-10) THEN
DO itype = 1, atoms%ntype
r = atoms%rmt(itype)
if (olap%l_real) THEN
......@@ -46,7 +46,7 @@ CONTAINS
DO ineq = 1, atoms%neq(itype)
icent = icent + 1
if (olap%l_real) THEN
olap%data_r(i, j) = olap%data_r(i, j) - fgr*exp(img*tpi_const*dot_product(dg, atoms%taual(:, icent)))
olap%data_r(i, j) = real(olap%data_r(i, j) - fgr*exp(img*tpi_const*dot_product(dg, atoms%taual(:, icent))))
else
olap%data_c(i, j) = olap%data_c(i, j) - fgr*exp(img*tpi_const*dot_product(dg, atoms%taual(:, icent)))
endif
......@@ -99,7 +99,7 @@ CONTAINS
dg = gpt(:, i) - gpt(:, j)
g = gptnorm(dg, cell%bmat)
olap_r(k) = 0
IF (g == 0) THEN
IF (abs(g) < 1e-10) THEN
DO itype = 1, atoms%ntype
r = atoms%rmt(itype)
olap_r(k) = olap_r(k) - atoms%neq(itype)*fpi_const*r**3/3/cell%omtil
......
......@@ -1473,7 +1473,7 @@ CONTAINS
SUBROUTINE ket_trafo1(vecout, vecin, ikpt0, isym, lreal, lsymmetrize, &
atoms, kpts, sym, mpbasis, hybrid, &
cell, maxindxm, nindxm, nbasm, ngptmall, nbasp, lcutm, maxlcutm)
maxindxm, nindxm, nbasm, ngptmall, nbasp, lcutm, maxlcutm)
USE m_constants
USE m_util, ONLY: modulo1
......@@ -1484,7 +1484,6 @@ CONTAINS
TYPE(t_mpbasis), intent(in) :: mpbasis
TYPE(t_hybrid), INTENT(IN) :: hybrid
TYPE(t_sym), INTENT(IN) :: sym
TYPE(t_cell), INTENT(IN) :: cell
TYPE(t_kpts), INTENT(IN) :: kpts
TYPE(t_atoms), INTENT(IN) :: atoms
......@@ -1702,9 +1701,9 @@ CONTAINS
COMPLEX, INTENT(OUT) :: vecout(nbasm(ikpt0))
! - private scalars -
INTEGER :: itype, ieq, ic, l, n, i, nn, i1, i2, j1, j2, m1
INTEGER :: m2, igptm, igptm2, igptp, iiop, isym
INTEGER :: ikpt1, isymi, rcent
INTEGER :: itype, ieq, ic, l, n, i, nn, i1, i2, j1, j2
INTEGER :: igptm, igptm2, igptp, iiop, isym
INTEGER :: ikpt1, rcent
LOGICAL :: trs
COMPLEX, PARAMETER :: img = (0.0, 1.0)
COMPLEX :: cexp, cdum
......
......@@ -30,7 +30,7 @@ contains
implicit none
integer, intent(in) :: n
integer :: identity(n, n)
integer :: i, j
integer :: i
identity = 0
do i = 1, n
identity(i, i) = 1
......@@ -344,7 +344,7 @@ contains
real, intent(in) :: mat1(:), mat2(:, :)
real :: matmat_dpdm(size(mat2, 1), size(mat2, 1))
real, allocatable :: vec(:), vec2(:)
integer :: nn, n, k1, i
integer :: nn, n, i
n = size(mat2, 1); nn = n*(n + 1)/2; allocate(vec(n), vec2(n))
if (size(mat2, 2) /= n) &
call juDFT_error( 'matmat_dpdm: dimensions of second input array differ.')
......@@ -364,7 +364,7 @@ contains
real, intent(in) :: mat1(:, :), mat2(:)
real :: matmat_dmdp(size(mat1, 1), size(mat1, 1))
real, allocatable :: vec(:), vec2(:)
integer :: nn, n, k1, i
integer :: nn, n, i
n = size(mat1, 1); nn = n*(n + 1)/2; allocate(vec(n), vec2(n))
if (size(mat1, 2) /= n) &
call juDFT_error( 'matmat_dmdp: dimensions of first input array differ.')
......@@ -402,7 +402,7 @@ contains
complex :: matmat_dpzm(size(mat2, 1), size(mat2, 1))
real, allocatable :: vecr(:), veci(:)
complex, allocatable :: vec2(:)
integer :: nn, n, k1, i
integer :: nn, n, i
n = size(mat2, 1)
nn = n*(n + 1)/2; allocate(vecr(n), veci(n), vec2(n))
if (size(mat2, 2) /= n) &
......@@ -425,7 +425,7 @@ contains
complex, intent(in) :: mat2(:)
complex :: matmat_dmzp(size(mat1, 1), size(mat1, 1))
complex, allocatable :: vec1(:), vec(:)
integer :: nn, n, k1, i
integer :: nn, n, i
n = size(mat1, 1); nn = n*(n + 1)/2; allocate(vec(n), vec1(n))
if (size(mat1, 2) /= n) &
call juDFT_error( 'matmat_dmzp: dimensions of first input array differ.')
......@@ -465,7 +465,7 @@ contains
real, intent(in) :: mat2(:, :)
complex :: matmat_zpdm(size(mat2, 1), size(mat2, 1))
complex, allocatable :: vec(:), vec2(:)
integer :: nn, n, k1, i
integer :: nn, n, i
n = size(mat2, 1); nn = n*(n + 1)/2; allocate(vec(n), vec2(n))
if (size(mat2, 2) /= n) &
call juDFT_error( 'matmat_zpdm: dimensions of second input array differ.')
......@@ -487,7 +487,7 @@ contains
complex :: matmat_zmdp(size(mat1, 1), size(mat1, 1))
complex, allocatable :: vec1(:)
real, allocatable :: vecr(:), veci(:)
integer :: nn, n, k1, i
integer :: nn, n, i
n = size(mat1, 1); nn = n*(n + 1)/2
allocate(vecr(n), veci(n), vec1(n))
if (size(mat1, 2) /= n) &
......@@ -528,7 +528,7 @@ contains
complex, intent(in) :: mat1(:), mat2(:, :)
complex :: matmat_zpzm(size(mat2, 1), size(mat2, 2))
complex, allocatable :: vec(:), vec2(:)
integer :: nn, n, k1, i, n2
integer :: nn, n, i, n2
n = size(mat2, 1); nn = n*(n + 1)/2; allocate(vec(n), vec2(n))
n2 = size(mat2, 2)
if (size(mat1) /= nn) &
......@@ -547,7 +547,7 @@ contains
complex, intent(in) :: mat1(:, :), mat2(:)
complex :: matmat_zmzp(size(mat1, 1), size(mat1, 1))
complex, allocatable :: vec(:), vec2(:)
integer :: nn, n, k1, i
integer :: nn, n, i
n = size(mat1, 1); nn = n*(n + 1)/2; allocate(vec(n), vec2(n))
if (size(mat1, 2) /= n) &
call juDFT_error( 'matmat_zmzp: dimensions of first input array differ.')
......
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