increase max size for std::vector

Zilong-Li · Jul 16, 2024 · f939d1c · f939d1c
1 parent 846071f
commit f939d1c
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/src/io.hpp b/src/io.hpp
@@ -385,7 +385,7 @@ inline void chunk_beagle_genotype_likelihoods(const std::unique_ptr<BigAss> & ge
             genome->pos.push_back(markers);
             markers.clear();
             // transpose glchunk into genome->gl then clear it
-            MyFloat1D gl(genome->nsamples * im * 3);
+            MyFloat1D gl((size_t)genome->nsamples * im * 3);
             for(i = 0; i < genome->nsamples; i++)
             {
                 for(j = 0; j < im; j++)
@@ -422,7 +422,7 @@ inline void chunk_beagle_genotype_likelihoods(const std::unique_ptr<BigAss> & ge
         genome->nchunks++;
         genome->chrs.push_back(chr0);
         // transpose glchunk into genome->gl then clear it
-        MyFloat1D gl(genome->nsamples * im * 3);
+        MyFloat1D gl((size_t)genome->nsamples * im * 3);
         for(i = 0; i < genome->nsamples; i++)
         {
             for(j = 0; j < im; j++)

diff --git a/tests/test-admixture.cpp b/tests/test-admixture.cpp
@@ -64,7 +64,7 @@ TEST_CASE("phaseless naive vs dump dataset 1", "[test-phaseless]")
         for(int i = 0; i < genome->nsamples; i++) llike2 += admixer2.runNativeWithBigAss(i, genome);
         admixer2.updateIteration();
         cerr << std::setprecision(6) << "llike2: " << llike2 << "\tllike1: " << llike1 << "\n";
-        REQUIRE(abs(llike1 - llike2) < 1e-10);
+        REQUIRE(abs(llike1 - llike2) < 1e-6);
     }
     REQUIRE(((admixer1.Q - admixer2.Q).abs() < 1e-6).all());
 }