all_close.hpp 3.63 KB
Newer Older
Scott Cyphers's avatar
Scott Cyphers committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
// ----------------------------------------------------------------------------
// Copyright 2017 Nervana Systems Inc.
// 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
// ----------------------------------------------------------------------------

#pragma once

17
#include <cmath>
Scott Cyphers's avatar
Scott Cyphers committed
18
#include <memory>
19
#include <vector>
Scott Cyphers's avatar
Scott Cyphers committed
20 21 22 23 24 25 26 27

#include "ngraph/runtime/parameterized_tensor_view.hpp"
#include "ngraph/types/element_type.hpp"

namespace ngraph
{
    namespace test
    {
28
        /// @brief Same as numpy.allclose
Robert Kimball's avatar
Robert Kimball committed
29 30
        /// @param a First tensor to compare
        /// @param b Second tensor to compare
31 32
        /// @param rtol Relative tolerance
        /// @param atol Absolute tolerance
Robert Kimball's avatar
Robert Kimball committed
33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49
        /// @returns true if shapes match and for all elements, |a_i-b_i| <= atol + rtol*|b_i|.
        template <typename T>
        bool all_close(const std::vector<T>& a,
                       const std::vector<T>& b,
                       T rtol = 1e-5f,
                       T atol = 1e-8f)
        {
            assert(a.size() == b.size());
            for (size_t i = 0; i < a.size(); ++i)
            {
                if (std::abs(a[i] - b[i]) > atol + rtol * std::abs(b[i]))
                {
                    return false;
                }
            }
            return true;
        }
50

Scott Cyphers's avatar
Scott Cyphers committed
51 52 53 54 55 56
        /// @brief Same as numpy.allclose
        /// @param a First tensor to compare
        /// @param b Second tensor to compare
        /// @param rtol Relative tolerance
        /// @param atol Absolute tolerance
        /// Returns true if shapes match and for all elements, |a_i-b_i| <= atol + rtol*|b_i|.
57 58 59 60 61
        template <typename T>
        bool all_close(const std::shared_ptr<ngraph::runtime::TensorView>& a,
                       const std::shared_ptr<ngraph::runtime::TensorView>& b,
                       T rtol = 1e-5f,
                       T atol = 1e-8f)
Robert Kimball's avatar
Robert Kimball committed
62 63 64 65 66 67
        {
            // Check that the layouts are compatible
            if (*a->get_tensor_view_layout() != *b->get_tensor_view_layout())
            {
                throw ngraph_error("Cannot compare tensors with different layouts");
            }
Scott Cyphers's avatar
Scott Cyphers committed
68

Robert Kimball's avatar
Robert Kimball committed
69 70
            if (a->get_shape() != b->get_shape())
                return false;
Scott Cyphers's avatar
Scott Cyphers committed
71

72
            return all_close(a->get_vector<T>(), b->get_vector<T>(), rtol, atol);
Robert Kimball's avatar
Robert Kimball committed
73
        }
Scott Cyphers's avatar
Scott Cyphers committed
74 75

        /// @brief Same as numpy.allclose
Robert Kimball's avatar
Robert Kimball committed
76 77
        /// @param as First tensors to compare
        /// @param bs Second tensors to compare
Scott Cyphers's avatar
Scott Cyphers committed
78 79
        /// @param rtol Relative tolerance
        /// @param atol Absolute tolerance
Robert Kimball's avatar
Robert Kimball committed
80
        /// Returns true if shapes match and for all elements, |a_i-b_i| <= atol + rtol*|b_i|.
81 82 83 84 85
        template <typename T>
        bool all_close(const std::vector<std::shared_ptr<ngraph::runtime::TensorView>>& as,
                       const std::vector<std::shared_ptr<ngraph::runtime::TensorView>>& bs,
                       T rtol,
                       T atol)
Robert Kimball's avatar
Robert Kimball committed
86 87 88 89 90 91 92 93 94 95 96 97 98 99
        {
            if (as.size() != bs.size())
            {
                return false;
            }
            for (size_t i = 0; i < as.size(); ++i)
            {
                if (!all_close(as[i], bs[i], rtol, atol))
                {
                    return false;
                }
            }
            return true;
        }
Scott Cyphers's avatar
Scott Cyphers committed
100 101
    }
}