MetaI: set new Dsigma after new sigma_max

chen-jialan · Jun 27, 2020 · 2ab5a63 · 2ab5a63
1 parent df5d276
commit 2ab5a63
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 2 deletions.
diff --git a/src/isdb/Metainference.cpp b/src/isdb/Metainference.cpp
@@ -1631,7 +1631,10 @@ void Metainference::get_sigma_mean(const unsigned iselect, const double fact, co
         if(optimized_step_==N_optimized_step_) {
           double isteps = 1./static_cast<double>(optimized_step_);
           sigmamax_opt_done_=true;
-          for(unsigned i=0; i<sigma_max_.size(); i++) sigma_max_[i]=sigma_max_est_[i]*isteps*sqrt(dnrep);
+          for(unsigned i=0; i<sigma_max_.size(); i++) {
+            sigma_max_[i]=sigma_max_est_[i]*isteps*sqrt(dnrep);
+            Dsigma_[i] = 0.05*(sigma_max_[i] - sigma_min_[i]);
+          }
         }
       }
       optimized_step_++;

diff --git a/src/isdb/MetainferenceBase.cpp b/src/isdb/MetainferenceBase.cpp
@@ -1398,7 +1398,10 @@ void MetainferenceBase::get_sigma_mean(const double fact, const double var_fact,
         if(optimized_step_==N_optimized_step_) {
           double isteps = 1./static_cast<double>(optimized_step_);
           sigmamax_opt_done_=true;
-          for(unsigned i=0; i<sigma_max_.size(); i++) sigma_max_[i]=sigma_max_est_[i]*isteps*sqrt(dnrep);
+          for(unsigned i=0; i<sigma_max_.size(); i++) {
+            sigma_max_[i]=sigma_max_est_[i]*isteps*sqrt(dnrep);
+            Dsigma_[i] = 0.05*(sigma_max_[i] - sigma_min_[i]);
+          }
         }
       }
       optimized_step_++;