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