BUG: initialize X to zero in lsmr Solve

Since x wasn't initialized, behavior was not well defined.
tvercaut · Jul 7, 2016 · 83bd9f7 · 83bd9f7
1 parent 5d83c91
commit 83bd9f7
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/Source/lsmrBase.cxx b/Source/lsmrBase.cxx
@@ -321,6 +321,8 @@ Solve( unsigned int m, unsigned int n, const double * b, double * x )
   std::copy( b, b+m, u );
   std::fill( v, v+n, zero);
   std::fill( w, v+n, zero);
+  std::fill( x, x+n, zero);
+
   this->Scale( m, (-1.0), u );
   this->Aprod1( m, n, x, u );
   this->Scale( m, (-1.0), u );