CompareVectorsByComponent.h

Go to the documentation of this file.

/* This file is part of MoDELib, the Mechanics Of Defects Evolution Library.
 *
 *
 * MoDELib is distributed without any warranty under the 
 * GNU General Public License (GPL) v2 <http://www.gnu.org/licenses/>.
 */
 
#ifndef model_COMPAREVECTORSBYCOMPONENT_H_
#define model_COMPAREVECTORSBYCOMPONENT_H_
 
#include <math.h>       /* fabs */
#include <Eigen/Dense>
#include <type_traits>
#include <NumericPrecision.h>
 
namespace model
{
 
    template<typename T, unsigned int N, typename PrecisionType=T>
    struct CompareVectorsByComponent
    {
        bool operator() (const Eigen::Matrix<T,N,1>& lhs, const Eigen::Matrix<T,N,1>& rhs) const
        {
//            std::cout<<PrecisionType(lhs(0)-rhs(0))<<" "<<NumericPrecision<PrecisionType>::epsilon<<std::endl;
            return ( fabs(PrecisionType(lhs(0)-rhs(0)))<=NumericPrecision<PrecisionType>::epsilon )? CompareVectorsByComponent<T,N-1,PrecisionType>()(lhs.template segment<N-1>(1),rhs.template segment<N-1>(1)) : (lhs(0)<rhs(0));
        }
    };
    
    template<typename T,typename PrecisionType>
    struct CompareVectorsByComponent<T,1,PrecisionType>
    {
        bool operator() (const Eigen::Matrix<T,1,1>& lhs, const Eigen::Matrix<T,1,1>& rhs) const
        {
//            std::cout<<PrecisionType(lhs(0)-rhs(0))<<" "<<NumericPrecision<PrecisionType>::epsilon<<std::endl;
            return ( fabs(PrecisionType(lhs(0)-rhs(0)))<=NumericPrecision<PrecisionType>::epsilon )? false : (lhs(0)<rhs(0));
        }
    };
 
    
//      template<typename T, unsigned int N, typename CastType=T>
//      struct CompareVectorsByComponent
//    {
//              bool operator() (const Eigen::Matrix<T,N,1>& lhs, const Eigen::Matrix<T,N,1>& rhs) const
//        {
//                      return ( ((CastType) lhs(0))==((CastType) rhs(0)) )? CompareVectorsByComponent<T,N-1,CastType>()(lhs.template segment<N-1>(1),rhs.template segment<N-1>(1)) : (lhs(0)<rhs(0));
//              }
//      };
//      
//      template<typename T,typename CastType>
//      struct CompareVectorsByComponent<T,1,CastType>
//    {
//              bool operator() (const Eigen::Matrix<T,1,1>& lhs, const Eigen::Matrix<T,1,1>& rhs) const
//        {
//                      return ((CastType) lhs(0))<((CastType) rhs(0));
//              }
//      };
                
} // namespace model
#endif