diff --git a/src/samplers/algorithm/algo_aitken.h b/src/samplers/algorithm/algo_aitken.h
index c5cac8e..8a37ac9 100644
--- a/src/samplers/algorithm/algo_aitken.h
+++ b/src/samplers/algorithm/algo_aitken.h
@@ -172,10 +172,6 @@ template<typename CONFIG=default_config> class algo_aitken {
                         assert((((present_iter->first.first - present_res_iter->first.first) < std::numeric_limits<time_type>::epsilon()) &&
                                 (present_iter->first.second == present_res_iter->first.second)) || pts_time_res_.empty());
 
-                        if (!pts_time_res_.empty()) {
-                            assert(!residual_l2_norm_.empty());
-                        }
-
                         auto pts_relx_value_iter = std::find_if(present_iter->second.begin(),
                             present_iter->second.end(), [focus](std::pair<point_type, REAL> b) {
                         return normsq(focus - b.first) < std::numeric_limits<REAL>::epsilon();
diff --git a/src/samplers/algorithm/algo_fixed_relaxation.h b/src/samplers/algorithm/algo_fixed_relaxation.h
index f6d94ef..9abf5bc 100644
--- a/src/samplers/algorithm/algo_fixed_relaxation.h
+++ b/src/samplers/algorithm/algo_fixed_relaxation.h
@@ -157,10 +157,6 @@ template<typename CONFIG=default_config> class algo_fixed_relaxation {
                     assert((((present_iter->first.first - present_res_iter->first.first) < std::numeric_limits<time_type>::epsilon()) &&
                             (present_iter->first.second == present_res_iter->first.second)) || pts_time_res_.empty());
 
-                    if (!pts_time_res_.empty()) {
-                        assert(!residual_l2_norm_.empty());
-                    }
-
                     auto pts_relx_val_iter = std::find_if(present_iter->second.begin(),
                         present_iter->second.end(), [focus](std::pair<point_type, REAL> b) {
                     return normsq(focus - b.first) < std::numeric_limits<REAL>::epsilon();