From c29d3575c05334722c82a710839f0e94fbbb775f Mon Sep 17 00:00:00 2001
From: Yuichi Motoyama <y-motoyama@issp.u-tokyo.ac.jp>
Date: Fri, 10 Nov 2023 17:49:15 +0900
Subject: [PATCH] update

---
 src/iTPS/core/ctm_single.cpp    | 4 ++++
 src/iTPS/finite_temperature.cpp | 2 ++
 2 files changed, 6 insertions(+)

diff --git a/src/iTPS/core/ctm_single.cpp b/src/iTPS/core/ctm_single.cpp
index 18be8c88..99a6d857 100644
--- 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
index dc6fcaa5..25e32b94 100644
--- 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_");