Merge pull request #16 from d-tk/public

Separate complex and real wrappers for ev method.
michael-lehn · Nov 25, 2015 · 0c855f6 · 0c855f6
2 parents 0f8dbe9 + a46d077
commit 0c855f6
Show file tree

Hide file tree

Showing 2 changed files with 23 additions and 6 deletions.
diff --git a/flens/lapack/ge/ev.h b/flens/lapack/ge/ev.h
@@ -160,13 +160,20 @@ template <typename MA, typename VW, typename MVL, typename MVR>
 //
 template <typename MA, typename VW>
     typename RestrictTo<(IsComplexGeMatrix<MA>::value
-                      && IsComplexDenseVector<VW>::value)
-             ||         (IsRealGeMatrix<MA>::value
-                      && IsRealDenseVector<VW>::value),
+                      && IsComplexDenseVector<VW>::value),
              typename RemoveRef<MA>::Type::IndexType>::Type
     ev(MA       &&A,
        VW       &&w);
 
+template <typename MA, typename VWR, typename VWI>
+    typename RestrictTo<(IsRealGeMatrix<MA>::value
+                      && IsRealDenseVector<VWR>::value
+                      && IsRealDenseVector<VWI>::value),
+             typename RemoveRef<MA>::Type::IndexType>::Type
+    ev(MA       &&A,
+       VWR      &&wr,
+       VWI      &&wi);
+
 } } // namespace lapack, flens
 
 #endif // FLENS_LAPACK_GE_EV_H
diff --git a/flens/lapack/ge/ev.tcc b/flens/lapack/ge/ev.tcc
@@ -1360,16 +1360,26 @@ ev(bool     computeVL,
 
 template <typename MA, typename VW>
 typename RestrictTo<(IsComplexGeMatrix<MA>::value
-                  && IsComplexDenseVector<VW>::value)
-         ||         (IsRealGeMatrix<MA>::value
-                  && IsRealDenseVector<VW>::value),
+                  && IsComplexDenseVector<VW>::value),
          typename RemoveRef<MA>::Type::IndexType>::Type
 ev(MA       &&A,
    VW       &&w)
 {
     return ev(false, false, A, w, A, A);
 }
 
+template <typename MA, typename VWR, typename VWI>
+typename RestrictTo<(IsRealGeMatrix<MA>::value
+                  && IsRealDenseVector<VWR>::value
+                  && IsRealDenseVector<VWI>::value),
+         typename RemoveRef<MA>::Type::IndexType>::Type
+ev(MA       &&A,
+   VWR      &&wr,
+   VWI      &&wi)
+{
+    return ev(false, false, A, wr, wi, A, A);
+}
+
 } } // namespace lapack, flens
 
 #endif // FLENS_LAPACK_GE_EV_TCC