[ITensors] Use QR in QN factorize

ITensor · Sep 13, 2023 · 0c59616 · 0c59616
1 parent c49da9e
commit 0c59616
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/src/tensor_operations/matrix_decomposition.jl b/src/tensor_operations/matrix_decomposition.jl
@@ -667,7 +667,7 @@ function factorize(A::ITensor, Linds...; kwargs...)
   might_truncate = !isnothing(cutoff) || maxdim < min(dL, dR)
 
   if isnothing(which_decomp)
-    if !might_truncate && !hasqns(A) && ortho != "none"
+    if !might_truncate && ortho != "none"
       which_decomp = "qr"
     elseif isnothing(cutoff) || cutoff ≤ automatic_cutoff
       which_decomp = "svd"
@@ -685,7 +685,6 @@ function factorize(A::ITensor, Linds...; kwargs...)
   elseif which_decomp == "eigen"
     L, R, spec = factorize_eigen(A, Linds...; kwargs...)
   elseif which_decomp == "qr"
-    hasqns(A) && error("QR factorization of an ITensor with QNs is not yet supported.")
     L, R = factorize_qr(A, Linds...; kwargs...)
     spec = Spectrum(nothing, 0.0)
   else