Promote to at least Float64 in init for now

dahong67 · Mar 5, 2024 · bafb1bb · bafb1bb
1 parent 774dc56
commit bafb1bb
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/src/GCPDecompositions.jl b/src/GCPDecompositions.jl
@@ -112,6 +112,7 @@ default_init(X, r, loss, constraints, algorithm) =
 function default_init(rng, X, r, loss, constraints, algorithm)
     # Generate CPD with random factors
     T, N = nonmissingtype(eltype(X)), ndims(X)
+    T = promote_type(T, Float64)
     M = CPD(ones(T, r), rand.(rng, T, size(X), r))
 
     # Normalize