where VK = vec(K), VW=[vec(W1),⋯,vec(WM)]∈ℝ(n×n)×M.