remove pointless initialisation of integration step size to max 0.01 ms

C.A.P. Linssen · C.A.P. Linssen · commit 6471863c3a12 · 2024-12-09T11:03:04.000+01:00
diff --git a/models/aeif_cond_alpha.cpp b/models/aeif_cond_alpha.cpp
@@ -384,9 +384,7 @@ nest::aeif_cond_alpha::init_buffers_()
   B_.logger_.reset();
 
   B_.step_ = Time::get_resolution().get_ms();
-
-  // We must integrate this model with high-precision to obtain decent results
-  B_.IntegrationStep_ = std::min( 0.01, B_.step_ );
+  B_.IntegrationStep_ = B_.step_;    // reasonable initial value for numerical integrator step size; this will anyway be overwritten by gsl_odeiv_evolve_apply(), but it might confuse the integrator if it contains uninitialised data
 
   if ( not B_.s_ )
   {
diff --git a/models/aeif_cond_alpha_astro.cpp b/models/aeif_cond_alpha_astro.cpp
@@ -386,9 +386,7 @@ nest::aeif_cond_alpha_astro::init_buffers_()
   B_.logger_.reset();
 
   B_.step_ = Time::get_resolution().get_ms();
-
-  // We must integrate this model with high-precision to obtain decent results
-  B_.IntegrationStep_ = std::min( 0.01, B_.step_ );
+  B_.IntegrationStep_ = B_.step_;    // reasonable initial value for numerical integrator step size; this will anyway be overwritten by gsl_odeiv_evolve_apply(), but it might confuse the integrator if it contains uninitialised data
 
   if ( not B_.s_ )
   {
diff --git a/models/aeif_cond_alpha_multisynapse.cpp b/models/aeif_cond_alpha_multisynapse.cpp
@@ -347,7 +347,7 @@ aeif_cond_alpha_multisynapse::Buffers_::Buffers_( aeif_cond_alpha_multisynapse&
   , c_( nullptr )
   , e_( nullptr )
   , step_( Time::get_resolution().get_ms() )
-  , IntegrationStep_( std::min( 0.01, step_ ) )
+  , IntegrationStep_( step_ )
   , I_stim_( 0.0 )
 {
 }
@@ -416,9 +416,7 @@ aeif_cond_alpha_multisynapse::init_buffers_()
   B_.logger_.reset();
 
   B_.step_ = Time::get_resolution().get_ms();
-
-  // We must integrate this model with high-precision to obtain decent results
-  B_.IntegrationStep_ = std::min( 0.01, B_.step_ );
+  B_.IntegrationStep_ = B_.step_;    // reasonable initial value for numerical integrator step size; this will anyway be overwritten by gsl_odeiv_evolve_apply(), but it might confuse the integrator if it contains uninitialised data
 
   if ( not B_.c_ )
   {
diff --git a/models/aeif_cond_beta_multisynapse.cpp b/models/aeif_cond_beta_multisynapse.cpp
@@ -355,7 +355,7 @@ aeif_cond_beta_multisynapse::Buffers_::Buffers_( aeif_cond_beta_multisynapse& n
   , c_( nullptr )
   , e_( nullptr )
   , step_( Time::get_resolution().get_ms() )
-  , IntegrationStep_( std::min( 0.01, step_ ) )
+  , IntegrationStep_( step_ )
   , I_stim_( 0.0 )
 {
 }
@@ -424,9 +424,7 @@ aeif_cond_beta_multisynapse::init_buffers_()
   B_.logger_.reset();
 
   B_.step_ = Time::get_resolution().get_ms();
-
-  // We must integrate this model with high-precision to obtain decent results
-  B_.IntegrationStep_ = std::min( 0.01, B_.step_ );
+  B_.IntegrationStep_ = B_.step_;    // reasonable initial value for numerical integrator step size; this will anyway be overwritten by gsl_odeiv_evolve_apply(), but it might confuse the integrator if it contains uninitialised data
 
   if ( not B_.c_ )
   {
diff --git a/models/aeif_cond_exp.cpp b/models/aeif_cond_exp.cpp
@@ -379,9 +379,7 @@ nest::aeif_cond_exp::init_buffers_()
   B_.logger_.reset();
 
   B_.step_ = Time::get_resolution().get_ms();
-
-  // We must integrate this model with high-precision to obtain decent results
-  B_.IntegrationStep_ = std::min( 0.01, B_.step_ );
+  B_.IntegrationStep_ = B_.step_;    // reasonable initial value for numerical integrator step size; this will anyway be overwritten by gsl_odeiv_evolve_apply(), but it might confuse the integrator if it contains uninitialised data
 
   if ( not B_.s_ )
   {
diff --git a/models/aeif_psc_alpha.cpp b/models/aeif_psc_alpha.cpp
@@ -374,9 +374,7 @@ nest::aeif_psc_alpha::init_buffers_()
   B_.logger_.reset();
 
   B_.step_ = Time::get_resolution().get_ms();
-
-  // We must integrate this model with high-precision to obtain decent results
-  B_.IntegrationStep_ = std::min( 0.01, B_.step_ );
+  B_.IntegrationStep_ = B_.step_;    // reasonable initial value for numerical integrator step size; this will anyway be overwritten by gsl_odeiv_evolve_apply(), but it might confuse the integrator if it contains uninitialised data
 
   if ( not B_.s_ )
   {
diff --git a/models/aeif_psc_delta.cpp b/models/aeif_psc_delta.cpp
@@ -351,9 +351,7 @@ nest::aeif_psc_delta::init_buffers_()
   B_.logger_.reset();
 
   B_.step_ = Time::get_resolution().get_ms();
-
-  // We must integrate this model with high-precision to obtain decent results
-  B_.IntegrationStep_ = std::min( 0.01, B_.step_ );
+  B_.IntegrationStep_ = B_.step_;    // reasonable initial value for numerical integrator step size; this will anyway be overwritten by gsl_odeiv_evolve_apply(), but it might confuse the integrator if it contains uninitialised data
 
   if ( not B_.s_ )
   {
diff --git a/models/aeif_psc_delta_clopath.cpp b/models/aeif_psc_delta_clopath.cpp
@@ -417,9 +417,7 @@ nest::aeif_psc_delta_clopath::init_buffers_()
   B_.logger_.reset();
 
   B_.step_ = Time::get_resolution().get_ms();
-
-  // We must integrate this model with high-precision to obtain decent results
-  B_.IntegrationStep_ = std::min( 0.01, B_.step_ );
+  B_.IntegrationStep_ = B_.step_;    // reasonable initial value for numerical integrator step size; this will anyway be overwritten by gsl_odeiv_evolve_apply(), but it might confuse the integrator if it contains uninitialised data
 
   if ( not B_.s_ )
   {
diff --git a/models/aeif_psc_exp.cpp b/models/aeif_psc_exp.cpp
@@ -369,9 +369,7 @@ nest::aeif_psc_exp::init_buffers_()
   B_.logger_.reset();
 
   B_.step_ = Time::get_resolution().get_ms();
-
-  // We must integrate this model with high-precision to obtain decent results
-  B_.IntegrationStep_ = std::min( 0.01, B_.step_ );
+  B_.IntegrationStep_ = B_.step_;    // reasonable initial value for numerical integrator step size; this will anyway be overwritten by gsl_odeiv_evolve_apply(), but it might confuse the integrator if it contains uninitialised data
 
   if ( not B_.s_ )
   {
diff --git a/models/glif_cond.cpp b/models/glif_cond.cpp
@@ -465,7 +465,7 @@ nest::glif_cond::Buffers_::Buffers_( glif_cond& n )
   , c_( nullptr )
   , e_( nullptr )
   , step_( Time::get_resolution().get_ms() )
-  , IntegrationStep_( std::min( 0.01, step_ ) )
+  , IntegrationStep_( step_ )
   , I_( 0.0 )
 {
 }
@@ -534,8 +534,6 @@ nest::glif_cond::init_buffers_()
   B_.logger_.reset();   // includes resize
 
   B_.step_ = Time::get_resolution().get_ms();
-  // We must integrate this model with high-precision to obtain decent results
-  B_.IntegrationStep_ = std::min( 0.01, B_.step_ );
 
   if ( not B_.c_ )
   {

Original file line number	Diff line number	Diff line change
`@@ -384,9 +384,7 @@ nest::aeif_cond_alpha::init_buffers_()`
`384`	`384`	`B_.logger_.reset();`
`385`	`385`
`386`	`386`	`B_.step_ = Time::get_resolution().get_ms();`
`387`		`-`
`388`		`- // We must integrate this model with high-precision to obtain decent results`
`389`		`- B_.IntegrationStep_ = std::min( 0.01, B_.step_ );`
	`387`	`+ B_.IntegrationStep_ = B_.step_; // reasonable initial value for numerical integrator step size; this will anyway be overwritten by gsl_odeiv_evolve_apply(), but it might confuse the integrator if it contains uninitialised data`
`390`	`388`
`391`	`389`	`if ( not B_.s_ )`
`392`	`390`	`{`
Original file line number	Diff line number	Diff line change
`@@ -386,9 +386,7 @@ nest::aeif_cond_alpha_astro::init_buffers_()`
`386`	`386`	`B_.logger_.reset();`
`387`	`387`
`388`	`388`	`B_.step_ = Time::get_resolution().get_ms();`
`389`		`-`
`390`		`- // We must integrate this model with high-precision to obtain decent results`
`391`		`- B_.IntegrationStep_ = std::min( 0.01, B_.step_ );`
	`389`	`+ B_.IntegrationStep_ = B_.step_; // reasonable initial value for numerical integrator step size; this will anyway be overwritten by gsl_odeiv_evolve_apply(), but it might confuse the integrator if it contains uninitialised data`
`392`	`390`
`393`	`391`	`if ( not B_.s_ )`
`394`	`392`	`{`
Original file line number	Diff line number	Diff line change
`@@ -347,7 +347,7 @@ aeif_cond_alpha_multisynapse::Buffers_::Buffers_( aeif_cond_alpha_multisynapse&`
`347`	`347`	`, c_( nullptr )`
`348`	`348`	`, e_( nullptr )`
`349`	`349`	`, step_( Time::get_resolution().get_ms() )`
`350`		`- , IntegrationStep_( std::min( 0.01, step_ ) )`
	`350`	`+ , IntegrationStep_( step_ )`
`351`	`351`	`, I_stim_( 0.0 )`
`352`	`352`	`{`
`353`	`353`	`}`
`@@ -416,9 +416,7 @@ aeif_cond_alpha_multisynapse::init_buffers_()`
`416`	`416`	`B_.logger_.reset();`
`417`	`417`
`418`	`418`	`B_.step_ = Time::get_resolution().get_ms();`
`419`		`-`
`420`		`- // We must integrate this model with high-precision to obtain decent results`
`421`		`- B_.IntegrationStep_ = std::min( 0.01, B_.step_ );`
	`419`	`+ B_.IntegrationStep_ = B_.step_; // reasonable initial value for numerical integrator step size; this will anyway be overwritten by gsl_odeiv_evolve_apply(), but it might confuse the integrator if it contains uninitialised data`
`422`	`420`
`423`	`421`	`if ( not B_.c_ )`
`424`	`422`	`{`
Original file line number	Diff line number	Diff line change
`@@ -355,7 +355,7 @@ aeif_cond_beta_multisynapse::Buffers_::Buffers_( aeif_cond_beta_multisynapse& n`
`355`	`355`	`, c_( nullptr )`
`356`	`356`	`, e_( nullptr )`
`357`	`357`	`, step_( Time::get_resolution().get_ms() )`
`358`		`- , IntegrationStep_( std::min( 0.01, step_ ) )`
	`358`	`+ , IntegrationStep_( step_ )`
`359`	`359`	`, I_stim_( 0.0 )`
`360`	`360`	`{`
`361`	`361`	`}`
`@@ -424,9 +424,7 @@ aeif_cond_beta_multisynapse::init_buffers_()`
`424`	`424`	`B_.logger_.reset();`
`425`	`425`
`426`	`426`	`B_.step_ = Time::get_resolution().get_ms();`
`427`		`-`
`428`		`- // We must integrate this model with high-precision to obtain decent results`
`429`		`- B_.IntegrationStep_ = std::min( 0.01, B_.step_ );`
	`427`	`+ B_.IntegrationStep_ = B_.step_; // reasonable initial value for numerical integrator step size; this will anyway be overwritten by gsl_odeiv_evolve_apply(), but it might confuse the integrator if it contains uninitialised data`
`430`	`428`
`431`	`429`	`if ( not B_.c_ )`
`432`	`430`	`{`
Original file line number	Diff line number	Diff line change
`@@ -379,9 +379,7 @@ nest::aeif_cond_exp::init_buffers_()`
`379`	`379`	`B_.logger_.reset();`
`380`	`380`
`381`	`381`	`B_.step_ = Time::get_resolution().get_ms();`
`382`		`-`
`383`		`- // We must integrate this model with high-precision to obtain decent results`
`384`		`- B_.IntegrationStep_ = std::min( 0.01, B_.step_ );`
	`382`	`+ B_.IntegrationStep_ = B_.step_; // reasonable initial value for numerical integrator step size; this will anyway be overwritten by gsl_odeiv_evolve_apply(), but it might confuse the integrator if it contains uninitialised data`
`385`	`383`
`386`	`384`	`if ( not B_.s_ )`
`387`	`385`	`{`
Original file line number	Diff line number	Diff line change
`@@ -374,9 +374,7 @@ nest::aeif_psc_alpha::init_buffers_()`
`374`	`374`	`B_.logger_.reset();`
`375`	`375`
`376`	`376`	`B_.step_ = Time::get_resolution().get_ms();`
`377`		`-`
`378`		`- // We must integrate this model with high-precision to obtain decent results`
`379`		`- B_.IntegrationStep_ = std::min( 0.01, B_.step_ );`
	`377`	`+ B_.IntegrationStep_ = B_.step_; // reasonable initial value for numerical integrator step size; this will anyway be overwritten by gsl_odeiv_evolve_apply(), but it might confuse the integrator if it contains uninitialised data`
`380`	`378`
`381`	`379`	`if ( not B_.s_ )`
`382`	`380`	`{`
Original file line number	Diff line number	Diff line change
`@@ -351,9 +351,7 @@ nest::aeif_psc_delta::init_buffers_()`
`351`	`351`	`B_.logger_.reset();`
`352`	`352`
`353`	`353`	`B_.step_ = Time::get_resolution().get_ms();`
`354`		`-`
`355`		`- // We must integrate this model with high-precision to obtain decent results`
`356`		`- B_.IntegrationStep_ = std::min( 0.01, B_.step_ );`
	`354`	`+ B_.IntegrationStep_ = B_.step_; // reasonable initial value for numerical integrator step size; this will anyway be overwritten by gsl_odeiv_evolve_apply(), but it might confuse the integrator if it contains uninitialised data`
`357`	`355`
`358`	`356`	`if ( not B_.s_ )`
`359`	`357`	`{`
Original file line number	Diff line number	Diff line change
`@@ -417,9 +417,7 @@ nest::aeif_psc_delta_clopath::init_buffers_()`
`417`	`417`	`B_.logger_.reset();`
`418`	`418`
`419`	`419`	`B_.step_ = Time::get_resolution().get_ms();`
`420`		`-`
`421`		`- // We must integrate this model with high-precision to obtain decent results`
`422`		`- B_.IntegrationStep_ = std::min( 0.01, B_.step_ );`
	`420`	`+ B_.IntegrationStep_ = B_.step_; // reasonable initial value for numerical integrator step size; this will anyway be overwritten by gsl_odeiv_evolve_apply(), but it might confuse the integrator if it contains uninitialised data`
`423`	`421`
`424`	`422`	`if ( not B_.s_ )`
`425`	`423`	`{`
Original file line number	Diff line number	Diff line change
`@@ -369,9 +369,7 @@ nest::aeif_psc_exp::init_buffers_()`
`369`	`369`	`B_.logger_.reset();`
`370`	`370`
`371`	`371`	`B_.step_ = Time::get_resolution().get_ms();`
`372`		`-`
`373`		`- // We must integrate this model with high-precision to obtain decent results`
`374`		`- B_.IntegrationStep_ = std::min( 0.01, B_.step_ );`
	`372`	`+ B_.IntegrationStep_ = B_.step_; // reasonable initial value for numerical integrator step size; this will anyway be overwritten by gsl_odeiv_evolve_apply(), but it might confuse the integrator if it contains uninitialised data`
`375`	`373`
`376`	`374`	`if ( not B_.s_ )`
`377`	`375`	`{`
Original file line number	Diff line number	Diff line change
`@@ -465,7 +465,7 @@ nest::glif_cond::Buffers_::Buffers_( glif_cond& n )`
`465`	`465`	`, c_( nullptr )`
`466`	`466`	`, e_( nullptr )`
`467`	`467`	`, step_( Time::get_resolution().get_ms() )`
`468`		`- , IntegrationStep_( std::min( 0.01, step_ ) )`
	`468`	`+ , IntegrationStep_( step_ )`
`469`	`469`	`, I_( 0.0 )`
`470`	`470`	`{`
`471`	`471`	`}`
`@@ -534,8 +534,6 @@ nest::glif_cond::init_buffers_()`
`534`	`534`	`B_.logger_.reset(); // includes resize`
`535`	`535`
`536`	`536`	`B_.step_ = Time::get_resolution().get_ms();`
`537`		`- // We must integrate this model with high-precision to obtain decent results`
`538`		`- B_.IntegrationStep_ = std::min( 0.01, B_.step_ );`
`539`	`537`
`540`	`538`	`if ( not B_.c_ )`
`541`	`539`	`{`