clean benchmarks

Author: vpuri3

benchmarks/tensor.jl

@@ -1,43 +1,48 @@
 using SciMLOperators, LinearAlgebra, BenchmarkTools
 using SciMLOperators: IdentityOperator, ⊗
 
-Id = IdentityOperator{12}()
-A = rand(12,12)
-B = rand(12,12)
-C = rand(12,12)
+N = 12
+K = 100
+Id = IdentityOperator(N)
+A = rand(N,N)
+B = rand(N,N)
+C = rand(N,N)
 
 println("#===============================#")
 println("2D Tensor Products")
 println("#===============================#")
 
 println("⊗(A, B)")
 
-u = rand(12^2, 100)
-v = rand(12^2, 100)
+u = rand(N^2, K)
+v = rand(N^2, K)
 
 T = ⊗(A, B)
 T = cache_operator(T, u)
 
+@btime *($T, $u)
 @btime mul!($v, $T, $u)
 
 println("⊗(I, B)")
 
-u = rand(12^2, 100)
-v = rand(12^2, 100)
+u = rand(N^2, K)
+v = rand(N^2, K)
 
 T = ⊗(Id, B)
 T = cache_operator(T, u)
 
+@btime *($T, $u)
 @btime mul!($v, $T, $u)
 
 println("⊗(A, I)")
 
-u = rand(12^2, 100)
-v = rand(12^2, 100)
+u = rand(N^2, K)
+v = rand(N^2, K)
 
 T = ⊗(A, Id)
 T = cache_operator(T, u)
 
+@btime *($T, $u)
 @btime mul!($v, $T, $u)
 
 println("#===============================#")
@@ -46,25 +51,25 @@ println("#===============================#")
 
 println("⊗(⊗(A, B), C)")
 
-u = rand(12^3, 100)
-v = rand(12^3, 100)
+u = rand(N^3, K)
+v = rand(N^3, K)
 
 T = ⊗(⊗(A, B), C)
 T = cache_operator(T, u)
 
-mul!(v, T, u) # dunny
+@btime *($T, $u)
 @btime mul!($v, $T, $u); #
 
 println("⊗(A, ⊗(B, C))")
 
-u = rand(12^3, 100)
-v = rand(12^3, 100)
+u = rand(N^3, K)
+v = rand(N^3, K)
 
 T = ⊗(A, ⊗(B, C))
 T = cache_operator(T, u)
 
-mul!(v, T, u) # dunny
+@btime *($T, $u)
 @btime mul!($v, $T, $u); #
 
 println("#===============================#")
-nothing
+#