Commit f53fc1b0 authored by Henning Janssen's avatar Henning Janssen

made test slightly less sensitive

parent 00cb19af
......@@ -18,9 +18,9 @@ $result+=jt::test_fileexists("$workdir/Hubbard1/2_f/hloc.cfg");
#test for one eigval file
$result+=jt::test_fileexists("$workdir/Hubbard1/2_f/eigval7part.dat");
#test density matrix
$result+=jt::test_grepnumber("$workdir/out","nmmp occupation distance:",": *([^ ]*)",7.00000000256298,0.0005);
$result+=jt::test_grepnumber("$workdir/out","nmmp element distance:",": *([^ ]*)",0.999909071638895,0.00001);
$result+=jt::test_grep_lastnumber("$workdir/out","total energy=","= *([^ ]*)",-22560.5556136242,0.00001);
$result+=jt::test_grepnumber("$workdir/out","distance of spin densities for it=","3: *([^ ]*)",25.053445,0.0001);
$result+=jt::test_grepnumber("$workdir/out","nmmp occupation distance:",": *([^ ]*)",7.0000,0.0005);
$result+=jt::test_grepnumber("$workdir/out","nmmp element distance:",": *([^ ]*)",0.9999,0.0005);
$result+=jt::test_grep_lastnumber("$workdir/out","total energy=","= *([^ ]*)",-22560.55561,0.00005);
$result+=jt::test_grepnumber("$workdir/out","distance of spin densities for it=","3: *([^ ]*)",25.0533,0.0002);
jt::stageresult($workdir,$result,"1");
......@@ -68,7 +68,7 @@ MODULE m_tetrahedronInit
DO icorn = 1, ncorn
ikpt = kpts%ntetra(icorn,itet)
IF(ikpt.GT.kpts%nkpt) CYCLE
fac = REAL(MERGE(1,COUNT(kpts%bkp(:).EQ.ikpt),kpts%nkptf.EQ.0))
fac = REAL(COUNT(kpts%bkp(:).EQ.ikpt))
vol = kpts%voltet(itet)/kpts%ntet*fac
!$OMP parallel do default(none) &
!$OMP shared(itet,neig,ikpt,icorn,film,ncorn,k,vol,l_weights_pres,l_weightsum_pres) &
......@@ -81,7 +81,7 @@ MODULE m_tetrahedronInit
IF( ALL(etetra>efermi) ) CYCLE
IF( ALL(etetra<efermi) ) THEN
w = vol/ncorn
w = vol/REAL(ncorn)
ELSE
w = getWeightSingleBand([efermi],etetra,icorn,vol,film,.FALSE.)
ENDIF
......@@ -120,7 +120,7 @@ MODULE m_tetrahedronInit
REAL, ALLOCATABLE :: calc_weights(:,:)
REAL, ALLOCATABLE :: calc_eMesh(:)
fac = REAL(MERGE(1,COUNT(kpts%bkp(:).EQ.ikpt),kpts%nkptf.EQ.0))
fac = REAL(COUNT(kpts%bkp(:).EQ.ikpt))
l_dos = PRESENT(dos)
IF(PRESENT(dos))THEN
......
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