fix after merge

Author: maki49

include/RI/physics/Exx.h

@@ -29,7 +29,7 @@ class Exx
 	constexpr static std::size_t Npos = Ndim;		// tmp
 	using Tatom_pos = std::array<Tpos,Npos>;		// tmp
 
-	Exx() = default;
+	Exx(const std::string method = "loop3") { this->method = method; }
 	Exx(Exx&&) noexcept = default;
 	Exx& operator=(Exx&&) noexcept = default;
 

include/RI/physics/Exx.hpp

@@ -295,7 +295,6 @@ template<typename TA, typename Tcell, std::size_t Ndim, typename Tdata>
 void Exx<TA,Tcell,Ndim,Tdata>::cal_force(
 	const std::array<std::string,5> &save_names_suffix)						// "Cs","Vs","Ds","dCs","dVs"
 {
-	std::cout << "post_2D adress: " << &(this->post_2D) << std::endl;
 	assert(this->flag_finish.stru);
 	assert(this->flag_finish.Cs);
 	assert(this->flag_finish.Vs);

include/RI/physics/LR.h

@@ -9,18 +9,16 @@ template<typename TA, typename Tcell, std::size_t Ndim, typename Tdata>
 class  LR : public Exx<TA,Tcell,Ndim,Tdata>
 {
 public:
-	LR(Exx<TA, Tcell, Ndim, Tdata>&& exx) : Exx<TA, Tcell, Ndim, Tdata>::Exx(std::move(exx)) {}
 	LR(const std::string method = "cvc") { this->method = method; }	// Exx default mehtod: "loop3"
+	LR(Exx<TA, Tcell, Ndim, Tdata>&& exx) : Exx<TA, Tcell, Ndim, Tdata>::Exx(std::move(exx)) {}
 	using Exx<TA,Tcell,Ndim,Tdata>::cal_force;
 	// New function: cal_force with two different density matrices
 	void cal_force(const std::map<TA, std::map<std::pair<TA, std::array<Tcell, Ndim>>, Tensor<Tdata>>>& Ds_left,
 		const std::array<std::string, 5>& save_names_suffix = { "","","","","" })	// "Cs","Vs","Ds","dCs","dVs"
 	{
 		// The only difference from Exx::cal_force: save Ds_left if not empty
-		std::cout << "post_2D adress: " << &(this->post_2D) << std::endl;
 		if (!Ds_left.empty())
 			this->post_2D.saves["Ds_" + save_names_suffix[2]] = this->post_2D.set_tensors_map2(Ds_left);
-		std::cout << "set OK" << std::endl;
 		this->cal_force(save_names_suffix);
 	}
 };