Removed broadcast() from matrix and vector pointers

mkatliar · Oct 6, 2024 · 6114d2e · 6114d2e
1 parent 85d47a1
commit 6114d2e
Show file tree

Hide file tree

Showing 11 changed files with 6 additions and 50 deletions.
diff --git a/include/blast/math/RowColumnVectorPointer.hpp b/include/blast/math/RowColumnVectorPointer.hpp
@@ -68,12 +68,6 @@ namespace blast
         }
 
 
-        SimdVecType broadcast() const noexcept
-        {
-            return ptr_.broadcast();
-        }
-
-
         void store(SimdVecType const& val) const noexcept
         {
             ptr_.store(transposeFlag, val);

diff --git a/include/blast/math/dense/DynamicMatrixPointer.hpp b/include/blast/math/dense/DynamicMatrixPointer.hpp
@@ -81,12 +81,6 @@ namespace blast
         }
 
 
-        SimdVecType broadcast() const noexcept
-        {
-            return *ptr_;
-        }
-
-
         void store(SimdVecType const& val) const noexcept
         {
             val.store(ptr_, AF);

diff --git a/include/blast/math/dense/DynamicVectorPointer.hpp b/include/blast/math/dense/DynamicVectorPointer.hpp
@@ -73,12 +73,6 @@ namespace blast
         }
 
 
-        SimdVecType broadcast() const noexcept
-        {
-            return *ptr_;
-        }
-
-
         void store(IntrinsicType val) const noexcept
         {
             // Non-optimized

diff --git a/include/blast/math/dense/StaticMatrixPointer.hpp b/include/blast/math/dense/StaticMatrixPointer.hpp
@@ -82,12 +82,6 @@ namespace blast
         }
 
 
-        SimdVecType broadcast() const noexcept
-        {
-            return SimdVecType {*ptr_};
-        }
-
-
         void store(SimdVecType const& val) const noexcept
         {
             val.store(ptr_, AF);

diff --git a/include/blast/math/dense/StaticVectorPointer.hpp b/include/blast/math/dense/StaticVectorPointer.hpp
@@ -81,12 +81,6 @@ namespace blast
         }
 
 
-        SimdVecType broadcast() const noexcept
-        {
-            return *ptr_;
-        }
-
-
         void store(SimdVecType val) const noexcept
         {
             if constexpr (S == 1)

diff --git a/include/blast/math/panel/DynamicPanelMatrixPointer.hpp b/include/blast/math/panel/DynamicPanelMatrixPointer.hpp
@@ -95,12 +95,6 @@ namespace blast
         }
 
 
-        SimdVecType broadcast() const noexcept
-        {
-            return SimdVecType {*ptr_};
-        }
-
-
         void store(SimdVecType const& val) const noexcept
         {
             static_assert(AF, "store() implemented only for aligned pointers");

diff --git a/include/blast/math/panel/StaticPanelMatrixPointer.hpp b/include/blast/math/panel/StaticPanelMatrixPointer.hpp
@@ -88,12 +88,6 @@ namespace blast
         }
 
 
-        SimdVecType broadcast() const noexcept
-        {
-            return SimdVecType {*ptr_};
-        }
-
-
         void store(SimdVecType const& val) const noexcept
         {
             static_assert(AF, "store() implemented only for aligned pointers");

diff --git a/include/blast/math/register_matrix/DynamicRegisterMatrix.hpp b/include/blast/math/register_matrix/DynamicRegisterMatrix.hpp
@@ -407,7 +407,7 @@ namespace blast
         #pragma unroll
         for (size_t j = 0; j < N; ++j) if (j < n_)
         {
-            SimdVecType bx = b(0, j).broadcast();
+            SimdVecType const bx = b[0, j];
 
             #pragma unroll
             for (size_t i = 0; i < RM; ++i) // TODO: !!! check i against m
@@ -430,7 +430,7 @@ namespace blast
         #pragma unroll
         for (size_t j = 0; j < N; ++j) if (j < n_)
         {
-            SimdVecType bx = b(0, j).broadcast();
+            SimdVecType bx = b[0, j];
 
             #pragma unroll
             for (size_t i = 0; i < RM; ++i) // TODO: !!! check i against m

diff --git a/include/blast/math/register_matrix/RegisterMatrix.hpp b/include/blast/math/register_matrix/RegisterMatrix.hpp
@@ -571,14 +571,14 @@ namespace blast
                     #pragma unroll
                     for (size_t k = 0; k < j; ++k)
                     {
-                        SimdVecType const a_kj = (~A)(k, j).broadcast();
+                        SimdVecType const a_kj = A[k, j];
 
                         #pragma unroll
                         for (size_t i = 0; i < RM; ++i)
                             v_[i][j] = fnmadd(a_kj, v_[i][k], v_[i][j]);
                     }
 
-                    IntrinsicType const a_jj = (~A)(j, j).broadcast();
+                    SimdVecType const a_jj = A[j, j];
 
                     #pragma unroll
                     for (size_t i = 0; i < RM; ++i)
@@ -759,7 +759,7 @@ namespace blast
             #pragma unroll
             for (size_t j = 0; j < N; ++j)
             {
-                SimdVecType bx = bu(0, j).broadcast();
+                SimdVecType const bx = bu[0, j];
 
                 #pragma unroll
                 for (size_t i = 0; i < ii; ++i)
@@ -796,7 +796,7 @@ namespace blast
                 #pragma unroll
                 for (size_t j = 0; j <= k; ++j)
                 {
-                    SimdVecType ax = au(0, j).broadcast();
+                    SimdVecType const ax = au[0, j];
 
                     #pragma unroll
                     for (size_t i = 0; i < RM; ++i)

diff --git a/include/blast/math/typetraits/MatrixPointer.hpp b/include/blast/math/typetraits/MatrixPointer.hpp
@@ -17,7 +17,6 @@ namespace blast
         p(i, j);
         p[i, j];
         p.load();
-        p.broadcast();
         p.vmove(i);
         p.hmove(j);
         p.spacing();

diff --git a/include/blast/math/typetraits/VectorPointer.hpp b/include/blast/math/typetraits/VectorPointer.hpp
@@ -17,7 +17,6 @@ namespace blast
         p(i);
         p[i];
         p.load();
-        p.broadcast();
         p.spacing();
         p.trans();
         ~p;