all_close.hpp 5.25 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
//*****************************************************************************
// Copyright 2017-2018 Intel Corporation
//
// 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.
//*****************************************************************************
Scott Cyphers's avatar
Scott Cyphers committed
16 17 18

#pragma once

19
#include <cmath>
Scott Cyphers's avatar
Scott Cyphers committed
20
#include <memory>
21
#include <vector>
Scott Cyphers's avatar
Scott Cyphers committed
22

23
#include "ngraph/type/element_type.hpp"
24
#include "test_tools.hpp"
Scott Cyphers's avatar
Scott Cyphers committed
25 26 27 28 29

namespace ngraph
{
    namespace test
    {
30 31 32 33 34 35
        /// \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|.
Robert Kimball's avatar
Robert Kimball committed
36
        template <typename T>
37 38 39 40 41
        typename std::enable_if<std::is_signed<T>::value, bool>::type
            all_close(const std::vector<T>& a,
                      const std::vector<T>& b,
                      T rtol = static_cast<T>(1e-5),
                      T atol = static_cast<T>(1e-8))
Robert Kimball's avatar
Robert Kimball committed
42
        {
43
            bool rc = true;
Robert Kimball's avatar
Robert Kimball committed
44
            assert(a.size() == b.size());
45
            size_t count = 0;
Robert Kimball's avatar
Robert Kimball committed
46 47
            for (size_t i = 0; i < a.size(); ++i)
            {
48 49
                if (std::abs(a[i] - b[i]) > atol + rtol * std::abs(b[i]) || !std::isfinite(a[i]) ||
                    !std::isfinite(b[i]))
Robert Kimball's avatar
Robert Kimball committed
50
                {
51 52 53 54 55
                    if (count < 5)
                    {
                        NGRAPH_INFO << a[i] << " is not close to " << b[i] << " at index " << i;
                    }
                    count++;
56
                    rc = false;
57 58
                }
            }
59 60 61 62
            if (!rc)
            {
                NGRAPH_INFO << "diff count: " << count << " out of " << a.size();
            }
63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87
            return rc;
        }

        /// \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|.
        template <typename T>
        typename std::enable_if<std::is_unsigned<T>::value, bool>::type
            all_close(const std::vector<T>& a,
                      const std::vector<T>& b,
                      T rtol = static_cast<T>(1e-5),
                      T atol = static_cast<T>(1e-8))
        {
            bool rc = true;
            assert(a.size() == b.size());
            for (size_t i = 0; i < a.size(); ++i)
            {
                T abs_diff = (a[i] > b[i]) ? (a[i] - b[i]) : (b[i] - a[i]);
                if (abs_diff > atol + rtol * b[i])
                {
                    NGRAPH_INFO << a[i] << " is not close to " << b[i] << " at index " << i;
                    rc = false;
Robert Kimball's avatar
Robert Kimball committed
88 89
                }
            }
90
            return rc;
Robert Kimball's avatar
Robert Kimball committed
91
        }
92

93 94 95 96 97
        /// \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
Scott Cyphers's avatar
Scott Cyphers committed
98
        /// Returns true if shapes match and for all elements, |a_i-b_i| <= atol + rtol*|b_i|.
99
        template <typename T>
100 101
        bool all_close(const std::shared_ptr<ngraph::runtime::Tensor>& a,
                       const std::shared_ptr<ngraph::runtime::Tensor>& b,
102 103
                       T rtol = 1e-5f,
                       T atol = 1e-8f)
Robert Kimball's avatar
Robert Kimball committed
104 105
        {
            // Check that the layouts are compatible
Scott Cyphers's avatar
Scott Cyphers committed
106
            if (*a->get_tensor_layout() != *b->get_tensor_layout())
Robert Kimball's avatar
Robert Kimball committed
107 108 109
            {
                throw ngraph_error("Cannot compare tensors with different layouts");
            }
Scott Cyphers's avatar
Scott Cyphers committed
110

Robert Kimball's avatar
Robert Kimball committed
111
            if (a->get_shape() != b->get_shape())
112
            {
Robert Kimball's avatar
Robert Kimball committed
113
                return false;
114
            }
Scott Cyphers's avatar
Scott Cyphers committed
115

116
            return all_close(read_vector<T>(a), read_vector<T>(b), rtol, atol);
Robert Kimball's avatar
Robert Kimball committed
117
        }
Scott Cyphers's avatar
Scott Cyphers committed
118

119 120 121 122 123
        /// \brief Same as numpy.allclose
        /// \param as First tensors to compare
        /// \param bs Second tensors to compare
        /// \param rtol Relative tolerance
        /// \param atol Absolute tolerance
Robert Kimball's avatar
Robert Kimball committed
124
        /// Returns true if shapes match and for all elements, |a_i-b_i| <= atol + rtol*|b_i|.
125
        template <typename T>
126 127
        bool all_close(const std::vector<std::shared_ptr<ngraph::runtime::Tensor>>& as,
                       const std::vector<std::shared_ptr<ngraph::runtime::Tensor>>& bs,
128 129
                       T rtol,
                       T atol)
Robert Kimball's avatar
Robert Kimball committed
130 131 132 133 134 135 136 137 138 139 140 141 142 143
        {
            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
144 145
    }
}