Indicate no-aliasing when realizing RHS into a buffer.

LTLA · Jul 3, 2024 · 9619dfc · 9619dfc
1 parent 78e6a6c
commit 9619dfc
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 4 deletions.
diff --git a/include/irlba/parallel.hpp b/include/irlba/parallel.hpp
@@ -356,7 +356,7 @@ class ParallelSparseMatrix {
             if constexpr(std::is_same<Right_, EigenVector_>::value) {
                 return rhs;
             } else {
-                work.buffer = rhs;
+                work.buffer.noalias() = rhs;
                 return work.buffer;
             }
         }();
@@ -374,7 +374,7 @@ class ParallelSparseMatrix {
             if constexpr(std::is_same<Right_, EigenVector_>::value) {
                 return rhs;
             } else {
-                work.buffer = rhs;
+                work.buffer.noalias() = rhs;
                 return work.buffer;
             }
         }();

diff --git a/include/irlba/wrappers.hpp b/include/irlba/wrappers.hpp
@@ -247,7 +247,7 @@ class Centered {
             if constexpr(std::is_same<Right_, EigenVector_>::value) {
                 return rhs;
             } else {
-                work.buffer = rhs;
+                work.buffer.noalias() = rhs;
                 return work.buffer;
             }
         }();
@@ -266,7 +266,7 @@ class Centered {
             if constexpr(std::is_same<Right_, EigenVector_>::value) {
                 return rhs;
             } else {
-                work.buffer = rhs;
+                work.buffer.noalias() = rhs;
                 return work.buffer;
             }
         }();