Commit 02300a33 authored by Henning Janssen's avatar Henning Janssen

Unify basic Hubbard 1 test

- made a little less sensitive
parent d23949bc
This diff is collapsed.
$test_name="Fleur Gd Hubbard 1";
$test_name="Fleur Gd Hubbard 1 SOC";
$test_code="Fleur";
%test_requirements=("SOC",0);
%test_requirements=("SOC",1);
$test_stages=1;
$test_desc=<<EOF
Simple testthe Hubbard 1 method without SOC:
Simple testthe Hubbard 1 method with SOC:
1. Generate starting density, run 2 Iterations with one Hubbard iteration in between
for f-orbitals. Ensure that the density matrix is reasonable
EOF
......
......@@ -10,10 +10,12 @@ $result+=jt::test_grepexists("$workdir/out","Hubbard 1 it= 1 is completed");
#Test for the input file of the solver
$result+=jt::test_fileexists("$workdir/Hubbard1/hubbard1.cfg");
$result+=jt::test_fileexists("$workdir/Hubbard1/hloc.cfg");
#Test if the SOC parameter is correct
$result+=jt::test_grepnumber("$workdir/Hubbard1/hloc.cfg","xiSOC","xiSOC *([^ ]*)",0.21978,0.00001);
#test for one eigval file
$result+=jt::test_fileexists("$workdir/Hubbard1/eigval7part.dat");
#test density matrix
$result+=jt::test_grepnumber("$workdir/out","nmmp occupation distance:",": *([^ ]*)",7.00015383049998,0.00001);
$result+=jt::test_grepnumber("$workdir/out","nmmp element distance:",": *([^ ]*)",1.00023177097908,0.00001);
$result+=jt::test_grepnumber("$workdir/out","nmmp occupation distance:",": *([^ ]*)",6.99968080424677,0.0001);
$result+=jt::test_grepnumber("$workdir/out","nmmp element distance:",": *([^ ]*)",0.997526549354331,0.0001);
jt::stageresult($workdir,$result,"1");
<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<fleurInput fleurInputVersion="0.32">
<comment>
hcp Gadolinium
</comment>
<calculationSetup>
<cutoffs Kmax="3.00000000" Gmax="10.00000000" GmaxXC="10.00000000" numbands="125"/>
<scfLoop itmax="2" 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="2" l_noco="F" swsp="F" lflip="F" l_onlyMtStDen="F"/>
<soc theta=".00000000" phi=".00000000" l_soc="T" spav="F"/>
<expertModes gw="0" secvar="F"/>
<geometryOptimization l_f="F" forcealpha="1.00000000" forcemix="BFGS" epsdisp=".00001000" epsforce=".00001000"/>
<ldaU l_linMix="F" mixParam=".050000" spinf="1.000000"/>
<greensFunction l_sphavg="T">
<realAxis ne="5400" ellow="-1.0" elup="1.0"/>
<contourSemicircle n="128" eb="-0.5" et="0.0" alpha="1.0"/>
</greensFunction>
<bzIntegration valenceElectrons="36.00000000" mode="tetra" fermiSmearingEnergy=".00100000">
<kPointList name="default" count="10">
<kPoint weight=" 12.0000000000000">-0.4000000000000 -0.4000000000000 0.3333333333333</kPoint>
<kPoint weight=" 12.0000000000000">-0.2000000000000 -0.4000000000000 0.3333333333333</kPoint>
<kPoint weight=" 12.0000000000000">-0.4000000000000 -0.6000000000000 0.3333333333333</kPoint>
<kPoint weight=" 12.0000000000000">-0.2000000000000 -0.2000000000000 0.3333333333333</kPoint>
<kPoint weight=" 2.0000000000000"> 0.0000000000000 0.0000000000000 0.3333333333333</kPoint>
<kPoint weight=" 6.0000000000000">-0.4000000000000 -0.4000000000000 0.0000000000000</kPoint>
<kPoint weight=" 6.0000000000000">-0.2000000000000 -0.4000000000000 0.0000000000000</kPoint>
<kPoint weight=" 6.0000000000000">-0.4000000000000 -0.6000000000000 0.0000000000000</kPoint>
<kPoint weight=" 6.0000000000000">-0.2000000000000 -0.2000000000000 0.0000000000000</kPoint>
<kPoint weight=" 1.0000000000000"> 0.0000000000000 0.0000000000000 0.0000000000000</kPoint>
<tetraeder ntet="174">
<tet vol=" 0.3866666666667">10 34 26 30</tet>
<tet vol=" 0.3866666666667">5 22 26 30</tet>
<tet vol=" 0.3866666666667">10 5 26 30</tet>
<tet vol=" 0.3866666666667">42 10 24 5</tet>
<tet vol=" 0.3866666666667">38 20 24 5</tet>
<tet vol=" 0.3866666666667">42 38 24 5</tet>
<tet vol=" 0.3866666666667">10 26 24 5</tet>
<tet vol=" 0.3866666666667">26 22 24 5</tet>
<tet vol=" 0.3866666666667">20 22 24 5</tet>
<tet vol=" 0.3866666666667">33 8 15 29</tet>
<tet vol=" 0.3866666666667">8 3 15 29</tet>
<tet vol=" 0.3866666666667">11 3 15 29</tet>
<tet vol=" 0.3866666666667">33 24 8 20</tet>
<tet vol=" 0.3866666666667">29 3 8 20</tet>
<tet vol=" 0.3866666666667">33 29 8 20</tet>
<tet vol=" 0.3866666666667">24 25 8 20</tet>
<tet vol=" 0.3866666666667">25 21 8 20</tet>
<tet vol=" 0.3866666666667">3 21 8 20</tet>
<tet vol=" 0.3866666666667">15 8 17 3</tet>
<tet vol=" 0.3866666666667">11 13 17 3</tet>
<tet vol=" 0.3866666666667">15 11 17 3</tet>
<tet vol=" 0.3866666666667">8 6 17 3</tet>
<tet vol=" 0.3866666666667">6 1 17 3</tet>
<tet vol=" 0.3866666666667">13 1 17 3</tet>
<tet vol=" 0.3866666666667">8 25 6 21</tet>
<tet vol=" 0.3866666666667">3 1 6 21</tet>
<tet vol=" 0.3866666666667">8 3 6 21</tet>
<tet vol=" 0.3866666666667">25 7 6 21</tet>
<tet vol=" 0.3866666666667">7 2 6 21</tet>
<tet vol=" 0.3866666666667">1 2 6 21</tet>
<tet vol=" 0.3866666666667">25 23 7 19</tet>
<tet vol=" 0.3866666666667">21 2 7 19</tet>
<tet vol=" 0.3866666666667">25 21 7 19</tet>
<tet vol=" 0.3866666666667">23 47 7 19</tet>
<tet vol=" 0.3866666666667">47 43 7 19</tet>
<tet vol=" 0.3866666666667">2 43 7 19</tet>
<tet vol=" 0.3866666666667">17 6 41 1</tet>
<tet vol=" 0.3866666666667">13 37 41 1</tet>
<tet vol=" 0.3866666666667">17 13 41 1</tet>
<tet vol=" 0.3866666666667">6 40 41 1</tet>
<tet vol=" 0.3866666666667">40 36 41 1</tet>
<tet vol=" 0.3866666666667">37 36 41 1</tet>
<tet vol=" 0.3866666666667">6 7 40 2</tet>
<tet vol=" 0.3866666666667">1 36 40 2</tet>
<tet vol=" 0.3866666666667">6 1 40 2</tet>
<tet vol=" 0.3866666666667">7 9 40 2</tet>
<tet vol=" 0.3866666666667">9 4 40 2</tet>
<tet vol=" 0.3866666666667">36 4 40 2</tet>
<tet vol=" 0.3866666666667">7 47 9 43</tet>
<tet vol=" 0.3866666666667">2 4 9 43</tet>
<tet vol=" 0.3866666666667">7 2 9 43</tet>
<tet vol=" 0.3866666666667">47 50 9 43</tet>
<tet vol=" 0.3866666666667">50 46 9 43</tet>
<tet vol=" 0.3866666666667">4 46 9 43</tet>
<tet vol=" 0.3866666666667">50 48 10 44</tet>
<tet vol=" 0.3866666666667">46 5 10 44</tet>
<tet vol=" 0.3866666666667">50 46 10 44</tet>
<tet vol=" 0.3866666666667">48 34 10 44</tet>
<tet vol=" 0.3866666666667">34 30 10 44</tet>
<tet vol=" 0.3866666666667">5 30 10 44</tet>
<tet vol=" 0.3866666666667">40 9 39 4</tet>
<tet vol=" 0.3866666666667">36 35 39 4</tet>
<tet vol=" 0.3866666666667">40 36 39 4</tet>
<tet vol=" 0.3866666666667">9 42 39 4</tet>
<tet vol=" 0.3866666666667">42 38 39 4</tet>
<tet vol=" 0.3866666666667">35 38 39 4</tet>
<tet vol=" 0.3866666666667">9 50 42 46</tet>
<tet vol=" 0.3866666666667">4 38 42 46</tet>
<tet vol=" 0.3866666666667">9 4 42 46</tet>
<tet vol=" 0.3866666666667">50 10 42 46</tet>
<tet vol=" 0.3866666666667">10 5 42 46</tet>
<tet vol=" 0.3866666666667">38 5 42 46</tet>
<tet vol=" 0.3866666666667">5 30 22 67</tet>
<tet vol=" 0.3866666666667">5 55 22 67</tet>
<tet vol=" 0.3866666666667">38 5 20 55</tet>
<tet vol=" 0.3866666666667">5 22 20 55</tet>
<tet vol=" 0.3866666666667">29 3 11 66</tet>
<tet vol=" 0.3866666666667">3 53 11 66</tet>
<tet vol=" 0.3866666666667">29 20 3 61</tet>
<tet vol=" 0.3866666666667">66 53 3 61</tet>
<tet vol=" 0.3866666666667">29 66 3 61</tet>
<tet vol=" 0.3866666666667">20 21 3 61</tet>
<tet vol=" 0.3866666666667">21 62 3 61</tet>
<tet vol=" 0.3866666666667">53 62 3 61</tet>
<tet vol=" 0.3866666666667">11 3 13 53</tet>
<tet vol=" 0.3866666666667">3 1 13 53</tet>
<tet vol=" 0.3866666666667">1 51 13 53</tet>
<tet vol=" 0.3866666666667">3 21 1 62</tet>
<tet vol=" 0.3866666666667">53 51 1 62</tet>
<tet vol=" 0.3866666666667">3 53 1 62</tet>
<tet vol=" 0.3866666666667">21 2 1 62</tet>
<tet vol=" 0.3866666666667">2 52 1 62</tet>
<tet vol=" 0.3866666666667">51 52 1 62</tet>
<tet vol=" 0.3866666666667">21 19 2 60</tet>
<tet vol=" 0.3866666666667">62 52 2 60</tet>
<tet vol=" 0.3866666666667">21 62 2 60</tet>
<tet vol=" 0.3866666666667">19 43 2 60</tet>
<tet vol=" 0.3866666666667">43 72 2 60</tet>
<tet vol=" 0.3866666666667">52 72 2 60</tet>
<tet vol=" 0.3866666666667">13 1 37 51</tet>
<tet vol=" 0.3866666666667">1 36 37 51</tet>
<tet vol=" 0.3866666666667">1 2 36 52</tet>
<tet vol=" 0.3866666666667">1 51 36 52</tet>
<tet vol=" 0.3866666666667">2 4 36 52</tet>
<tet vol=" 0.3866666666667">4 54 36 52</tet>
<tet vol=" 0.3866666666667">2 43 4 72</tet>
<tet vol=" 0.3866666666667">52 54 4 72</tet>
<tet vol=" 0.3866666666667">2 52 4 72</tet>
<tet vol=" 0.3866666666667">43 46 4 72</tet>
<tet vol=" 0.3866666666667">46 75 4 72</tet>
<tet vol=" 0.3866666666667">54 75 4 72</tet>
<tet vol=" 0.3866666666667">46 44 5 73</tet>
<tet vol=" 0.3866666666667">75 55 5 73</tet>
<tet vol=" 0.3866666666667">46 75 5 73</tet>
<tet vol=" 0.3866666666667">44 30 5 73</tet>
<tet vol=" 0.3866666666667">30 67 5 73</tet>
<tet vol=" 0.3866666666667">55 67 5 73</tet>
<tet vol=" 0.3866666666667">36 4 35 54</tet>
<tet vol=" 0.3866666666667">4 38 35 54</tet>
<tet vol=" 0.3866666666667">4 46 38 75</tet>
<tet vol=" 0.3866666666667">4 54 38 75</tet>
<tet vol=" 0.3866666666667">46 5 38 75</tet>
<tet vol=" 0.3866666666667">5 55 38 75</tet>
<tet vol=" 0.3866666666667">10 26 63 34</tet>
<tet vol=" 0.3866666666667">55 10 63 34</tet>
<tet vol=" 0.3866666666667">71 55 61 10</tet>
<tet vol=" 0.3866666666667">42 24 61 10</tet>
<tet vol=" 0.3866666666667">71 42 61 10</tet>
<tet vol=" 0.3866666666667">55 63 61 10</tet>
<tet vol=" 0.3866666666667">63 26 61 10</tet>
<tet vol=" 0.3866666666667">24 26 61 10</tet>
<tet vol=" 0.3866666666667">53 8 56 33</tet>
<tet vol=" 0.3866666666667">15 8 56 33</tet>
<tet vol=" 0.3866666666667">33 8 53 24</tet>
<tet vol=" 0.3866666666667">8 25 53 24</tet>
<tet vol=" 0.3866666666667">56 53 58 8</tet>
<tet vol=" 0.3866666666667">15 17 58 8</tet>
<tet vol=" 0.3866666666667">56 15 58 8</tet>
<tet vol=" 0.3866666666667">53 51 58 8</tet>
<tet vol=" 0.3866666666667">51 6 58 8</tet>
<tet vol=" 0.3866666666667">17 6 58 8</tet>
<tet vol=" 0.3866666666667">8 6 51 25</tet>
<tet vol=" 0.3866666666667">53 8 51 25</tet>
<tet vol=" 0.3866666666667">52 7 51 25</tet>
<tet vol=" 0.3866666666667">6 7 51 25</tet>
<tet vol=" 0.3866666666667">25 7 52 23</tet>
<tet vol=" 0.3866666666667">7 47 52 23</tet>
<tet vol=" 0.3866666666667">58 51 70 6</tet>
<tet vol=" 0.3866666666667">17 41 70 6</tet>
<tet vol=" 0.3866666666667">58 17 70 6</tet>
<tet vol=" 0.3866666666667">51 69 70 6</tet>
<tet vol=" 0.3866666666667">69 40 70 6</tet>
<tet vol=" 0.3866666666667">41 40 70 6</tet>
<tet vol=" 0.3866666666667">51 52 69 7</tet>
<tet vol=" 0.3866666666667">6 40 69 7</tet>
<tet vol=" 0.3866666666667">51 6 69 7</tet>
<tet vol=" 0.3866666666667">52 54 69 7</tet>
<tet vol=" 0.3866666666667">54 9 69 7</tet>
<tet vol=" 0.3866666666667">40 9 69 7</tet>
<tet vol=" 0.3866666666667">7 9 54 47</tet>
<tet vol=" 0.3866666666667">52 7 54 47</tet>
<tet vol=" 0.3866666666667">9 50 54 47</tet>
<tet vol=" 0.3866666666667">50 10 55 48</tet>
<tet vol=" 0.3866666666667">10 34 55 48</tet>
<tet vol=" 0.3866666666667">69 54 68 9</tet>
<tet vol=" 0.3866666666667">40 39 68 9</tet>
<tet vol=" 0.3866666666667">69 40 68 9</tet>
<tet vol=" 0.3866666666667">54 71 68 9</tet>
<tet vol=" 0.3866666666667">71 42 68 9</tet>
<tet vol=" 0.3866666666667">39 42 68 9</tet>
<tet vol=" 0.3866666666667">9 42 71 50</tet>
<tet vol=" 0.3866666666667">54 9 71 50</tet>
<tet vol=" 0.3866666666667">55 10 71 50</tet>
<tet vol=" 0.3866666666667">42 10 71 50</tet>
</tetraeder>
</kPointList>
</bzIntegration>
<energyParameterLimits ellow="-1.80000000" elup="1.00000000"/>
<symmetryOperations>
<symOp>
<row-1>1 0 0 .0000000000</row-1>
<row-2>0 1 0 .0000000000</row-2>
<row-3>0 0 1 .0000000000</row-3>
</symOp>
<symOp>
<row-1>-1 0 0 .0000000000</row-1>
<row-2>0 -1 0 .0000000000</row-2>
<row-3>0 0 1 .5000000000</row-3>
</symOp>
<symOp>
<row-1>0 1 0 .0000000000</row-1>
<row-2>-1 -1 0 .0000000000</row-2>
<row-3>0 0 1 .0000000000</row-3>
</symOp>
<symOp>
<row-1>-1 -1 0 .0000000000</row-1>
<row-2>1 0 0 .0000000000</row-2>
<row-3>0 0 1 .0000000000</row-3>
</symOp>
<symOp>
<row-1>1 1 0 .0000000000</row-1>
<row-2>-1 0 0 .0000000000</row-2>
<row-3>0 0 1 .5000000000</row-3>
</symOp>
<symOp>
<row-1>0 -1 0 .0000000000</row-1>
<row-2>1 1 0 .0000000000</row-2>
<row-3>0 0 1 .5000000000</row-3>
</symOp>
<symOp>
<row-1>1 0 0 .0000000000</row-1>
<row-2>0 1 0 .0000000000</row-2>
<row-3>0 0 -1 .5000000000</row-3>
</symOp>
<symOp>
<row-1>-1 0 0 .0000000000</row-1>
<row-2>0 -1 0 .0000000000</row-2>
<row-3>0 0 -1 .0000000000</row-3>
</symOp>
<symOp>
<row-1>0 1 0 .0000000000</row-1>
<row-2>-1 -1 0 .0000000000</row-2>
<row-3>0 0 -1 .5000000000</row-3>
</symOp>
<symOp>
<row-1>-1 -1 0 .0000000000</row-1>
<row-2>1 0 0 .0000000000</row-2>
<row-3>0 0 -1 .5000000000</row-3>
</symOp>
<symOp>
<row-1>1 1 0 .0000000000</row-1>
<row-2>-1 0 0 .0000000000</row-2>
<row-3>0 0 -1 .0000000000</row-3>
</symOp>
<symOp>
<row-1>0 -1 0 .0000000000</row-1>
<row-2>1 1 0 .0000000000</row-2>
<row-3>0 0 -1 .0000000000</row-3>
</symOp>
</symmetryOperations>
</calculationSetup>
<cell>
<bulkLattice scale="1.0000000000">
<bravaisMatrix>
<row-1> -5.9406819509 3.4298543235 0.0000000000</row-1>
<row-2> -5.9406819509 -3.4298543235 0.0000000000</row-2>
<row-3> 0.0000000000 0.0000000000 10.9549358120</row-3>
</bravaisMatrix>
</bulkLattice>
</cell>
<xcFunctional name="pbe" relativisticCorrections="F"/>
<atomSpecies>
<species name="Gd-1" element="Gd" atomicNumber="64" flipSpinPhi=".00000000" flipSpinTheta=".00000000" flipSpinScale="F">
<mtSphere radius="2.80000000" gridPoints="925" logIncrement=".01400000"/>
<atomicCutoffs lmax="8" lnonsphr="6"/>
<electronConfig>
<coreConfig>(1s1/2) (2s1/2) (2p1/2) (2p3/2) (3s1/2) (3p1/2) (3p3/2) (4s1/2) (3d3/2) (3d5/2) (4p1/2) (4p3/2) (4d3/2) (4d5/2)</coreConfig>
<valenceConfig>(5s1/2) (5p1/2) (5p3/2) (6s1/2) (4f5/2) (4f7/2) (5d3/2) (5d5/2)</valenceConfig>
<stateOccupation state="(4f5/2)" spinUp="2.55000000" spinDown=".45000000"/>
<stateOccupation state="(4f7/2)" spinUp="4.00000000" spinDown=".00000000"/>
<stateOccupation state="(5d3/2)" spinUp=".40000000" spinDown=".00000000"/>
<stateOccupation state="(5d5/2)" spinUp=".60000000" spinDown=".00000000"/>
</electronConfig>
<energyParameters s="6" p="6" d="5" f="4"/>
<ldaHIA l="3" U="6.7" J="0.7" l_amf="F" init_occ="7.0">
<exc l="2" J="0.0952" init_mom="0.5"/>
</ldaHIA>
<lo type="SCLO" l="0" n="5" eDeriv="0"/>
<lo type="SCLO" l="1" n="5" eDeriv="0"/>
</species>
</atomSpecies>
<atomGroups>
<atomGroup species="Gd-1">
<relPos label=" 1">1.000/3.000 1.000/3.000 1.000/4.000</relPos>
<relPos label=" 2">-1.000/3.000 -1.000/3.000 -1.000/4.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"/>
<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>
$test_name="Fleur Gd Hubbard 1 SOC";
$test_code="Fleur";
%test_requirements=("SOC",1);
$test_stages=1;
$test_desc=<<EOF
Simple testthe Hubbard 1 method with SOC:
1. Generate starting density, run 2 Iterations with one Hubbard iteration in between
for f-orbitals. Ensure that the density matrix is reasonable
EOF
;
#juDFT Testscript
jt::copyfile("files/inp.xml",$workdir);
jt::testrun("$executable ",$workdir);
#now test output
$result=jt::test_grepexists("$workdir/out","it= 1 is completed");
$result+=jt::test_grepexists("$workdir/out","Hubbard 1 it= 1 is completed");
#Test for the input file of the solver
$result+=jt::test_fileexists("$workdir/Hubbard1/hubbard1.cfg");
$result+=jt::test_fileexists("$workdir/Hubbard1/hloc.cfg");
#Test if the SOC parameter is correct
$result+=jt::test_grepnumber("$workdir/Hubbard1/hloc.cfg","xiSOC","xiSOC *([^ ]*)",0.21978,0.00001);
#test for one eigval file
$result+=jt::test_fileexists("$workdir/Hubbard1/eigval7part.dat");
#test density matrix
$result+=jt::test_grepnumber("$workdir/out","nmmp occupation distance:",": *([^ ]*)",6.99968080424677,0.00001);
$result+=jt::test_grepnumber("$workdir/out","nmmp element distance:",": *([^ ]*)",0.997526549354331,0.00001);
jt::stageresult($workdir,$result,"1");
......@@ -53,7 +53,7 @@ endif()
#Tests for EDsolver
if (FLEUR_USE_EDSOLVER)
set(SerialParallelTests ${SerialParallelTests} Gd_Hubbard1 Gd_Hubbard1_SOC Gd_Hubbard1_noSYM)
set(SerialParallelTests ${SerialParallelTests} Gd_Hubbard1 Gd_Hubbard1_noSYM)
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