CLHEP-154 initialize sigmas before calling prepareUsigmas

ChangeLog

+
+2019-06-11  Lynn Garren <garren@fnal.gov>
+
+  * RandomObjects/src/RandMultiGauss.cc:
+    initialize sigmas before calling prepareUsigmas
+  * for CLHEP-154
+
 ==============================
 11.06.18 Release CLHEP-2.4.1.0
 ==============================

RandomObjects/ChangeLog

+
+2019-06-11  Lynn Garren <garren@fnal.gov>
+
+  * RandomObjects/src/RandMultiGauss.cc:
+    initialize sigmas before calling prepareUsigmas
+  * for CLHEP-154
+
 ==============================
 11.06.18 Release CLHEP-2.4.1.0
 ==============================

RandomObjects/src/RandMultiGauss.cc

@@ -220,7 +220,7 @@ HepVector RandMultiGauss::fire() {
 HepVector RandMultiGauss::fire( const HepVector& mu, const HepSymMatrix& S ) {
 
   HepMatrix U;  
-  HepVector sigmas;
+  HepVector sigmas(mu.num_row());
 
   if (mu.num_row() == S.num_row()) {
     prepareUsigmas ( S, U, sigmas );
@@ -260,7 +260,7 @@ void RandMultiGauss::fireArray( const int size, HepVector* array,
   // on that U and sigmas.
 
   HepMatrix U;  
-  HepVector sigmas;
+  HepVector sigmas(mu.num_row());
   HepVector mu_ (mu);
 
   if (mu.num_row() == S.num_row()) {