1313
1414#include " Particle/LongRange/KContainer.h"
1515#include " Particle/Lattice/CrystalLattice.h"
16- #include " Particle/tests/ MinimalParticlePool.h"
16+ #include < MinimalParticlePool.h>
1717#include " Utilities/for_testing/checkVector.hpp"
1818#include " Utilities/for_testing/NativeInitializerPrint.hpp"
1919
@@ -39,8 +39,7 @@ TEST_CASE("kcontainer at gamma in 3D", "[longrange]")
3939 for (int ik = 0 ; ik < kcs.size (); ik++)
4040 {
4141 const double kc = kcs[ik] + 1e-6 ;
42- klists.updateKLists (lattice, kc,
43- ndim);
42+ klists.updateKLists (lattice, kc, ndim);
4443 CHECK (klists.getKpts ().size () == nks[ik]);
4544 }
4645 const int mxk = klists.getKpts ().size ();
@@ -76,14 +75,12 @@ TEST_CASE("kcontainer at twist in 3D", "[longrange]")
7675
7776 PosType twist;
7877 twist[0 ] = 0.1 ;
79- klists.updateKLists (lattice, kc,
80- ndim, twist);
78+ klists.updateKLists (lattice, kc, ndim, twist);
8179 CHECK (klists.getKpts ().size () == 1 );
8280 CHECK (klists.getKptsCartWorking ()[0 ][0 ] == Approx (blat * (twist[0 ] - 1 )));
8381
8482 twist = {-0.5 , 0 , 0.5 };
85- klists.updateKLists (lattice, kc,
86- ndim, twist);
83+ klists.updateKLists (lattice, kc, ndim, twist);
8784 int gvecs[3 ][3 ] = {{0 , 0 , -1 }, {1 , 0 , -1 }, {1 , 0 , 0 }};
8885 CHECK (klists.getKpts ().size () == 3 );
8986 for (int ik = 0 ; ik < klists.getKpts ().size (); ik++)
@@ -109,8 +106,7 @@ TEST_CASE("kcontainer at gamma in 2D", "[longrange]")
109106 for (int ik = 0 ; ik < kcs.size (); ik++)
110107 {
111108 const double kc = kcs[ik] + 1e-6 ;
112- klists.updateKLists (lattice, kc,
113- ndim);
109+ klists.updateKLists (lattice, kc, ndim);
114110 CHECK (klists.getKpts ().size () == nks[ik]);
115111 }
116112 const int mxk = klists.getKpts ().size ();
@@ -146,23 +142,20 @@ TEST_CASE("kcontainer at twist in 2D", "[longrange]")
146142
147143 PosType twist;
148144 twist[0 ] = 0.1 ;
149- klists.updateKLists (lattice, kc,
150- ndim, twist);
145+ klists.updateKLists (lattice, kc, ndim, twist);
151146 CHECK (klists.getKpts ().size () == 1 );
152147 CHECK (klists.getKptsCartWorking ()[0 ][0 ] == Approx (blat * (twist[0 ] - 1 )));
153148
154149 twist[1 ] = 0.1 ;
155- klists.updateKLists (lattice, kc,
156- ndim, twist);
150+ klists.updateKLists (lattice, kc, ndim, twist);
157151 CHECK (klists.getKpts ().size () == 2 );
158152 CHECK (klists.getKptsCartWorking ()[0 ][0 ] == Approx (blat * (twist[0 ] - 1 )));
159153 CHECK (klists.getKptsCartWorking ()[0 ][1 ] == Approx (blat * twist[1 ]));
160154 CHECK (klists.getKptsCartWorking ()[1 ][0 ] == Approx (blat * (twist[0 ])));
161155 CHECK (klists.getKptsCartWorking ()[1 ][1 ] == Approx (blat * (twist[1 ] - 1 )));
162156
163157 twist = {-0.5 , 0.5 , 0 };
164- klists.updateKLists (lattice, kc,
165- ndim, twist);
158+ klists.updateKLists (lattice, kc, ndim, twist);
166159 CHECK (klists.getKpts ().size () == 3 );
167160 // for (int ik=0;ik<3;ik++)
168161 // app_log() << klists.getKptsCartWorking()[ik] << std::endl;
@@ -189,7 +182,7 @@ TEST_CASE("kcontainer for diamond", "[longrange]")
189182 lattice.LR_dim_cutoff = 15 ;
190183 SimulationCell cell (lattice);
191184
192- const KContainer& klists = cell.getKLists ();
185+ const KContainer& klists = cell.getKLists ();
193186 std::remove_cv_t <std::remove_reference_t <decltype (KContainer ().getKpts ())>> kpoint_lists = {
194187 {-1 , -1 , -1 }, {-1 , 0 , 0 }, {0 , -1 , 0 }, {0 , 0 , -1 }, {0 , 0 , 1 }, {0 , 1 , 0 }, {1 , 0 , 0 }, {1 , 1 , 1 },
195188 {-1 , -1 , 0 }, {-1 , 0 , -1 }, {0 , -1 , -1 }, {0 , 1 , 1 }, {1 , 0 , 1 }, {1 , 1 , 0 }, {-2 , -1 , -1 }, {-1 , -2 , -1 },
0 commit comments