Commit c5a94e1d authored by Stefan Rost's avatar Stefan Rost

added gw=1 and gw=2 interface test

parent 4f8b3f45
Start testing: Dec 05 14:59 CET
----------------------------------------------------------
End testing: Dec 05 14:59 CET
<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<fleurInput fleurInputVersion="0.31">
<comment>
Si bulk
</comment>
<calculationSetup>
<cutoffs Kmax="3.40000000" Gmax="11.10000000" GmaxXC="9.20000000" numbands="0"/>
<scfLoop itmax="1" minDistance=".00001000" maxIterBroyd="99" imix="Anderson" alpha=".05000000" precondParam="0.0" spinf="2.00000000"/>
<coreElectrons ctail="T" frcor="F" kcrel="0" coretail_lmax="0"/>
<magnetism jspins="1" l_noco="F" swsp="F" lflip="F"/>
<soc theta=".00000000" phi=".00000000" l_soc="F" spav="F"/>
<expertModes gw="1" secvar="F"/>
<geometryOptimization l_f="F" forcealpha="1.00000000" forcemix="BFGS" epsdisp=".00001000" epsforce=".00001000"/>
<ldaU l_linMix="F" mixParam=".050000" spinf="1.000000"/>
<bzIntegration valenceElectrons="8.00000000" mode="hist" fermiSmearingEnergy=".00100000">
<kPointCount count="6" gamma="F"/>
<altKPointSet purpose="bands">
<kPointCount count=" 240" gamma="F"/>
</altKPointSet>
</bzIntegration>
<energyParameterLimits ellow="-.80000000" elup="1.00000000"/>
</calculationSetup>
<cell>
<symmetryFile filename="sym.out"/>
<bulkLattice scale="1.0000000000" latnam="any">
<bravaisMatrix>
<row-1>.000000000000000 5.130608533500000 5.130608533500000</row-1>
<row-2>5.130608533500000 .000000000000000 5.130608533500000</row-2>
<row-3>5.130608533500000 5.130608533500000 .000000000000000</row-3>
</bravaisMatrix>
</bulkLattice>
</cell>
<xcFunctional name="pbe" relativisticCorrections="F"/>
<atomSpecies>
<species name="Si-1" element="Si" atomicNumber="14" coreStates="4" magMom=".00000000" flipSpin="T">
<mtSphere radius="2.17000000" gridPoints="717" logIncrement=".01600000"/>
<atomicCutoffs lmax="8" lnonsphr="6"/>
<energyParameters s="3" p="3" d="3" f="4"/>
</species>
</atomSpecies>
<atomGroups>
<atomGroup species="Si-1">
<relPos label=" 1">1.000/8.000 1.000/8.000 1.000/8.000</relPos>
<relPos label=" 2">-1.000/8.000 -1.000/8.000 -1.000/8.000</relPos>
<force calculate="T" relaxXYZ="TTT"/>
</atomGroup>
</atomGroups>
<output dos="F" band="F" vacdos="F" slice="F" mcd="F">
<checks vchk="F" cdinf="F"/>
<densityOfStates ndir="0" minEnergy="-.50000000" maxEnergy=".50000000" sigma=".01500000"/>
<vacuumDOS layers="0" integ="F" star="F" nstars="0" locx1=".00000" locy1=".00000" locx2=".00000" locy2=".00000" nstm="0" tworkf=".00000"/>
<unfoldingBand unfoldBand="F" supercellX="1" supercellY="1" supercellZ="1"/>
<plotting iplot="0" score="F" plplot="F"/>
<chargeDensitySlicing numkpt="0" minEigenval=".00000000" maxEigenval=".00000000" nnne="0" pallst="F"/>
<specialOutput eonly="F" bmt="F"/>
<magneticCircularDichroism energyLo="-10.00000000" energyUp=".00000000"/>
</output>
<!-- We include the file relax.inp here to enable relaxations (see documentation) -->
<xi:include xmlns:xi="http://www.w3.org/2001/XInclude" href="relax.xml"> <xi:fallback/> </xi:include>
</fleurInput>
48 48 F ! nop,nop2,symor
! 1
1 0 0 0.00000
0 1 0 0.00000
0 0 1 0.00000
! 2
0 1 0 0.00000
1 0 0 0.00000
0 0 1 0.00000
! 3
-1 0 0 0.00000
0 -1 0 0.00000
1 1 1 0.50000
! 4
0 -1 0 0.00000
-1 0 0 0.00000
1 1 1 0.50000
! 5
0 1 0 0.00000
1 0 0 0.00000
-1 -1 -1 0.50000
! 6
1 0 0 0.00000
0 1 0 0.00000
-1 -1 -1 0.50000
! 7
-1 0 0 0.00000
0 -1 0 0.00000
0 0 -1 0.00000
! 8
0 -1 0 0.00000
-1 0 0 0.00000
0 0 -1 0.00000
! 9
0 0 1 0.00000
1 0 0 0.00000
-1 -1 -1 0.50000
! 10
0 0 1 0.00000
0 1 0 0.00000
-1 -1 -1 0.50000
! 11
1 1 1 0.50000
0 -1 0 0.00000
0 0 -1 0.00000
! 12
1 1 1 0.50000
-1 0 0 0.00000
0 0 -1 0.00000
! 13
-1 0 0 0.00000
1 1 1 0.50000
0 0 -1 0.00000
! 14
1 0 0 0.00000
0 0 1 0.00000
-1 -1 -1 0.50000
! 15
0 1 0 0.00000
0 0 1 0.00000
-1 -1 -1 0.50000
! 16
0 -1 0 0.00000
1 1 1 0.50000
0 0 -1 0.00000
! 17
-1 0 0 0.00000
0 0 -1 0.00000
0 -1 0 0.00000
! 18
0 1 0 0.00000
-1 -1 -1 0.50000
1 0 0 0.00000
! 19
0 -1 0 0.00000
0 0 -1 0.00000
-1 0 0 0.00000
! 20
1 0 0 0.00000
-1 -1 -1 0.50000
0 1 0 0.00000
! 21
1 1 1 0.50000
0 0 -1 0.00000
0 -1 0 0.00000
! 22
0 0 1 0.00000
-1 -1 -1 0.50000
1 0 0 0.00000
! 23
1 1 1 0.50000
0 0 -1 0.00000
-1 0 0 0.00000
! 24
0 0 1 0.00000
-1 -1 -1 0.50000
0 1 0 0.00000
! 25
0 0 -1 0.00000
-1 0 0 0.00000
0 -1 0 0.00000
! 26
0 0 -1 0.00000
0 -1 0 0.00000
-1 0 0 0.00000
! 27
-1 -1 -1 0.50000
0 1 0 0.00000
1 0 0 0.00000
! 28
-1 -1 -1 0.50000
1 0 0 0.00000
0 1 0 0.00000
! 29
0 0 1 0.00000
0 1 0 0.00000
1 0 0 0.00000
! 30
1 1 1 0.50000
-1 0 0 0.00000
0 -1 0 0.00000
! 31
1 1 1 0.50000
0 -1 0 0.00000
-1 0 0 0.00000
! 32
0 0 1 0.00000
1 0 0 0.00000
0 1 0 0.00000
! 33
0 0 -1 0.00000
1 1 1 0.50000
-1 0 0 0.00000
! 34
0 0 -1 0.00000
1 1 1 0.50000
0 -1 0 0.00000
! 35
-1 -1 -1 0.50000
0 0 1 0.00000
1 0 0 0.00000
! 36
-1 -1 -1 0.50000
0 0 1 0.00000
0 1 0 0.00000
! 37
-1 0 0 0.00000
1 1 1 0.50000
0 -1 0 0.00000
! 38
0 1 0 0.00000
0 0 1 0.00000
1 0 0 0.00000
! 39
0 -1 0 0.00000
1 1 1 0.50000
-1 0 0 0.00000
! 40
1 0 0 0.00000
0 0 1 0.00000
0 1 0 0.00000
! 41
1 0 0 0.00000
-1 -1 -1 0.50000
0 0 1 0.00000
! 42
0 1 0 0.00000
-1 -1 -1 0.50000
0 0 1 0.00000
! 43
0 -1 0 0.00000
0 0 -1 0.00000
1 1 1 0.50000
! 44
-1 0 0 0.00000
0 0 -1 0.00000
1 1 1 0.50000
! 45
-1 -1 -1 0.50000
1 0 0 0.00000
0 0 1 0.00000
! 46
-1 -1 -1 0.50000
0 1 0 0.00000
0 0 1 0.00000
! 47
0 0 -1 0.00000
0 -1 0 0.00000
1 1 1 0.50000
! 48
0 0 -1 0.00000
-1 0 0 0.00000
1 1 1 0.50000
$test_name="GW=1 switch for Si";
$test_code="Fleur";
$test_stages=1;
$test_desc=<<EOF
Simple test of Fleur with one step:
1.Generate a starting density and run 1 iteration with gw=1 and compare several quantities written out
EOF
;
#juDFT Testscript
jt::copyfile("files/inp.xml",$workdir);
jt::copyfile("files/sym.out",$workdir);
jt::testrun("$executable ",$workdir);
#now test output
$result=jt::test_grepexists("$workdir/out","it= 1 is completed");
$result=jt::test_fileexists("$workdir/basis.hdf");
$result=jt::test_fileexists("$workdir/pot.hdf");
$result=jt::test_fileexists("$workdir/ecore");
$result=jt::test_fileexists("$workdir/sym.out");
jt::stageresult($workdir,$result,"1");
Start testing: Dec 05 14:59 CET
----------------------------------------------------------
End testing: Dec 05 14:59 CET
<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<fleurInput fleurInputVersion="0.31">
<comment>
Si bulk
</comment>
<calculationSetup>
<cutoffs Kmax="3.70000000" Gmax="11.10000000" GmaxXC="9.20000000" numbands="0"/>
<scfLoop itmax="1" minDistance=".00001000" maxIterBroyd="99" imix="Anderson" alpha=".05000000" precondParam="0.0" spinf="2.00000000"/>
<coreElectrons ctail="T" frcor="F" kcrel="0" coretail_lmax="0"/>
<magnetism jspins="1" l_noco="F" swsp="F" lflip="F"/>
<soc theta=".00000000" phi=".00000000" l_soc="F" spav="F"/>
<expertModes gw="2" secvar="F"/>
<geometryOptimization l_f="F" forcealpha="1.00000000" forcemix="BFGS" epsdisp=".00001000" epsforce=".00001000"/>
<ldaU l_linMix="F" mixParam=".050000" spinf="1.000000"/>
<bzIntegration valenceElectrons="8.00000000" mode="hist" fermiSmearingEnergy=".00100000">
<kPointCount count="17" gamma="F"/>
<altKPointSet purpose="bands">
<kPointCount count=" 240" gamma="F"/>
</altKPointSet>
<altKPointSet purpose="GW">
<kPointListFile filename="kpts_gw"/>
</altKPointSet>
</bzIntegration>
<energyParameterLimits ellow="-.80000000" elup="1.00000000"/>
</calculationSetup>
<cell>
<symmetryFile filename="sym.out"/>
<bulkLattice scale="1.0000000000" latnam="any">
<bravaisMatrix>
<row-1>.000000000000000 5.130608533500000 5.130608533500000</row-1>
<row-2>5.130608533500000 .000000000000000 5.130608533500000</row-2>
<row-3>5.130608533500000 5.130608533500000 .000000000000000</row-3>
</bravaisMatrix>
</bulkLattice>
</cell>
<xcFunctional name="pbe" relativisticCorrections="F"/>
<atomSpecies>
<species name="Si-1" element="Si" atomicNumber="14" coreStates="4" magMom=".00000000" flipSpin="T">
<mtSphere radius="2.17000000" gridPoints="717" logIncrement=".01600000"/>
<atomicCutoffs lmax="8" lnonsphr="6"/>
<energyParameters s="3" p="3" d="3" f="4"/>
</species>
</atomSpecies>
<atomGroups>
<atomGroup species="Si-1">
<relPos label=" 1">1.000/8.000 1.000/8.000 1.000/8.000</relPos>
<relPos label=" 2">-1.000/8.000 -1.000/8.000 -1.000/8.000</relPos>
<force calculate="T" relaxXYZ="TTT"/>
</atomGroup>
</atomGroups>
<output dos="F" band="F" vacdos="F" slice="F" mcd="F">
<checks vchk="F" cdinf="F"/>
<densityOfStates ndir="0" minEnergy="-.50000000" maxEnergy=".50000000" sigma=".01500000"/>
<vacuumDOS layers="0" integ="F" star="F" nstars="0" locx1=".00000" locy1=".00000" locx2=".00000" locy2=".00000" nstm="0" tworkf=".00000"/>
<unfoldingBand unfoldBand="F" supercellX="1" supercellY="1" supercellZ="1"/>
<plotting iplot="0" score="F" plplot="F"/>
<chargeDensitySlicing numkpt="0" minEigenval=".00000000" maxEigenval=".00000000" nnne="0" pallst="F"/>
<specialOutput eonly="F" bmt="F"/>
<magneticCircularDichroism energyLo="-10.00000000" energyUp=".00000000"/>
</output>
<!-- We include the file relax.inp here to enable relaxations (see documentation) -->
<xi:include xmlns:xi="http://www.w3.org/2001/XInclude" href="relax.xml"> <xi:fallback/> </xi:include>
</fleurInput>
3 2.0000000000
0.00000 0.00000 0.00000 1.00000
0.00000 0.00000 1.00000 4.00000
0.00000 1.00000 1.00000 3.00000
48 48 F ! nop,nop2,symor
! 1
1 0 0 0.00000
0 1 0 0.00000
0 0 1 0.00000
! 2
0 1 0 0.00000
1 0 0 0.00000
0 0 1 0.00000
! 3
-1 0 0 0.00000
0 -1 0 0.00000
1 1 1 0.50000
! 4
0 -1 0 0.00000
-1 0 0 0.00000
1 1 1 0.50000
! 5
0 1 0 0.00000
1 0 0 0.00000
-1 -1 -1 0.50000
! 6
1 0 0 0.00000
0 1 0 0.00000
-1 -1 -1 0.50000
! 7
-1 0 0 0.00000
0 -1 0 0.00000
0 0 -1 0.00000
! 8
0 -1 0 0.00000
-1 0 0 0.00000
0 0 -1 0.00000
! 9
0 0 1 0.00000
1 0 0 0.00000
-1 -1 -1 0.50000
! 10
0 0 1 0.00000
0 1 0 0.00000
-1 -1 -1 0.50000
! 11
1 1 1 0.50000
0 -1 0 0.00000
0 0 -1 0.00000
! 12
1 1 1 0.50000
-1 0 0 0.00000
0 0 -1 0.00000
! 13
-1 0 0 0.00000
1 1 1 0.50000
0 0 -1 0.00000
! 14
1 0 0 0.00000
0 0 1 0.00000
-1 -1 -1 0.50000
! 15
0 1 0 0.00000
0 0 1 0.00000
-1 -1 -1 0.50000
! 16
0 -1 0 0.00000
1 1 1 0.50000
0 0 -1 0.00000
! 17
-1 0 0 0.00000
0 0 -1 0.00000
0 -1 0 0.00000
! 18
0 1 0 0.00000
-1 -1 -1 0.50000
1 0 0 0.00000
! 19
0 -1 0 0.00000
0 0 -1 0.00000
-1 0 0 0.00000
! 20
1 0 0 0.00000
-1 -1 -1 0.50000
0 1 0 0.00000
! 21
1 1 1 0.50000
0 0 -1 0.00000
0 -1 0 0.00000
! 22
0 0 1 0.00000
-1 -1 -1 0.50000
1 0 0 0.00000
! 23
1 1 1 0.50000
0 0 -1 0.00000
-1 0 0 0.00000
! 24
0 0 1 0.00000
-1 -1 -1 0.50000
0 1 0 0.00000
! 25
0 0 -1 0.00000
-1 0 0 0.00000
0 -1 0 0.00000
! 26
0 0 -1 0.00000
0 -1 0 0.00000
-1 0 0 0.00000
! 27
-1 -1 -1 0.50000
0 1 0 0.00000
1 0 0 0.00000
! 28
-1 -1 -1 0.50000
1 0 0 0.00000
0 1 0 0.00000
! 29
0 0 1 0.00000
0 1 0 0.00000
1 0 0 0.00000
! 30
1 1 1 0.50000
-1 0 0 0.00000
0 -1 0 0.00000
! 31
1 1 1 0.50000
0 -1 0 0.00000
-1 0 0 0.00000
! 32
0 0 1 0.00000
1 0 0 0.00000
0 1 0 0.00000
! 33
0 0 -1 0.00000
1 1 1 0.50000
-1 0 0 0.00000
! 34
0 0 -1 0.00000
1 1 1 0.50000
0 -1 0 0.00000
! 35
-1 -1 -1 0.50000
0 0 1 0.00000
1 0 0 0.00000
! 36
-1 -1 -1 0.50000
0 0 1 0.00000
0 1 0 0.00000
! 37
-1 0 0 0.00000
1 1 1 0.50000
0 -1 0 0.00000
! 38
0 1 0 0.00000
0 0 1 0.00000
1 0 0 0.00000
! 39
0 -1 0 0.00000
1 1 1 0.50000
-1 0 0 0.00000
! 40
1 0 0 0.00000
0 0 1 0.00000
0 1 0 0.00000
! 41
1 0 0 0.00000
-1 -1 -1 0.50000
0 0 1 0.00000
! 42
0 1 0 0.00000
-1 -1 -1 0.50000
0 0 1 0.00000
! 43
0 -1 0 0.00000
0 0 -1 0.00000
1 1 1 0.50000
! 44
-1 0 0 0.00000
0 0 -1 0.00000
1 1 1 0.50000
! 45
-1 -1 -1 0.50000
1 0 0 0.00000
0 0 1 0.00000
! 46
-1 -1 -1 0.50000
0 1 0 0.00000
0 0 1 0.00000
! 47
0 0 -1 0.00000
0 -1 0 0.00000
1 1 1 0.50000
! 48
0 0 -1 0.00000
-1 0 0 0.00000
1 1 1 0.50000
$test_name="GW=2 switch for Si";
$test_code="Fleur";
$test_stages=1;
$test_desc=<<EOF
Simple test of Fleur with one step:
1.Generate a starting density and run 1 iteration with gw=1 and compare several quantities written out
EOF
;
#juDFT Testscript
jt::copyfile("files/inp.xml",$workdir);
jt::copyfile("files/sym.out",$workdir);
jt::testrun("$executable ",$workdir);
#now test output
$result+=jt::test_grepexists("$workdir/out","0.50000.*0.37500");
$result=jt::test_fileexists("$workdir/eig_gw.hdf");
jt::stageresult($workdir,$result,"1");
......@@ -4,7 +4,7 @@ set(SerialParallelTests CuBulk CuBulkXML SiLOXML Fe_1l Fe_1lXML Fe-Atom CuBand C
Fe_bct Fe_bctXML PTO PTOXML Fe_1l_SOCXML PTO-SOC PTO-SOCXML Fe_bct_SOC Fe_bct_SOCXML Fe_fccXML
GaAsMultiUForceXML SiFilmPlotXML SiFilmSlicePlotXML CoMCDXML Fe_Kerker Fe_bct_LOXML )#SiHybridGamma)
set(SerialOnlyTests KClHybridPBE0 GaAsHybridPBE0 FeHybridPBE0 Fe_bct_LO Fe_fcc CoUnfold)# TiO2eels TiO2eelsXML)
set(SerialOnlyTests KClHybridPBE0 GaAsHybridPBE0 FeHybridPBE0 Fe_bct_LO Fe_fcc CoUnfold gw1Interface gw2Iterface)# TiO2eels TiO2eelsXML)
set(InpgenTests Si_plain Si_plain_explicit Si_full_para)# Si_kpt Si_kden Si_round_trip)
set(Testdirs ${SerialParallelTests} ${SerialOnlyTests})
......
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