@@ -105,7 +105,9 @@ void DownhillSimplexOptimizer::init()
{
matrix_type simplex(m_numberOfParameters,m_numberOfParameters+1);
- // random initialization of the simplex
+ // previously, we did a random initialization
+ // right now, we rely on our first guess and move it in one dimension by m_scales[i]
+ // as suggested in Numerical Recipes
for(int i = 0; i < static_cast<int>(m_numberOfParameters); i++)
for(int j = 0; j < static_cast<int>(m_numberOfParameters+1); j++)