1 files changed, 142 insertions, 0 deletions
diff --git a/src/core/NEON/kernels/assembly/NEGEMMInterleavedMatrixMultiplyWrapper.cpp b/src/core/NEON/kernels/assembly/NEGEMMInterleavedMatrixMultiplyWrapper.cpp
new file mode 100644
index 0000000000..3d42f8a51f
--- /dev/null
+++ b/src/core/NEON/kernels/assembly/NEGEMMInterleavedMatrixMultiplyWrapper.cpp
@@ -0,0 +1,142 @@
+/*
+ * Copyright (c) 2018 ARM Limited.
+ *
+ * SPDX-License-Identifier: MIT
+ *
+ * Permission is hereby granted, free of charge, to any person obtaining a copy
+ * of this software and associated documentation files (the "Software"), to
+ * deal in the Software without restriction, including without limitation the
+ * rights to use, copy, modify, merge, publish, distribute, sublicense, and/or
+ * sell copies of the Software, and to permit persons to whom the Software is
+ * furnished to do so, subject to the following conditions:
+ *
+ * The above copyright notice and this permission notice shall be included in all
+ * copies or substantial portions of the Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+ * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+ * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+ * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+ * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+ * SOFTWARE.
+ */
+
+#include "arm_compute/core/NEON/kernels/assembly/NEGEMMInterleavedMatrixMultiplyWrapper.h"
+
+#include "NEGEMMInterleavedStrategies.h"
+#include "arm_compute/core/Helpers.h"
+#include "arm_compute/core/ITensor.h"
+#include "arm_compute/core/Utils.h"
+#include "arm_compute/core/Validate.h"
+#include "arm_compute/core/WindowIterator.h"
+
+namespace arm_compute
+{
+template <typename To, typename Tr, bool use_dot>
+void NEGEMMInterleavedMatrixMultiplyWrapperTemplate<To, Tr, use_dot>::configure(const ITensor *prepared_a, const ITensor *transformed_b, ITensor *tmp_c, ITensor *c, const Window &block_walker,
+                                                                                const BlockSizes &block_sizes, const INEGEMMWrapperKernel::Params &params, bool b_is_pretransposed, float alpha, float beta, unsigned int max_num_threads)
+{
+    using strategy = typename Kernel<To>::strategy;
+
+    _prepared_a         = prepared_a;
+    _transformed_b      = transformed_b;
+    _tmp_c              = tmp_c;
+    _c                  = c;
+    _block_walker       = block_walker;
+    _block_sizes        = block_sizes;
+    _params             = params;
+    _b_is_pretransposed = b_is_pretransposed;
+    _alpha              = alpha;
+    _beta               = beta;
+
+    auto_init_if_empty(*_tmp_c->info(), c->info()->clone()->set_tensor_shape(TensorShape{ _block_sizes.x_block * strategy::out_height(), max_num_threads }));
+}
+
+template <typename To, typename Tr, bool use_dot>
+void NEGEMMInterleavedMatrixMultiplyWrapperTemplate<To, Tr, use_dot>::transform(const MatrixMultiplyWorkload &wl, const ThreadInfo &info, const Window &batch_window, const Coordinates &start_offset,
+                                                                                const Coordinates &end_offset)
+{
+    using strategy = typename Kernel<To>::strategy;
+
+    strategy           strat(info.cpu_info);
+    TensorAccessor<To> prepared_a(*_prepared_a);
+    TensorAccessor<To> transformed_b(*_transformed_b);
+    TensorAccessor<Tr> c(*_c);
+    TensorAccessor<Tr> tmp_c(*_tmp_c);
+
+    int  prev_batch      = -1;
+    To *a_ptr           = nullptr;
+    auto window_iterator = arm_compute::create_window_iterator(batch_window, start_offset, end_offset, [&](const Coordinates & id)
+    {
+        const unsigned int y     = id.x();
+        const unsigned int batch = id.y();
+        const unsigned int ymax  = std::min(_params.M, y + strategy::out_height());
+
+        // If it's the first block of a new batch then reset the pointer to A.
+        if(prev_batch != static_cast<int>(batch))
+        {
+            const unsigned int first_m = id.x();
+            a_ptr                      = prepared_a(0, first_m, batch);
+            prev_batch                 = batch;
+        }
+
+        // Call matrix multiply assembly routine to process the block:
+        strat.kernel(a_ptr, transformed_b(wl._offset_transformed_b), tmp_c(0, info.thread_id), 1, wl._bblocks, wl._kern_k);
+        a_ptr += strategy::out_height() * wl._kern_k;
+
+        // Merge the result with the other blocks' results:
+        strat.transforms.Merge(c(0, 0, batch, wl._multi), tmp_c(0, info.thread_id), c.stride(1), y, ymax, wl._x0, wl._xmax, _alpha, (wl._k0 == 0 ? _beta : static_cast<Tr>(1)));
+    });
+    auto on_new_row_size = [&](unsigned int start, unsigned int end)
+    {
+        //Nothing to do
+    };
+    window_iterator.iterate_2D(on_new_row_size);
+}
+
+template <typename To, typename Tr, bool use_dot>
+void NEGEMMInterleavedMatrixMultiplyWrapperTemplate<To, Tr, use_dot>::create_workloads(std::vector<MatrixMultiplyWorkload> &workloads)
+{
+    using strategy = typename Kernel<To>::strategy;
+
+    unsigned int offset_transformed_b = 0;
+    execute_window_loop(_block_walker, [&](const Coordinates & id)
+    {
+        const unsigned int x0    = id.x();
+        const unsigned int k0    = id.y();
+        const unsigned int multi = id.z();
+
+        const unsigned int xmax = std::min(x0 + _block_walker.x().step(), _params.N);
+        const unsigned int kmax = std::min(k0 + _block_walker.y().step(), _params.K);
+
+        // Figure out how many "K" the kernel will actually process.
+        const int kern_k  = ceil_to_multiple(kmax - k0, strategy::k_unroll());
+        const int bblocks = DIV_CEIL(xmax - x0, strategy::out_width());
+
+        workloads.push_back(MatrixMultiplyWorkload(offset_transformed_b, x0, xmax, k0, kmax, multi, kern_k, bblocks));
+
+        if(_b_is_pretransposed)
+        {
+            offset_transformed_b += bblocks * strategy::out_width() * kern_k;
+        }
+        else
+        {
+            ARM_COMPUTE_ERROR("Not supported");
+        }
+    });
+}
+
+//TODO: regroup somewhere ?
+template class NEGEMMInterleavedMatrixMultiplyWrapperTemplate<float, float>;
+#ifdef __aarch64__
+template class NEGEMMInterleavedMatrixMultiplyWrapperTemplate<uint8_t, uint32_t>;
+template class NEGEMMInterleavedMatrixMultiplyWrapperTemplate<int8_t, int32_t>;
+template class NEGEMMInterleavedMatrixMultiplyWrapperTemplate<uint8_t, uint32_t, true>;
+template class NEGEMMInterleavedMatrixMultiplyWrapperTemplate<int8_t, int32_t, true>;
+#endif /* __aarch64__ */
+
+#ifdef __ARM_FEATURE_FP16_VECTOR_ARITHMETIC
+template class NEGEMMInterleavedMatrixMultiplyWrapperTemplate<float16_t, float16_t>;
+#endif /* __ARM_FEATURE_FP16_VECTOR_ARITHMETIC */
+} // namespace arm_compute