update

issp-center-dev · Nov 10, 2023 · c29d357 · c29d357
1 parent d9966ba
commit c29d357
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 0 deletions.
diff --git a/src/iTPS/core/ctm_single.cpp b/src/iTPS/core/ctm_single.cpp
@@ -220,6 +220,8 @@ void Calc_projector_left_block_single(const tensor &C1, const tensor &C4,
       int cut = std::min(std::min(std::min(peps_parameters.CHI, e78 * t41), e12 * t12), e56 * t34);
       s_c.resize(cut);
       for (int i = 0; i < cut; ++i) {
+        // DEBUG
+        std::cerr << "DEBUG s " << i << " " << s[i] / denom << " " << peps_parameters.Inverse_projector_cut << std::endl;
         if (s[i] / denom > peps_parameters.Inverse_projector_cut) {
           s_c[i] = 1.0 / sqrt(s[i]);
         } else {
@@ -357,6 +359,8 @@ void Calc_projector_updown_blocks_single(
       s_c.resize(cut);
 
       for (int i = 0; i < s.size(); ++i) {
+        // DEBUG
+        std::cerr << "DEBUG s " << i << " " << s[i] / denom << " " << peps_parameters.Inverse_projector_cut << std::endl;
         if (s[i] / denom > peps_parameters.Inverse_projector_cut) {
           s_c[i] = 1.0 / sqrt(s[i]);
         } else {

diff --git a/src/iTPS/finite_temperature.cpp b/src/iTPS/finite_temperature.cpp
@@ -21,6 +21,8 @@ namespace itps {
 
 template <class tensor>
 void iTPS<tensor>::finite_temperature() {
+  // DEBUG
+  std::cerr << std::scientific << std::setprecision(std::numeric_limits<double>::digits10);
   double beta = 0.0;
   measure_density(beta, "FT_");