Heat trace for Laplace type operators with non-scalar symbols