a*trans(b) for NxN * PxN matrices without transposing b N size of an simd vector