Add nocopy twists in tensorproduct

Author: lkdvos

ext/TensorKitChainRulesCoreExt/linalg.jl

@@ -42,15 +42,15 @@ function ChainRulesCore.rrule(::typeof(⊗), A::AbstractTensorMap, B::AbstractTe
             ipA = (codomainind(A), domainind(A))
             pB = (allind(B), ())
             dA = zerovector(A, promote_contract(scalartype(ΔC), scalartype(B)))
-            tB = twist(B, filter(x -> isdual(space(B, x)), allind(B)))
+            tB = _twist_nocopy(B, filter(x -> isdual(space(B, x)), allind(B)))
             dA = tensorcontract!(dA, ΔC, pΔC, false, tB, pB, true, ipA)
             return projectA(dA)
         end
         dB_ = @thunk let
             ipB = (codomainind(B), domainind(B))
             pA = ((), allind(A))
             dB = zerovector(B, promote_contract(scalartype(ΔC), scalartype(A)))
-            tA = twist(A, filter(x -> isdual(space(A, x)), allind(A)))
+            tA = _twist_nocopy(A, filter(x -> isdual(space(A, x)), allind(A)))
             dB = tensorcontract!(dB, tA, pA, true, ΔC, pΔC, false, ipB)
             return projectB(dB)
         end

ext/TensorKitChainRulesCoreExt/tensoroperations.jl

@@ -191,7 +191,7 @@ end
 
 # temporary function to avoid copies when not needed
 # TODO: remove once `twist(t; copy=false)` is defined
-function _twist_nocopy(t, inds)
+function _twist_nocopy(t, inds; kwargs...)
     (BraidingStyle(sectortype(t)) isa Fermionic && !isempty(inds)) || return t
-    return twist(t, inds)
+    return twist(t, inds; kwargs...)
 end