source: issm/oecreview/Archive/25834-26739/ISSM-26579-26580.diff

../trunk-jpl/src/c/solutionsequences/solutionsequence_sampling.cpp

@@ -17 +17 @@
         double      rdnumber;
 
         /*Define seed*/
-        normal_distribution_rnd distribution;
+        //normal_distribution_rnd distribution;
         if(seed<0){
                 std::random_device rd;
                 seed = rd();
@@ -28 +28 @@
                 ISSM_MPI_Comm_rank(ISSM_MPI_COMM_WORLD,&my_rank);
                 seed = seed + 783728*my_rank; // change default seed for parallel simulations (by considering an arbitrary shif based on the rank number)
         }
-        distribution.seed(seed);
+        //distribution.seed(seed);
 
         int        *local_indices = NULL;
         IssmDouble *local_vector = NULL;
@@ -37 +37 @@
   int M;
         ppf->GetLocalSize(&M);
         for(int i=0;i<M;i++){
-                rdnumber = distribution.generator();
+                rdnumber = 1.0;//distribution.generator();
                 ppf->SetValue(local_indices[i],rdnumber,INS_VAL);
         }
         ppf->Assemble();