1// This file is part of Eigen, a lightweight C++ template library
2// for linear algebra.
3//
4// Copyright (C) 2008-2014 Gael Guennebaud <gael.guennebaud@inria.fr>
5//
6// This Source Code Form is subject to the terms of the Mozilla
7// Public License v. 2.0. If a copy of the MPL was not distributed
8// with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
9
10#ifndef EIGEN_SPARSE_FUZZY_H
11#define EIGEN_SPARSE_FUZZY_H
12
13namespace Eigen {
14
15template<typename Derived>
16template<typename OtherDerived>
17bool SparseMatrixBase<Derived>::isApprox(const SparseMatrixBase<OtherDerived>& other, const RealScalar &prec) const
18{
19  const typename internal::nested_eval<Derived,2,PlainObject>::type actualA(derived());
20  typename internal::conditional<bool(IsRowMajor)==bool(OtherDerived::IsRowMajor),
21    const typename internal::nested_eval<OtherDerived,2,PlainObject>::type,
22    const PlainObject>::type actualB(other.derived());
23
24  return (actualA - actualB).squaredNorm() <= prec * prec * numext::mini(actualA.squaredNorm(), actualB.squaredNorm());
25}
26
27} // end namespace Eigen
28
29#endif // EIGEN_SPARSE_FUZZY_H
30