diff --git a/src/api/ExaGeoStat.cpp b/src/api/ExaGeoStat.cpp
index f739ef83..7abc2b46 100644
--- a/src/api/ExaGeoStat.cpp
+++ b/src/api/ExaGeoStat.cpp
@@ -27,6 +27,8 @@ using namespace exageostat::configurations;
 template<typename T>
 void ExaGeoStat<T>::ExaGeoStatLoadData(Configurations &aConfigurations, std::unique_ptr<ExaGeoStatData<T>> &aData) {
 
+    int seed = 0;
+    std::srand(seed);
     aConfigurations.PrintSummary();
     LOGGER("** ExaGeoStat data generation/loading **")
     // Register and create a kernel object
diff --git a/src/linear-algebra-solvers/concrete/tlr/HicmaImplementation.cpp b/src/linear-algebra-solvers/concrete/tlr/HicmaImplementation.cpp
index 03e8ca2b..da7e378e 100644
--- a/src/linear-algebra-solvers/concrete/tlr/HicmaImplementation.cpp
+++ b/src/linear-algebra-solvers/concrete/tlr/HicmaImplementation.cpp
@@ -233,10 +233,6 @@ T HicmaImplementation<T>::ExaGeoStatMLETile(std::unique_ptr<ExaGeoStatData<T>> &
     VERBOSE("LR: Cholesky factorization of Sigma...")
     START_TIMING(time_facto);
 
-#ifdef USE_R
-    ((double *)HICMA_descCrk->mat)[1] = 0;
-#endif
-
     this->ExaGeoStatPotrfTile(EXAGEOSTAT_LOWER, HICMA_descCUV, 0, HICMA_descCD, HICMA_descCrk, max_rank,
                               pow(10, -1.0 * acc));