Fix factorizations for CPU

Author: Katharine Hyatt

src/factorizations/matrixalgebrakit.jl

@@ -226,15 +226,15 @@ end
 function MAK.initialize_output(::typeof(eigh_full!), t::AbstractTensorMap, ::AbstractAlgorithm)
     V_D = fuse(domain(t))
     T = real(scalartype(t))
-    D = DiagonalTensorMap{T, typeof(V_D), storagetype(t)}(undef, V_D)
+    D = DiagonalTensorMap{T, typeof(V_D), similarstoragetype(t, T)}(undef, V_D)
     V = similar(t, codomain(t) ← V_D)
     return D, V
 end
 
 function MAK.initialize_output(::typeof(eig_full!), t::AbstractTensorMap, ::AbstractAlgorithm)
     V_D = fuse(domain(t))
     Tc = complex(scalartype(t))
-    D = DiagonalTensorMap{Tc, typeof(V_D), storagetype(t)}(undef, V_D)
+    D = DiagonalTensorMap{Tc, typeof(V_D), similarstoragetype(t, Tc)}(undef, V_D)
     V = similar(t, Tc, codomain(t) ← V_D)
     return D, V
 end
@@ -248,7 +248,7 @@ end
 function MAK.initialize_output(::typeof(eig_vals!), t::AbstractTensorMap, alg::AbstractAlgorithm)
     V_D = fuse(domain(t))
     Tc = complex(scalartype(t))
-    return D = DiagonalTensorMap{Tc, typeof(V_D), storagetype(t)}(undef, V_D)
+    return D = DiagonalTensorMap{Tc, typeof(V_D), similarstoragetype(t, Tc)}(undef, V_D)
 end
 
 # QR decomposition