From 3b7de792d5c3628bf9f0c4df64962a76a05cf5ac Mon Sep 17 00:00:00 2001
From: paboyle <paboyle@ph.ed.ac.uk>
Date: Tue, 18 Apr 2017 13:28:04 +0100
Subject: [PATCH] Type comparison in the traits work

---
 lib/simd/Grid_vector_types.h | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/lib/simd/Grid_vector_types.h b/lib/simd/Grid_vector_types.h
index 7087d724..1bafe12e 100644
--- a/lib/simd/Grid_vector_types.h
+++ b/lib/simd/Grid_vector_types.h
@@ -91,10 +91,12 @@ template <> struct is_complex<std::complex<float> > : public std::true_type {};
 template <typename T> using IfReal       = Invoke<std::enable_if<std::is_floating_point<T>::value, int> >;
 template <typename T> using IfComplex    = Invoke<std::enable_if<is_complex<T>::value, int> >;
 template <typename T> using IfInteger    = Invoke<std::enable_if<std::is_integral<T>::value, int> >;
+template <typename T1,typename T2> using IfSame = Invoke<std::enable_if<is_same<T1,T2>::value, int> >;
 
 template <typename T> using IfNotReal    = Invoke<std::enable_if<!std::is_floating_point<T>::value, int> >;
 template <typename T> using IfNotComplex = Invoke<std::enable_if<!is_complex<T>::value, int> >;
 template <typename T> using IfNotInteger = Invoke<std::enable_if<!std::is_integral<T>::value, int> >;
+template <typename T1,typename T2> using IfNotSame = Invoke<std::enable_if<!is_same<T1,T2>::value, int> >;
 
 ////////////////////////////////////////////////////////
 // Define the operation templates functors