add matrix dct performance test
Y = A * (X * A_transpose) seems to be a little faster
matrix dct: performance improvement: computing A and A_transpose only once!