Loading include/dca/phys/dca_analysis/bse_solver/bse_lattice_solver.hpp +6 −4 Original line number Diff line number Diff line Loading @@ -143,7 +143,8 @@ private: DcaDataType& MOMS; func::function<std::complex<ScalarType>, HOST_matrix_dmn_t> Gamma_lattice; func::function<std::complex<ScalarType>, func::dmn_variadic<b_b, b_b, k_HOST_VERTEX, WVertexDmn>> chi_0_lattice; func::function<std::complex<ScalarType>, func::dmn_variadic<b_b, b_b, k_HOST_VERTEX, WVertexDmn>> chi_0_lattice; // Matrix in \vec{k} and \omega_n with the diagonal = chi_0_lattice. func::function<std::complex<ScalarType>, HOST_matrix_dmn_t> chi_0_lattice_matrix; Loading Loading @@ -310,7 +311,7 @@ void BseLatticeSolver<ParametersType, DcaDataType, ScalarType>::computeChi0Latti clustermapping::coarsegraining_tp<ParametersType, k_HOST_VERTEX> coarsegraining_tp(parameters); // DCA+/DCA with post-interpolation: Compute \chi_0 from continuous lattice self-energy. if (parameters.do_dca_plus() || parameters.do_post_interpolation()) { if (parameters.do_dca_plus() || parameters.doPostInterpolation()) { latticemapping::lattice_mapping_sp<ParametersType, k_DCA, k_HOST> lattice_map_sp(parameters); MOMS.Sigma_lattice_interpolated = 0.; Loading Loading @@ -373,7 +374,7 @@ void BseLatticeSolver<ParametersType, DcaDataType, ScalarType>::computeGammaLatt std::cout << "\n" << __FUNCTION__ << std::endl; // DCA+/DCA with post-interpolation: Compute Gamma_lattice with continuous momentum dependence. if (parameters.do_dca_plus() || parameters.do_post_interpolation()) { if (parameters.do_dca_plus() || parameters.doPostInterpolation()) { latticemapping::lattice_mapping_tp<ParametersType, k_DCA, k_HOST_VERTEX> lattice_map_tp( parameters); lattice_map_tp.execute(Gamma_cluster, Gamma_lattice); Loading Loading @@ -814,7 +815,8 @@ void BseLatticeSolver<ParametersType, DcaDataType, ScalarType>::printOnShell() { std::complex<ScalarType> norm = 0; for (int j = 0; j < k_HOST_VERTEX::dmn_size(); j++) { scal_prod += conj(psi_k(j, l)) * leading_eigenvectors(i, 0, 0, j, WVertexDmn::dmn_size() / 2); scal_prod += conj(psi_k(j, l)) * leading_eigenvectors(i, 0, 0, j, WVertexDmn::dmn_size() / 2); norm += conj(leading_eigenvectors(i, 0, 0, j, WVertexDmn::dmn_size() / 2)) * leading_eigenvectors(i, 0, 0, j, WVertexDmn::dmn_size() / 2); } Loading include/dca/phys/dca_data/dca_data.hpp +1 −1 Original line number Diff line number Diff line Loading @@ -452,7 +452,7 @@ void DcaData<Parameters>::write(Writer& writer) { if (parameters_.dump_lattice_self_energy()) { if (parameters_.do_dca_plus()) writer.execute(Sigma_lattice); else if (parameters_.do_post_interpolation()) else if (parameters_.doPostInterpolation()) writer.execute(Sigma_lattice_interpolated); } Loading include/dca/phys/dca_loop/dca_loop.hpp +1 −1 Original line number Diff line number Diff line Loading @@ -314,7 +314,7 @@ void DcaLoop<ParametersType, DcaDataType, MCIntegratorType>::perform_lattice_map if (concurrency.id() == concurrency.first()) std::cout << "\n\t\t lattice-mapping " << dca::util::print_time(); if (parameters.do_dca_plus() || parameters.do_post_interpolation()) { if (parameters.do_dca_plus() || parameters.doPostInterpolation()) { if (parameters.hts_approximation()) { DcaDataType MOMS_HTS(parameters); Loading include/dca/phys/parameters/dca_parameters.hpp +1 −1 Original line number Diff line number Diff line Loading @@ -84,7 +84,7 @@ public: return interacting_orbitals_; } bool do_post_interpolation() const { bool doPostInterpolation() const { return do_post_interpolation_; } Loading include/dca/phys/parameters/parameters.hpp +1 −1 Original line number Diff line number Diff line Loading @@ -292,7 +292,7 @@ void Parameters<Concurrency, Threading, Profiler, Model, RandomNumberGenerator, KSpHostDmn::parameter_type::print(std::cout); // Host grid for two-particle functions (tp-lattice) if (do_dca_plus() || do_post_interpolation()) { if (do_dca_plus() || doPostInterpolation()) { domains::cluster_domain_initializer<RTpHostDmn>::execute(Model::get_r_DCA_basis(), DomainsParameters::get_tp_host()); } Loading Loading
include/dca/phys/dca_analysis/bse_solver/bse_lattice_solver.hpp +6 −4 Original line number Diff line number Diff line Loading @@ -143,7 +143,8 @@ private: DcaDataType& MOMS; func::function<std::complex<ScalarType>, HOST_matrix_dmn_t> Gamma_lattice; func::function<std::complex<ScalarType>, func::dmn_variadic<b_b, b_b, k_HOST_VERTEX, WVertexDmn>> chi_0_lattice; func::function<std::complex<ScalarType>, func::dmn_variadic<b_b, b_b, k_HOST_VERTEX, WVertexDmn>> chi_0_lattice; // Matrix in \vec{k} and \omega_n with the diagonal = chi_0_lattice. func::function<std::complex<ScalarType>, HOST_matrix_dmn_t> chi_0_lattice_matrix; Loading Loading @@ -310,7 +311,7 @@ void BseLatticeSolver<ParametersType, DcaDataType, ScalarType>::computeChi0Latti clustermapping::coarsegraining_tp<ParametersType, k_HOST_VERTEX> coarsegraining_tp(parameters); // DCA+/DCA with post-interpolation: Compute \chi_0 from continuous lattice self-energy. if (parameters.do_dca_plus() || parameters.do_post_interpolation()) { if (parameters.do_dca_plus() || parameters.doPostInterpolation()) { latticemapping::lattice_mapping_sp<ParametersType, k_DCA, k_HOST> lattice_map_sp(parameters); MOMS.Sigma_lattice_interpolated = 0.; Loading Loading @@ -373,7 +374,7 @@ void BseLatticeSolver<ParametersType, DcaDataType, ScalarType>::computeGammaLatt std::cout << "\n" << __FUNCTION__ << std::endl; // DCA+/DCA with post-interpolation: Compute Gamma_lattice with continuous momentum dependence. if (parameters.do_dca_plus() || parameters.do_post_interpolation()) { if (parameters.do_dca_plus() || parameters.doPostInterpolation()) { latticemapping::lattice_mapping_tp<ParametersType, k_DCA, k_HOST_VERTEX> lattice_map_tp( parameters); lattice_map_tp.execute(Gamma_cluster, Gamma_lattice); Loading Loading @@ -814,7 +815,8 @@ void BseLatticeSolver<ParametersType, DcaDataType, ScalarType>::printOnShell() { std::complex<ScalarType> norm = 0; for (int j = 0; j < k_HOST_VERTEX::dmn_size(); j++) { scal_prod += conj(psi_k(j, l)) * leading_eigenvectors(i, 0, 0, j, WVertexDmn::dmn_size() / 2); scal_prod += conj(psi_k(j, l)) * leading_eigenvectors(i, 0, 0, j, WVertexDmn::dmn_size() / 2); norm += conj(leading_eigenvectors(i, 0, 0, j, WVertexDmn::dmn_size() / 2)) * leading_eigenvectors(i, 0, 0, j, WVertexDmn::dmn_size() / 2); } Loading
include/dca/phys/dca_data/dca_data.hpp +1 −1 Original line number Diff line number Diff line Loading @@ -452,7 +452,7 @@ void DcaData<Parameters>::write(Writer& writer) { if (parameters_.dump_lattice_self_energy()) { if (parameters_.do_dca_plus()) writer.execute(Sigma_lattice); else if (parameters_.do_post_interpolation()) else if (parameters_.doPostInterpolation()) writer.execute(Sigma_lattice_interpolated); } Loading
include/dca/phys/dca_loop/dca_loop.hpp +1 −1 Original line number Diff line number Diff line Loading @@ -314,7 +314,7 @@ void DcaLoop<ParametersType, DcaDataType, MCIntegratorType>::perform_lattice_map if (concurrency.id() == concurrency.first()) std::cout << "\n\t\t lattice-mapping " << dca::util::print_time(); if (parameters.do_dca_plus() || parameters.do_post_interpolation()) { if (parameters.do_dca_plus() || parameters.doPostInterpolation()) { if (parameters.hts_approximation()) { DcaDataType MOMS_HTS(parameters); Loading
include/dca/phys/parameters/dca_parameters.hpp +1 −1 Original line number Diff line number Diff line Loading @@ -84,7 +84,7 @@ public: return interacting_orbitals_; } bool do_post_interpolation() const { bool doPostInterpolation() const { return do_post_interpolation_; } Loading
include/dca/phys/parameters/parameters.hpp +1 −1 Original line number Diff line number Diff line Loading @@ -292,7 +292,7 @@ void Parameters<Concurrency, Threading, Profiler, Model, RandomNumberGenerator, KSpHostDmn::parameter_type::print(std::cout); // Host grid for two-particle functions (tp-lattice) if (do_dca_plus() || do_post_interpolation()) { if (do_dca_plus() || doPostInterpolation()) { domains::cluster_domain_initializer<RTpHostDmn>::execute(Model::get_r_DCA_basis(), DomainsParameters::get_tp_host()); } Loading