mul_compute.cc 2.7 KB
Newer Older
S
superjomn 已提交
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37
// Copyright (c) 2019 PaddlePaddle Authors. All Rights Reserved.
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
//
//     http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
// See the License for the specific language governing permissions and
// limitations under the License.

#include <Eigen/Core>
#include "paddle/fluid/lite/core/kernel.h"
#include "paddle/fluid/lite/core/op_registry.h"
#include "paddle/fluid/lite/core/types.h"

namespace paddle {
namespace lite {
namespace kernels {
namespace host {

template <typename T>
void mul_compute_eigen(const T* x, int x_h, int x_w, const T* y, int y_h,
                       int y_w, T* out) {
  using matrix_t =
      Eigen::Matrix<T, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor>;

  Eigen::Map<const matrix_t> X(x, x_h, x_w);
  Eigen::Map<const matrix_t> Y(y, y_h, y_w);
  Eigen::Map<matrix_t> Out(out, x_h, y_w);

  Out = X * Y;
}

S
superjomn 已提交
38
class MulCompute : public KernelLite<TARGET(kHost), PRECISION(kFloat)> {
S
superjomn 已提交
39 40 41 42
 public:
  using param_t = operators::MulParam;

  void Run() override {
43
    auto& param = Param<operators::MulParam>();
44 45 46 47 48 49 50 51 52 53 54 55 56 57
    core::dim2 x_shape(
        {static_cast<int>(
             param.x->dims().Slice(0, param.x_num_col_dims).production()),
         static_cast<int>(
             param.x->dims()
                 .Slice(param.x_num_col_dims, param.x->dims().size())
                 .production())});
    core::dim2 y_shape(
        {static_cast<int>(
             param.y->dims().Slice(0, param.y_num_col_dims).production()),
         static_cast<int>(
             param.y->dims()
                 .Slice(param.y_num_col_dims, param.y->dims().size())
                 .production())});
S
superjomn 已提交
58

59 60
    mul_compute_eigen(param.x->data<float>(), x_shape.x, x_shape.y,  //
                      param.y->data<float>(), y_shape.x, y_shape.y,  //
61
                      param.output->mutable_data<float>());
62 63 64
    LOG(INFO) << "MUL x " << *param.x;
    LOG(INFO) << "MUL W " << *param.y;
    LOG(INFO) << "MUL out " << *param.output;
S
superjomn 已提交
65 66 67 68 69 70 71 72 73 74
  }

  virtual ~MulCompute() = default;
};

}  // namespace host
}  // namespace kernels
}  // namespace lite
}  // namespace paddle

S
superjomn 已提交
75
REGISTER_LITE_KERNEL(mul, kHost, kFloat, kNCHW,
76
                     paddle::lite::kernels::host::MulCompute, def)
S
Superjomn 已提交
77 78 79
    .BindInput("X", {LiteType::GetTensorTy(TARGET(kHost))})
    .BindInput("Y", {LiteType::GetTensorTy(TARGET(kHost))})
    .BindOutput("Out", {LiteType::GetTensorTy(TARGET(kHost))})
S
superjomn 已提交
80
    .Finalize();