From 8c39827061c12598130e21129f162b6273a2a22b Mon Sep 17 00:00:00 2001
From: Alejandro Gallo <aamsgallo@gmail.com>
Date: Thu, 27 Jan 2022 20:39:00 +0100
Subject: [PATCH] Add type traits

---
 atrip.org | 12 +++++++++++-
 1 file changed, 11 insertions(+), 1 deletion(-)

diff --git a/atrip.org b/atrip.org
index 30be0ca..43a7f04 100644
--- a/atrip.org
+++ b/atrip.org
@@ -8,7 +8,6 @@ The algorithm uses two main data types, the =Slice= and the
 
 ** The slice
 
-
 #+begin_src c++ :tangle (atrip-slice-h)
 #pragma once
 #include <iostream>
@@ -18,9 +17,20 @@ The algorithm uses two main data types, the =Slice= and the
 
 #include <atrip/Tuples.hpp>
 #include <atrip/Utils.hpp>
+#include <atrip/Blas.hpp>
 
 namespace atrip {
 
+namespace traits {
+  template <typename FF> bool isComplex() { return false; };
+  template <> bool isComplex<Complex>() { return true; };
+namespace mpi {
+  template <typename FF> MPI_Datatype datatypeOf(void);
+  template <> MPI_Datatype datatypeOf<double>() { return MPI_DOUBLE; }
+  template <> MPI_Datatype datatypeOf<Complex>() { return MPI_DOUBLE_COMPLEX; }
+}
+}
+
 
 struct Slice {