[NDTensorsCUDAExt] Remove and test for scalar indexing

NDTensors/src/Unwrap/test/runtests.jl

@@ -174,4 +174,13 @@ include("../../../test/device_list.jl")
   Cp = zero(C)
   Cpt = NDTensors.mul!!(Cp', A', B, true, false)
   @test Cpt' ≈ C
+
+  ##################################
+  ### Add test for transpose(reshape(adjoint )) failure in CUDA
+  A = dev(transpose(reshape(randn(elt, 2, 12)', (12, 2))))
+  B = dev(randn(elt, 2, 2))
+  C = dev(zeros(elt, (2, 12)))
+  NDTensors.mul!!(C, B, A, true, false)
+  Cp = B * copy(A)
+  @test @allowscalar C ≈ Cp
 end

NDTensors/test/diag.jl

@@ -51,12 +51,6 @@ using GPUArraysCore: @allowscalar
     S1 = contract(S, (2, -1), V, (3, 4, -1))
     S2 = contract(dense(S), (2, -1), copy(V), (3, 4, -1))
     @test @allowscalar S1 ≈ S2
-
-    S1M = dev(Matrix{elt}(undef, (2, 12)))
-    SM = array(S)
-    VM = transpose(reshape(array(V), (12, 2)))
-    NDTensors.mul!!(S1M, SM, VM, 1.0, 0.0)
-    @test @allowscalar reshape(S1M, (2, 3, 4)) ≈ array(S1)
   end
 end
 @testset "DiagTensor contractions" begin