allclose_op.h

// Copyright (c) 2018 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.

#pragma once

#include "paddle/fluid/framework/eigen.h"
#include "paddle/fluid/framework/operator.h"
#include "paddle/fluid/platform/place.h"

namespace paddle {
namespace operators {
using Tensor = framework::Tensor;

template <typename DeviceContext, typename T>
class AllcloseKernel : public framework::OpKernel<T> {
 public:
  void Compute(const framework::ExecutionContext& ctx) const override {
    // get attrs
    float rtol = ctx.Attr<float>("rtol");
    float atol = ctx.Attr<float>("atol");
    bool equal_nan = ctx.Attr<bool>("equal_nan");
    // get input/output
    auto* input = ctx.Input<Tensor>("Input");
    auto* other = ctx.Input<Tensor>("Other");
    auto* out = ctx.Output<Tensor>("Out");
    out->mutable_data<bool>(ctx.GetPlace());
    // get place
    auto& place = *ctx.template device_context<DeviceContext>().eigen_device();

    auto input_v = framework::EigenVector<T>::Flatten(*input);
    auto other_v = framework::EigenVector<T>::Flatten(*other);
    auto out_v = framework::EigenScalar<bool>::From(*out);

    auto left = (input_v - other_v).abs();
    auto right = static_cast<T>(atol) + static_cast<T>(rtol) * other_v.abs();
    auto compare_res = left <= right;

    if (equal_nan) {
      auto input_nan = input_v.isnan();
      auto other_nan = other_v.isnan();
      out_v.device(place) =
          (input_nan == other_nan).all() && (compare_res != input_nan).all();
    } else {
      out_v.device(place) = compare_res.all();
    }
  }
};

}  // namespace operators
}  // namespace paddle