softmax_op.h 2.0 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
/* Copyright (c) 2016 PaddlePaddle Authors. All Rights Reserve.

   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. */

#pragma once

17
#include "paddle/operators/type_alias.h"
18 19 20 21

namespace paddle {
namespace operators {

Q
qijun 已提交
22
template <typename Place, typename T>
23
class SoftmaxKernel : public OpKernel {
24
public:
25 26 27
  void Compute(const ExecutionContext& context) const override {
    auto input = context.Input<Tensor>(0);
    auto output = context.Output<Tensor>(0);
Q
qijun 已提交
28
    output->mutable_data<T>(context.GetPlace());
Q
qijun 已提交
29

30
    auto logits = EigenMatrix<T>::From(*input);
31
    auto softmax = EigenMatrix<T>::From(*output);
Q
qijun 已提交
32 33 34 35 36 37 38 39 40 41 42

    const int kBatchDim = 0;
    const int kClassDim = 1;

    const int batch_size = logits.dimension(kBatchDim);
    const int num_classes = logits.dimension(kClassDim);

    Eigen::DSizes<int, 1> along_class(kClassDim);
    Eigen::DSizes<int, 2> batch_by_one(batch_size, 1);
    Eigen::DSizes<int, 2> one_by_class(1, num_classes);

Q
qijun 已提交
43 44 45 46 47
    auto shifted_logits = (logits -
                           logits.maximum(along_class)
                               .eval()
                               .reshape(batch_by_one)
                               .broadcast(one_by_class));
Q
qijun 已提交
48

49
    softmax.device(context.GetEigenDevice<Place>()) = shifted_logits.exp();
Q
qijun 已提交
50

51
    softmax.device(context.GetEigenDevice<Place>()) =
Q
qijun 已提交
52 53 54 55 56 57
        (softmax *
         softmax.sum(along_class)
             .inverse()
             .eval()
             .reshape(batch_by_one)
             .broadcast(one_by_class));
58 59 60 61
  }
};
}  // namespace operators
}  // namespace paddle