Merge pull request #109 from SpM-lab/terasaki/fix-fit-part2

Use std::copy
SpM-lab · Feb 6, 2025 · 68f8563 · 68f8563
2 parents 36b6ff2 + 68b5d67
commit 68f8563
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 4 deletions.
diff --git a/include/sparseir/sampling.hpp b/include/sparseir/sampling.hpp
@@ -240,7 +240,7 @@ inline Eigen::Tensor<T, N>& div_noalloc_inplace(
 
         // Remove the extra argument '0'
         ldiv_noalloc_inplace<T, N>(flatbuffer, svd, flatarr, workarr);
-
+        std::copy(flatbuffer.data(), flatbuffer.data() + flatbuffer.size(), buffer_perm.data());
         auto inv_perm = getperm<N>(0, dim);
         buffer = buffer_perm.shuffle(inv_perm).eval();
         return buffer;

diff --git a/test/sampling.cxx b/test/sampling.cxx
@@ -263,7 +263,6 @@ TEST_CASE("Two-dimensional TauSampling test", "[sampling]") {
     }
 }
 
-/*
 TEST_CASE("Sampling Tests") {
     double beta = 1.0;
     vector<double> lambdas = {10.0, 42.0};
@@ -333,7 +332,7 @@ TEST_CASE("Sampling Tests") {
                 //std::cout << "gl(0, 0, 0, 0) " << gl(0, 0, 0, 0) << std::endl;
                 //std::cout << "gl(0, 0, 0, 1) " << gl(0, 0, 0, 1) << std::endl;
 
-                //REQUIRE(sparseir::tensorIsApprox(gl_from_tau, gl, 1e-10));
+                REQUIRE(sparseir::tensorIsApprox(gl_from_tau, gl, 1e-10));
             }
 
             // Test evaluate and fit
@@ -346,7 +345,6 @@ TEST_CASE("Sampling Tests") {
         }
     }
 }
-*/
 
 TEST_CASE("Matsubara Sampling Tests") {
     double beta = 1.0;